Explore the vast range of titles available.

Reynolds-averaged Navier–Stokes (RANS) simulations with turbulence closure models continue to play important roles in industrial flow simulations. However, the commonly used linear eddy-viscosity models are intrinsically unable to handle flows with non-equilibrium turbulence (e.g. flows with massive separation). Reynolds stress models, on the other hand, are plagued by their lack of robustness. Recent studies in plane channel flows found that even substituting Reynolds stresses with errors below 0.5 % from direct numerical simulation databases into RANS equations leads to velocities with large errors (up to 35 %). While such an observation may have only marginal relevance to traditional Reynolds stress models, it is disturbing for the recently emerging data-driven models that treat the Reynolds stress as an explicit source term in the RANS equations, as it suggests that the RANS equations with such models can be ill-conditioned. So far, a rigorous analysis of the condition of such models is still lacking. As such, in this work we propose a metric based on local condition number function for a priori evaluation of the conditioning of the RANS equations. We further show that the ill-conditioning cannot be explained by the global matrix condition number of the discretized RANS equations. Comprehensive numerical tests are performed on turbulent channel flows at various Reynolds numbers and additionally on two complex flows, i.e. flow over periodic hills, and flow in a square duct. Results suggest that the proposed metric can adequately explain observations in previous studies, i.e. deteriorated model conditioning with increasing Reynolds number and better conditioning of the implicit treatment of the Reynolds stress compared to the explicit treatment. This metric can play critical roles in the future development of data-driven turbulence models by enforcing the conditioning as a requirement on these models.

In this work, we propose using an ensemble Kalman method to learn a nonlinear eddy viscosity model, represented as a tensor basis neural network, from velocity data. Data-driven turbulence models have emerged as a promising alternative to traditional models for providing closure mapping from the mean velocities to Reynolds stresses. Most data-driven models in this category need full-field Reynolds stress data for training, which not only places stringent demand on the data generation but also makes the trained model ill-conditioned and lacks robustness. This difficulty can be alleviated by incorporating the Reynolds-averaged Navier–Stokes (RANS) solver in the training process. However, this would necessitate developing adjoint solvers of the RANS model, which requires extra effort in code development and maintenance. Given this difficulty, we present an ensemble Kalman method with an adaptive step size to train a neural-network-based turbulence model by using indirect observation data. To our knowledge, this is the first such attempt in turbulence modelling. The ensemble method is first verified on the flow in a square duct, where it correctly learns the underlying turbulence models from velocity data. Then the generalizability of the learned model is evaluated on a family of separated flows over periodic hills. It is demonstrated that the turbulence model learned in one flow can predict flows in similar configurations with varying slopes.

Characteristics of turbulence anisotropy in flow over two-dimensional rigid dunes are analysed. The Reynolds stress anisotropy is envisaged from the perspective of the stress ellipsoid shape. The spatial evolutions of the anisotropic invariant map (AIM), anisotropic invariant function, eigenvalues of the scaled Reynolds stress tensor and eccentricities of the stress ellipsoid are investigated at various streamwise distances along the vertical. The data plots reveal that the oblate spheroid axisymmetric turbulence appears near the top of the crest, whereas the prolate spheroid axisymmetric turbulence dominates near the free surface. At the dune trough, the axisymmetric contraction to the oblate spheroid diminishes, as the vertical distance below the crest increases. At the reattachment point and one-third of the stoss-side, the oblate spheroid axisymmetric turbulence formed below the crest appears to be more contracted, as the vertical distance increases. The AIMs suggest that the turbulence anisotropy up to edge of the boundary layer follows a looping pattern. As the streamwise distance increases, the turbulence anisotropy at the edge of the boundary layer approaches the plane-strain limit up to two-thirds of the stoss-side, intersecting the plane-strain limit at the top of the crest and thereafter moving towards the oblate spheroid axisymmetric turbulence.

