2. Partitioning framework

The VGT, $\boldsymbol{\mathsf{A}} = \boldsymbol {\nabla } \boldsymbol {u}$, can be expressed in its principal reference frame, denoted by $({\cdot })^*$, as $\boldsymbol{\mathsf{A}}^* = \boldsymbol{\mathsf{Q}}\boldsymbol{\mathsf{A}}\boldsymbol{\mathsf{Q}}^{\rm T}$, where $\boldsymbol{\mathsf{Q}}$ is unitary and $({\cdot })^{\rm T}$ represents the transpose. In this frame, the VGT is quasi-triangular and it can be decomposed as

(2.1) \begin{equation} \boldsymbol{\mathsf{A}}^* = \underbrace{\begin{bmatrix} \dot{\epsilon}_1^* & 0 & 0 \\ 0 & \dot{\epsilon}_2^* & 0 \\ 0 & 0 & \dot{\epsilon}_3^* \end{bmatrix}}_{\textstyle\boldsymbol{\mathsf{A}}_\epsilon^*} + \underbrace{\begin{bmatrix} 0 & 0 & 0 \\ 0 & 0 & \dot{\varphi}_1^* \\ 0 & -\dot{\varphi}_1^* & 0 \end{bmatrix}}_{\textstyle\boldsymbol{\mathsf{A}}_\varphi^*} + \underbrace{\begin{bmatrix} 0 & \dot{\gamma}_3^* & \dot{\gamma}_2^* \\ 0 & 0 & \dot{\gamma}_1^* \\ 0 & 0 & 0 \end{bmatrix}}_{\textstyle\boldsymbol{\mathsf{A}}_\gamma^*}, \end{equation}

where $\boldsymbol{\mathsf{A}}_\epsilon ^*$, $\boldsymbol{\mathsf{A}}_\varphi ^*$ and $\boldsymbol{\mathsf{A}}_\gamma ^*$ denote the normal straining, rigid rotation and pure shearing tensors, respectively. These tensors can be determined and transformed to the original coordinate system using the ordered real Schur decomposition of $\boldsymbol{\mathsf{A}}$ (Kronborg & Hoffman Reference Kronborg and Hoffman2023).

Correspondingly, the strength of the velocity gradients can be expressed as

(2.2) \begin{equation} A^2 = {\rm tr}(\boldsymbol{\mathsf{A}}^{\rm T} \boldsymbol{\mathsf{A}} ) = \underbrace{{\rm tr}( \boldsymbol{\mathsf{A}}_\epsilon^{\rm T} \boldsymbol{\mathsf{A}}_\epsilon^{})}_{\textstyle A_\epsilon^2} + \underbrace{{\rm tr}( \boldsymbol{\mathsf{A}}_\varphi^{\rm T} \boldsymbol{\mathsf{A}}_\varphi^{})}_{\textstyle A_\varphi^2} + \underbrace{{\rm tr}( \boldsymbol{\mathsf{A}}_\gamma^{\rm T} \boldsymbol{\mathsf{A}}_\gamma^{})}_{\textstyle A_\gamma^2} + \underbrace{2 {\rm tr}( \boldsymbol{\mathsf{A}}_\varphi^{\rm T} \boldsymbol{\mathsf{A}}_\gamma^{})}_{\textstyle A_{\varphi\gamma}^2}, \end{equation}

where ${\rm tr}({\cdot })$ represents the trace. The first three terms represent the strengths of the constituents in (2.1) and the last term represents the interaction between shearing and rigid rotation. The velocity gradient partitioning is defined in terms of the relative contributions of these constituents to $A^2$. In normalized form, these contributions are bounded as $A_\epsilon ^2 / A^2 \in [0,1]$, $A_\varphi ^2 / A^2 \in [0,1]$, $A_\gamma ^2 / A^2 \in [0,1]$ and $A_{\varphi \gamma }^2 / A^2 \in [0,(\sqrt {2}+1)^{-1}]$ (Das & Girimaji Reference Das and Girimaji2020).

Ensemble averaging can be used to characterize the statistical relevance of the contributions in this partitioning. We define the averaged partitioning as

(2.3) \begin{equation} \langle {A_\zeta^2} \rangle_{A^2} = \langle A_\zeta^2 \rangle/ \langle A^2 \rangle, \quad \zeta \in \{\epsilon,\varphi,\gamma,\varphi\gamma\}, \end{equation}

