Applied Mathematics and Mechanics (English Edition) Non-equilibrium turbulent phenomena in the flow over a backward-facing ramp

Transcription

1 Appl. Math. Mech. -Engl. Ed., 40(2), (2019) Applied Mathematics and Mechanics (English Edition) Non-equilibrium turbulent phenomena in the flow over a backward-facing ramp Le FANG 1, Hongkai ZHAO 2, Weidan NI 2,3, Jian FANG 3, Lipeng LU 2, 1. Laboratoire de Mecanique Physique (LMP), Ecole Centrale de Pékin, Beihang University, Beijing , China; 2. National Key Laboratory of Science and Technology on Aero-Engine Aero-Thermodynamics, School of Energy and Power Engineering, Beihang University, Beijing , China; 3. Scientific Computing Department, Science and Technology Facilities Council (STFC), Daresbury Laboratory, Warrington WA4 4AD, U. K. (Received Sept. 5, 2018 / Revised Oct. 25, 2018) Abstract Non-equilibrium turbulence phenomena have raised great interests in recent years. Significant efforts have been devoted to non-equilibrium turbulence properties in canonical flows, e.g., grid turbulence, turbulent wakes, and homogeneous isotropic turbulence (HIT). The non-equilibrium turbulence in non-canonical flows, however, has rarely been studied due to the complexity of the flows. In the present contribution, a directnumerical simulation (DNS) database of a turbulent flow is analyzed over a backwardfacing ramp, the flow near the boundary is demonstrated, and the non-equilibrium turbulent properties of the flow in the wake of the ramp are presented by using the characteristic parameters such as the dissipation coefficient C ɛ and the skewness of longitudinal velocity gradient S k, but with opposite underlying turbulent energy transfer properties. The equation of Lagrangian velocity gradient correlation is examined, and the results show that non-equilibrium turbulence is the result of phase de-coherence phenomena, which is not taken into account in the modeling of non-equilibrium turbulence. These findings are expected to inspire deeper investigation of different non-equilibrium turbulence phenomena in different flow conditions and the improvement of turbulence modeling. Key words non-equilibrium turbulence, backward ramp, direct numerical simulation (DNS) Chinese Library Classification O Mathematics Subject Classification 76F02 Citation: FANG, L., ZHAO, H. K., NI, W. D., FANG, J., and LU, L. P. Non-equilibrium turbulent phenomena in the flow over a backward-facing ramp. Applied Mathematics and Mechanics (English Edition), 40(2), (2019) Corresponding author, lulp@buaa.edu.cn Project supported by the National Natural Science Foundation of China (Nos , , and ), the National Basic Research Program of China (No. 2014CB046405), and the U. K. Engineering and Physical Sciences Research Council (EPSRC) (Nos. EP/K024574/1 and EP/L000261/1) c Shanghai University and Springer-Verlag GmbH Germany, part of Springer Nature 2019

2 216 Le FANG, Hongkai ZHAO, Weidan NI, Jian FANG, and Lipeng LU 1 Introduction As the foundation of turbulence research, the Richardson-Kolmogorov theory [1 4] leads to two conclusions in the inertial range of high-reynolds number turbulence, i.e., the scaling law of k 5/3 where k is the wavenumber for the energy distribution at an instant and the equilibrium energy transfer process which leads to the equilibrium dissipation scaling ɛ = C ɛ U 3 /L, (1) in which the dissipation coefficient C ɛ is constant, ɛ is the dissipation rate, U 2 is the turbulent kinetic energy, and L is the integral length scale. The same resultant dissipation scaling has also been suggested by simplifying the generic Kármán-Howarth equation (without assumptions of homogeneity and isotropy) [5 8]. In fact, in some situations, the assumptions involved in the above approaches are not satisfied, and this might result in non-constant C ɛ, which has been considered as a consequence of the non-equilibrium turbulent kinetic energy transfer process. This anomalous dissipation scaling has been found in many experiments and numerical simulations for fractal grid turbulence [9 13], turbulent wakes [14 16], homogeneous isotropic turbulence (HIT) with unsteady forcing [17 18], transitional channel flows [19], etc. In general, the non-equilibrium dissipation scaling can be described as C ɛ Re m I /Ren L or C ɛ Re m/2 I /Re n λ with m 1 n, wherere λ is the local Reynolds number defined as Re λ = Uλ/ν in which λ is the Taylor scale defined as λ 2 = νu 2 /ɛ, ν is the viscosity, Re L is the integral scale Reynolds number defined as Re L = UL/ν, andre I the inlet Reynolds number. In recent years, the skewness of longitudinal velocity derivative [20 21] defined by S k = ( u x ) 3 ( u x ) 2 3 2, (2) where u is the fluctuating velocity component in the x-direction and is the ensemble average, is suggested to indicate the similar non-equilibrium phenomena as C ɛ. In the experiments of Hearst and Lavoie [22], the consistent revolution tendency between C ɛ and S k along the downstream of the grid was observed. This definition has also been extended to the resolvedscale skewness to enable the investigation of multi-scale non-equilibrium phenomena [23]. Though the two quantities described above, i.e., C ɛ and S k, can describe non-equilibrium turbulence phenomena, they are limited in simple geometries and are not applicable for complex turbulent flows. Indeed, the calculation of C ɛ needs to calculate the integral scale involving spacial integration, which by definition is only valid in homogeneous turbulence. Meanwhile, S k is defined by only the local variables, and can be used in inhomogeneous turbulence, which cannot be simply extended to anisotropic turbulence. In order to extend these definitions to complex flows, Fang et al. [24] used the equation of Lagrangian velocity gradient correlation to estimate the non-equilibrium phenomena. Specifically, the equation is Φ ijij := d A ij A ij dt = Q ijij +Π ijij + V ijij, (3) where the superscript denotes fluctuation, A ij = ui x j is the fluctuating velocity gradient tensor, d/dt represents the Lagrangian derivative [24 25],and Q ijij = 2 A ij A kj A ik, Π ijij = 2 ρ V A ij 2 p, x i x j ijij νa =2 2 A ij ij x k x k (4)