The three-dimensional vortex clusters, and the structures based on the quadrant classification of the intense tangential Reynolds stress (Qs), are studied in direct numerical simulations of statistically stationary homogeneous shear turbulence (HST) at Taylor microscale Reynolds number $Re_{\\unicode[STIX]{x1D706}}\\approx 50{-}250$ , with emphasis on comparisons with turbulent channels (CHs). The Qs and vortex clusters in HST are found to be versions of the corresponding detached (in the sense of del Álamo et al. (J. Fluid Mech., vol. 561 (2006), pp. 329–358)) structures in CHs, although statistically symmetrised with respect to the substitution of sweeps by ejections and vice versa. In turn, these are more symmetric versions of the corresponding attached Qs and clusters. In both flows, only co-gradient sweeps and ejections larger than the local Corrsin scale are found to couple with the shear. They are oriented anisotropically, and are responsible for carrying most of the total Reynolds stress. Most large eddies in CHs are attached to the wall, but it is shown that this is probably a geometric consequence of their size, rather than the reason for their dynamical significance. Most small Q structures associated with different quadrants are far from each other in comparison to their size, but those that are close to each other tend to form quasi-streamwise trains of groups of a sweep and an ejection paired side by side in the spanwise direction, with a vortex cluster in between, generalising to three dimensions the corresponding arrangement of attached eddies in CHs. These pairs are organised around an inclined large-scale conditional vortex ‘roller’, and it is shown that the composite structure tends to be located at the interface between high- and low-velocity streaks, as well as in strong ‘co-gradient’ shear layers that separate streaks of either sign in which velocity is more uniform. It is further found that the conditional rollers are terminated by ‘hooks’ reminiscent of hairpins, both upright and inverted. The inverted hook weakens as the structures approach the wall, while the upright one changes little. At the same time, the inclination of the roller with respect to the mean velocity decreases from $45^{\\circ }$ in HST to quasi-streamwise for wall-attached eddies. Many of these observations are generalised to intense Reynolds stresses formed with different pairs of velocity components, and it is shown that most properties of the small structures can be traced to their definitions, rather than to their dynamics. It is concluded that the larger Reynolds-stress structures are associated with shear turbulence, rather than with the presence of a wall, while the smaller ones are generic to turbulence in general, whether sheared or not.

Direct numerical simulation (DNS) of a turbulent boundary layer over the Gaussian (Boeing) bump is performed. This boundary layer exhibits a series of adverse and favourable pressure gradients and convex and concave curvature effects before separating. These effects on turbulent boundary layers are characterised and compared with a lower-Reynolds-number flow over the same geometry. The momentum budgets are analysed in the streamline-aligned coordinate system upstream of the separation region. These momentum budgets allow the simplification of equations to facilitate an integral analysis. Integral-analysis-based approximations for Reynolds stresses in the inner and outer regions of the boundary layer are also formulated. The shear and wall-normal Reynolds stress profiles normalised by these approximations exhibit a better collapse compared with friction velocity and Zagarola–Smits normalisations in the strong favourable pressure gradient region and in the mild adverse pressure region that precedes it in this flow. Simplification of these Reynolds stress approximations along with results from the DNS are used to obtain semi-empirical approximations that are able to provide stress closure in terms of wall solution fields for the turbulent boundary layer under consideration.

Turbulent flows are investigated upstream of a bar rack system that is recommended as optimum in recent literature from tests with several fish species of different morphology, swimming ability, and behavior. Both two‐dimensional two‐component and two‐dimensional three‐component state‐of‐the‐art particle image velocimetry were used to quantify and analyze hydrodynamic metrics important for downstream migrating species. The inclination angles of the bar and rack were 45° and 30°, respectively, and the thickness of the bottom overlay was 13% of the water depth. The two Reynolds numbers investigated, based on incoming velocity and bar thickness, were 4,000 and 6,000. The statistical and structural characteristics of turbulent flows in the streamwise‐spanwise plane at 5% water depth, and the streamwise‐vertical plane at channel mid‐span are discussed. Upstream of the bottom overlay, the mean flow is deflected and accelerated toward the bypass, leading to an increase in the Reynolds stresses, while the turbulence eddies become smaller. For effective fish guidance, it is recommended that sweeping velocity (Vp) be larger than normal velocity (Vn), with Vp parallel and Vn perpendicular to the bar rack and bottom overlay. In the downstream half of the bar rack, Vn may increase sufficiently to surpass Vp near the bypass, possibly reducing effective guidance for some species and sizes. Upstream of the bars, the levels of streamwise mean velocity vary abruptly, which may deter fish from contacting the bars. Although inferences on passage effectiveness are made based on previous studies, tests with different species and sizes are needed to confirm fish responses. Key Points Turbulent flows upstream of a bar rack system with bottom overlay are characterized using particle image velocimetry Upstream, levels of the mean velocity and flow acceleration are significant and vary abruptly, which may deter fish from contacting the bars Flow conditions upstream of the bottom overlay may disorient fish and may affect guidance toward the bypass for safer downstream passage

Pore-resolved direct numerical simulations are performed to investigate the interactions between streamflow turbulence and groundwater flow through a randomly packed porous sediment bed for three permeability Reynolds numbers, $Re_K=2.56$, 5.17 and 8.94, representative of natural stream or river systems. Time–space averaging is used to quantify the Reynolds stress, form-induced stress, mean flow and shear penetration depths, and mixing length at the sediment–water interface (SWI). The mean flow and shear penetration depths increase with $Re_K$ and are found to be nonlinear functions of non-dimensional permeability. The peaks and significant values of the Reynolds stresses, form-induced stresses, and pressure variations are shown to occur in the top layer of the bed, which is also confirmed by conducting simulations of just the top layer as roughness elements over an impermeable wall. The probability distribution functions (p.d.f.s) of normalized local bed stress are found to collapse for all Reynolds numbers, and their root-mean-square fluctuations are assumed to follow logarithmic correlations. The fluctuations in local bed stress and resultant drag and lift forces on sediment grains are mainly a result of the top layer; their p.d.f.s are symmetric with heavy tails, and can be well represented by a non-Gaussian model fit. The bed stress statistics and the pressure data at the SWI potentially can be used in providing better boundary conditions in modelling of incipient motion and reach-scale transport in the hyporheic zone.

