skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on December 10, 2025

Title: Informative and non-informative decomposition of turbulent flow fields
Not all the information in a turbulent field is relevant for understanding particular regions or variables in the flow. Here, we present a method for decomposing a source field into its informative$$\boldsymbol {\varPhi }_{I}(\boldsymbol {x},t)$$and residual$$\boldsymbol {\varPhi }_{R}(\boldsymbol {x},t)$$components relative to another target field. The method is referred to as informative and non-informative decomposition (IND). All the necessary information for physical understanding, reduced-order modelling and control of the target variable is contained in$$\boldsymbol {\varPhi }_{I}(\boldsymbol {x},t)$$, whereas$$\boldsymbol {\varPhi }_{R}(\boldsymbol {x},t)$$offers no substantial utility in these contexts. The decomposition is formulated as an optimisation problem that seeks to maximise the time-lagged mutual information of the informative component with the target variable while minimising the mutual information with the residual component. The method is applied to extract the informative and residual components of the velocity field in a turbulent channel flow, using the wall shear stress as the target variable. We demonstrate the utility of IND in three scenarios: (i) physical insight into the effect of the velocity fluctuations on the wall shear stress; (ii) prediction of the wall shear stress using velocities far from the wall; and (iii) development of control strategies for drag reduction in a turbulent channel flow using opposition control. In case (i), IND reveals that the informative velocity related to wall shear stress consists of wall-attached high- and low-velocity streaks, collocated with regions of vertical motions and weak spanwise velocity. This informative structure is embedded within a larger-scale streak–roll structure of residual velocity, which bears no information about the wall shear stress. In case (ii), the best-performing model for predicting wall shear stress is a convolutional neural network that uses the informative component of the velocity as input, while the residual velocity component provides no predictive capabilities. Finally, in case (iii), we demonstrate that the informative component of the wall-normal velocity is closely linked to the observability of the target variable and holds the essential information needed to develop successful control strategies.  more » « less
Award ID(s):
2140775
PAR ID:
10637468
Author(s) / Creator(s):
;
Publisher / Repository:
Cambridge University
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
1000
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Let$$\varphi _1,\ldots ,\varphi _r\in {\mathbb Z}[z_1,\ldots z_k]$$be integral linear combinations of elementary symmetric polynomials with$$\text {deg}(\varphi _j)=k_j\ (1\le j\le r)$$, where$$1\le k_1<\cdots . Subject to the condition$$k_1+\cdots +k_r\ge \tfrac {1}{2}k(k-~1)+2$$, we show that there is a paucity of nondiagonal solutions to the Diophantine system$$\varphi _j({\mathbf x})=\varphi _j({\mathbf y})\ (1\le j\le r)$$. 
    more » « less
  2. Well-resolved direct numerical simulations (DNS) have been performed of the flow in a smooth circular pipe of radius$$R$$and axial length$$10{\rm \pi} R$$at friction Reynolds numbers up to$$Re_\tau =5200$$using the pseudo-spectral code OPENPIPEFLOW. Various turbulence statistics are documented and compared with other DNS and experimental data in pipes as well as channels. Small but distinct differences between various datasets are identified. The friction factor$$\lambda$$overshoots by$$2\,\%$$and undershoots by$$0.6\,\%$$the Prandtl friction law at low and high$$Re$$ranges, respectively. In addition,$$\lambda$$in our results is slightly higher than in Pirozzoliet al.(J. Fluid Mech., vol. 926, 2021, A28), but matches well the experiments in Furuichiet al.(Phys. Fluids, vol. 27, issue 9, 2015, 095108). The log-law indicator function, which is nearly indistinguishable between pipe and channel up to$$y^+=250$$, has not yet developed a plateau farther away from the wall in the pipes even for the$$Re_\tau =5200$$cases. The wall shear stress fluctuations and the inner peak of the axial turbulence intensity – which grow monotonically with$$Re_\tau$$– are lower in the pipe than in the channel, but the difference decreases with increasing$$Re_\tau$$. While the wall value is slightly lower in the channel than in the pipe at the same$$Re_\tau$$, the inner peak of the pressure fluctuation shows negligible differences between them. The Reynolds number scaling of all these quantities agrees with both the logarithmic and defect-power laws if the coefficients are properly chosen. The one-dimensional spectrum of the axial velocity fluctuation exhibits a$$k^{-1}$$dependence at an intermediate distance from the wall – also seen in the channel. In summary, these high-fidelity data enable us to provide better insights into the flow physics in the pipes as well as the similarity/difference among different types of wall turbulence. 
    more » « less
  3. Abstract Cohesive powersof computable structures are effective analogs of ultrapowers, where cohesive sets play the role of ultrafilters. Let$$\omega $$,$$\zeta $$, and$$\eta $$denote the respective order-types of the natural numbers, the integers, and the rationals when thought of as linear orders. We investigate the cohesive powers of computable linear orders, with special emphasis on computable copies of$$\omega $$. If$$\mathcal {L}$$is a computable copy of$$\omega $$that is computably isomorphic to the usual presentation of$$\omega $$, then every cohesive power of$$\mathcal {L}$$has order-type$$\omega + \zeta \eta $$. However, there are computable copies of$$\omega $$, necessarily not computably isomorphic to the usual presentation, having cohesive powers not elementarily equivalent to$$\omega + \zeta \eta $$. For example, we show that there is a computable copy of$$\omega $$with a cohesive power of order-type$$\omega + \eta $$. Our most general result is that if$$X \subseteq \mathbb {N} \setminus \{0\}$$is a Boolean combination of$$\Sigma _2$$sets, thought of as a set of finite order-types, then there is a computable copy of$$\omega $$with a cohesive power of order-type$$\omega + \boldsymbol {\sigma }(X \cup \{\omega + \zeta \eta + \omega ^*\})$$, where$$\boldsymbol {\sigma }(X \cup \{\omega + \zeta \eta + \omega ^*\})$$denotes the shuffle of the order-types inXand the order-type$$\omega + \zeta \eta + \omega ^*$$. Furthermore, ifXis finite and non-empty, then there is a computable copy of$$\omega $$with a cohesive power of order-type$$\omega + \boldsymbol {\sigma }(X)$$. 
    more » « less
  4. The rheological behaviour of dense suspensions of ideally conductive particles in the presence of both electric field and shear flow is studied using large-scale numerical simulations. Under the action of an electric field, these particles are known to undergo dipolophoresis (DIP), which is the combination of two nonlinear electrokinetic phenomena: induced-charge electrophoresis (ICEP) and dielectrophoresis (DEP). For ideally conductive particles, ICEP is predominant over DEP, resulting in transient pairing dynamics. The shear viscosity and first and second normal stress differences$$N_1$$and$$N_2$$of such suspensions are examined over a range of volume fractions$$15\,\% \leq \phi \leq 50\,\%$$as a function of Mason number$$Mn$$, which measures the relative importance of viscous shear stress over electrokinetic-driven stress. For$$Mn < 1$$or low shear rates, the DIP is shown to dominate the dynamics, resulting in a relatively low-viscosity state. The positive$$N_1$$and negative$$N_2$$are observed at$$\phi < 30\,\%$$, which is similar to Brownian suspensions, while their signs are reversed at$$\phi \ge 30\,\%$$. For$$Mn \ge 1$$, the shear thickening starts to arise at$$\phi \ge 30\,\%$$, and an almost five-fold increase in viscosity occurs at$$\phi = 50\,\%$$. Both$$N_1$$and$$N_2$$are negative for$$Mn \gg 1$$at all volume fractions considered. We illuminate the transition in rheological behaviours from DIP to shear dominance around$$Mn = 1$$in connection to suspension microstructure and dynamics. Lastly, our findings reveal the potential use of nonlinear electrokinetics as a means of active rheology control for such suspensions. 
    more » « less
  5. Abstract Letfbe an$$L^2$$-normalized holomorphic newform of weightkon$$\Gamma _0(N) \backslash \mathbb {H}$$withNsquarefree or, more generally, on any hyperbolic surface$$\Gamma \backslash \mathbb {H}$$attached to an Eichler order of squarefree level in an indefinite quaternion algebra over$$\mathbb {Q}$$. Denote byVthe hyperbolic volume of said surface. We prove the sup-norm estimate$$\begin{align*}\| \Im(\cdot)^{\frac{k}{2}} f \|_{\infty} \ll_{\varepsilon} (k V)^{\frac{1}{4}+\varepsilon} \end{align*}$$ with absolute implied constant. For a cuspidal Maaß newform$$\varphi $$of eigenvalue$$\lambda $$on such a surface, we prove that$$\begin{align*}\|\varphi \|_{\infty} \ll_{\lambda,\varepsilon} V^{\frac{1}{4}+\varepsilon}. \end{align*}$$ We establish analogous estimates in the setting of definite quaternion algebras. 
    more » « less