where $\langle ({\cdot }) \rangle$ denotes averaging over homogeneous spatial directions and time. The present study focuses on this averaged partitioning, which can be used to characterize contributions to enstrophy and dissipation (Das & Girimaji Reference Das and Girimaji2020; Arun & Colonius Reference Arun and Colonius2024).

The practical relevance of the partitioning is tied, in part, to its ability to inform modelling efforts, including for non-canonical flows. Some modelling paradigms, such as large-eddy simulations (LES) and Lagrangian approaches, operate on the basis of the total velocity. Other paradigms, such as the Reynolds-averaged Navier–Stokes (RANS) equations and input-output analyses, operate on the basis of the velocity fluctuations. In these contexts, it is therefore important to distinguish the partitioning of the total velocity gradients from the partitioning of the velocity gradient fluctuations, which we denote as $\langle {A_\zeta ^2} \rangle _{A^2}$ and $\langle {A_\zeta ^2} \rangle _{A^2}'$, respectively. This distinction is particularly relevant since the normality-based triple decomposition does not generally commute with filtering or averaging operations. Physically, the modes of deformation for the total partitioning reflect what a fluid parcel would actually experience. By contrast, the modes for the fluctuation partitioning reflect what it would experience if advected only by the velocity fluctuations. We consider both the total partitioning and the fluctuation partitioning for the wall-bounded flows in the present study.

3. Turbulence datasets

As summarized in table 1, we analyse the partitioning in several well-validated turbulence datasets obtained from direct numerical simulations. These datasets include forced isotropic turbulence (FIT315 and FIT610) and wall-bounded turbulence (Ch0186, Ch1000, BL0729 and BL1024) over a broad range of Reynolds numbers. Cases FIT610 and Ch1000 are obtained from the Johns Hopkins Turbulence Databases (Li et al. Reference Li, Perlman, Wan, Yang, Meneveau, Burns, Chen, Szalay and Eyink2008) and Ch0186 represents a minimal flow unit for near-wall turbulence (Jiménez & Moin Reference Jiménez and Moin1991). The references in table 1 provide further computational details and validation for each dataset.

Table 1.

Turbulence datasets considered in the present analyses. Taylor-scale and friction Reynolds numbers are denoted by $Re_\lambda$ and $Re_\tau$, respectively, and $N_t$ denotes the number of snapshots. The grid sizes correspond to the streamwise ($x$), wall-normal ($y$) and spanwise ($z$) directions, respectively.

We establish the isotropic partitioning using FIT610 and test its sensitivity to $Re_\lambda$ using FIT315. The snapshots for FIT315 are spaced roughly one integral time unit apart and $N_t$ is selected to produce a similar number of samples to FIT610. We use Ch1000 to characterize the partitioning in a turbulent channel at a moderate $Re_\tau$. Its snapshots are spaced roughly 0.65 eddy turnover time units apart and they span roughly one flow-through time unit (Graham et al. Reference Graham2016). Case Ch0186 allows us to investigate how the wall-bounded partitioning changes when $Re_\tau$ is barely large enough to sustain turbulence. Its snapshots span roughly 160 eddy turnover time units. We use BL0729 and BL1024 to further characterize the partitioning for wall-bounded turbulence subject to mild spatial development. Their broad ranges of $Re_\tau$ allow us to characterize how the partitioning evolves as the flows become increasingly turbulent. The snapshots for BL0729 and BL1024 span more than 20 and 7 eddy turnover time units, respectively (Towne et al. Reference Towne2023).

When computing the VGT, we adopt the differentiation techniques employed in the original simulations where possible. We further require that all elements of the VGT be collocated prior to partitioning. For FIT315 and FIT610, we use a spectral method to compute all velocity gradients. For Ch0186, we use the second-order accurate staggered finite differences employed in the original simulation and subsequently shift staggered quantities to cell centres. This shifting is performed by adjusting the phases of the Fourier modes in $x$ and $z$ and by averaging adjacent values in $y$ (Arun, Bae & McKeon Reference Arun, Bae and McKeon2023). For Ch1000, we use a spectral method in $x$ and $z$ and collocated finite differences with a stencil size of $N_s = 7$ in $y$. The original BL0729 and BL1024 simulations employed second-order accurate staggered finite differences; however, the published datasets are collocated and subsampled by a factor of two in $y$ and $z$ (Towne et al. Reference Towne2023). As a result, we employ collocated finite differences with $N_s = 3$ in $x$ and $y$. Since the $z$ direction has periodic boundary conditions, we compute spanwise derivatives using a spectral method with the modified wavenumbers associated with the original staggered finite difference scheme.

