About the influences of compressibility, heat transfer and pressure gradients in compressible turbulent boundary layers

Abstract This paper presents a comprehensive analysis of the momentum and energy transfer in compressible turbulent boundary layers based on integral identities. By considering data obtained from direct numerical simulations for a wide parameter range, the superordinate influences of compressibility, wall heat transfer and pressure gradient on the terms of the governing equations are identified and visualized. This allows us both to determine to what degree cases corresponding to different Mach number, heat transfer and pressure-gradient conditions have physically comparable behaviour and to design turbulent boundary-layer cases with specific sought-after behaviour. To this end, newly formulated identities for the skin-friction coefficient $c_f$ and the specific heat-transfer coefficient $c_h$ from wall-normal integrals based on the non-dimensional compressible momentum and total-enthalpy equations are derived and evaluated. As the individual terms of the resulting identities stay formally close to the terms of the governing equations, the integral analysis further allows the evaluation of common arguments derived from the ‘strong’ Reynolds analogy from an integral perspective. A particular formulation of the Eckert number $Ec$ is proposed as a similarity parameter, mapping cases with different Mach numbers and wall heat transfer conditions.


Introduction
In compressible turbulent boundary layers (TBLs), the interaction of compressibility, wall heat transfer and pressure gradient results in a complex parameter space. Thus an intuitive statement about how the single terms of the governing equations behave for different conditions is difficult to formulate. As only one example, it is well known that some aspects of our simplified understanding of the analogy between momentum and energy † Email address for correspondence: wenzel@iag.uni-stuttgart.de transfer in TBLs have to be modified for 'strongly' cooled cases; see e.g. Zhang, Duan & Choudhari (2018). However, if cases at various Mach numbers are considered, it is not at all obvious what 'strongly' actually means for each Mach-number case. Hence, as all the influences of compressibility, wall heat transfer and pressure gradient lead to a shift in the relevance of the single terms in the overall transfer processes with respect to each other, the question of comparability between certain cases always arises.
For describing the complex interrelation between the momentum and energy transfer within the TBL, the 'strong' Reynolds analogy (SRA) provides an essential building block, giving us a reasonably good understanding of how the compressible TBL behaves. Nearly all theoretical descriptions of compressible TBLs are related in some way to the idea of the SRA -for example, Morkovin's hypothesis that compressible zero-pressure-gradient (ZPG) TBLs behave similarly to 'strangely' heated incompressible ZPG TBLs for large regions of the boundary layer (Smits & Dussauge 2006). Essentially, the main idea of the SRA is simple. If the Prandtl number Pr = 1, the (dimensionless) x-momentum and total-enthalpy boundary-layer equations can be rewritten in a similar form, Hereinafter, for an arbitrary variable f Reynolds (ensemble) averages are denoted by an overbar,f , density-weighted Favre averages are denoted by a tilde,f = ρf /ρ, and values that do not have a fluctuation component are denoted without an overbar/tilde; fluctuations around the Reynolds and Favre averages are denoted by single and double primes, f = f −f and f = f −f , respectively. The similarity between the momentum and total-enthalpy equations leads to the assumption that individual terms of the momentum and total-enthalpy equations can be directly coupled with each other under a 'strong' assumption of direct proportionality of total-enthalpy and velocity fluctuations, h 0 = Ec k 1 u , where k 1 = −1/(Pr Ec)(q w /τ w ) is a proportionality constant; see Morkovin (1961), Gaviglio (1987) and Zhang et al. (2014). Here,τ w is the mean wall shear stress andq w the wall heat flux in the wall-normal direction y. The assumed proportionality between h 0 and u can be interpreted as a turbulence model for the energy equation, which allows calculation of the mean-temperature profile based on the mean-velocity profile as a linear relation of both profiles; see e.g. Walz (1959). An obvious weakness of the SRA becomes apparent when considering adiabatic walls (q y,w = 0), for which k 1 and thus h 0 becomes zero. This implication has been proven incorrect since the total-temperature fluctuations h 0 are of the same order of magnitude as the temperature fluctuations; see e.g. Guarini et al. (2000). Nevertheless, the errors made in some predictions by the SRA are often surprisingly small but can also become large for selected parameter ranges. In recent years, several improvements have been proposed, which have proven to be very powerful for various flow conditions; see e.g.  and Zhang et al. (2014). For the isolated influence of the (non-hypersonic) Mach number under adiabatic conditions, our present knowledge can describe the most fundamental behaviour reasonably well. According to Morkovin (1961), the effects of compressibility are essentially passive, allowing a large number of compressible fluid-flow quantities to be transformed to incompressible reference cases. As a consequence, Mach-number influences on the SRA are only minor, and the success of the Walz relation in determining the temperature profile from a given velocity profile is almost Mach-number independent.
For compressible TBLs, the most complete studies that work out the effect of the Mach number under adiabatic wall conditions can be found in Guarini et al. (2000), Duan, Beekman & Martin (2011), Bernardini & Pirozzoli (2011) and Wenzel et al. (2018), to mention a few; all those are based on direct numerical simulation (DNS) data.
From a theoretical perspective, the influence of wall heat transfer can imitate a temperature profile that (for large regions of the boundary layer) is comparable to that caused by a change of the Mach number, and vice versa. In contrast to Mach-number influences, however, the effect of wall heat transfer can affect the TBL behaviour to a much greater extent (strong heating/cooling), resulting in conditions that are difficult to predict. Since the temperature gradient changes its sign along the wall-normal direction inside the TBL for cooled cases (negative at the wall, positive further outside), the signs of the conductive and turbulent heat transfer also change within the boundary layer. Consequently, certain characteristics such as the correlation between temperature and velocity fluctuations can behave differently from adiabatic cases, see Zhang et al. (2014), thus challenging the basic assumption made in the SRA of a direct proportionality between velocity and total-enthalpy fluctuations.
In contrast to the effects of Mach number and wall heat transfer, our knowledge about the influence of streamwise pressure gradients on the SRA's validity is scarce. This does not mean that we expect our general understanding of the SRA to not apply in these cases; see Gaviglio (1987). However, since the understanding has to be based on evidence, one cannot simply estimate to what extent it still applies to pressure-gradient cases without specific evidence/validation. For general cases with the combined influence of (strong) pressure gradients and heat transfer, Fernholz & Finley (1980) and Gaviglio (1987) have pointed out that it appears difficult to incorporate both influencing factors adequately. For the special class of self-similar cases with 'moderate' pressure-gradient strength, it was shown by Wenzel et al. (2019) that the basic assumptions of Morkovin (and hence implicitly the SRA) still work in good approximation. This behaviour is supported by the self-similarity analysis in Gibis et al. (2019), which suggests that many conclusions found for adiabatic self-similar cases translate to non-adiabatic cases if the heat-transfer distribution is appropriately designed. These findings are further evidenced in Wenzel et al. (2021), where the Reynolds analogy factor s = 2c h /c f was shown to be only weakly influenced by the local Mach number for the weakly non-adiabatic self-similar cases considered. Here, c f ≡τ w /(ρ e u 2 e /2) is the skin-friction coefficient and c h ≡q w /(ρ e u e c p (T w −T r )) the specific heat-transfer coefficient, i.e. the Stanton number. As the Reynolds analogy factor can be interpreted as an integral measure of the SRA, see § 3, the only weak Mach-number dependence of s was interpreted as an indication that the main principles of the SRA will work comparably to the ZPG cases, at least for the near-adiabatic self-similar pressure-gradient cases considered.