3 Non-equilibrium turbulent phenomena in the flow over a backward-facing ramp 217 are the production, pressure, and viscosity terms, respectively. This equation can also be extended to the large-eddy simulation (LES) databases containing only the information of the resolved scale, by adding a subgrid-scale term according to the expression of the subgrid stress [24]. In equilibrium HIT, the constant value of S k implies that Φ ijij = 0, illustrating that Φ ijij can be considered as a new quantity for estimating non-equilibrium turbulence[24]. In Subsection 3.4 of Ref. [24], the consistency was discussed among the criteria of C ɛ, S k,and Φ ijij in an interior flow in a compressor, but the results were still quite preliminary. In order to study the non-equilibrium turbulence under complex geometries and to find more different mechanisms of non-equilibrium turbulence, in the present contribution, we focus on a direct-numerical simulation (DNS) database of a backward-facing ramp to resolve the whole range of energetic eddies of turbulence. Such numerical experiments have been widely adopted in the studies in turbulence control [26 30] and turbulence physics [31 34]. In this study, the non-equilibrium phenomena raised by two typical flow conditions, i.e., the boundary layer flow and the wake with a separation vortex downstream the ramp, are investigated, respectively, showing different underlying mechanisms. In Section 2, we will introduce the basic settings of the DNS database and its validation. In Section 3, we will show the consistency among the three indicators of non-equilibrium turbulence mentioned above, explain the mechanisms of energy and transfer spectra evolutions, and quantitatively show the non-equilibrium turbulent phenomena at different filtering scales. Concluding remarks are given in Section 4 finally. 2 Descriptions of backward-facing ramp and DNS settings In this section, the settings of DNS simulation are described in Subsection 2.1 in three aspects: the governing equations, the numerical method, and the computation setup, respectively. The simulation validations will be found in Subsection 2.2. The averaged flow condition will be described in Subsection Methodology Governing equations The three-dimensional (3-D) unsteady compressible Navier-Stokes (N-S) equations in a general time-invariant system are solved numerically. This set of N-S equations can be written in a strong conservation non-dimensional form as follows: Q t + E i x i = F i x i, (5) where Q =(ρ, ρu, ρv, ρw, E) T is the solution vector, in which ρ is the density, u, v, andw are components of the velocity, and E is the total energy. The static temperature T and the pressure P are related to the density via the equation of state of the ideal gas as P = ρrt. The convection terms E i and the diffusion terms F i in Eq. (5) are, respectively, expressed as E i =(ρu i,ρu 1 u i + δ 1i P, ρu 2 u i + δ 2i P, ρu 3 u i + δ 3i P, (E + P )u i ) T, (6) F i =(0,τ 1i,τ 2i,τ 3i,b i ) T with the Einstein summation convention. The indices i =1, 2, 3 correspond to the streamwise, wall-normal, and spanwise directions, respectively. The notation x i =(x 1,x 2,x 3 ) is adopted to represent (x, y, z), and u i =(u 1,u 2,u 3 )=(u, v, w) corresponds to the streamwise, wall-normal, and spanwise velocity components, respectively. δ ij is the standard Kronecker delta. The total energy E is expressed as follows: E = 1 2 (ρu iu i )+ ρt γ(γ 1)M 2. (7)

4 218 Le FANG, Hongkai ZHAO, Weidan NI, Jian FANG, and Lipeng LU The stress tensor and the heat flux vector are expressed as follows: τ ij = μ Re b i = u j τ ij + ( ui + u j 2 x j x i 3 δ ij u ) k, (8) x k μ PrRe(γ 1)M 2 T x i. (9) The dynamic viscosity coefficient μ is calculated via the Sutherland law as follows: μ = T 1.5 T S/T ref +1 T + T S /T ref, (10) where T ref = K, and T S = 110.4K. The N-S equations are non-dimensionalized with the reference density ρ ref, the reference velocity u ref, the temperature at wall T ref, and the dynamic viscosity μ ref as well as the height of the ramp H. The resulting dimensionless parameters are the Reynolds number Re = ρ ref u ref H/μ ref and the Mach number M = u ref / γrt ref. A constant Prandtl number Pr = μc p /k =0.72 is used, where c p = γr/(γ 1) is the specific heat capacity of gas at constant pressure, and k is the thermal conductivity. R and γ are the gas constant and the specific heat capacity ratio, which are set to be R = 287.1J kg 1 K 1 and γ =1.4, respectively Numerical method The in-house code ASTR, which has been applied in DNS [35 36] and LES [37] of wall bounded turbulent flows, is used to solve Eq. (5). The N-S equations are projected to the computational coordinate system and solved with the classic sixth-order compact central scheme [38].Toremove the small-scale wiggles due to the aliasing errors resulting from the discrete evaluation of the nonlinear convection terms, the 10th-order compact filter is incorporated, which limits the filter being implemented only at high wave numbers [39]. After all the spatial terms are solved, the 3rd-order three-step total variation diminishing Runge-Kutta method [40] is used for the time integration Computation setup The considered geometry corresponds to the configuration of Lardeau and Leschziner [41] and Bentaleb et al. [42] (see Fig. 1). The Reynolds number, based on the height of the ramp H and the incoming free-stream velocity, is Re = The corresponding friction Reynolds number is Re τ = 755, which is based on the wall friction velocity u τ at the inlet plane x = The freestream Mach number is M =0.2, defining a weakly compressible flow. The size of the computational domain L x L y L z ranges from to55 9 5, which are discretized with a mesh of points. Fig. 1 Sketch of the computation domain and reduced numerical meshes, where the meshes are plotted every 5th grid lines in both the x- andy-directions, respectively As shown in Fig. 2, the mesh is refined above the ramp along the streamwise direction, stretched towards the walls in the wall-normal direction, and evenly distributed in the spanwise