4. Partitioning in nearly isotropic turbulence

The isotropic velocity gradient partitioning characterizes the contributions of $\boldsymbol{\mathsf{A}}_\epsilon$, $\boldsymbol{\mathsf{A}}_\varphi$ and $\boldsymbol{\mathsf{A}}_\gamma$ in the idealized setting of forced isotropic turbulence. It has been established previously for $Re_\lambda \approx 1\unicode{x2013}588$ (Das & Girimaji Reference Das and Girimaji2020) and is roughly invariant for $Re_\lambda \gtrsim 200$. Here, we use FIT610 to confirm this isotropic partitioning, which we denote by $\langle {A_\zeta ^2} \rangle _{A^2}^{iso}$. We characterize deviations from this partitioning using the following metric:

(4.1) \begin{equation} \varDelta_{iso} = \left.\sum_{\zeta \in \{\epsilon,\varphi,\gamma,\varphi\gamma\}} \left\lvert \langle A_\zeta^2 \rangle_{A^2} - \langle A_\zeta^2 \rangle_{A^2}^{iso} \right\rvert \right/ \sum_{\zeta \in \{\epsilon,\varphi,\gamma,\varphi\gamma\}} \left\lvert \langle A_\zeta^2 \rangle_{A^2}^{iso} \right\rvert,\end{equation}

in which the denominator sums to unity. An analogous metric, $\varDelta _{iso}'$, is defined for the fluctuation partitioning by replacing $\langle {A_\zeta ^2} \rangle _{A^2}$ with $\langle {A_\zeta ^2} \rangle _{A^2}'$ in (4.1). One advantage of these metrics is that they are not affected by further decomposing $\boldsymbol{\mathsf{A}}_{\gamma }$ into its symmetric and antisymmetric parts. In the present study, they produce results similar to those produced by relative root-mean-square deviations, which are more commonplace.

Table 2 shows the partitioning alongside the deviation metrics for each dataset. Consistent with previous results (Das & Girimaji Reference Das and Girimaji2020), the FIT315 partitioning is nearly identical to the FIT610 partitioning. For the wall-bounded flows, the partitioning is taken from regions where the mean shearing is relatively weak. Further, for the boundary layers, it is taken sufficiently far from the boundary-layer thickness to mitigate the imprint of the exterior potential flow.

Table 2.

Velocity gradient partitioning for each flow and the corresponding deviation metrics. The partitioning is reported at the channel centreline for Ch0186 and Ch1000 and at $(Re_\tau,y^+) \approx (729,159)$ and $(1000,155)$ for BL0729 and BL1024, respectively. The column shadings reflect our partitioning colour scheme.

The partitioning in the selected regions of Ch1000, BL0729 and BL1024 is remarkably similar to the isotropic partitioning, with deviations of 1.2 % or less. Since the mean flow has a minimal imprint on the velocity gradients in these regions, the total partitioning is quite similar to the fluctuation partitioning. These results highlight that the isotropic partitioning is broadly applicable in appropriate regions of inhomogeneous turbulent flows.

5. Partitioning in wall-bounded turbulence

5.1. Effect of mean shearing

The velocity gradient partitioning in wall-bounded turbulence is heavily influenced by the mean shearing imposed by the wall. Figure 1 shows the partitioning profiles as a function of both wall-normal location (in inner units) and mean shearing strength. Consistent with table 2, the partitioning approaches the isotropic values far from the wall, starting near the top of the log-law region. However, the boundary-layer partitioning begins to diverge from these values near the edge of the boundary layer, reflecting the transition to a potential flow. Beyond the boundary-layer thickness (not shown), this transition is associated with monotonic enhancement of normal straining and decay of the other constituents. As observed in table 2, the Ch0186 partitioning does not converge as well onto the isotopic values since it is barely turbulent. The partitioning throughout this minimal channel is generally associated with enhanced contributions from shearing.

Figure 1.

Total (a,c,e) and fluctuation (b,d,f) partitioning profiles for the channels and boundary layers in terms of wall-normal location in inner units (a,b) and mean shearing strength (c–f). The vertical lines represent the isotropic values. The BL0729 and BL1024 profiles are shown for $Re_\tau \approx 729$ and $1000$, respectively, and the top boundary of the log-law region represents Ch1000. The dashed white lines in (a,b) represent the locations of the partitioning values reported in table 2 for BL0729 and BL1024. In (a,b), the markers are used to distinguish between the profiles and, in (c–f), they represent actual data points. In (c–f), the mean shearing axis is reversed and the dashed lines represent comparable linear-log trends for each dataset, with the partitioning as the dependent variable.

