A new theoretical approach for turbulent flows based on Lie-group analysis is presented. It unifies a large set of ‘solutions’ for the mean velocity of stationary parallel turbulent shear flows. These results are not solutions in the classical sense but instead are defined by the maximum number of possible symmetries, only restricted by the flow geometry and other external constraints. The approach is derived from the Reynolds-averaged Navier–Stokes equations, the fluctuation equations, and the velocity product equations, which are the dyad product of the velocity fluctuations with the equations for the velocity fluctuations. The results include the logarithmic law of the wall, an algebraic law, the viscous sublayer, the linear region in the centre of a Couette flow and in the centre of a rotating channel flow, and a new exponential mean velocity profile not previously reported that is found in the mid-wake region of high Reynolds number flat-plate boundary layers. The algebraic scaling law is confirmed in both the centre and the near-wall regions in both experimental and DNS data of turbulent channel flows. In the case of the logarithmic law of the wall, the scaling with the distance from the wall arises as a result of the analysis and has not been assumed in the derivation. All solutions are consistent with the similarity of the velocity product equations to arbitrary order. A method to derive the mean velocity profiles directly from the two-point correlation equations is shown.

A new direct numerical simulation of a Poiseuille channel flow has been conducted for a friction Reynolds number of $10000$, using the pseudospectral code LISO. The mean streamwise velocity presents a long logarithmic layer, extending from 400 to 2500 wall units, longer than it was thought. The maximum of the intensity of the streamwise velocity increases with the Reynolds number, as expected. Also, the elusive second maximum of this intensity has not appeared yet. In case it exists, its location will be around $y^{+}\approx 120$, for a friction Reynolds number extrapolated to approximately 13 500. The small differences in the near-wall gradient of this intensity for several Reynolds numbers are related to the scaling failure of the dissipation, confirming this hypothesis. The scaling of the turbulent budgets in the center of the channel is almost perfect above 1000 wall units. Finally, the peak of the pressure intensity grows with the Reynolds number and does not scale in wall units. If the pressure at the wall is modeled as an inverse quadratic power of ${\text{Re}}_{\tau }$, then $p_{\infty }^{\prime +}\approx 4.7$ at the limit of infinite Reynolds number.

A new set of numerical simulations of turbulent plane Couette flow in a large box of dimension ( $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}20\pi h,\, 2h,\, 6\pi h$ ) at Reynolds number $(\mathit{Re}_{\tau }) =125$ , 180, 250 and 550 is described and compared with simulations at lower Reynolds numbers, Poiseuille flows and experiments. The simulations present a logarithmic near-wall layer and are used to verify and revise previously known results. It is confirmed that the fluctuation intensities in the streamwise and spanwise directions do not scale well in wall units. The scaling failure occurs both near to and away from the wall. On the contrary, the wall-normal intensity scales in inner units in the near-wall region and in outer units in the core region. The spectral ridge found by Hoyas & Jiménez (Phys. Fluids, vol. 18, 2003, 011702) for the turbulent Poiseuille flow can also be seen in the present flow. Away from the wall, very large-scale motions are found spanning through all the length of the channel. The statistics of these simulations can be downloaded from the webpage of the Chair of Fluid Dynamics.