This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

causal analysis of inner and outer motions in near-wall turbulent flow

Jingxuan Zhang Center for Particle-Laden Turbulence
Lanzhou University
Lanzhou, China
   Zhengping Zhu Zhejiang Laboratory
Hangzhou, China
   Ricardo Vinuesa FLOW, Engineering Mechanics
KTH Royal Institute of Technology
Stockholm, Sweden
Email: rvinuesa@mech.kth.se
   Ruifeng Hu Center for Particle-Laden Turbulence
Lanzhou University
Lanzhou, China
Email: hurf@lzu.edu.cn

ABSTRACT

In this work, we study the causality of near-wall inner and outer turbulent motions. Here we define the inner motions as the self-sustained near-wall cycle and the outer motions as those living in the logarithmic layer exhibiting a footprint on the near-wall region. We perform causal analysis using two different methods: one is the transfer entropy, based on the information theory, and the other one is the Liang–Kleeman information-flow theory. The causal-analysis methods are applied to several scenarios, including a linear and a non-linear problem, a low-dimensional model of the near-wall cycle of turbulence, as well as the interaction between inner and outer turbulent motions in a channel at a friction Reynolds number of Reτ=1000Re_{\tau}=1000. We find that both methods can well predict the causal links in the linear problem, and the information flow can identify more of the nonlinear problem. Despite richer causalities revealed by the transfer entropy for turbulent-flow problems, both methods can successfully identify the streak-vortex regeneration mechanism that majorly sustains the near-wall turbulence. It is also indicated that both bottom-up and top-down influences of inner and outer motions may coexist in addition to the multiscale self-sustaining mechanism. Lastly, we mention that the computation of the information flow is much more efficient than the transfer entropy. The present study suggests that the information flow can have great potential in causal inference for turbulent-flow problems besides the transfer entropy.

INTRODUCTION

It is well known that wall-bounded turbulent flows are populated with energy-containing coherent structures spanning a wide range of spatial and temporal scales (Robinson, 1991; Adrian, 2007; Smits et al., 2011; Jiménez, 2018), e.g. streamwise streaks, hairpin vortices, large-scale motions, very-large-scale motions and superstructures among others. An open question is what are the causal links of these turbulent motions, especially in the context of inner and outer motions, either bottom-up (Adrian et al., 2000), top-down (Hunt & Morrison, 2000), co-supporting (Toh & Itano, 2005; Zhou et al., 2022) or independently self-sustained at all scales (Cossu & Hwang, 2017). For this purpose, We resort to the two widely adopted metrics which have been applied in physics, atmospheric sciences as well as fluid dynamics research recently, i.e. the transfer entropy (Schreiber, 2000) and the Liang–Kleeman information flow (Liang & Kleeman, 2005; Liang, 2013, 2014, 2016). Here we define the inner motions as the universal near-wall cycle (Hamilton et al., 1995; Waleffe, 1997), and the outer motions as those living in the logarithmic layer exhibiting a footprint on the near-wall region (Hutchins & Marusic, 2007; Marusic et al., 2010).

METHODOLOGY

Inner-outer decomposition of near-wall turbulent motions

The turbulent-flow velocity components of inner and outer motions are decomposed by a scaling-improved inner-outer decomposition method Wang et al. (2021a), which is generally based on the predictive inner-outer model (PIOM) Marusic et al. (2010); Mathis et al. (2011); Baars et al. (2016). In the PIOM, the near-wall turbulent-flow velocities can be decomposed by ui=ui,S+ui,Lu_{i}=u_{i,S}+u_{i,L}, where uiu_{i} is the turbulent-flow velocity in the iith direction, ui,Su_{i,S} is the small-scale velocity of the inner motions or near-wall cycle modulated by large-scale outer footprint velocity ui,Lu_{i,L}.

The large-scale outer footprint velocity is calculated by:

ui,L(x,y,z,t)=Fx1{Hi,L(λx,y)Fx[ui(x,yO,z,t)]},u_{i,L}(x,y,z,t)=\\ F_{x}^{-1}\left\{H_{i,L}(\lambda_{x},y)F_{x}\left[u_{i}(x,y_{O},z,t)\right]\right\}, (1)