5 Non-equilibrium turbulent phenomena in the flow over a backward-facing ramp 219 direction. A non-slip isothermal boundary condition is applied at the wall. The digital filter method proposed by Touber and Sandham [43] is used to generate the synthetic inflow turbulence, and a transitional region of 12δ 0 is incorporated to let synthetic fluctuations evolve into fully developed turbulence. The generated artificial turbulent fluctuations are super-imposed onto the turbulent mean velocity and temperature profiles. The inflow condition is then adopted to prescribe the flow variables, and the pressure is extrapolated from the inner grid points. At the upper and outlet planes, the generalized non-reflecting boundary conditions [44 45] are used. The periodic boundary condition is prescribed in the spanwise direction. Due to the use of periodic boundary condition, it is important to examine the influence of the computational domain size on the simulation results [31]. There are different criteria for choosing the computational domain size. The most strict criterion is to contain all energetic eddies at different scales [32,46], which usually requires a very large computational domain size and requires a large amount of computer sources. This criterion is usually used when there exist very large scale coherent structures in the flow [47 51]. A more practical criterion is to ensure that the statistical results are insensitive to the change in the computational domain size. This criterion is used in the present study. It is found that reducing the spanwise computational domain size by 1/3 results in the change of less than 5% in the mean velocity. Based on this test, it is believed that the results presented in this paper are not affected by the periodic boundary condition. The mesh settings at x = 10.0, which is located in the fully developed turbulent boundary layer upstream of the ramp, for the present simulation are as follows: Re τ = 650, Δx = , Δx + =13.2, Δy 1 = , Δy 1 + =0.662, Δy δ = , Δy + δ =9.05, Δz = , Δz+ =5.16, where Δy δ refers to the mesh spacing in the wall-normal direction at y = δ, andδ is the local boundary layer thickness at x = Fig. 2 Mesh resolution at three representative streamwise locations, where Δ is the effective mesh spacing, and η is the Kolmogorov scale The mesh resolution scaled with the local Kolmogorov length scale [52] is checked (see Fig. 2). Δ = 3 ΔxΔyΔz is the effective mesh spacing, where Δx, Δy, and Δz refertothelocal grid spacing in the streamwise, wall-normal, and spanwise directions, respectively. η is the Kolmogorov scale calculated by ( μ z ) 1 4 η =, ρ z ε z where ε is the local dissipation rate of the turbulence kinetic energy [53]. The superscript stands for the time-averaged operator, and z stands for the spanwise-averaged operator. The

6 220 Le FANG, Hongkai ZHAO, Weidan NI, Jian FANG, and Lipeng LU ratio of the effective mesh spacing to the local Kolmogorov length scale is less than 4.2 throughout the computational zone, indicating that the mesh matches the resolution requirement of DNS. After the transient period, the flow becomes fully turbulent and reaches a statistically steady state after about 133 time units. The data samples are collected for about 182 time units for the post-process and analysis. 2.2 Validation The mean velocity profile in the equilibrium zone (x = 10.0) is compared with the classic law of wall and the incompressible DNS data at Re τ in Ref. [54] in Fig. 3. An agreement in both the linear sub-layer and log-law layer is seen to be close and the difference in the wake layer is attributed to the Reynolds number effects [36]. The root mean square (RMS) velocity fluctuations u i,rms u = i, z u i, z (i =1, 2, 3) and Reynolds shear stress at the same streamwise position are compared with the DNS data in Refs. [54] and [55] in Fig. 4. Good agreements for RMS velocity fluctuations and Reynolds shear stress are obtained in the nearwall region. Fig. 3 Mean velocity profile in inner scaling, where the von Kármán constant is κ =0.41, the data from Ref. [54] are at Re θ = 1 000, and the present data are at x = x start and Re θ = ' ' ' ' ' ' ' ' ' ' Fig. 4 RMS velocity fluctuations and Reynolds shear stress in the inner and outer scaling, where the data from Ref. [54] are at Re θ = 1 000, the data from Ref. [55] are at Re θ = 1 100, and the present data are at x = x start and Re θ = 1 120

7 Non-equilibrium turbulent phenomena in the flow over a backward-facing ramp Averaged flow condition As shown in Fig. 5, in front of the ramp, a steady flat plate boundary flow condition is established with a thickness of δ/h 1. In the wake of the ramp, the flow separates around x/h 1 and reattaches around x/h 5 within which a recirculating vortex appears. Meanwhile, streamlines are bent toward the lower wall to refill the area after the expansion after the ramp. This expansion causes a larger low speed region and a lower main flow velocity. Fig. 5 Averaged streamlines with contours of u 1 (color online) 3 Analysis of non-equilibrium turbulence In this section, we first demonstrate the consistency among C ɛ, S k,andφ ijij by the analysis of their distributions in space, along a streamline and their statistical correlations, respectively. Then, we show the distributions of energy and transfer at different scales by presenting energy spectra and transfer spectra at different points. At the end of this section, we quantitatively show the non-equilibrium phenomena at different filtering scales by applying filter operations to the DNS database. 3.1 Consistency among C ɛ, S k,andφ ijij It is well-known that under homogeneity and isotropy conditions, C ɛ, S k,andφ ijij are consistent with each other in the description of non-equilibrium phenomena. Here, we provide evidence by comparing the distributions in space, along a streamline and the statistical correlation, respectively. In the present sample, as the flow is homogeneous in the spanwise direction and in time, the ensemble average involved in the definitions of all the three quantities represents the average over the spanwise direction and over time. The integral scale L in Eq. (1) is defined as L w = R ww (r)dr, 0 where R ww (r) = w(z)w(z + r), along the spanwise direction, and U in Eq. (1) is adopted as the RMS of the spanwise fluctuating velocity w. The streamwise components are used to calculate the values of S k. The values of Φ ijij are normalized by 1/ τ K 3 s in which τ K is the Kolmogorov time scale defined by ( ν ) 1/2, τ K = (11) ɛ where s denotes a mass-flow-rate-weighted average operator [24]. We select a streamwise crosssection x/h = 10 to represent the incoming flow with a fully developed flat plate boundary and without being disturbed by the ramp to perform the normalization process Consistency in space As shown in Fig. 6, in the regions far from the wall and the ramp, the flow is under equilibrium condition with a typical value of C ɛ 0.6. This is in agreement with the traditional results in HIT [56]. The non-equilibrium turbulent phenomena concentrate in the boundary, the recirculation zone, and a long distance downstream of the ramp, where C ɛ departs from the equilibrium value. In the boundary layer upstream of the ramp, non-equilibrium turbulent phenomena stem from the wall and spread to about 0.2 times the boundary layer thickness, which