For all wall-normal profiles in figure 1, the near-wall partitioning is dominated by shearing. This feature reflects the strong imprint of the mean shearing on the near-wall velocity gradients. The imprint of the mean shearing is also responsible for the enhanced contribution of shearing (at the expense of the other constituents) to the total partitioning relative to the fluctuation partitioning for $y^+ \lesssim 20$. Between the near-wall regime and the nearly isotropic regime far from the wall, the collapse of the partitioning profiles is particularly striking for BL0729, BL1024 and (to a lesser extent) Ch1000. The most significant differences between the channel and boundary-layer profiles occur in the buffer layer. For the channels, the contributions of shearing and normal straining in this region are enhanced and reduced, respectively, relative to their contributions in the boundary layers. While beyond the scope of the present work, characterizing how the flow structures in the buffer layer reflect these differences (e.g. through the imprint of the exterior potential flow) would be interesting future work.

The mean shearing strength profiles in figure 1 provide a complementary view to the wall-normal profiles. The mean shearing strength parameter is normalized as $\bar {A}_\gamma ^2 / \bar {A}_{\gamma,w}^2$, where $\bar {A}_{\gamma,w}^2$ represents the (maximum) value at the wall. This parameter quantifies the effect of the wall in terms of velocity gradients and can be determined directly from the mean flow. Further, for the channels and boundary layers we consider, it can be well-approximated using the wall-normal gradient of the mean streamwise velocity, $\partial \bar {u} / \partial y$. For these profiles, we focus primarily on the buffer layer and log-law region since they capture the majority of the evolution from the near-wall regime to the nearly isotropic regime. As observed for the wall-normal profiles, the partitioning in these regions collapses very well as a function of the mean shearing strength for BL0729, BL1024 and (to a lesser extent) Ch1000.

The partitioning in wall-bounded flows becomes similar to the isotropic partitioning when $\bar {A}_\gamma ^2 / \bar {A}_{\gamma,w}^2 \lesssim 10^{-3}$. The dashed lines in figure 1 illustrate the similarity of the partitioning profiles to log-linear variations with the mean shearing strength for $10^{-3} \lesssim \bar {A}_\gamma ^2 / \bar {A}_{\gamma,w}^2 \lesssim 10^{-1}$. The slopes of these variations for the total partitioning are slightly steeper than those for the fluctuation partitioning due to the enhanced contribution from the mean shearing for $\bar {A}_\gamma ^2 / \bar {A}_{\gamma,w}^2 \gtrsim 10^{-1}$. The slopes for shearing and normal straining tend to be similar in magnitude and steeper than those of rigid rotation and shear-rotation interactions. This feature mirrors the relative contributions in the isotropic partitioning, for which $\langle {A_\epsilon ^2} \rangle _{A^2}^{iso} \sim \tfrac {1}{2}\langle {A_\gamma ^2} \rangle _{A^2}^{iso} \sim \langle {A_\varphi ^2} \rangle _{A^2}^{iso} + \langle {A_{\varphi \gamma }^2} \rangle _{A^2}^{iso}$.

More rigorous modelling approaches may help enable predictions of the partitioning profiles in terms of mean flow variables. While we do not propose an explicit model for these profiles in the present study, our analysis suggests that both $y^+$ and the mean shearing are appropriate for modelling the partitioning and sufficiently high $Re_\tau$. Beyond mean flow variables, characterizing how the strength of the velocity gradient fluctuations relative to the mean shearing strength impacts these profiles may provide further insight.

The pronounced, roughly monotonic variations of the normality-based partitioning profiles strikingly capture the development from the near-wall regime to the nearly isotropic regime far from the wall. By contrast, as depicted and discussed in the Appendix, the variations for the symmetry-based partitioning are non-monotonic and do not exceed $\pm 0.02$ of their isotropic values in the regimes of interest. Supplemented by previous findings (Arun & Colonius Reference Arun and Colonius2024), these results highlight that the expressivity of the normality-based partitioning provides a key advantage over considering symmetry alone.

5.2. Effect of friction Reynolds number