in which FxF_{x} and Fx1F_{x}^{-1} denote the Fourier transform and inverse Fourier transform, respectively, λx\lambda_{x} is the streamwise wavelength, (x,y,z)(x,y,z) are the coordinates in streamwise, wall-normal and spanwise directions, respectively, yO+=100y^{+}_{O}=100 is the outer reference wall-normal height (in inner units) for calculating near-wall footprint (Wang et al., 2021a), and Hi,LH_{i,L} is the scale-dependent complex-valued kernel function of the spectral linear stochastic estimation (Baars et al., 2016), which is:

Hi,L(λx,y)=u^i(λx,y,z,t)u^i(λx,yO,z,t)¯z,tu^i(λx,yO,z,t)u^i(λx,yO,z,t)¯z,t,H_{i,L}(\lambda_{x},y)=\frac{\langle\hat{u}_{i}(\lambda_{x},y,z,t)\overline{\hat{u}_{i}(\lambda_{x},y_{O},z,t)}\rangle_{z,t}}{\langle\hat{u}_{i}(\lambda_{x},y_{O},z,t)\overline{\hat{u}_{i}(\lambda_{x},y_{O},z,t)}\rangle_{z,t}}, (2)

where, ()^\hat{(\cdot)} represents variables in Fourier space, z,t\langle\cdot\rangle_{z,t} denotes averaging in spanwise and time.

Causal-analysis methods

Causal discovery or causal inference has recently become one of the most popular topics and a tool to discover the underlying causal structure of physical systems. There is an ever-growing number of methods designed to work under different assumptions since the seminal work of Granger (1969). There are several reviews dedicated to introducing different methods (Pearl, 2009; Liang, 2013; Camps-Valls et al., 2023; Runge et al., 2023). In this study, we choose two widely-used causal-inference methods, i.e. transfer entropy and information flow, and apply them to several scenarios including interactions of near-wall inner and outer turbulent motions.

Transfer entropy.

The framework of information theory (Shannon, 1948) can be employed to quantify causality among time signals of different variables. A key metric is the transfer entropy (Schreiber, 2000), which has been applied in turbulent-flow problems recently (Lozano-Durán et al., 2020; Wang et al., 2021b, 2022; Martínez-Sánchez et al., 2023).

Here we follow Martínez-Sánchez et al. (2023) to calculate the Shannon entropy S(X)S(X) of variable XX through estimations of the probability density function of each signal using the kk-nearest-neighbour entropy estimator (Kozachenko & Leonenko, 1987), as:

S(X)=ψ(N)ψ(k)+logcd+dNi=1Nlogϵ(i),{S}(X)=\psi(N)-\psi(k)+\log c_{d}+\frac{d}{N}\sum_{i=1}^{N}\log\epsilon(i), (3)

where XX is a discrete-valued variable, NN is the number of discrete data, ψ()\psi(\cdot) donates the digamma function, dd is the dimension of XX, cd=1c_{d}=1 for the LL_{\infty}-norm, and ϵ(i)\epsilon(i) is the distance from xix_{i} to its kkth neighbour.

For a multivariate problem, we use the following formula

Tji(Δt)=S(Vi(t)|𝐕(tΔt))S(Vi(t)|𝐕(tΔt)),T_{j\to i}(\Delta t)=S(V_{i}(t)|\mathbf{V}^{\not{j}}(t-\Delta t))-S(V_{i}(t)|\mathbf{V}(t-\Delta t)), (4)

to calculate the conditional transfer entropy (Lizier, 2014; Lozano-Durán et al., 2020), where Δt\Delta t is the time lag, S(Vi|𝐕)S(V_{i}|\mathbf{V}) is the conditional Shannon entropy and 𝐕\mathbf{V}^{\not{j}} is equivalent to 𝐕\mathbf{V} but excluding the component jj. The conditional Shannon entropy is defined by S(Y|X)=S(X,Y)S(X)S(Y|X)=S(X,Y)-S(X). Transfer entropy quantifies the amount of uncertainty in a future signal reduced by a past signal. A larger transfer entropy TjiT_{j\to i} represents a stronger causality from VjV_{j} to ViV_{i}.