8 222 Le FANG, Hongkai ZHAO, Weidan NI, Jian FANG, and Lipeng LU corresponds to the thickness of the inner layer of a flat plate boundary layer, with a relatively higher value of C ɛ than the equilibrium region. In the wake behind the ramp, a lower value of C ɛ is presented in the center of the recirculation vortex, where C ɛ shows a similar value within the equilibrium region. This can be explained as that the flow is slow, and the non-equilibrium turbulence has enough time to reestablish an equilibrium state. This reestablishment is a slow process as argued in our previous study [24]. Therefore, the results in a long distance of nonequilibrium turbulence phenomena are downstream the ramp till more than twenty times the ramp height. Fig. 6 Distribution of C ɛ (color online) According to traditional values in HIT, the typical value of S k in the equilibrium region should be around 0.45 [23]. In Fig. 7, it is shown that, in the main flow above the beginning point of the ramp, S k shows a higher value than that in equilibrium regions. We remark that, this is not the result of the non-equilibrium energy transfer process, since, in this area, the distributions of both C ɛ and Φ ijij indicate an equilibrium turbulence state (see Figs. 6 and 8(a)). Indeed, this phenomenon should be related to the curvature of streamlines. Since the isotropy condition is not satisfied in this area, the traditional values of S k in HIT are not satisfied. Similar to C ɛ, in the boundary layer upstream of the ramp, the flow shows non-equilibrium transfer properties with a smaller value of S k compared with that in the equilibrium region. Also, the effects of the non-equilibrium region continue till about 0.2 times the boundary layer thickness. Downstream the ramp, S k shows a higher value, and the non-equilibrium region is larger than that shown by C ɛ, which should be the combined effects of the non-equilibrium phenomenon and the streamline curvature. The non-equilibrium region lasts for a long distance similar to C ɛ. Fig. 7 Distribution of S k (color online) As shown in Fig. 8(c), Π ijij is much smaller than Q ijij and V ijij for about one decade. It only affects a small region at the beginning of the ramp with a diffusion effect, which is not obvious downstream the ramp. Φ ijij is mainly the balance between the generating Q ijij and the dissipating V ijij. In the boundary layer upstream the ramp, both Q ijij and V ijij affect till about 0.2 times the boundary layer thickness. By contrast, their resultant balance Φ ijij is only observed within 0.05 times the boundary layer, which corresponds to y + 15 and is located in the buffer layer. The results of Φ ijij are negative, indicating the decrease in the dissipation rate when following the fluid particles. Downstream the ramp, the influence regions of Q ijij and V ijij are smaller than those of C ɛ and S k. The influence region of their balance Φ ijij is even much smaller. According to the distributions of the three indicators shown above, the three indicators,

9 Non-equilibrium turbulent phenomena in the flow over a backward-facing ramp 223 ' ' ' ' Fig. 8 Distributions of Φ ijij and its budget in the equation of Lagrangian velocity gradient correlation (color online) i.e., C ɛ, S k,andφ ijij, show qualitative agreement, and thus can be used to estimate the nonequilibrium phenomena Consistency along a streamline The definition of Φ ijij is the Lagrangian derivatives of the pseudodissipation, which gives another way to think about the question. The changes of Φ ijij by following a fluid particle may reveal the non-equilibrium process it experienced. The advection term is originally expanded as follows: A ij U A ij k = U k A ij A ij + u A ij A ij k. (12) x k x k x k The two terms on the right-hand side of Eq. (12) are the advection effects induced by the average velocity and fluctuating velocity, respectively. We denote the latter one as T ijij to represent the turbulent advection of the pseudodissipation when following a fluid particle along a statistical streamline. The distribution is shown in Fig. 9. It is only evident in a small region at the beginning of the ramp, and acts as a diffusion effect similar to the pressure term. By contrast, in other regions, it is relatively small, which indicates the fact that the turbulent advection is not essential in the analysis of non-equilibrium phenomena and the investigation along the streamlines can be considered to use the statistical streamlines in the regions away from the beginning of the ramp. Though the turbulent advection cannot be neglected in some regions, we can still approximately analyze along a statistical streamline with the understanding that the values of the indicators for non-equilibrium phenomena at one single point in these regions, where the fluid particles pass through, represent only the statistical results at this point

10 224 Le FANG, Hongkai ZHAO, Weidan NI, Jian FANG, and Lipeng LU but not the particles experienced along a streamline, since the turbulent advection effect is too obvious to be ignored and the statistical streamline cannot be considered as the streamline along which the particles go through. Fig. 9 Distribution of T ijij (color online) We select a streamline passing through a point inside the inner layer upstream the ramp whose coordinate is x/h = 10, and y/h =1.1 (see Fig. 10). The variations of C ɛ, S k,and Φ ijij and the variations of Φ ijij and its components versus x/h are shown in Fig. 11. As shown in Fig. 11(a), the three quantities show good consistency with each other. In the inner layer upstream the ramp (x/h 0), the flow is under the non-equilibrium state, but is steady with a high level of C ɛ and a low level of S k ; similar results can also be found in Ref. [19]. This nonequilibrium phenomenon should be caused by the influence of the wall. When the flow passes the ramp (0 x/h 3), all the three quantities fluctuate greatly, indicating a severe nonequilibrium transfer phenomenon in this region. Note that in this region, the results along the streamline are not experienced by the fluid particles on this statistical streamline, but are the average results on each single point as explained before. Around about 20 times the ramp height downstream the ramp, the flow returns to the equilibrium state gradually. Though the transfer from the non-equilibrium state to the equilibrium state takes a long distance downstream the disturbance, the flow will always reestablish an equilibrium transfer state, unless it is disturbed all the time. We stress here that the values of the three quantities deviate to different sides of the equilibrium value in the boundary layer and the wake of the ramp, which reveals the different non-equilibrium phenomena in the boundary layer and the wake. Fig. 10 Characteristic points and selected streamlines with contours of C ɛ (color online) As shown in Fig. 11(b), the diffusing Π ijij has only a little effect at the separation point, and is one order smaller than Q ijij and V ijij. Φ ijij is mainly the balance between the generating Q ijij and the dissipating V ijij Consistency by statistical analysis As first found by Hearst and Lavoie [22], C ɛ and S k are consistent with each other in the wake of fractal grids under the non-equilibrium state. We calculate the joint probability distribution function (joint-pdf) between C ɛ and S k in three typical regions, i.e., the equilibrium region as reference (set as 20 <x/h< 10, 1.3 <y/h<1.7), the wake region downstream the ramp (set as 0 <x/h<20, y/h <1.2), and the inner layer region (set as 20 <x/h< 10, y/h<1.3). As shown in Fig. 12, C ɛ and S k obviously correlate. In the wake, C ɛ and the absolute value of S k are lower than those in the equilibrium region. The results in the wake and in the equilibrium region show good agreement with the results of Hearst and Lavoie [22].