Methodology and objectives of this study
From the current state of the literature, it can be concluded that the isolated influences of Mach number, heat transfer and, in some part, the pressure gradient on the momentum and energy transfer in TBLs are described to a reasonably good extent. If more than one parameter is changed simultaneously, however, it is not apparent how the governing equations (and hence the SRA's validity) are affected. This study's primary goal is to provide a broad view on the behaviour of compressible TBLs in the complex parameter space, visualizing the effects of several influencing parameters on the single terms of the governing equations in a more tangible fashion. This view enables both an evaluation of the comparability of different TBLs under different conditions and the design of TBL cases with specific sought-after behaviour. For this purpose, direct numerical simulations of both sub-and supersonic flat-plate TBLs are conducted at various wall-temperature conditions in both zero and adverse pressure-gradient conditions (self-similar cases). The obtained data are analysed using an integral method, where the momentum and total-enthalpy equations are integrated in the wall-normal direction. Consequently, the dependency of the individual terms on the wall-normal coordinate is eliminated, which in turn allows a scalar and thus compact perspective over the streamwise direction. This approach results in integral identities for c f and c h , which hence quantify the integral effects of the overall transfer processes occurring inside the boundary layer. Both identities rely on the same mathematical approach as the Fukagata, Iwamoto & Kasagi (2002), or FIK, identity, which provides a tool for analysing the c f distribution under both compressible and incompressible conditions. To the authors' best knowledge, the identity for c h is the first of this kind for spatially developing boundary layers. Furthermore, by utilizing the Reynolds analogy factor s = 2c h /c f , the two identities are coupled with each other, allowing insight into the analogy between momentum and energy transfer, and hence some aspects of the SRA.
The study is structured as follows: first, the DNS data are summarized in § 2. To allow for a clear understanding of the integral analysis used, its main idea and its derivation are introduced in § 3. To provide an understanding of the raw data used, the spatial evolutions of the integral identities of c f and c h are described qualitatively in § 4. With the goal of examining the isolated effects of compressibility, heat transfer and pressure gradients on the respective identities, the c f and c h identities are analysed in more detail in § § 5 and 6, respectively. Further remarks about resulting consequences and implications are given in § 7, and the conclusions are summarized in § 8.

Simulation details
This study is based on DNS results from self-similar compressible TBLs with adverse pressure gradients (APGs) and ZPGs for inflow Mach numbers of M ∞,0 = 0.3, 0.5, 0.95 and 2.0. All cases are computed on the DNS set-ups presented in Wenzel et al. (2019), one each for sub-and supersonic cases. To allow for a meaningful determination of c h ≡ q y,w /(ρ e u e c p (T w −T r )), all cases have been recomputed with near-adiabatic fixed wall temperatures T w ; some of the cases have been discussed already in Wenzel et al. (2021).
To enable reliable conclusions to be drawn about the influence of the Eckert number Ec = (γ − 1)M 2 e T e /(T r − T w ) on the investigations, the wall temperature T w was altered to 2K, 10K and 20K above and 10K below the adiabatic wall temperatureT r taken from the fully adiabatic reference cases in Wenzel et al. (2019). For the M ∞,0 = 0.3 and M ∞,0 = 0.95 cases, new adiabatic reference cases were computed. Note that the adiabatic wall temperature is a function of x for cases with pressure gradient, and so is the prescribed wall temperature. Besides, a strongly cooled ZPG case at M ∞ = 2.0 and (T w −T r ) = −308K has been simulated, which was specially designed to verify observed trends; see § 6.
The kinematic Rotta-Clauser parameter β K = (δ * K /τ w )(dp e /dx), constructed with the incompressible form of the displacement thickness δ * K = δ 99 0 (1 −ū/u e ) dy as a length scale, allows a comparison of pressure-gradient effects between compressible and incompressible cases; see Wenzel et al. (2019) and Gibis et al. (2019). Boundary-layer edge values (index e) are determined at the wall-normal position where the spanwise vorticitȳ ω z drops below 10 −6 of its wall value, similar to Spalart & Strelets (2000). Relevant boundary-layer properties, evaluated at the beginning and the end of the self-similar regions used, are summarized in table 1. Given parameters are the Mach number at the boundary-layer edge M e , the Eckert number Ec, the Reynolds numbers Re τ =ρ w u τ δ 99 /μ w and Re θ = ρ e u e θ/μ e , the wall-to-edge temperature ratio T w /T e and the wall-to-adiabatic temperature ratio T w /T r . The displacement and momentum thicknesses used are δ * = δ 99 0 (1 − (ρū)/(ρ e u e )) dy and θ = δ 99 0 (ρū)/(ρ e u e )(1 −ū/u e ) dy, respectively. Details about the spatial extent of the computational domains are identical to those in the fully adiabatic cases, see Wenzel et al. (2019), and the same holds for the numerical method.

Simulation parameters
As for the study in Wenzel et al. (2021), but in contrast to the adiabatic DNS results presented in Wenzel et al. (2019), the temperature at the solid wall is fixed according to a prescribed distribution T w = T w (x). Except for this condition the numerical set-up of the non-adiabatic cases is identical to that of the adiabatic ones; see Wenzel et al. (2019). The reference thermodynamic flow properties are the inflow far-field temperature T ∞,0 = 288.15 K, inflow far-field density ρ ∞,0 = 1.225 kg m −3 , Prandtl number Pr = 0.71, specific gas constant R = 287 J kg −1 K −1 and ratio of specific heats γ = 1.4 and are set equal for all cases. Data averaging is performed over both time and the spanwise direction and does not start before the flow has passed the whole domain at least twice. Time averages were computed over a flow-through time corresponding to at least 250 local boundary-layer thicknesses δ 99 ( t u e /δ 99 = 250, where t is the averaging period), which corresponds to at least 8 eddy-turnover times t u τ /δ 99 for the most restricting APG cases. Both the appropriate convergence of the statistics and the initial transients have been assessed as described in Wenzel et al. (2018) and Wenzel et al. (2019).

Integral analysis -introduction and derivation
Integral identities are valuable tools for the investigation of boundary-layer flows. By a wall-normal integration, mean-flow quantities varying in the wall-normal direction through the boundary layer are reduced to scalar quantities, allowing for a compact characterization of their streamwise evolution. Utilizing integral identities, momentum losses, for instance, can be assigned to different physical loss mechanisms, or they can help us understand and quantify flow-control strategies, among others. In the literature, integral relations have been obtained by integrating different types of balance laws such as momentum (von Kármán 1921;Fukagata et al. 2002) or different types of energy such as mean kinetic energy (Renard & Deck 2016), mean exergy (Arntz, Atinault & Merlen 2015) or mechanical energy (Drela 2009).
The main reason for the variety of identities is the difficulty in physically interpreting the results. In TBLs, for instance, theρ u v stress is the most important stress increasing the shear stress in comparison with laminar flows, and thus should play a dominant role in the overall c f composition. In one of the simplest c f relations, namely the von Kármán momentum integral, however, theρ u v term does not contribute at all. Since theρ u v term simply redistributes momentum through turbulent mixing and feeds the viscous dissipation mechanisms close to the wall, its wall-normal integral is equal to zero (Deck et al. 2014): This shows the difficulty of deriving integral relationships that actually reflect the relevance of important physical mechanisms. In this study, the behaviour of both the x-momentum and the total-enthalpy equations are analysed by using integral relations for c f and c h . To utilize essential aspects of the SRA, the contributing terms of the identities should remain formally very close to those of the x-momentum and the total-enthalpy equations. Moreover, a similar derivation procedure should be applicable to both equations. From this perspective, an integral approach comparable to the FIK identity is advantageous compared with other approaches such as the identity derived by Renard & Deck (2016), henceforth referred to as the RD identity. After briefly introducing the basic idea of the classical FIK identity in the following, the integral identities used are presented.

The FIK identity
The FIK identity (Fukagata et al. (2002)) is based on a threefold repeated integration of the Navier-Stokes x-momentum equation in the wall-normal direction. This introduces weight functions like ( y b − y) to terms such as theρ u v stress, where y b is the non-dimensionalized upper integration bound. These weights enhance the integrand near the wall, ( y b − y)| y=0 = y b , and reduce it in the outer region, ( y b − y)| y=y b = 0. According to its original formulation, the resulting terms of the FIK identity can be associated with four different skin-friction generating mechanisms c f = c L f + c T f + c I f + c t f , a laminar one c L f , a turbulent one c T f , an inhomogeneous one c I f and a transient one c t f , of which the second is the weighted integral of the Reynolds shear-stress 930 A1-6 distribution (3.2); see Fukagata et al. (2002). In recent years, the incompressible FIK identity has been further developed and modified, such as by Mehdi & White (2011), Mehdi et al. (2014, Peet & Sagaut (2009), Bannier, Garnier & Sagaut (2015, Yoon et al. (2016) and Modesti et al. (2018), and has become a common tool used in turbulent flow analysis and control; see e.g. Iwamoto et al. (2005), Kametani & Fukagata (2011), Kametani et al. (2015) and Stroh et al. (2016). In contrast to the incompressible regime, only a few detailed studies are available for compressible flat-plate TBLs. In Gomez, Flutet & Sagaut (2009), a possible compressible form of the FIK identity has been derived and corresponding effects on the mean skin-friction drag generation were studied, but without a full clarification of the compressibility effects on the various contributions as mentioned in Li et al. (2019).
It is a drawback of the FIK identity, however, that some of the contributing terms can easily be misinterpreted (especially for spatially evolving TBLs); see Deck et al. (2014), Renard & Deck (2016) and . One of the key controversies is that there is no simple justification for the threefold repeated integration and thus no physics-based explanation for the linearly weighted Reynolds shear stress (3.2). In the authors' opinion, the problem in justifying the threefold repeated integration is further evidenced by the fact that the threefold repeated integration is not necessary for deriving a valid integral identity for c f . According to Cauchy's formula for repeated integration, n repeated integrations of a continuous function can be rewritten into a single integral, This means that any number of repetitions greater than one would result in a FIK-like identity, since c f already appears as a variable in the equations after the first integration, as shown in Appendix A; variations in the number of repetitions essentially only affect the exponent of the respective weighting factors, as in (3.3).
In the authors' opinion, there are major advantages of using only a twofold repeated integration. For a clear argument, the resulting integral identities for c f and c h are introduced first, before this choice is justified in § 3.4.

The compressible integral identities for c f and c h
In this study, the same mix of Reynolds-averaged and Favre-averaged mean quantities is used as in Huang, Coleman & Bradshaw (1995). The Favre decomposition is applied to the convective terms; for all other terms, the Reynolds-averaged decomposition is applied. To allow for a unique interpretation of the effects of Reynolds number, Mach number and heat transfer, the integral identities are presented in dimensionless form with f = f /f n , wheref is the dimensional quantity. Reference values for non-dimensionalization of f n are defined as ρ n = ρ e , p n = p e , μ n =μ w , u n = u e , L n = δ e , T n =T r − T w and h n = h 0,n = c p (T r − T w ) to ensure the appropriate non-dimensionalization of c f ≡τ w /(ρ e u 2 e /2) and c h ≡q y,w /(ρ e u e c p (T w −T r )). Withh 0 =h + (Ec/2)ũ iũi + (Ec/2) u i u i , and with the Reynolds number, Mach number and Eckert number defined as the simplified two-dimensional (2-D) boundary-layer equations for the x-momentum and total-enthalpy can be written as 3.3. Integral identities Using the assumptions for quasi-2-D boundary-layer flow and a twofold repeated integration of (3.5), the integral identity for c f yields (3.7) The non-dimensionalized upper integration bound y b becomes equal to one if the integration is performed up to the boundary-layer edge δ e ; see also § 3.4. For the integral identity of the Stanton number c h , a twofold repeated integration of the total-enthalpy equation in (3.6) results in (3.8) A detailed derivation of the c f and c h identities for an arbitrary number n of repeated integrations is given in Appendix A. In contrast to the more complete forms of the identities in (A2) and (A4), both the viscous-stress fluctuation terms c VF f and c VF h are neglected in (3.7) and (3.8) for simplification, respectively. However, these terms have a noticeable share near the inflow regions of the simulation domain, where the flow is not yet fully turbulent. For the composition of c h , for instance, c VF h can contribute up to 30 %. For boundary layers rapidly changing in the streamwise direction, these terms therefore should be accounted for. For all figures shown hereafter, the inlet regions are cropped.

Remarks
The following remarks should be made about the benefits and limitations of the integral identities.

General remarks on the FIK identity
The first remark is related to the threefold repeated integration used in the original derivation of the FIK identity; see (3.2). As mentioned, the three integrations are not a requirement in deriving an integral condition for c f or c h , since an arbitrary number of repeated integrations n provides a valid decomposition of c f or c h . Exemplarily evaluated for n = 2, 3, 4 and 10 for the subsonic iZPG +10K M=0.5 case, the streamwise evolution of the integral identity (A2) is evaluated in Appendix B.2. Briefly, increasing the number of integrations n increases the weights of terms that are dominant near the wall. Nevertheless, even though n can be freely chosen from a mathematical point of view, not every choice of n allows a physically plausible interpretation. In the original formulation of the FIK identity, the motivation for the threefold repeated integration is based on the viscous mean stress term, which reads in the dimensional formulation as By an appropriate choice, the viscosity can be normalized to one by introducing the Reynolds number into the dimensionless representation, and the integrand simplifies to ∂/∂y(∂ū/∂y). According to Fukagata et al. (2002), the first integration then results in a force balance, the second leads to the mean velocity profile and the third is akin to obtaining the flow rate. As discussed in Renard & Deck (2016), however, this interpretation is only a consequence of the dimensionless formulation of the identity and therefore not intuitive, since the product of a force (first integral) and a length (due to the second integral) has the dimension of energy rather than velocity in dimensional representation. From the authors' point of view, there are several reasons why using a twofold repeated integration is more appropriate. (i) From a physical point of view, the first integration gives a force/energy balance between the wall and all locations within the boundary layer. The second integration represents its average in the wall-normal direction. In other words, both integral identities simply quantify the wall-normal average of how much momentum/energy is transported towards the wall (positive) or away from the wall (negative) by each term. (ii) Compared with higher-fold integrations, the result is easier to interpret since all terms in the boundary-layer equations appearing with a derivative in the wall-normal direction are integrated without being weighted by a wall-distance-dependent term as in the classical FIK identity. In the present formulation chosen, this holds both for the c T f and the c L f terms in (3.7). The resulting values for c L f and c T f can be therefore simply associated with the area enclosed by theμ∂ū/∂y andρ u v distributions, respectively. Note that the main results derived below concerning the effects of pressure gradients, heat transfer and compressibility are qualitatively identical for n = 2 and n = 3.

Remarks on the upper integration limit
Regardless of the upper integration limit y b , both identities yield a valid decomposition of c f and c h . However, the choice of y b greatly influences the portions of the individual terms with respect to each other. Usually, both the upper integration limit and the reference length for non-dimensionalization are set to δ 99 , which simplifies y d to one in (3.7) and (3.8). An exception is provided by the study of Mehdi et al. (2014), in which the influence of the lower and upper integration limits on the accuracy of the c f calculation is evaluated. For the compressible boundary layers in the present study, it seems reasonable to choose an upper integration bound of y b = δ e ≈ 1.3δ 99 to cover the thermal boundary layer. As this is quite an arbitrary choice, y b has been kept as a variable in the derivation of both identities; its influence on the resulting contributions is further discussed in Appendix B.3. In summary, a lower integration bound y b crops terms at y b , thus increasing the portion of near-wall dominant terms like the c L f term. A higher integration bound decreases the percentage share of near-wall dominant terms. Note that -as long as y b remains at approximately 1 -the particular choice of y b does not affect the qualitative conclusions drawn in the following.

Remarks on the splitting of terms
Variants of the original FIK identity are found in the literature, modified for specific applications. For the original purpose of understanding flow-control strategies and drag-reducing mechanisms, terms are often separated into different components of interest. In experimental studies, where often only the wall-normal profiles are available, also all x-derivatives are usually replaced with y-derivatives to allow for an accurate evaluation. For the formulations chosen in the present study, however, individual components of the momentum and total-enthalpy equations are intended to be easily recognizable in a consistent form. Hence, a further subdivision of the individual terms is omitted. Furthermore, the continuity equation has been used to reformulate the mean-convection terms to avoid large, counteracting components that essentially sum to zero.

Streamwise evolution of c f and c h
In the following sections the integral identities for c f and c h are evaluated. Since all results can be simultaneously affected by compressibility, heat transfer and the pressure gradient, meaningful conclusions can only be derived if the trends observed are isolated and clearly associated with the variation of only one of the influencing parameters. To allow for a better understanding of the raw data used, the spatial evolutions of the integral identities    for c f and c h are described qualitatively in figure 1 first. Here the focus is on working out relevant trends, on which the detailed investigations in § § 5 and 6 are built. Depicted are the subsonic, nearly incompressible iZPG +10K M=0.5 and iAPG +10K β K =0.6 cases on the left-hand side and the supersonic cZPG +20K M=2.0 and cAPG +20K β K =0.6 cases on the right-hand side. The lower abscissa gives the spatial position in the domain normalized with the boundary-layer thickness at the inlet of the domain x/δ 99,0 (δ 99,0 is equal for all cases), and the upper abscissa gives the corresponding Reynolds number Re τ . All contributing terms c L  ; c L f only provides a small portion of the overall c f (especially for high Re; see e.g. Fan et al. 2020). For both the sub-and supersonic ZPG cases in figure 1(a,c), respectively, the spatial-development term c D f constitutes the largest portion. For the suband supersonic APG cases in figure 1(b,d), the c D f distributions are influenced in a different manner by the pressure gradient in regions of approximated self-similarity, increased in panel (b) and strongly decreased in panel (d). The mean-convection term c M f essentially balances the changes in c D f . A comparison of the turbulence terms c T f (blue) reveals almost 930 A1-11 identical behaviour for the sub-and supersonic ZPG cases in panels (a,c) and the sub-and supersonic APG cases in panels (b,d), respectively. The percentage contribution of c T f is therefore expected to be predominantly influenced by the pressure-gradient parameter β K , and only marginally by the local Mach number.

Streamwise evolution of c h
According to (3.8), the results of the c h distributions are depicted in figure 1(e-h). For two reasons, the results are more complex to describe than for the c f distributions. First, the total enthalpy describes the overall change in enthalpy, mean kinetic energy and mean turbulent kinetic energy. Hence, more terms have to be taken into account to capture a complete representation of c h . Second, the c h distributions are strongly influenced by the pressure-gradient parameter β K , the local Mach number and heat transfer. For the subsonic iZPG +10K M=0.5 case in figure 1(e), comparable trends can be found to those for the c f contributions in panel (a); only the relative proportions of c T h and c D h differ. A comparison between the subsonic iZPG +10K M=0.5 and iAPG +10K β K =0.6 cases in panels (e,f ), respectively, shows only a small influence of the streamwise pressure gradient. For the supersonic cases in panels (g,h), the blue c T h contributions are noticeably increased in comparison with the subsonic cases. For the supersonic cZPG +20K M=2.0 case in panel (g), c T h is increased from approximately 60 % of c h to approximately 120 % of c h at Re τ = 450, and for the supersonic cAPG +20K β K =0.6 case in (h) to approximately 140 % of c h . The increase in c T h is mainly balanced by a decrease of the spatial evolution term c D h to negative values. Note that, in addition to the overall terms c L h , c M h , c T h and c D h shown here, the portions of the individual terms also show a complex behaviour; see figure 14.

Detailed evaluation of the c f identity
This section's primary goal is to clearly assign the previously observed variations in the c f identity to the physical effects of compressibility, heat transfer and pressure gradient. In figure 1(a-d) it was shown that the pressure-gradient parameter β K is the decisive influencing parameter for the c f identity. Consequently, all results are plotted against their local β K value in this section. This approach isolates the main influencing factor and allows an evaluation of how much the other factors, compressibility and heat transfer, affect the portions of each component in the c f identity.
Normalized with c f and 2c h , respectively, the ratios of the single components c x f /c f are depicted in figure 2(a) and the ratios of c x f /(2c h ) in figure 2(b). The corresponding c f and c h distributions are depicted in Appendix B.1 in figure 10. In figure 2(a,b) results are plotted as dots with transparent colour for every 20th data point in the streamwise direction in regions of approximate self-similarity. Regions near the inlet of the computational domain, where the streamwise pressure gradient is slowly imposed and the resulting flow is not yet self-similar, are depicted as thick, transparent lines for each case. Since the mean-convection term c M f and the spatial evolution term c D f exhibit counteracting behaviour, the two will be considered as one term  behaviour for the blue subsonic and the red supersonic cases. Hence, to allow for a simple argument, the influence of the pressure gradient is discussed based on the subsonic cases first, and then the implications of the Mach number/heat transfer are examined.
All in all, the blue subsonic results depicted in panel (a) reflect the trends expected from figure 1. The boundary-layer term c L f is only slightly influenced by the pressure gradient's strength and contributes only approximately 2 %-3 % to c f . For ZPG cases at β K = 0, the turbulent-convection term c T f contributes approximately 44 % to c f . For increasing APG strength this portion increases steeply, while an extrapolation for favourable-pressure-gradient cases (β K < 0) suggests a decrease. The sum of c D f and c M f constitutes the remaining part and thus essentially behaves in a contrary manner to c T f with respect to β K . In principle, the explanation of the trends shown is straightforward: the normalization of the individual components c x f with c f essentially corresponds to a normalization with the wall shear stressτ w (inner scaling). Hence, the share of c L f /c f corresponds to the average value of (μ∂ū/∂y)/τ w between 0 ≤ y ≤ y b , and the share of c T f /c f to the average value of (ρ u v )/τ w . As known, the turbulence term (ρ u v )/τ w increases considerably for APG cases compared with ZPG cases, explaining the increasing share of c T f /c f for increasing β K . As shown in Wenzel et al. (2019), the distributions of (μ∂ū/∂y)/τ w are almost unaffected by β K .
Next, the effects of Mach number and heat transfer are assessed. Indicated by the clustering of both sub-and supersonic cases, only a small influence of the local Mach number is visible in panel (a). According to Morkovin's transformation, the distributions ofρ u v /τ w -and hence also its integral version c T f /c f -are Mach-number invariant; compare the good agreement of the sub-and supersonic cases (coloured blue and red) in panels (a) and (a * ) for cases with the same β K . To assess the influence of heat transfer, a comparison between the supersonic near-adiabatic cZPG +20K M=2.0 case (in red) and the strongly cooled supersonic ZPG −308K M=2.0 case (in orange) does not reveal a recognizable influence of the wall temperature. The same holds for the near-adiabatic pressure-gradient cases; however, it cannot be excluded that this changes for strongly heated or cooled conditions. In summary, the influences of Mach number and heat transfer can be seen to act comparably on the individual c x f terms for the cases considered. This causes their relative contributions to the overall c f always to remain approximately the same, even if c f changes significantly; compare the c f and c h distributions given in figure 10.
For the remaining part of this study, it is useful to approximate the shown trends as a function of the local pressure gradient β K . As it is simpler than a direct dependence on β K , an empirical approximation as a function of the Reynolds analogy factor s = f (β K , . . .) is used: see the solid black lines in panel (a). As introduced in more detail in Wenzel et al. (2021), the distribution for s is determined from the analytically derived relation by So (1994), which is based on classical, incompressible self-similarity assumptions; see also panel (b): where Re δ * K,w = 16 000, κ θ = 0.41/Pr t , B = 4.9 and B θ = 3.8. Both A(β K ) and A θ (β K , Pr t ) are obtained from the solutions of differential equations for the wall-normal velocity and temperature profiles; see Mellor & Gibson (1966), So (1994) and Wenzel et al. (2021) for details. The depicted approximations further allow a qualitative assessment of the non-self-similar regions of flow represented by the transparent lines in figure 2(a). Both the c T f and (c D f + c M f ) terms deviate notably from the approximation curves for self-similar regions, for both sub-and supersonic cases. In APG regions, the turbulent contribution c T f has a delayed reaction to the pressure gradient compared with a self-similar case at the same β K ; the sum of (c D f + c M f ) reacts prematurely. A closer look into the history effects in non-self-similar incompressible TBLs can be found in Tanarro, Vinuesa & Schlatter (2020).

Evaluation of c x
f /(2c h ) Next, the results in figure 2 essentially represents a scaling of the c x f /c f distributions depicted in panel (a) with s, and hence the approximations (5.1a-c) become Finally, if the non-self-similar flow regions are assessed for the subsonic cases, the blue transparent lines for 1/s closely follow the black dotted reference curve for self-similar cases. This indicates a synchronous response of both c f and c h on the imposed pressure gradient also in regions of non-self-similar flows, where the individual contributions c x f /c f clearly have not yet converged to a self-similar state. For the supersonic cases (red transparent lines), in contrast, the values for 1/s are higher than the settled values in regions of non-self-similar flows, implying an asynchronous response of c f and c h to the imposed pressure gradient; see also Wenzel et al. (2021).

Detailed evaluation of the c h identity
In contrast to the c x f components of the c f identity, the c x h components of the c h identity are affected by all three influencing parameters; see § 4.2. The following analysis is based on the following assumption/consideration: in the non-dimensionalized form of the total-enthalpy equation chosen, the Eckert number Ec = u 2 e /(c p (T r − T w )) = (γ − 1)M 2 e T e /(T r − T w ) appears as a relevant parameter accounting for temperature effects. Note that this definition differs from the very basic one (in gas dynamics) where just T e is used in the denominator (and thus Ec = (γ − 1)M 2 e ), and from the one commonly used for heat-transfer investigations, where (T w − T e ) is used in the denominator. Since the Ec number incorporates both the Mach number and the wall temperature in the specific dimensionless form chosen, it appears to be a suitable parameter to commonly account for both effects. Its use essentially implies that the same integral behaviour is expected for flow cases with the same Ec, regardless of whether a variation of Ec is caused by a change To get an idea of the Ec-number behaviour in the parameter space, its dependence on both the local Mach number M e and the wall-to-recovery temperature ratio T w /T r is depicted in figure 3(a,b) for cooled and heated cases, respectively, if both γ and r are assumed constant. Both the heated and the cooled curves show the same behaviour, only differing in the Ec number's sign. As depicted, both increased heating and cooling as well as a reduction of M e reduces the absolute Ec number. Note that it becomes impossible to reach arbitrarily low positive Ec numbers for cooled cases and all Mach numbers in panel (a), e.g. reaching Ec 1 at M e = 2. For very high Mach numbers, Ec becomes merely a function of T w /T r (and variations in r if pressure gradients are applied; see Wenzel et al. 2019). Further, with increasing Mach number M e , an increasingly strong heating or cooling (T r − T w ) is necessary to achieve a low Ec number.

Turbulent contribution c T h
For understanding of the integral composition of c h , see (3.8), the turbulent-convection term c T h is of central importance and will be considered first: The enthalpy fluctuations are denoted by c T,1 h and the fluctuations of the mean kinetic and turbulent kinetic energy (triple fluctuations) by c T,2 h and c T,3 h , respectively. Following the idea that the SRA essentially describes the analogy between enthalpy transfer and kinetic-energy transfer, both kinetic-energy terms c T,2 h and c T,3 h are mainly considered as only a single term c T,2+3 h in the following. 6.1.1. Implications from the SRA Since the individual terms of the derived integral identities stay formally close to the terms of the momentum and total-enthalpy equations, implications derived from the SRA (Morkovin relations) can be directly transferred to the turbulence terms of the integral identity of c h . Thus, by comparison with the DNS data, the validity of selected implications of the SRA can be evaluated from an integral perspective, depending on where a selected case is located in the parameter space considered.
As introduced in § 1, the SRA assumes a direct proportionality between the streamwise velocity fluctuations u and the total-enthalpy fluctuations h 0 according to where k 1 = −1/(Pr Ec)(q w /τ w ) takes into account wall heat fluxes in the case of non-adiabatic wall conditions. Multiplication byρv , time averaging and rearrangement yieldsρ which finally predicts the ratio between enthalpy-associated and kinetic-energy-associated fluctuations to be In terms of the integral identity, (6.3) becomes Here, C 1 (β K ) denotes a constant mainly depending on the pressure gradient's strength applied. With the same argument for the integrand of c SRA h and the additional assumption that the ratioq w /τ w essentially resembles the Ec-number-independent Reynolds analogy factor s, the course of c SRA h /c h is expected to behave as a constant C 2 (β K ) only depending on the local pressure-gradient strength, Using (6.5) and the modelling approaches (6.7) and (6.8), the c T,1 h /c h term should depend linearly on the Ec number: Finally, with (6.7) and (6.9), the ratio c in panel (c); the contribution of c T,3 h will be discussed afterwards. Contrary to figure 2, ZPG cases are plotted in red, APG cases in blue, and the strongly cooled supersonic ZPG case in orange. The red and orange coloured ZPG cases are discussed first, the blue coloured pressure-gradient cases afterward.
ZPG cases. Six ZPG cases are considered at M e = 0.3, 0.5, 0.95 and 2.0 with various relative wall temperatures of (T r − T w ) = +2K, +10K, +20K and −308K. As previewed M=0.5 c Z P G (6.14) its distribution is depicted as a black solid line in panel (c). While there are some differences in the constants between (6.13) and (6.14), the characteristic behaviour of the resulting hyperbolas is the same: for Ec → ±∞, the ratio tends towards a constant value; with −1.1059 for the DNS data, this value is approximately 10 % lower than predicted by the SRA. For these cases, both the turbulent heat transfer and the turbulent kinetic-energy transfer are of the same order and act in the opposite wall-normal direction. For Ec → ±0, both distributions exhibit a singularity, when the turbulent heat transfer greatly exceeds the turbulent kinetic-energy transfer. The different constants in the numerators of the Ec-number-dependent terms (1.6235 in (6.13) and 1.4882 in (6.14)) only have a minor influence on the slope of the respective curves. Particular attention should be paid to the Ec-number range 0 < Ec 1.5. Due to the offset of c T,1 h at Ec = 0, both terms c T,1 h and c T,2+3 h are positive in the considered range in panel (b), causing also a positive ratio in panel (c). Hence, for cases with 0 < Ec 1.5, both the turbulent heat transfer and the turbulent kinetic-energy transfer act in the same wall-normal direction; see § 7.1 for further discussion.
Pressure-gradient cases. Next, the blue pressure-gradient cases in figure 4 are considered. Depicted are cases with β K = 0.2, 0.6 and 1.0. For β K = 0.2 and 1.0 only a single subsonic case is shown. For β K = 0.6 four cases are shown, three subsonic ones with (T r − T w ) = +2K, +10K and −10K, and a supersonic one with (T r − T w ) = +20K. Implied by their equal β K value, the four β K = 0.6 cases are expected to exhibit similar behaviour; they are approximated by a single set of linear fits depicted as blue lines in figure 4. In panel (a) (6.18) see the black dashed line in panel (c). As shown in figure 4(b), the slope of the approximations (6.16a,b) differs distinctly between the pressure-gradient and ZPG cases; for the APG cases with β K = 0.6, both slopes are approximately 27 % steeper than in the ZPG cases; compare (6.12a,b) and (6.16a,b). Considering Ec = 0, the zero offset of the c T,1 h /c h contribution (≈0.48) is approximately the same for β K = 0.6 as for the ZPG cases. In contrast to the ZPG cases discussed before, this offset is not correctly balanced by the heat-transfer modelling term c SRA h by the SRA in (6.15) (≈−0.72). Due to the large schematic similarity between the APG and ZPG cases in panel (b), however, the ratio c T at Ec = 0, implying that pressure-gradient effects are not adequately incorporated in the present formulation of the SRA chosen; see § 7 for a more detailed discussion. Nevertheless, from an integral point of view, the underlying schematic, predicting the trends depicted in panel (c), can be assessed to apply to both the ZPG and APG cases considered. It is noteworthy that this holds in good approximation for both regions with self-similar flow (symbols) and regions with non-self-similar flow (thick blue lines).
Briefly summarizing the main result for both the ZPG cases and the APG cases, figure 4 provides a schematic to evaluate how the turbulent transfer processes within the total-enthalpy equation are related to each other for a wide parameter range as a function of the local Ec number. In § 7 some further remarks are made on how far the behaviour shown for the integral progressions of c T,1 h /c T,2(+3) h depicted can be transferred to the local profiles of its integrandsρ v h and Ecũρ u v .  contrast, greatly deviates from the approximation shown. However, it is emphasized that c T h is the sum of two large opposite terms, which makes its distribution very sensitive. Although Mach-number influences have been almost completely eliminated in previous representations by scaling with c f and/or c h , they seem to be of relevance for the quantities considered here, at least for the supersonic pressure-gradient case.

Evaluation of the overall c h identity
In the following, the overall integral identity of c h (3.8) is evaluated, putting the turbulent-convection term c T h in its overall context. All quantities evaluated are normalized with c h in figure 6(a) and with 1/2c f in panel (b).
First, figure 6(a) is discussed. Using (6.19a,b) for c T h /c h , the remaining ZPG approximations depicted are fitted by The boundary-layer term c L h is negligibly small for small absolute Ec numbers, but increases notably with decreasing Ec number; for example, it contributes with approximately 12 % at Ec = −25, a considerable part of c h . Hence, the dominant behaviour occurring with varying Ec number is the variation of are virtually unaffected by the strength of the pressure gradient for the subsonic cases; compare the blue dots and red lines. With represents a multiplication of the c x h /c h distributions depicted in panel (a) by s; the same holds for the approximations (6.19a,b) and (6.20a,b). As the sum of all contributions results in the Reynolds analogy factor s = 2c h /c f , see Wenzel et al. (2021), the individual portions 2c x h /c f quantify their relative contributions to s. While in panel (a) all distributions add up to 100 per cent and thus allow a simple comparison between the different cases, the pressure-gradient and ZPG cases are fanned out in panel (b) due to the strong pressure-gradient dependence of s. This results in a separate set of curves for each β K , adding up to the respective s value; again, only the sets for β K = 0 and β K = 0.6 are approximated by red and blue lines, respectively. Since the  . Symbols represent results evaluated at every 20th grid point in the streamwise direction in regions of self-similarity. ZPG cases are coloured red, APG cases blue and the strongly cooled supersonic case yellow. Red and blue lines linearly approximate ZPG and β K = 0.6 cases.
Reynolds analogy factor is only weakly influenced by the local Ec number, the sum of the respective curves is also virtually invariant to the Ec number for cases with the same β K in near-adiabatic conditions. To what extent the course of the pressure-gradient cases would be affected under strongly non-adiabatic conditions cannot be evaluated based on the available data, however. At least for the strongly cooled supersonic ZPG case, no influence of the wall temperature can be determined.

Discussion
Having presented the integral results in the previous section, some of the resulting consequences and implications are examined in more detail in the following. In § 7.1, trends obtained from the integral analyses are related to local, wall-normal distributions. In § 7.2, results from the integral c f and c h analysis are related to each other. Following the intuition given by Prandtl's mixing-length model, the signs of Ecũρ u v andρ v h depend on the signs of ∂ū/∂y and ∂T/∂y, respectively, making the interpretation of the different cases straightforward. For heated cases (Ec < 0), ∂T/∂y and ∂ū/∂y are of opposite signs throughout the entire boundary layer, implying a negative sign ofρ v h /(Ecũρ u v ). For Ec → −∞ (near adiabatic), the ratio is close to −1, and for Ec → −0, the ratio is reduced; see cases (I) → (III). For cooled cases (Ec > 0), ∂T/∂y and ∂ū/∂y are of the same sign in the near-wall region, also causing regions with a positive sign ofρ v h /(Ecũρ u v ) near the wall. For Ec → +∞ (near adiabatic), the negative regions dominate the positive regions, resulting in a ratio close to −1, whereas for cases with Ec → +0, the positive regions increasingly dominate the negative regions; see cases (IV) → (VI). From an integral point of view, an Ec number of approximately 1.5 represents cases where the wall-proximal positive and the wall-distant negative contributions equalize each other. According to figure 3, this is only (at least theoretically) possible for cases with M e 3.4, where Ec(T w /T r = 0) 1.5. It is mentioned that this switch of the integral terms does not signal an abrupt change in the local behaviour. Note that the topological differences due to the change of the sign of ∂T/∂y between cooled and heated cases are hidden in the integral evaluation. In the integral sense, cooling and heating appear to behave in the same way, and only the sign of the additional energy flux introduced is altered.
To assess to what extent the integral trends found for c T,1 h /c T,2+3 h are transferable to the actual ratio of its integrandsρ v h / (Ecũρ u v ), the wall-normal distributions For cases with high absolute Ec numbers, the Ec-number-dependent term in (7.1) becomes negligible and the resulting distributions stay close to −1 for y/δ 99 0.6. Hence, for high absolute Ec numbers, the integral trends closely resemble the distributions of ρ v h /(Ecũρ u v ). For cases with low absolute Ec numbers, e.g. the iZPG +10K M=0.3 case or the cZPG −308K M=2.0 case, the Ec-number-dependent term in (7.1) 'rounds' the resulting distribution in the near-wall region for y/δ 99 0.2 with a wall-distance-dependent portion depending on 1/ũ; see the red and orange shaded regions in panel (a). Hence, for small absolute Ec numbers, the integral values only estimate the magnitude of the local profiles.
To evaluate (7.1) for pressure gradients, the sum ofρ v h /(Ecũρ u v ) and 1/(PrEc)(q w /τ w )(1/ũ) is depicted in figure 8(b). As all ZPG cases resemble a constant trend close to −1 independent of the Ec number considered, the Morkovin relation (7.1) can be assumed to incorporate effects of Ec-number variations accurately. For the pressure-gradient cases depicted in blue, however, the distributions differ distinctly from the expected value of −1, suggesting that the pressure gradient must be accounted for by an additional term in the wall shear stress in (7.1). In the integral analysis, this portion would increase the value of the constant C 1 (β K ) depicted in figure 4(a).

About the connection between the integral identities of c f and c h
In the previous section the two integral identities for c f and c h were evaluated almost independently of each other. In this section the definition of the Reynolds analogy factor is utilized to link the two identities. In this way, additional insight is gained regarding the weighting of the respective terms of the momentum and total-enthalpy equations to each other from an integral perspective.

About the influence of Ec and β K on c f , c h and s
As previously discussed and shown in Wenzel et al. (2021), the value of s = 2c h /c f strongly depends on the local pressure-gradient strength and only weakly on the local Ec number. Altering the Ec number thus alters c h and c f in a proportional manner, leaving s nearly unaffected. An alteration of the local pressure-gradient strength, conversely, alters c h and c f in a non-proportional manner, so that s varies as a function of the local pressure-gradient strength. As shown in this study, the percentage share of c x f in the total c f primarily depends on the pressure gradient and only weakly on the Ec number. The percentage share of c x h in the total c h , conversely, primarily depends on the Ec number and only weakly on the pressure-gradient strength. As a consequence, the following three scenarios can be distinguished for the cases considered here.

About the relation between c T h and c T f
As introduced in § 1, it is the main assumption of the SRA that the total-enthalpy fluctuations h 0 and the velocity fluctuations u can be directly coupled by the 'strong' assumption of a direct proportionality between h 0 and u , k 1 ≈ h 0 /(Ec u ). Multiplication byρv and time averaging yields or written in terms of the integral identities, if k 1 is assumed to be a constant for all wall-normal positions within the boundary layer. Hence, by evaluating the ratio of the turbulent-convection terms of the c h and the c f M=0.5 c A P G + 2 0 K β , both plotted as a function of the Eckert number Ec = (γ − 1)M 2 e T e /(T r − T w ). Symbols represent results evaluated at every 20th numerical grid point in the streamwise direction in regions of self-similarity. ZPG cases are coloured red, APG cases blue and the strongly cooled supersonic case yellow. Red and blue lines are linear best fits of the ZPG and β K = 0.6 cases, in (b) fitted to the data and in (a) according to (7.5).
identities, (1/Ec)(2c T h /c T f ), a nominal value for the proportionality constant k 1 can be obtained. The value (1/Ec)(2c T h /c T f ) thus represents the value that the proportionality constant of the strong SRA k 1 would have to take, to satisfy the relation of turbulent correlations in the integral sense; see (7.4). Hence, by comparison of this nominal, integral value for k 1 and the directly computed value for k 1 = −1/(Pr Ec)(q w /τ w ) proposed by the SRA, the uncertainty of the 'strong' assumption k 1 ≈ h 0 /(Ec u ) can be assessed from an integral perspective in the considered Ec-number range. As both k 1 in (7.3) and (1/Ec)(2c T h /c T f ) depend on 1/Ec, both distributions are hyperbolas with the same qualitative behaviour as depicted in figure 7. To allow for a better analysis of their respective shapes, both distributions are plotted multiplied by Ec in figure 9, the directly computed values for Ec k 1 in panel (a) and the integral ratio of 2c T h /c T f in panel (b). Using the approximations obtained from the DNS data (5.3a-c) and (6.19a,b) (a), all distributions are constant; for ZPG cases, this value is around 1.18, for APG cases, it is increased. This is in line with the discussion of figure 4(a). In panel (b), the trends for the ZPG and APG cases both exhibit a clear Ec-number dependence, increasing for negative Ec numbers and decreasing for positive Ec numbers. As the c T f distributions in the denominator have been found to be almost Ec-number independent in figure 2, this trend essentially resembles the qualitative trends found for the c T h distributions in figure 5. For ZPG cases at Ec = 0, the value is comparable to panel (a) (black cross); however, for APG cases, it is reduced in contrast to panel (a). This essentially resembles the same systematic error made by the SRA already seen from the discussion of figures 4(c) and 8(b).
Based on the curve progressions of Ec k 1 and 2c T h /c T f shown, the progressions of the hyperbolas of k 1 and (1/Ec)(2c T h /c T f ) are discussed; only ZPG cases are considered in the following. For the progression of k 1 , panel (a) implies a hyperbola of the form k 1 ≈ 1.18/Ec, and thus a hyperbola tending towards ±0 in its limit for Ec → ±∞. This limit is the direct consequence of the assumption by the SRA that total-enthalpy fluctuations (or in our caseρ h 0 v ) become zero for near-adiabatic cases. For the progression of (1/Ec .408/(Ec s) − 0.0944/s. This hyperbola has a comparable shape to that of k 1 , but is slightly shifted downwards by an offset of −0.0944/s over the entire Ec-number range compared with the hyperbola of k 1 . This Ec-number-independent offset ultimately quantifies the total-enthalpy fluctuations integrally accounted for in c T h , which are non-zero for near-adiabatic cases. Thus, this offset can be understood as a quantification of the integral error made by the SRA through its assumption of vanishing total-enthalpy fluctuations for adiabatic conditions. Whereas this modelling error becomes dominant for large absolute Ec numbers for near-adiabatic cases, its influence decreases successively for smaller Ec numbers and thus cases with strong heating or cooling.

Conclusions
To gain deeper insight into the influences of Mach number, heat transfer and pressure gradient on the compressible momentum and total-enthalpy equations, an integral study has been conducted. To this end, newly formulated integral identities for c f and c h have been derived and evaluated based on the dimensionless compressible momentum and total-enthalpy equations, which can be linked by the Reynolds analogy factor s. Both identities use the same mathematical approach as the known FIK identity, see Fukagata et al. (2002), but with a twofold repeated integration instead of a threefold one, making its interpretation physically more sound. From a physical point of view, the first integration essentially gives a force/energy balance between the wall and all locations within the integration domain; the second integration represents its average in the wall-normal direction. For both identities, three contributions are distinguished: a turbulent one; a viscous or conductive one; and the one of the sum of the mean convection and spatial development. As the individual terms of the resulting identities stay formally close to the terms of the governing equations for momentum and total enthalpy, the integral identities directly reflect the behaviour of both equations' terms.
Characterizing the behaviour of the momentum equation, the c f identity essentially condenses the findings already discussed in previous studies in a compact manner; see e.g. Wenzel et al. (2019). With increasing APG strength the turbulent portion of c f increases and the mean-convective and spatial-development portions decrease. Furthermore, no significant dependence on the Mach number or wall temperature is observed (at least for the near-adiabatic cases considered here) if the c f identity is normalized by c f , which essentially corresponds to a normalization with wall quantities and hence Morkovin's transformation for the turbulent terms. As a result, all portions of the c f identity representing the portions of the relevant transfer mechanisms in the momentum equation are (almost) only a function of the pressure-gradient parameter β K = (δ * K /τ w )(dp e /dx). Characterizing the behaviour of the total-enthalpy equation, the interpretation of the c h identity is more complex: first, as it is very sensitive and requires the consideration of small, often neglected terms; and second, as the shares of the c h identity are affected by the local Mach number, heat transfer and, in some part, pressure gradients if normalized by c h . It is a main conclusion of this study that the particular definition of the Eckert number Ec = u 2 e /(c p (T r − T w )) = (γ − 1)M 2 e T e /(T r − T w ) allows us to account for the effects of both Mach number and wall temperature on the boundary layer for a wide parameter range. This essentially implies that flow cases with the same Ec number exhibit the same integral behaviour and are thus physically comparable to each other, regardless of how this Ec number is achieved by adjusting the Mach number or the wall temperature. For hypersonic Mach numbers, mostly the wall-to-recovery temperature determines the Ec number.
Implied by the data considered, all shares of the c h identity (and hence the total-enthalpy equation) show a linear Ec-number dependence, which can be considered a direct consequence of both the SRA and the Morkovin hypotheses; pressure-gradient effects are small. Applying arguments of the SRA to the individual terms of the c h identity, comprehensive predictions can be made about the integral behaviour of the turbulence terms as a function of the local Ec number for a fixed β K . Conversely, this allows us to estimate the accuracy of the approximation of the SRA to the integral behaviour of the boundary layer for a wide parameter range consisting of Mach number, wall heat transfer and pressure gradients. The main conclusions from the c h identity are the following. (i) Plotted over the Ec number, the integral energy transfer caused by the turbulence terms is reflected in an easy-to-understand manner. For adiabatic conditions (Ec → ∓∞), the turbulent contributions of heat and kinetic-energy transfer are of the same order and act in the opposite wall-normal direction. For strongly heated or cooled cases (Ec → ∓0), the ratio tends to ∓∞, implying that the turbulent contribution of heat transfer tends to largely exceed that of the kinetic-energy transfer. (ii) This behaviour is qualitatively well predicted by the SRA for ZPG cases, but only if the often neglected heat-flux term incorporating non-adiabatic wall conditions is taken into account. Quantitatively, the constants found by the SRA deviate slightly from the DNS data. (iii) Due to the counteracting nature of the turbulent contributions of heat and kinetic-energy transfer for adiabatic cases, the triple-correlation term is non-negligible for the complete turbulence term; its significance increases with increasing absolute Ec number. (iv) Up to moderate supersonic Mach numbers, the Ec number and thus the relative portions of the terms of the total-enthalpy equation are influenced by both the wall-to-recovery temperature ratio T w /T r and the local Mach number. For hypersonic Mach numbers, the influence of the Mach number becomes negligible. (v) For pressure-gradient cases, the qualitative trends are not appropriately predicted by the SRA, implying that a modification for pressure-gradient influences needs to be constructed.
In summary, the integral analysis conducted allows us to represent the relevant transfer mechanisms in a simplified and compact form, clearly visualizing superordinate mechanisms and making them intuitively tangible. Besides providing hints for improved turbulence modelling, it forms a schematic that allows us to easily determine the relative importance of the terms of the governing equations, if the local Eckert number Ec = (γ − 1)M 2 e T e /(T r − T w ) and the pressure-gradient parameter β K = (δ * K /τ w )(dp e /dx) are known.

B.2. Remarks on the number of repeated integrations
Exemplarily evaluated for n = 2, 3, 4 and 10 for the subsonic iZPG +10K M=0.5 case, the streamwise evolution of the integral identity (A2) is evaluated in figure 11 for c f as a function of Re τ . According to Cauchy's formula for repeated integration in (3.3), variations in n essentially only affect the exponent of the weighting terms ( y b − y) n−1 . Thus, the greater the number of repeated integrations is, the more the near-wall area of the individual terms is weighted. Consequently, the proportions of the c T f and the c L f terms become increasingly dominant for larger n, both of which have their maximum in the near-wall region.

B.3. Remarks on the upper integration limit
To estimate the influence of the upper integration limit on the resulting contributions, the integral identity for c f (3.7) is evaluated for the iZPG +10K M=0.5 case with y d = 0.2δ 99 , 1.0δ 99 , 1.3δ 99 and 2.0δ 99 . Results are depicted in figure 12. As can be seen by comparison of figure 12(b,c), the choice of y d = 1.0δ 99 or of approximately 1.3δ 99 does not influence the qualitative trends. 300 400 500