Information flow.

Information flow refers to the transfer of information between two entities in a dynamical system. In information-flow theory, the quantity that quantifies causalities is the rate of information flow. Liang & Kleeman (2005) argued that, as a two-dimensional system ({x1,x2}\{x_{1},x_{2}\}) steers a state forward, the marginal entropy of x1x_{1} is replenished from two different sources: one is from x1x_{1} itself, and another one from x2x_{2}. The latter is through the very mechanism namely information flow. This gives a decomposition of the marginal entropy increase according to the underlying mechanisms: dS1/dt=dS1/dt+L21{\rm d}S_{1}/{\rm d}t={\rm d}S_{1}^{*}/{\rm d}t+L_{2\to 1}, the term dS1{\rm d}S_{1}^{*} is the x1x_{1} own contribution, the other term L21L_{2\to 1} is the rate of information flow from x2x_{2} to x1x_{1}. Liang (2013) firstly established a rigorous formalism of information flow for deterministic and stochastic systems. After that, Liang (2014) used maximum-likelihood estimation to derive an estimation formula for information flow consisting only of correlation. A preliminary application of information flow in turbulence can be found in Liang & Lozano-Durán (2016).

We evaluate multivariate information flow according to Liang (2021), as

Lji=dSidtdSidt1det(𝐂)m=1dΔjmCm,diCijCii,L_{j\to i}=\frac{{\rm d}S_{i}}{{\rm d}t}-\frac{{\rm d}S_{i\not j}}{{\rm d}t}\approx\frac{1}{{\rm det}(\mathbf{C})}\cdot\sum\limits_{m=1}^{d}\Delta_{jm}C_{m,di}\cdot\frac{C_{ij}}{C_{ii}}, (5)

where CijC_{ij} is the covariance between ViV_{i} and VjV_{j}, and Ci,djC_{i,dj} is the covariance between ViV_{i} and V˙j\dot{V}_{j}, V˙j\dot{V}_{j} is the time derivative of VjV_{j}, 𝐂\mathbf{C} is the covariance matrix, det(){\rm det}(\cdot) is the determinant of a matrix, and dd is the number of variables. The Δij\Delta_{ij} is the cofactors of the covariance matrix. dSi/dt{\rm d}S_{i}/dt is the temporal variation rate of S(Vi)S(V_{i}), and dSi/dt{\rm d}S_{i\not j}/dt is the evolution of S(𝐕i)S(\mathbf{V}_{i}) with the effect of VjV_{j} excluded. So, LjiL_{j\to i} denotes the rate of information flow from VjV_{j} to ViV_{i}.

In the following, the self-causalities are set to zero. The largest causality is used for normalization so that the causalities are between 0 and 1. The absolute value of the normalized metrics is preserved only if greater than 0.01, otherwise it is set to zero.

RESULTS AND DISCUSSION

Linear problem.

First, we study the causalities of a linear problem, which is a six-dimensional vector autoregressive process (Liang, 2021):

𝐗(n+1)=α+𝐀𝐗(n)+𝐁𝐞(n+1),\displaystyle{\mathbf{X}(n+1)=\mathbf{\alpha}+\mathbf{AX}(n)+\mathbf{Be}(n+1)}, (6)
𝐀=(000.60000.500000.800.700000000.70.400000.200.7000000.5),\displaystyle{\mathbf{A}=\left(\begin{array}[]{cccccc}0&0&-0.6&0&0&0\\ -0.5&0&0&0&0&0.8\\ 0&0.7&0&0&0&0\\ 0&0&0&0.7&0.4&0\\ 0&0&0&0.2&0&0.7\\ 0&0&0&0&0&-0.5\end{array}\right)},
α=(0.1,0.7,0.5,0.2,0.8,0.3)T,\displaystyle{\mathbf{\alpha}=(0.1,0.7,0.5,0.2,0.8,0.3)^{T}},