11 Non-equilibrium turbulent phenomena in the flow over a backward-facing ramp 225 ' ijij 'ijij 'ijij ' ijij 'ijij 'ijij ' ijij ' ijij ' ijij ' ijij Fig. 11 Variations of C ɛ, S k,andφ ijij and the variations of Φ ijij and its components along the streamline: (a) variations of C ɛ, S k,andφ ijij versus the streamwise distance normalized by the ramp height; (b) variations of Φ ijij and its components versus the streamwise distance normalized by the ramp height By contrast, in the inner layer, C ɛ and the absolute value of S k are higher than those in the equilibrium region, which seems to be a new phenomenon which, to our knowledge, has not been reported in the literature. This high value appears in the region near the wall, showing great similarity to the result in the named equilibrium region of a transitional channel flow obtained in Ref. [19]. The distribution on the phase space may indicate a different non-equilibrium energy transfer mechanism. Indeed, in the boundary layer, the turbulence structures grow from the wall and form large structures like the streaks under non-linear interaction, which is a procedure of generating large scales from small scales. As a result, the portion of the energy transfer from large scales to small scales should be lower than that in the equilibrium region. However, in the wake, the process is on the contrary, and the large-scale flow structures break down to small vortices, indicating a process of generation of small scales from large scales. Fig. 12 Distributions of the joint-pdf between C ɛ and S k in three typical regions (color online) From Fig. 13, similar consistency between C ɛ and Φ ijij can be found. In the inner layer region, high values of C ɛ correspond to negative values of Φ ijij. In the wake, the phenomenon is on the contrary. We calculate the joint-pdf between C ɛ and Re 1 λ in three typical regions, respectively. The results are shown in Fig. 14. As shown in Fig. 14(b), inside the equilibrium region, the values of C ɛ and Re 1 λ remain nearly unchanged. For the third region, i.e., the inner layer shown in Fig. 14(c), the value of C ɛ varies in a large range from a typical value under the

12 226 Le FANG, Hongkai ZHAO, Weidan NI, Jian FANG, and Lipeng LU 'ijij ' ijij ' ijij Fig. 13 Distributions of the joint-pdf between C ɛ and Φ ijij in three typical regions (color online) equilibrium state to higher values. We observe an obvious proportional correlation between C ɛ and Re 1 λ. However, in the wake, the present calculation domain is not long enough for the flow to recover a steady equilibrium state with constant C ɛ and S k. Therefore, the results in Fig. 14(a) only concentrate in the lower left quarter without evolution toward equilibrium typical values. However, a similar proportional correlation should be expected according to the results found by Seoud and Vassilicos [9]. In short, C ɛ is inversely proportional to Re λ in the wake and inner layer. According to the Richardson-Kolmogorov theory, in turbulence with equilibrium energy transfer, a higher Reynolds number implies a wider range of scales that the turbulent energy has to travel through from large to small scales, which results in a larger ratio of L to λ. Also, from the definitions of C ɛ and λ, one can write immediately L λ C ɛre λ. (13) We can get the same conclusion from Eq. (13) that if the flow is under the equilibrium state with a constant C ɛ, a higher Reynolds number results in a larger ratio of L to λ [56].Ascanbe seen from Fig. 13, both in the boundary layer and the wake, the local Reynolds number Re λ is not constant in the whole region. From Eq. (13), we can see that inconstant Re λ should lead to an inconstant ratio of L to λ under the equilibrium energy transfer condition. On the contrary, according to Eq. (13), the inversely proportional correlation between C ɛ and Re 1 λ indicates a constant ratio of L to λ. This conflict reveals a turbulent energy transfer phenomenon different from the equilibrium state described by the Richardson-Kolmogorov theory, and should imply a non-equilibrium energy transfer phenomenon. λ λ λ Fig. 14 Distributions of the joint-pdf between C ɛ and Re 1 λ in three typical regions (color online)

13 Non-equilibrium turbulent phenomena in the flow over a backward-facing ramp Distributions of energy and transfer at different scales In order to investigate more about the underlying mechanism of the non-equilibrium transfer process, in this subsection, we examine the scale-dependent energy and transfer behaviors. The energy spectrum is defined as E(k) = û i (k)û i(k)da(k), (14) where k = k in which k is the wave number vector, ˆ denotes the spectral quantities, the superscript denotes the conjugate of a complex quantity, and da(k) is the area element of the sphere whose radius equals k in the spectrum space. Also, the transfer spectrum is defined as follows: ( T (k) = i ) k p (û p (k k )û i (k )û i (k) û p (k k )û i (k )û i (k)) da(k), (15) k where i is the imaginary unit. Here, we denote the convection term as follows: c i (x) =u p (x) u i(x) x p, (16) where x is the position vector in the physic space. Then, its spectral expression is ĉ i (k) =i k k p û p (k k )û i (k ). (17) Combining Eq. (17) with Eq. (15), we can obtain T (k) = ( ĉ i (k)û i(k) ĉ i (k)û i (k))da(k). (18) Because of the existence of convolution calculation in the original definition of transfer spectra, which is caused by the convection term, we use Eq. (18) instead of Eq. (15) to calculate the transfer spectra, in which the convection term is calculated in the physical space with a pseudo-spectral method. In the present case, homogeneity is only valid in the spanwise direction. Therefore, the results are presented as one-dimension spectra in the z-direction. E(k) is normalized by (ɛν 5 ) 1/4, T (k) is normalized by ɛ/k, and the wave number k is normalized by 1/η where η is the dissipation scale. We calculate the energy spectra and transfer spectra at six points, respectively. The wavenumber is normalized by the local dissipation scale, the energy spectra are normalized by the local dissipation rate and viscosity, while the transfer spectra are normalized by the wavenumber and local dissipation rate. Four of the points are selected downstream the ramp along the streamline discussed in Subsection and marked in Fig. 10. The other two points are selected inside the boundary layer and in the equilibrium region upstream the ramp. Point 1 is located near the separation point. Point 2 is the middle point between the non-equilibrium state of the boundary layer and the non-equilibrium state of the wake of the ramp reflected by C ɛ. Point 3 is around the middle region of the wake, while Point 4 is located near the reattachment point. We select (x/h = 10, y/h = 1.5) as the reference equilibrium point and (x/h = 10, y/h=1.1) as the inner layer point. The energy spectra at these points are shown in Fig. 15(a). All spectra overlap with each other in the inertial range, indicating that even for non-equilibrium turbulence, the energy distribution is similar to those of equilibrium cases. Within the region 0.1 <kη<0.2, the energy distribution approximately coincides with Kolmogorov s 3/5 law and can be considered as the inertial range. Meanwhile, this region is the same as the interval between the wave numbers