Beyond wall-normal variations, the streamwise spatial development in the boundary layers is associated with increasing $Re_\tau$. Figure 2 shows how the region where the isotropic partitioning is applicable evolves as a function of $Re_\tau$. For BL0729, this region grows appreciably in size as $Re_\tau$ increases. The same is true to a lesser extent for BL1024, where the turbulence is more well-developed throughout the domain. The $\varDelta _{iso}$ and $\varDelta _{iso}'$ contours identify regions where the isotropic partitioning is most applicable. Throughout BL1024 and at the downstream end of BL0729, it is most applicable around $y^+ \approx 170$ and $y^+ \approx 140$ for the total and fluctuating velocity gradients, respectively. The domain of applicability of the isotropic partitioning is slightly larger for the velocity fluctuations than for the total velocity.

Figure 2.

Streamwise development of BL0729 (a,b) and BL1024 (c,d) in terms of $Re_\tau$, where the colour axis represents $\bar {A}_\gamma ^2 / \bar {A}_{\gamma,w}^2$. The white, grey and black contours represent $\varDelta _{iso} = 1\,\%$, $\varDelta _{iso} = 2\,\%$ and $\varDelta _{iso} = 5\,\%$, respectively, for the total partitioning (a,c) and $\varDelta _{iso}' = 1\,\%$, $\varDelta _{iso}' = 2\,\%$ and $\varDelta _{iso}' = 5\,\%$, respectively, for the fluctuation partitioning (b,d). The dashed and dash-dotted black lines represent the top of the viscous sublayer $(y^+ = 5)$ and the top of the buffer layer $(y^+ = 30)$, respectively, and the dotted black lines represent $y^+ = 100$ and $y^+ = 150$. The black circles represent the locations of the partitioning values reported in table 2 and $\delta$ represents the boundary-layer thickness.

Except at the upstream end of BL0729, $\varDelta _{iso}$ and $\varDelta _{iso}'$ are less than roughly 5 % for $y^+ \gtrsim 100$. Alongside figure 1, these results complement the findings of Johnson et al. (Reference Johnson, Hamilton, Burns and Meneveau2017), which suggest that velocity gradient statistics associated with vortex stretching approach their isotropic values for $y^+ \gtrsim 100$. Therefore, our results further support their conclusion that this collapse provides support for the local isotropy hypothesis sufficiently far from the wall.

The inner unit scaling of this collapse contrasts with the boundary-layer thickness, which scales in outer units and captures the divergence from the isotropic partitioning near the free stream. Our results therefore suggest that, while inner unit scalings can be used to determine when the partitioning approaches the isotropic values, outer unit scalings may more appropriately capture the divergence to a potential flow in the boundary layers. They further suggest that $Re_\tau \gtrsim 700$ is a reasonable regime in which to expect a collapsed partitioning. This estimate is consistent with the collapse of the partitioning profiles for BL0729, BL1024 and (to a lesser extent) Ch1000 as well as the lack of collapse for Ch0186 in figure 1.

6. Concluding remarks

We have analysed the normality-based partitioning of velocity gradients in several canonical turbulent flows. The partitioning we compute for forced isotropic turbulence agrees well with previous results (Das & Girimaji Reference Das and Girimaji2020). Moreover, we show that the isotropic partitioning also applies to velocity gradients near and beyond the top of the log-law region in wall-bounded flows over a broad range of $Re_\tau$. The broad applicability of the isotropic partitioning for $y^+ \gtrsim 100$ complements previous results (Johnson et al. Reference Johnson, Hamilton, Burns and Meneveau2017), thereby providing further support for the local isotropy hypothesis for well-developed turbulence far from solid boundaries. Our results suggest that $Re_\tau \gtrsim 700$ is sufficiently high to expect the partitioning profiles to collapse as they transition from the shearing-dominated near-wall regime to the nearly isotropic regime. Further, the mean shearing provides a reasonable mean flow parameter for modelling their variations in the buffer layer and the log-law region. Altogether, our results highlight expressivity as a key advantage of the normality-based partitioning over symmetry-based approaches.

Moving forward, analysing the partitioning profiles at higher $Re_\tau$ would help further characterize their sensitivity and collapse for wall-bounded flows. Developing more rigorous models for the partitioning in terms of mean flow variables would also be useful, especially models that do not depend strongly on the flow configuration. The partitioning may aid turbulence modelling efforts in RANS, LES, Lagrangian or input-output settings, e.g. by directly informing closure models or by providing an evaluation metric for models of interest. Finally, connecting the statistical features we report to the turbulence structures that produce them would provide an enhanced view of the roles of the partitioning constituents.