where 𝐁\mathbf{B} is a diagonal matrix with Bii=100B_{ii}=100 (i=1,,6)(i=1,...,6), and the errors eiN(0,1)e_{i}\sim N(0,1) are independent variables. We generated 500 independent sets of data consisting of six series with 10000 steps (randomly initialized). The time lag for causal analysis is set to Δt\Delta t.

The results are shown in figure 1. It shows that the two methods yield the same qualitative causalities in the linear problem, despite some numerical differences. The identified non-zero causalities of x1x2x_{1}\to x_{2}, x2x3x_{2}\to x_{3}, x3x1x_{3}\to x_{1}, x4x5x_{4}\to x_{5}, x5x4x_{5}\to x_{4}, x6x2x_{6}\to x_{2} and x6x5x_{6}\to x_{5} are consistent with equation (6). Therefore, both the transfer entropy and information flow are effective in evaluating the causalities of the linear problem.

Refer to caption
Refer to caption
Figure 1: Causal maps of the linear problem: (a) transfer entropy; (b) information flow. Red boxes: causalities identified by both transfer entropy and information flow.
Refer to caption
Refer to caption
Figure 2: Causal maps of the nonlinear problem: (a) transfer entropy; (b) information flow. Red boxes: causalities identified by both transfer entropy and information flow. Blue boxes: causalities identified by transfer entropy but information flow. Green boxes: causalities identified by information flow but transfer entropy.

Nonlinear problem.

The second case is the causalities of a nonlinear problem. We choose the coupled system investigated by Paluš et al. (2018):