14 228 Le FANG, Hongkai ZHAO, Weidan NI, Jian FANG, and Lipeng LU of the integral scale and the Taylor scale (see Table 1), indicating that the energy distribution between the integral scale and the Taylor scale is coherent with the distribution of the inertial region described by Kolmogorov s law. As shown in Fig. 15(b), the transfer spectra at each point show a similar tendency that negative values are at small wave numbers while positive values are at large wave numbers, which indicates the fact that turbulent kinetic energy is transferred from large scales to small scales. However, the transfer spectra do not show exactly the similar distributions like energy spectra, which implies the different presences of the distribution of energy and its transfer under different non-equilibrium states. This can be explained since the non-equilibrium turbulence can stem from the phase de-coherence, which is sensitive to the different perturbations raised from complex geometries, leading to the non-injection facts from energy to transfer [19]. Table 1 The corresponding wave numbers of integral scale and Taylor scale at each characteristic point 1 Point L η 1 λ η Equilibrium point Boundary layer point Point Point Point Point υ Fig. 15 Spectra of the turbulent kinetic energy and transfer versus the wave number It is interesting to see the detailed evolutions of energy and transfer in the boundary layer and the wake. At the inner layer point, the energy is relatively lower in large scales as well as small scales than that at the equilibrium point. For the energy transfer, the starting point and the peak of energy output move toward smaller scales with a larger output peak value. Meanwhile, energy is also gained at smaller scales compared with that at the equilibrium point but with a similar input peak value. Because the physical process inside the boundary is the generation of turbulent flow from small structures to large structures with limitations on the size of the largest structures, it is natural that energy is lower in large scales while the energy output moves to smaller scales. Also, the portion of dissipated energy is relatively smaller than that under the equilibrium state for the generating process. This leads to a lower energy distribution in small scales. At Point 1, i.e., the separation point, the energy distribution remains nearly the same as that in the inner layer point. However, the transfer spectrum

15 Non-equilibrium turbulent phenomena in the flow over a backward-facing ramp 229 changes a lot. At large scales, it loses energy at a smaller scale than that in the boundary, while the output peak value is still low, which results in the remaining of energy in large scales. However, at small scales, the transfer spectrum presents a similar distribution at the boundary layer point, at which the energy transfer process is changed at large scales while remains nearly the same in small scales, indicating the scale effect of non-equilibrium disturbance. This means that the non-equilibrium disturbance firstly affects certain scales, which are usually the large scales, since the disturbance originates mostly from the complex geometry, e.g., the ramp, the backward step, the grid, and the cascade, whose scale is relatively larger than that of the turbulence. Then, the non-equilibrium effect will spread to smaller scales with the evolution of the turbulence. From Point 2 to Point 3, the energy at large scales increases gradually and the energy at small scales keeps nearly the same inside the wake. While from Point 3 to Point 4, the large scale energy rests unchanged, while increases at small scales. For the energy transfer from Point 2 to Point 4, the energy output and input move toward large scales. This is because of the gradual increase in the energy at large scales. In fact, the physical process inside the wake is a gradual breaking down of large scale structures formed inside the boundary. The breaking down process gives large scales more energy from the main flow. Then, it will be transferred toward small scales with the turbulence evolution. 3.3 Non-equilibrium turbulence at different filtering scales In order to analyze the non-equilibrium turbulence phenomena at different scales, we perform filter operations on the DNS database with different cutoff wave numbers (see Figs ). We remark that, in the spectral space, the filter operators are easy to be implemented, but in physical space, it is more difficult. Indeed, the sharp cutoff filters correspond to the convolution calculations, and they are implemented by applying a template operation to take the average of several layers of points around the selected one and traverse all the mesh points in the computational domain. With the filter operator, the evolution equation of Lagrangian velocity gradient correlation should be modified as follows: ΦF (n)ijij := d A < ij A < ij dt = QF (n)ijij +ΠF (n)ijij + VF (n)ijij +ΓF (n)ijij, (19) where the superscript < is the resolved scale, F (n) is the results obtained in the filtering process of n layers of points around the center point, and the subgrid scale term is ΓF (n)ijij =2 A < ij 2 τ ik < x k x j with the subgrid stress τ ij < = u < i u< j (u iu j ) <. After the filter operations, the value of Φ ijij increases greatly and the non-equilibrium turbulence regions in large scales are larger than those in small scales. With the increases in the filtering layers, in the wake region, the non-equilibrium turbulence region is separated into two parts gradually, one is located above the wake, and the other is located inside the wake and extends toward downstream. Inside the inner layer, Φ ijij changes from negative to positive, and affects the outer region of the boundary. The distributions of the production term and viscosity term remain similar to those under the non-filtered situation upstream the ramp and in the wake, but the absolute values decrease greatly with the increase in the filtering level. However, downstream the wake, the non-equilibrium turbulence region is further extended, which indicates a slower evolution speed of the non-equilibrium turbulence at large scales. Similarly, the absolute value of the pressure

16 230 Le FANG, Hongkai ZHAO, Weidan NI, Jian FANG, and Lipeng LU ' ' ' ' ' ' Fig. 16 Distributions of Φ ijij and its filtered results (color online) term also decreases. The affecting regions of pressure term extend toward the main flow above the wake. By contrast, Φ ijij is the balance between the production and the viscosity term before the filter operation. After the filter operation, it is dominated by the subgrid term, reflecting the fact that non-equilibrium turbulence phenomena in large scales originate from small scales. Downstream the ramp, near the boundary around 2 <x/h<3, there is a region under a relatively equilibrium state before the filter operation. However, with the increase in the filter level, this region diminishes gradually and becomes invisible in the end. As we have mentioned before, the non-equilibrium disturbance should be transferred from scale to scale. As a result, the reason why it shows an equilibrium state before the filter operation is that the