Two-phase bubbly flows are found in many industrial applications. These flows involve complex local phenomena that are still poorly understood. For instance, two-phase turbulence modelling is still commonly based on single-phase flow analyses. A direct numerical simulation (DNS) database is described here to improve the understanding of two-phase turbulent channel flow at a parietal Reynolds number of 127. Based on DNS results, a physical interpretation of the Reynolds stress and momentum budgets is proposed. First, surface tension is found to be the strongest force in the direction of migration so that budgets of the momentum equations suggest a significant impact of surface tension in the migration process, whereas most modelling used in industrial application does not include it. Besides, the suitability of the design of our cases to study the interaction between bubble-induced fluctuations (BIF) and single-phase turbulence (SPT) is shown. Budgets of the Reynolds stress transport equation computed from DNS reveal an interaction between SPT and BIF, revealing weaknesses in the classical way in which pseudoturbulence and perturbations to standard single-phase turbulence are modelled. An SPT reduction is shown due to changes in the diffusion because of the presence of bubbles. An increase of the redistribution leading to a more isotropic SPT has been observed as well. BIF is comprised of a turbulent (wake-induced turbulence, WIT) and a non-turbulent (wake-induced fluctuations, WIF) part which are statistically independent. WIF is related to averaged wake and potential flow, whereas WIT appears when wakes become unstable or interact with each other for high-velocity bubbles. In the present low gravity conditions, BIF is reduced to WIF only. A thorough analysis of the transport equations of the Reynolds stresses is performed in order to propose an algebraic closure for the WIF towards an innovative two-phase turbulence model.

This study presents an approach to investigate the role of eddy viscosity in linearized mean-field analysis of broadband turbulent flows. The procedure is based on spectral proper orthogonal decomposition (SPOD), resolvent analysis and the energy budget of coherent structures and is demonstrated using the example of a turbulent jet. The focus is on the coherent component of the Reynolds stresses, the nonlinear interaction term of the fluctuating velocity component in frequency space, which appears as an unknown in the derivation of the linearized Navier–Stokes equations and which is the quantity modelled by the Boussinesq approach. For the considered jet the coherent Reynolds stresses are found to have a mostly dissipative effect on the energy budget of the dominant coherent structures. Comparison of the energy budgets of SPOD and resolvent modes demonstrates that dissipation caused by nonlinear energy transfer must be explicitly considered within the linear operator to achieve satisfactory results with resolvent analysis. Non-modelled dissipation distorts the energy balance of the resolvent modes and is not, as often assumed, compensated for by the resolvent forcing vector. A comprehensive analysis, considering different predictive and data-driven eddy viscosities, demonstrates that the Boussinesq model is highly suitable for modelling the dissipation caused by nonlinear energy transfer for the considered flow. Suitable eddy viscosities are analysed with regard to their frequency, azimuthal wavenumber and spatial dependence. In conclusion, the energetic considerations reveal that the role of eddy viscosity is to ensure that the energy the structures receive from the mean-field is dissipated.

Trip-resolved large-eddy simulations of the DARPA SUBOFF are performed to investigate the development of turbulent boundary layers (TBLs) in model-scale studies. The primary consideration of the study is the extent to which the details of tripping affect statistics in large-eddy simulations of complex geometries, which are presently limited to moderate Reynolds number TBLs. Two trip wire configurations are considered, along with a simple numerical trip (wall-normal blowing), which serves as an exemplar of artificial computational tripping methods often used in practice. When the trip wire height exceeds the laminar boundary layer thickness, shedding from the trip wire initiates transition, and the near field is characterized by an elevation of the wall-normal Reynolds stress and a modification of the turbulence anisotropy and mean momentum balance. This trip wire also induces a large jump in the boundary layer thickness, which affects the way in which the TBL responds to the pressure gradients and streamwise curvature of the hull. The trip-induced turbulence decays along the edge of the TBL as a wake component that sits on top of the underlying TBL structure, which dictates the evolution of the momentum and displacement thicknesses. In contrast, for a trip wire height shorter than the laminar boundary layer thickness, transition is initiated at the reattachment point of the trip-induced recirculation bubble, and the artificial trip reasonably replicates the resolved trip wire behaviour relatively shortly downstream of the trip location. For each case, the inner layer collapses rapidly in terms of the mean profile, Reynolds stresses and mean momentum balance, which is followed by the collapse of the Reynolds stresses in coordinates normalized by the local momentum thickness, and finally against the 99 % thickness. By this point, the lasting impact of the trip is the offset in boundary layer thickness due to the trip itself, which becomes a diminishing fraction of the total boundary layer thickness as the TBL grows. The importance of tripping the model appendages is also highlighted due to their lower Reynolds numbers and susceptibility to laminar separations.