{dx1/dt=ω1x2(t)x3(t),dx2/dt=ω1x1(t)+0.15x2(t),dx3/dt=0.2+x3(t)[x1(t)10],dy1/dt=ω2y2(t)y3(t)+ϵ[x1(t)y1(t)],dy2/dt=ω2y1(t)+0.15y2(t),dy3/dt=0.2+y3(t)[y1(t)10],\left\{\begin{aligned} &{\rm d}x_{1}/{\rm d}t=-\omega_{1}x_{2}(t)-x_{3}(t),\\ &{\rm d}x_{2}/{\rm d}t=\omega_{1}x_{1}(t)+0.15x_{2}(t),\\ &{\rm d}x_{3}/{\rm d}t=0.2+x_{3}(t)\left[x_{1}(t)-10\right],\\ &{\rm d}y_{1}/{\rm d}t=-\omega_{2}y_{2}(t)-y_{3}(t)+\epsilon\left[x_{1}(t)-y_{1}(t)\right],\\ &{\rm d}y_{2}/{\rm d}t=\omega_{2}y_{1}(t)+0.15y_{2}(t),\\ &{\rm d}y_{3}/{\rm d}t=0.2+y_{3}(t)\left[y_{1}(t)-10\right],\end{aligned}\right. (7)

in which ω1=1.015\omega_{1}=1.015 and ω2=0.985\omega_{2}=0.985, and ϵ\epsilon is set to 0.250.25. The above nonlinear system is solved using the fourth-order Runge-Kutta scheme with a time step Δt=0.001\Delta t=0.001. Initialized with random numbers, the equations are integrated forward for N=50000N=50000 steps, in which the first 1000010000 steps are discarded to eliminate the effects of initial conditions. A total of 500 independent sets of data are generated finally. The time lag for causal analysis is set to 2Δt2\Delta t.

The results are displayed in figure 2, It is seen that the causalities identified by the two metrics are not the same. Both of the two methods can identify the causalities of x2x1x_{2}\to x_{1}, x1x2x_{1}\to x_{2}, y2y1y_{2}\to y_{1}, x1y1x_{1}\to y_{1} and y1y2y_{1}\to y_{2}, which exist in equation (7). However, the causalities of x3x1x_{3}\to x_{1}, x1x3x_{1}\to x_{3}, y3y1y_{3}\to y_{1} and y1y3y_{1}\to y_{3} are missing in the results of transfer entropy, which are preserved in the results of information flow. Therefore, in this nonlinear problem, the information flow is more effective for identifying causalities than the transfer entropy. It is also noted that some identified causalities are not explicitly found in equation (7) which may be attributed to the nonlinear effects.

Refer to caption
Refer to caption
Figure 3: Causal maps of the low-dimensional model of self-sustaining near-wall turbulence: (a) transfer entropy; (b) information flow. Red boxes: causalities identified by both transfer entropy and information flow. Blue boxes: causalities identified by transfer entropy but information flow. Green boxes: causalities identified by information flow but transfer entropy.

Low-dimensional model of the near-wall cycle of turbulence.

Next, we analyze the causal relations present in a low-dimensional model of the plane Couette flow with a sinusoidal body force developed by Moehlis et al. (2004). This model presents an improvement of the eight-mode model of Waleffe (1997). The self-sustaining process of near-wall turbulence, i.e. the streak-vortex cycle, can be well captured by the low-dimensional model. There are nine modes of the model, which are: the basic profile mode (1), the streak mode (2), the downstream vortex mode (3), the spanwise flow modes (4, 5), the normal vortex modes (6, 7), a fully three-dimensional mode (8) and the mode of the modification of the basic profile (9). The first eight modes are the same as those in Waleffe (1997). By applying the Galerkin projection to the Navier-Stokes equations without the pressure term, one can obtain a set of nine ordinary differential equations for the temporal amplitude coefficient ai(t)a_{i}(t), featuring quadratic nonlinearities. In the present study, we follow Martínez-Sánchez et al. (2023) to use the same parameter setup. The ordinary differential equations were numerically solved to obtain 500 sets of solutions, each spanning 4000 time units and the time step is 0.01 time units (400 000 time steps). These data are generated by introducing a random perturbation to a4a_{4} (Martínez-Sánchez et al., 2023). The time lag for the causal analysis is set to one time step.

The results are given in figure 3. It is seen that there are quite a lot of differences between the results of the two approaches. First, we focus on the core sustaining mechanism of near-wall turbulence, i.e. the generation of streaks (mode 2). It has been well recognized that the streaks are generated by the streamwise vortex acting on the mean flow, that is the so-called lift-up mechanism. In the low-order model, it is the causality from mode 1 and 3 to mode 2. We can see from figure 3 that the information flow can successfully capture the causal links of a1a2a_{1}\to a_{2} and a3a2a_{3}\to a_{2}, while the transfer entropy only can predict a3a2a_{3}\to a_{2}. For the other nonlinear interaction mechanism for a2a_{2}, both transfer entropy and information flow can identify the causal links which are present in the dynamical equation of a2a_{2}, except that a4a2a_{4}\to a_{2} is only detected by the transfer entropy. Second, we look into the generation of streamwise vortex (mode 3). Both the transfer entropy and information flow can identify a2a3a_{2}\to a_{3} (streak to streamwise vortex), a6a3a_{6}\to a_{3} (normal vortex to streamwise vortex), and a9a3a_{9}\to a_{3}. However, the causal links of a4a3a_{4}\to a_{3} (spanwise flow to streamwise vortex) and a8a3a_{8}\to a_{3} (three-dimensional mode to streamwise vortex), which exist in the dynamical equation of a3a_{3}, can be found in the result of the transfer entropy analysis, but that of the information flow. Lastly, if further surveying other modes, one can find that the transfer entropy can identify more causal links in the dynamical equations than the information flow.

Refer to caption
Refer to caption
Figure 4: Causal maps of the inner and outer motions of near-wall turbulence at y+=15y^{+}=15: (a) transfer entropy; (b) information flow. Red boxes: causalities identified by both transfer entropy and information flow. Blue boxes: causalities identified by transfer entropy but information flow. Green boxes: causalities identified by information flow but transfer entropy.
Refer to caption
Refer to caption
Figure 5: Causal maps of inner and outer motions of near-wall turbulence at y+=70y^{+}=70: (a) transfer entropy; (b) information flow. Red boxes: causalities identified by both transfer entropy and information flow. Blue boxes: causalities identified by transfer entropy but information flow. Green boxes: causalities identified by information flow but transfer entropy.

Near-wall inner and outer turbulent motions.

Finally, we apply the transfer entropy and information flow to analyze the causalities of inner and outer motions of near-wall turbulence. A time-resolved direct numerical simulation (DNS) of turbulent channel flow at a friction Reynolds number of Reτ=1000Re_{\tau}=1000 is carried out in this work. The DNS adopts a fourth-order accurate compact difference scheme in the homogeneous directions and a second-order accurate central difference scheme in the wall-normal direction for the discretization of the incompressible Navier-Stokes equations on a staggered grid (Hu et al., 2018). A series of low-Reynolds-number channel DNS (up to Reτ=600Re_{\tau}=600) were conducted using the code (Hu & Zheng, 2018; Wang et al., 2021a) and the results were well validated against Lee & Moser (2015). In the present DNS, the flow time span is 19.8δ/uτ\delta/u_{\tau} (δ\delta is channel half-height and uτu_{\tau} is the friction velocity) and the DNS data are stored with a time step of Δt=0.03\Delta t=0.03 (20 000 snapshots in total). The transfer entropy in equation (4) is estimated using a time lag Δt=0.03\Delta t=0.03 and the nearest-neighbour parameter k=4k=4 following Martínez-Sánchez et al. (2023).

Figure 4 shows the causal maps of local velocity fluctuations of inner and outer turbulent motions and pressure fluctuations, i.e. 𝐕=(uS,vS,wS,uL,vL,wL,p)\mathbf{V}=\left(u_{S},v_{S},w_{S},u_{L},v_{L},w_{L},p\right), where uu, vv and ww are the streamwise, wall-normal and spanwise velocity fluctuations, respectively, and pp is pressure fluctuation. The subscripts ”SS” and ”LL” indicate inner and outer motions, respectively. The time series of local velocity and pressure fluctuations at 10 000 random locations at y+=15y^{+}=15 are adopted to calculate the causalities. It is seen that there are many more causalities identified by the transfer entropy including all of those by the information flow. On the other hand, the compactness of the causal links identified by the information flow could provide a concise picture that helps to find the most crucial causalities. For example, one can find the causal links of uSvSu_{S}\leftrightarrow v_{S} and uLvLu_{L}\leftrightarrow v_{L}, which are the streak-vortex self-sustaining cycles of inner and outer motions, respectively. It is also seen that there exists a non-zero causality of uSvLu_{S}\to v_{L}, indicating a possible bottom-up generation mechanism for the wall-normal outer motions. The pressure fluctuation is also found to be an active component in the causal cycles with inner and outer motions.

Figure 5 shows the causal maps of the local velocity and pressure fluctuations of inner and outer turbulent motions at y+=70y^{+}=70. Compared to figure 4, more causal links can be identified as the outer motions are stronger. Similar to figure 4, the transfer entropy can predict more causalities than the information flow. However, we still discuss the results focusing on the causal links that can be identified by both the transfer entropy and information flow. Similarly, the inner and outer self-sustaining causalities of the streak-vortex cycle can also be well observed at this height as well as the bottom-up causality of uSvLu_{S}\to v_{L} and the pressure-velocity causalities. In addition, two new groups of links with non-zero causalities can be found. The first one includes uSwSu_{S}\to w_{S}, uLwLu_{L}\to w_{L} and wSvSw_{S}\to v_{S}, which may be related to the streak instabilities and generation of streamwise vortex. The other one reflects the top-down influence, like (uL,vL,wL)(uS,vS,wS)(u_{L},v_{L},w_{L})\to(u_{S},v_{S},w_{S}). These results indicate that the inner and outer turbulent motions are not only self-sustained but also engage in complex interaction mechanisms that we may not fully understand at the current stage.

SUMMARY

In this work, we have applied two popular causal-analysis methods, i.e. transfer entropy and information flow, to several problems from simple linear and nonlinear dynamical systems to complex turbulent flows. The two methods can predict the same results for the linear problem, but some differences exist for the other nonlinear problems and the complex turbulent-flow problems. To be more specific, more causal links that exist in the dynamical system can be revealed by the information flow in the nonlinear problem. However, for the low-dimensional model of the near-wall cycle of turbulence, it seems that richer causalities can be identified by the transfer entropy. Despite this, both methods can capture the major self-sustaining mechanism of near-wall turbulence, i.e. the streak-vortex regeneration cycle. Some insight into the inner and outer near-wall turbulent motions has been achieved through the present causal analysis. The inner and outer motions may be not only self-sustained but also exhibit non-negligible influence on each other, beyond the common opinion of either bottom-up or top-down, and more details on these will be investigated in future work.

Lastly, we need to mention that the computation is time-consuming using the transfer entropy for high-dimensional problems, while it is much more efficient with the information flow (nearly a hundred of speed-up in our experience). Therefore, we believe that the information flow may have great potential in applications of causal inference for turbulent-flow problems besides the transfer entropy.

ACKNOWLEDGEMENT

The authors are grateful to Prof. X. San Liang for very helpful discussions.

References

  • Adrian (2007) Adrian, RJ 2007 Hairpin vortex organization in wall turbulence. Physics of Fluids 19 (4), 041301.
  • Adrian et al. (2000) Adrian, RJ, Meinhart, CD & Tomkins, CD 2000 Vortex organization in the outer region of the turbulent boundary layer. Journal of Fluid Mechanics 422, 1–54.
  • Baars et al. (2016) Baars, WJ, Hutchins, N & Marusic, I 2016 Spectral stochastic estimation of high-Reynolds-number wall-bounded turbulence for a refined inner-outer interaction model. Physical Review Fluids 1 (5), 054406.
  • Camps-Valls et al. (2023) Camps-Valls, G et al. 2023 Discovering causal relations and equations from data. Physics Reports 1044, 1–68.
  • Cossu & Hwang (2017) Cossu, C & Hwang, Y 2017 Self-sustaining processes at all scales in wall-bounded turbulent shear flows. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences 375 (2089).
  • Granger (1969) Granger, CWJ 1969 Investigating causal relations by econometric models and cross-spectral methods. Econometrica 37 (3), 424–438.
  • Hamilton et al. (1995) Hamilton, JM, Kim, J & Waleffe, F 1995 Regeneration mechanisms of near-wall turbulence structures. Journal of Fluid Mechanics 287, 317–348.
  • Hu et al. (2018) Hu, R, Wang, L, Wang, P, Wang, Y & Zheng, X 2018 Application of high-order compact difference scheme in the computation of incompressible wall-bounded turbulent flows. Computation 6 (2), 31.
  • Hu & Zheng (2018) Hu, R & Zheng, X 2018 Energy contributions by inner and outer motions in turbulent channel flows. Physical Review Fluids 3 (8), 084607.
  • Hunt & Morrison (2000) Hunt, JCR & Morrison, JF 2000 Eddy structure in turbulent boundary layers. European Journal of Mechanics - B/Fluids 19 (5), 673–694.
  • Hutchins & Marusic (2007) Hutchins, N & Marusic, I 2007 Evidence of very long meandering features in the logarithmic region of turbulent boundary layers. Journal of Fluid Mechanics 579, 1–28.
  • Jiménez (2018) Jiménez, J 2018 Coherent structures in wall-bounded turbulence. Journal of Fluid Mechanics 842, P1.
  • Kozachenko & Leonenko (1987) Kozachenko, LF & Leonenko, NN 1987 Sample estimate of the entropy of a random vector. Problemy Peredachi Informatsii 23 (2), 9–16.
  • Lee & Moser (2015) Lee, MK & Moser, RD 2015 Direct numerical simulation of turbulent channel flow up to Reτ5200{R}e_{\tau}\approx 5200. Journal of Fluid Mechanics 774, 395–415.
  • Liang (2013) Liang, XS 2013 The Liang-Kleeman information flow: Theory and applications. Entropy 15 (1), 327–360.
  • Liang (2014) Liang, XS 2014 Unraveling the cause-effect relation between time series. Physical Review E 90 (5), 052150.
  • Liang (2016) Liang, XS 2016 Information flow and causality as rigorous notions ab initio. Physical Review E 94 (5), 052201.
  • Liang (2021) Liang, XS 2021 Normalized multivariate time series causality analysis and causal graph reconstruction. Entropy 23 (6), 679.
  • Liang & Kleeman (2005) Liang, XS & Kleeman, R 2005 Information transfer between dynamical system components. Physical Review Letters 95 (24), 244101.
  • Liang & Lozano-Durán (2016) Liang, XS & Lozano-Durán, A 2016 A preliminary study of the causal structure in fully developed near-wall turbulence. CTR Proc. of the Summer Prog p. 233.
  • Lizier (2014) Lizier, JT 2014 JIDT: An information-theoretic toolkit for studying the dynamics of complex systems. Frontiers in Robotics and AI 1, 11.
  • Lozano-Durán et al. (2020) Lozano-Durán, A, Bae, HJ & Encinar, MP 2020 Causality of energy-containing eddies in wall turbulence. Journal of Fluid Mechanics 882, A2.
  • Martínez-Sánchez et al. (2023) Martínez-Sánchez, Á, López, E, Le Clainche, S, Lozano-Durán, A, Srivastava, A & Vinuesa, R 2023 Causality analysis of large-scale structures in the flow around a wall-mounted square cylinder. Journal of Fluid Mechanics 967, A1.
  • Marusic et al. (2010) Marusic, I, Mathis, R & Hutchins, N 2010 Predictive model for wall-bounded turbulent flow. Science 329 (5988), 193–196.
  • Mathis et al. (2011) Mathis, R, Hutchins, N & Marusic, I 2011 A predictive inner-outer model for streamwise turbulence statistics in wall-bounded flows. Journal Fluid Mechanics 681, 537–566.
  • Moehlis et al. (2004) Moehlis, J, Faisst, H & Eckhardt, B 2004 A low-dimensional model for turbulent shear flows. New Journal of Physics 6 (1), 56.
  • Paluš et al. (2018) Paluš, M, Krakovská, A, Jakubík, J & Chvosteková, M 2018 Causality, dynamical systems and the arrow of time. Chaos: An Interdisciplinary Journal of Nonlinear Science 28 (7).
  • Pearl (2009) Pearl, J 2009 Causal inference in statistics: An overview. Statistics Surveys 3, 96–146.
  • Robinson (1991) Robinson, SK 1991 Coherent motions in the turbulent boundary layer. Annual Review of Fluid Mechanics 23, 601–639.
  • Runge et al. (2023) Runge, J et al. 2023 Causal inference for time series. Nature Reviews Earth & Environment 4 (7), 487–505.
  • Schreiber (2000) Schreiber, T 2000 Measuring information transfer. Physical Review Letters 85 (2), 461.
  • Shannon (1948) Shannon, CE 1948 A mathematical theory of communication. The Bell System Technical Journal 27 (3), 379–423.
  • Smits et al. (2011) Smits, AJ, McKeon, BJ & Marusic, I 2011 High-Reynolds number wall turbulence. Annual Review of Fluid Mechanics 43, 353–375.
  • Toh & Itano (2005) Toh, S & Itano, T 2005 Interaction between a large-scale structure and near-wall structures in channel flow. Journal of Fluid Mechanics 524 (2005), 249–262.
  • Waleffe (1997) Waleffe, F 1997 On a self-sustaining process in shear flows. Physics of Fluids 9 (4), 883–900.
  • Wang et al. (2021a) Wang, L, Hu, R & Zheng, X 2021a A scaling improved inner–outer decomposition of near-wall turbulent motions. Physics of Fluids 33 (4), 045120.
  • Wang et al. (2022) Wang, W, Lozano-Durán, A, Helmig, R & Xu, C 2022 Spatial and spectral characteristics of information flux between turbulent boundary layers and porous media. Journal of Fluid Mechanics 949, A16.
  • Wang et al. (2021b) Wang, W, Xu, C, Lozano-Durán, A, Helmig, R & Weigand, B 2021b Information transfer between turbulent boundary layers and porous media. Journal of Fluid Mechanics 920, A21.
  • Zhou et al. (2022) Zhou, Z, Xu, C-X & Jiménez, J 2022 Interaction between near-wall streaks and large-scale motions in turbulent channel flows. Journal of Fluid Mechanics 940, A23.