17 Non-equilibrium turbulent phenomena in the flow over a backward-facing ramp 231 ' ' ' ' ' ' Fig. 17 Distributions of Q ijij and its filtered results (color online) non-equilibrium disturbance at large scales have not reached the mesh scale. 4 Conclusions In this paper, we present the non-equilibrium turbulence phenomena in the boundary layer and the wake of a backward ramp flow. The contributions are as follows: (i) Based on the distributions of C ɛ, S k,andφ ijij in the physical space and their joint distributions, the consistency among them is verified in the wake and in the boundary layer. C ɛ is inversely proportional to the local Reynolds number, indicating the non-equilibrium nature

18 232 Le FANG, Hongkai ZHAO, Weidan NI, Jian FANG, and Lipeng LU ' ' ' ' ' ' Fig. 18 Distributions of Π ijij and its filtered results (color online) in these regions. It is further shown that the mechanisms of non-equilibrium transfer differ between the wake and the boundary layer. (ii) From the energy spectra and transfer spectra, it is shown that the energy transfer at each scale varies greatly, indicating that the underlying mechanism of non-equilibrium turbulence might be a phase de-coherence phenomenon. (iii) By analyzing the filtered results, it is shown that, at different scales, the non-equilibrium turbulence phenomena are quantitatively different, indicating the inter-scale spread of nonequilibrium disturbances.

19 Non-equilibrium turbulent phenomena in the flow over a backward-facing ramp 233 ' ' ' ' ' ' Fig. 19 Distributions of V ijij and its filtered results (color online) References [1] RICHARDSON, L. F. Weather Prediction by Numerical Process, Cambridge University Press, Cambridge (2007) [2] KOLMOGOROV, A. N. The local structure of turbulence in incompressible viscous fluid for very large Reynolds number. Proceedings: Mathematical and Physical Sciences, 30, 301 (1941) [3] KOLMOGOROV, A. N. On the degeneration of isotropic turbulence in an incompressible viscous fluid. Doklady Akademii Nauk SSSR, 31, (1941) [4] KOLMOGOROV, A. N. Dissipation of energy in locally isotropic turbulence. Proceedings: Mathematical and Physical Sciences, 434, (1991)

20 234 Le FANG, Hongkai ZHAO, Weidan NI, Jian FANG, and Lipeng LU ' ' ' ' ' Fig. 20 Distributions of the subgrid terms (color online) [5] HILL, R. Exact second-order structure function relationship. Journal of Fluid Mechanics, 468, (2002) [6] MARATI, N., CASCIOLA, C., and PIVA, R. Energy cascade and spatial fluxes in wall turbulence. Journal of Fluid Mechanics, 521, (2004) [7] DANAILA, L., KRAWCZYNSKI, J., THIESSET, F., and RENOU, B. Yaglom-like equation in axisymmetric anisotropic turbulence. Physica D: Nonlinear Phenomena, 241, (2012) [8] VALENTE, P. and VASSILICOS, J. The energy cascade in grid-generated non-equilibrium decaying turbulence. Physics of Fluids, 27, (2015) [9] SEOUD, R. and VASSILICOS, J. Dissipation and decay of fractal-generated turbulence. Physics of Fluids, 19, (2007) [10] MAZELLIER, N. and VASSILICOS, J. Turbulence without Richardson-Kolmogorov cascade. Physics of Fluids, 22, (2010) [11] NAGATA, K., SAKAI, Y., INABA, T., SUZUKI, H., TERASHIMA, O., and SUZUKI, H. Turbulence structure and turbulence kinetic energy transport in multiscale/fractal-generated turbulence. Physics of Fluids, 25, (2013)

21 Non-equilibrium turbulent phenomena in the flow over a backward-facing ramp 235 [12] GOMES-FERNANDES, R., GANAPATHISUBRAMANI, B., and VASSILICOS, J. Particle image velocimetry study of fractal-generated turbulence. Journal of Fluid Mechanics, 711, (2012) [13] VALENTE, P. and VASSILICOS, J. Universal dissipation scaling for nonequilibrium turbulence. Physical Review Letters, 108, (2012) [14] NEDIC, J., VASSILICOS, J., and GANAPATHISUBRAMANI, B. Axisymmetric turbulent wakes with new nonequilibrium similarity scalings. Physical Review Letters, 111, (2013) [15] DAIRAY, T., OBLIGADO, M., and VASSILICOS, J. Non-equilibrium scaling laws in axisymmetric turbulent wakes. Journal of Fluid Mechanics, 781, (2015) [16] OBLIGADO, M., DAIRAY, T., and VASSILICOS, J. Nonequilibrium scalings of turbulent wakes. Physical Review Fluids, 1, (2016) [17] VALENTE, P., ONISHI, R., and DA SILVA, C. Origin of the imbalance between energy cascade and dissipation in turbulence. Physical Review E, 90, (2014) [18] GOTO, S. and VASSILICOS, J. Energy dissipation and flux laws for unsteady turbulence. Physics Letters A, 379, (2015) [19] LIU, F., LU, L., and FANG, L. Non-equilibrium turbulent phenomena in transitional channel flows. Journal of Turbulence, 19, (2018) [20] ISAZA, J., SALAZAR, R., and WARHAFT, Z. On grid-generated turbulence in the near- and far-field regions. Journal of Fluid Mechanics, 753, (2014) [21] AYYALASOMAYAJULA, S. and WARHAFT, Z. Nonlinear interactions in strained axisymmetric high-reynolds-number turbulence. Journal of Fluid Mechanics, 566, (2006) [22] HEARST, R. and LAVOIE, P. Velocity derivative skewness in fractal-generated, non-equilibrium grid turbulence. Physics of Fluids, 27, (2015) [23] FANG, L., ZHU, Y., LIU, Y., and LU, L. Spectral non-equilibrium property in homogeneous isotropic turbulence and its implication in subgrid-scale modeling. Physics Letters A, 379, (2015) [24] FANG, L., ZHAO, H. K., LU, L. P., LIU, Y. W., and YAN, H. Quantitative description of nonequilibrium turbulent phenomena in compressors. Aerospace Science and Technology, 71, (2017) [25] FANG, L., BOS, W., and JIN, G. Short-time evolution of Lagrangian velocity gradient correlations in isotropic turbulence. Physics of Fluids, 27, (2015) [26] CHOI, H., MOIN, P., and KIM, J. Direct numerical simulation of turbulent flow over riblets. Journal of Fluid Mechanics, 255, (1993) [27] DENG, B. and XU, C. Influence of active control on STG-based generation of streamwise vortices in near-wall turbulence. Journal of Fluid Mechanics, 710, (2012) [28] DENG, B., XU, C., HUANG, W., and CUI, G. Strengthened opposition control for skin-friction reduction in wall-bounded turbulent flows. Journal of Turbulence, 15, (2014) [29] GE, M., XU, C., and CUI, G. Active control of turbulence for drag reduction based on the detection of near-wall streamwise vortices by wall information. Acta Mechanica Sinica, 31, (2015) [30] DENG, B., HUANG, W., and XU, C. Origin of effectiveness degradation in active drag reduction control of turbulent channel flow at Re λ = Journal of Turbulence, 17, (2016) [31] MOSER, R., KIM, J., and MANSOUR, N. Direct numerical simulation of turbulent channel flow up to Re τ = 590. Physics of Fluids, 11, (1999) [32] ALAMO, J. D. and JIMENEZ, J. Spectra of the very large anisotropic scales in turbulent channels. Physics of Fluids, 15, L41 L44 (2003) [33] GE, M., ZUO, Y., DENG, Y., and LI, Y. Spatial relation between fluctuating wall pressure and near-wall streamwise vortices in wall bounded turbulent flow. Applied Mathematics and Mechanics (English Edition), 36(6), (2015) [34] GE, M., TIAN, D., and LIU, Y. Dynamic evolution process of turbulent channel flow after opposition control. Fluid Dynamics Research, 49, (2017) [35] FANG, J., YAO, Y., LI, Z., and LU, L. Investigation of low-dissipation monotonicity-preserving scheme for direct numerical simulation of compressible turbulent flows. Computers and Fluids, 104, (2014)

22 236 Le FANG, Hongkai ZHAO, Weidan NI, Jian FANG, and Lipeng LU [36] FANG, J., YAO, Y., ZHELTOVODOV, A. A., LI, Z., and LU, L. Direct numerical simulation of supersonic turbulent flows around a tandem expansion-compression corner. Physics of Fluids, 27, (2015) [37] FANG, J., YAO, Y., ZHELTOVODOV, A. A., and LU, L. Investigation of three-dimensional shock wave/turbulent-boundary-layer interaction initiated by a single fin. AIAA Journal, 55, (2016) [38] LELE, S. Compact finite difference schemes with spectral-like resolution. Journal of Computational Physics, 103, (1992) [39] GAITONDE, D. V. and VISBAL, M. R. Pade-type higher-order boundary filters for the Navier- Stokes equations. AIAA Journal, 38, (2000) [40] GOTTLIEB, S. and SHU, C. W. Total variation diminishing Runge-Kutta schemes. Mathematics of Computation of the American Mathematical Society, 67, (1998) [41] LARDEAU, S. and LESCHZINER, M. The interaction of round synthetic jets with a turbulent boundary layer separating from a rounded ramp. Journal of Fluid Mechanics, 683, (2011) [42] BENTALEB, Y., LARDEAU, S., and LESCHZINER, M. A. Large-eddy simulation of turbulent boundary layer separation from a rounded step. Journal of Turbulence, 13, 1 4 (2012) [43] TOUBER, E. and SANDHAM, N. D. Large-eddy simulation of low-frequency unsteadiness in a turbulent shock-induced separation bubble. Theoretical and Computational Fluid Dynamics, 23, (2009) [44] KIM, J. W. and LEE, D. J. Generalized characteristic boundary conditions for computational aeroacoustics. AIAA Journal, 38, (2000) [45] KIM, J. W. and LEE, D. J. Generalized characteristic boundary conditions for computational aeroacoustics, part 2. AIAA Journal, 42, (2004) [46] YANG, Z. and WANG, B. C. Capturing Taylor-Görtler vortices in a streamwise-rotating channel at very high rotation numbers. Journal of Fluid Mechanics, 838, (2018) [47] YANG, Z., CUI, G., ZHANG, Z., and XU, C. A modified nonlinear sub-grid scale model for large eddy simulation with application to rotating turbulent channel flows. Physics of Fluids, 24, (2012) [48] YANG, Z., CUI, G., and ZHANG, Z. Large eddy simulation of rotating turbulent channel flow with a new dynamic global-coefficient nonlinear subgrid stress model. Journal of Turbulence, 13, N48 (2012) [49] AVSARKISOV, V., HOYAS, S., OBERLACK, M., and GARCIA-GALLACHE, J. Turbulent plane Couette flow at moderately high Reynolds number. Journal of Fluid Mechanics, 751, R1 (2014) [50] YANG, Z., DENG, B., WANG, B., and SHEN, L. Letter: the effects of streamwise system rotation on pressure fluctuations in a turbulent channel flow. Physics of Fluids, 30, (2018) [51] GANDIA-BARBERA, S., HOYAS, S., OBERLACK, M., and KRAHEBERGER, S. Letter: the link between the Reynolds shear stress and the large structures of turbulent Couette-Poiseuille flow. Physics of Fluids, 30, (2018) [52] YANG, Z., DENG, B. Q., and SHEN, L. Direct numerical simulation of wind turbulence over breaking waves. Journal of Fluid Mechanics, 850, (2018) [53] YANG, Z. and WANG, B. C. On the topology of the eigenframe of the subgrid-scale stress tensor. Journal of Fluid Mechanics, 798, (2016) [54] SCHLATTER, P. and ORLU, R. Assessment of direct numerical simulation data of turbulent boundary layers. Journal of Fluid Mechanics, 659, (2010) [55] JIMENEZ-SENDIN, J., HOYAS, S., SIMENS, M. P., and MIZUNO, Y. Turbulent boundary layers and channels at moderate Reynolds numbers. Journal of Fluid Mechanics, 657, (2010) [56] VASSILICOS, J. Dissipation in turbulent flows. Annual Review of Fluid Mechanics, 47, (2015)