Direct Numerical Simulations (DNSs) of high Reynolds number turbulent flows, encountered in engineering, earth sciences, and astrophysics, are not tractable because of the curse of dimensionality associated with the number of degrees of freedom required to resolve all the dynamically significant spatio-temporal scales. Designing efficient and accurate Machine Learning (ML)-based reduced models of fluid

A numerical solution of the Navier–Stokes equations describing the instability of a two-dimensional subsonic flow about a flat plate is presented. It is obtained with the 16th-order multioperators-based scheme without introducing artificial exciters. It shows how small numerical disturbances produced by the scheme generate downstream from the leading edge a succession of wave packets formed by the

The constant growth of computational resources allows performing Large-Eddy Simulation (LES) on realistic flow configurations. In this context, it is important to properly model boundary conditions and particularly inflow turbulence. This study addresses wind tunnel applications where the object under investigation is downstream of a turbulence grid. This grid aims at generating a highly turbulent

The effects of Reynolds number on flow and mixing characteristics of a merged single swinging jet induced by a V-shaped fluidic oscillator were studied experimentally. The jet Reynolds number was varied from 80 to 3000. A high-speed digital camera was used to capture the instantaneous flow evolution processes using the laser light-sheet-assisted flow visualisation method. Long-exposure images together

The Discrete Element Method (DEM) for modelling of flow over rough walls is revised in the framework of the DANS (Double Averaged Navier-Stokes) equations, with a special focus on the drag term appearing in the mean momentum equation as a key to the robustness of the model. A set of 14 Direct Numerical Simulations (DNS) of channel flows with systematically varying roughness topographies is considered

We performed direct numerical simulations of a turbulent channel flow in viscoelastic fluids to assess various spatial scaling methods and to find a scaling method that can visualise the spatial features regardless of fluid properties. The Giesekus and Oldroyd-B models were used as constitutive equations to model the viscoelastic fluids alongside a Newtonian fluid. The initial scaling method that can

The effect of wall expansion on the structural and statistical characteristics of wall-shear stress (WSS) fluctuations was investigated by direct numerical simulations of a supersonic turbulent boundary layer over a sharp expansion corner with various deflection angles (β = 00, 20, 50 and 100). It is found that the two-dimensional fields of WSS are characterised as streamwise-elongated streaky structures

A closure for the effective relaxation time of the Boltzmann–BGK kinetic equation for fluid turbulence is presented, based on a double-averaging procedure over both kinetic and turbulent fluctuations. The resulting effective relaxation time appears to agree with values obtained via a renormalisation group treatment of the Navier–Stokes equation only at low values of k/T, the ratio of turbulent kinetic

In the present investigation, the LES WALE method along with the Zwart cavitation model is utilised to predict the cavitating flow around a Delft Twist-11 hydrofoil and the numerical results show a reasonable agreement with the experimental data. A novel elucidation is provided to the formation of the U-type structure and the Lumley method is introduced into the numerical simulation which can offer

Thermal convection is ubiquitous in nature as well as in many industrial applications. The identification of effective control strategies to, e.g. suppress or enhance the convective heat exchange under fixed external thermal gradients is an outstanding fundamental and technological issue. In this work, we explore a novel approach, based on a state-of-the-art Reinforcement Learning (RL) algorithm, which

Direct numerical simulations were performed of a free shear layer developing downstream from a no-slip splitter plate. The simulation results show that three-dimensional streaky perturbations developing in the boundary layer on the high-speed side of the splitter plate are convected downstream from the trailing edge and influence the development of three-dimensional perturbations in the planar free

We carry out Direct Numerical Simulation (DNS) of flows in closed rectangular ducts with several aspect ratios. The Navier–Stokes equations are discretised through a second-order finite difference scheme, with non-uniform grids in two directions. The duct cross-sectional area is maintained constant as well as the flow rate, which allows to investigate which is the appropriate length scale in the Reynolds

Temporally developing direct numerical simulations (T-DNS) are performed for free-stream turbulence induced bypass transition flow over a flat-plate under zero-pressure gradient, and results are validated using experimental data and spatially developing DNS (S-DNS) results. The temporal simulations predict the growth of near-wall Klebanoff modes in the pre-transition regime and their subsequent breakdown

The physical complexity and the large number of degrees of freedom that can be resolved today by direct numerical simulations of turbulent flows, and by the most sophisticated experimental techniques, require new strategies to reduce and analyse the data so generated, and to model the turbulent behaviour. We discuss a few concrete examples for which the turbulence data have been analysed by machine

We address the important point of the proportionality between the longitudinal integral lengthscale (L) and the characteristic mean flow width (δ) using experimental data of an axisymmetric wake and a turbulent planar jet. This is a fundamental hypothesis when deriving the self-similar scaling laws in free shear flow. We show that L/δ is indeed constant, at least in a range of streamwise distances

It is illustrated that a sharply truncated initial kinetic energy spectrum evolves to a staircase-shaped spectrum at short times. This effect is directly associated with the triadic nature of the energy transfer. A rigorous analysis leads to predictions on the time-dependence of this effect, and these predictions are verified by both direct numerical simulations and eddy-damped quasi-normal Markovian

Three-dimensional inert and reactive shock-interface interactions are simulated by solving NS equation with ninth-order WENO scheme. An initial planar shock wave with Mach number 1.7 accelerates a single mode and finite thickness interface with initial wavelength 4 mm and amplitude 0.4 mm from the reactant composed of C2H4+3O2+4N2. Then the transmitting shock branch is reflected at the end boundary

Direct Numerical Simulations (DNS) of forced homogeneous isotropic turbulence in a dense gas (FC-70), accurately described by a complex EoS, are computed for a turbulent Mach number of 0.8. In a numerical experiment, results are compared to the ones obtained when considering the fluid as a perfect gas. It is found that the dense gas displays a deeply modified shocklets' structure. The amplitude of

In a double-pipe heat exchanger, heat transfer enhancement can be carried out by active or passive methods. In the present study, heat transfer characteristics are studied by the addition of hemispherical turbulators placed in the annulus region in a double-pipe heat exchanger. The turbulators are placed in the annulus region on the inner pipe separated at 90° to each other around the circumference

Uniform blowing in wall bounded shear flows is well known for its drag reducing effects and has long been investigated ever since. However, many contemporary and former research on this topic has confirmed the drag reducing effect but very less is known regarding how blowing is effecting the Reynolds stresses at high Reynolds number. Therefore, effect of uniform blowing has been experimentally investigated

This paper explores how far the scientific discovery process can be automated. Using the identification of causally significant flow structures in two-dimensional turbulence as an example, it probes how far the usual procedure of planning experiments to test hypotheses can be substituted by ‘blind’ randomised experiments and notes that the increased efficiency of computers is beginning to make such

There is a rapidly growing interest in using general-purpose CFD codes based on second-order finite volume methods for Large-Eddy Simulation (LES) in a wide range of applications, and in many cases involving wall-bounded flows. However, such codes are strongly affected by numerical dissipation and the accuracy obtained for typical LES resolutions is often poor. In the present study, we approach the

Turbulence structures in flow over three types of wall roughness: sand-grain, cube roughness and a realistic, multi-scale turbine-blade roughness, are compared to structures observed in flow over a smooth wall in open channel flow at Reτ=1000, using direct numerical simulations. Two-point velocity correlations, length scales, inclination angles, and velocity spectra are analysed, and the applicability

Reynolds-averaged Navier-Stokes (RANS) equations are presently one of the most popular models for simulating turbulence. Performing RANS simulation requires additional modelling for the anisotropic Reynolds stress tensor, but traditional Reynolds stress closure models lead to only partially reliable predictions. Recently, data-driven turbulence models for the Reynolds anisotropy tensor involving novel

Transpired boundary layers are of major interest for many industrial applications. Although well described, there is no turbulence model specifically dedicated to the prediction of boundary layers for both blowing and suction configurations. Revisiting closure relations of turbulence models, a general strategy was established to recover Stevenson's law of the wall that described the behaviour of transpired

Scalar iso-surfaces have been used to describe many interface problems in turbulence, such as flame surfaces in turbulent reacting flows and turbulent/non-turbulent interfaces. In this paper we report on direct numerical simulations of the behaviour of iso-surfaces of two differently initialised conserved scalars evolving in isotropic turbulence. The terms in the equation for the iso-surface area density

Spalart-Allmaras (SA) model is a low Reynolds number (Re) model, which means that the first off-wall grid point should be placed in the viscous sub-layer with y+≃1. This restriction of placing the first off-wall grid point so close to the wall leads to an increase in the mesh size, and thus the computation. The wall function approach is an alternative to this problem. The standard wall function method

A novel numerical method to perform direct numerical simulations (DNS) of turbulent flow was established to predict the characteristics of compressible gas–liquid two-phase flow. A cavitation model was introduced to reproduce the phase change in the cavitation phenomena, and a DNS of wall turbulence of water with cavitation was performed. Vortex cavitation was observed in the wall turbulence. The growing

Adverse-pressure-gradient turbulent boundary layer flow was inspected at Reynolds number based on momentum thickness, Reθ≈1200, using particle image velocimetry in a refractive-index-matching flume. Proper orthogonal decomposition was used to quantify the effect of large-scale motions on the Reynolds stresses at the onset of separation and within the separated flow. Results show that approximately

We investigate the properties of the velocity gradient tensor for spatially evolving turbulent flows (a near-wake, two axisymmetric jets and a planar mixing layer). Emphasis is placed on the study of the normal and non-normal parts of the tensor. Non-normality plays a greater role in the dynamics than is observed for HIT and does so for all spatial locations examined. This implies a greater role for

This paper investigated the integral length scales of turbulence in a low-roughness atmospheric surface layer (ASL), characterised by very smooth terrain in the Utah desert during near-neutral conditions, and evaluated the Engineering Sciences Data Unit (ESDU) 85020 and 86010 predictions for the turbulence length scales in a low-roughness ASL. The correlation integral method was used to estimate the

The incompressible temporally developing turbulent boundary layer (TBL) is analysed using the map-based stochastic one-dimensional turbulence (ODT) model. The TBL is a canonical flow problem, which is, in the present study, formed by a planar moving wall and a free stream at rest. An understanding of this idealised flow is of fundamental relevance for the numerical analysis of turbulent boundary-layer-type

In the present investigation, the evolution of stratified turbulence submitted to a non-vertical shear is studied using second-order models. Different cases of turbulent flows are considered. Firstly, the case of a purely horizontal shear is considered ( θ=π/2). In the second case, we study a purely vertical shear θ=0. Finally, we studied an inclined shear corresponding to θ=π/6,θ=π/4,θ=π/3,π/8,3π/8

The effects of non-uniformity of flows and initial disturbance intensity of the incident shock wave on Richtmyer–Meshkov instability (RMI) when a sinusoidal shock wave with Ma = 1.25 impinging an unperturbed interface are numerically investigated. The interface morphology, turbulent mixing zone width (TMZW), Y-integrated vorticity, circulation, and turbulent kinetic energy (TKE) are qualitatively and

Numerical studies into the dynamics of non-spherical particles in turbulent channel flow have, until now, been mostly confined to low Reynolds numbers. In this paper, we investigate the dynamics of tracer non-spherical particles in a channel flow at Reτ≈1000 for the first time. Tracer spheroidal particles suspended in a turbulent channel flow are computed by means of direct numerical simulation coupled

An efficient algorithm for computing the terms appearing in the Generalised Kolmogorov Equation (GKE) written for the indefinite plane channel flow is presented. The algorithm, which features three distinct strategies for parallel computing, is designed such that CPU and memory requirements are kept to a minimum, so that high-Re wall-bounded flows can be afforded. Computational efficiency is mainly

We develop a theory for the cascade mixing terms in a moment closure approach to binary active scalar mixing in variable-density turbulence. To address the variable-density complications we apply, as a principle and constraint, the conservation of the probability density function (PDF) through a Fokker–Planck equation with bounded sample space whose attractor is the beta PDF with skewness. Mixing is

Recently Brouwers [Dissipation equals production in the log layer of wall-induced turbulence. Phys Fluids. 2007;19:101702] carried out an asymptotic analysis using the RANS based turbulence energy transport equation and showed that the energy dissipation equals its production in the inertial layer of wall-induced turbulence. Assuming log-law profile to the mean velocity, pressure, viscous and energy

In wall-bounded time-periodic flows, nonlinearity, associated with higher harmonic term(s) in velocity and/or acceleration outside the boundary layer, can significantly change the wall turbulence compared with that in the linear Stokes Boundary Layer. A significant feature of a nonlinear wall-bounded turbulent time-periodic flow is the formation of a net current which has not yet been mechanistically

The partially averaged Navier–Stokes (PANS) model, proposed in Girimaji [Partially-averaged Navier-Stokes model for turbulence: a Reynolds-averaged Navier-Stokes to direct numerical simulation bridging method. ASME J Appl Mech. 2006;73(3):413–421], can be used to simulate turbulent flows either as a RANS, LES or DNS. The PANS model includes fk which denotes the ratio of modelled to total kinetic energy

The flow characteristics of the trailing edge of vertical vanes installed at the intersection of a T-junction duct were experimentally investigated using particle image velocimetry. The measured velocity field in the branch duct with/without single circular cylinder was studied under different cross velocities and velocity ratios. Additionally, the effect of the locations of cylinder on the flow field

The topological evolution near the turbulent/non-turbulent interface (TNTI) in turbulent mixing layer is studied by means of statistical analysis of the invariants of velocity gradient tensor (VGT) based on direct numerical simulation data. The dynamics of topological evolution is investigated in terms of the source terms of the evolution equations for the invariants, including the pressure effect

In light of recent data from hot-wire anemometry and laser Doppler velocimetry, this article explores experimentally the momentum balance and kinetic energy production in fully developed turbulent pipe flow for shear Reynolds numbers in the range 8⋅102≤R+≤4⋅104 from two pipe facilities. It has become common practice to indirectly deduce the Reynolds shear stress via the mean flow data and the mean-momentum

We use direct numerical simulations to study the dynamics of incompressible homogeneous turbulence subjected to a uniform magnetic field B in a rotating frame with rotation vector Ω. We consider two cases: Ω∥B and Ω⊥B. The initial state is homogeneous isotropic hydrodynamic turbulence with Reynolds number Re=ull/ν≃170. The magnetic Prandtl number Pm=ν/η=1 and the Elsasser number Λ=B2/(2Ωη)=0.5, 0.9

In the present work, the use of cylindrical turbulators in a double pipe heat exchanger has been investigated. Cylindrical fin type of turbulators has been placed circumferentially separated by 90° on the outer side of an inner pipe at a regular pitch. Experimental studies were undertaken for different air flow rates in a turbulent regime whose Reynolds number range between 2500 and 10000. Heat transfer

This paper investigates the applicability of a hybrid data assimilation approach, namely ensemble-based variational method (EnVar), to optimise Reynolds Averaged Navier-Stokes (RANS) simulations in convergent-divergent channel from the perspective of Bayesian inference. Concretely, the ensemble-based variational method is applied to infer the inlet velocity and turbulence model corrections by assimilating

Particle Image Velocimetry, PIV, is employed to characterise the near-homogeneous, isotropic, turbulence generated inside a large spherical vessel by four rotating fans. Spatial and temporal distributions of mean and root mean square, rms, velocity fluctuations are investigated, as well as integral length scales, L, Taylor microscales, λ, and Kolmogorov length scales, η, in the fan speed range, 1000–6000 rpm

Multi-scale statistics are used to analyse the flow structure of wake flow in the boundary layer of a wind turbine array. Experimentally, a 3×3 wind turbine array is tested with X-type hot-wire anemometry, providing a velocity signal at discrete locations downstream of the array along the centreline of the centre turbine. Based on the Markov property, the turbulent cascade can be taken as a stochastic

Jets and plumes are shear flows created by momentum and buoyant sources. They can be classified as either pure jets, forced plumes, pure plumes, or lazy plumes. Lazy plumes are characterised as having a high source buoyancy flux relative to the momentum flux. We use direct numerical simulations to simulate lazy plumes to ascertain the effects of increasing plume ‘laziness’ on the higher order moments

Understanding of the structure of turbulent flows at extreme Reynolds numbers (Re) is relevant because of several reasons: almost all turbulence theories are only valid in the high Re limit, and most turbulent flows of practical relevance are characterized by very high Re. Specific questions about wall-bounded turbulent flows at extreme Re concern the asymptotic laws of the mean velocity and turbulence

The structural attributes of turbulent flow over a complex roughness topography are explored using high-frame-rate stereo particle-image velocimetry measurements in the wall-normal–spanwise plane. The roughness under consideration was replicated from a turbine blade damaged by deposition of foreign materials and contains a broad range of topographical scales arranged in a highly irregular manner. Previous

(2019). Letter. Journal of Turbulence: Vol. 20, Special Issue on Wind Energy, pp. 1-1.

(2019). Correction. Journal of Turbulence: Vol. 20, No. 2, pp. i-i.

The accelerating growth of wind energy in recent years mandates improved understanding of wind turbine, wind farm and atmospheric turbulence interactions. Fluid turbulence plays a vital role in these interactions, motivating the present formulation of several pertinent questions for turbulence research. These questions touch upon the need for better analytical, synthetic and reduced order models of

(2019). Review of the monograph by Pierre Sagaut and Claude Cambon entitled Homogeneous Turbulence Dynamics. Journal of Turbulence: Vol. 20, No. 3, pp. 240-244.

An approach to devising a consistency formulation for Pk/εPk/ε (production-to-dissipation ratio) is proposed to obtain a non-singular CμCμ (coefficient of eddy-viscosity) embedded in the one-equation model based on the turbulent kinetic energy k. The dissipation rate ε is evaluated with an algebraically prescribed length scale having only one adjustable coefficient, accompanied by an anisotropic function

Wind energy has been growing steadily in the U.S. and worldwide in the past decades. As wind farms are projected to increase in size and number, however, concerns are rising about possible undesirable effects of wind turbines near the Earth's surface. The literature is highly divided about what these effects could be, including warming, cooling, both, or neither. Only one mechanism, however, has been

Large-eddy simulation (LES) has been used to model turbulent channel flow over urban-like, fractal topographies, constructed via iterated function system (IFS). By using the IFS approach, the topography fractal dimension, D, and a desired number of generations, Ng, can be nominated a priori. The topographies considered herein all featured the same central square-based block for generation one, while

The universality and mathematical physical structure of wall-bounded turbulent flows is a topic of discussions over many decades. There is no agreement about questions like what is the physical mean flow structure, how universal is it, and how universal are theoretical concepts for local and global flow variations. These questions are addressed by using latest direct numerical simulation (DNS) data

A description of the turbulence kinetic energy budget is investigated to understand the dynamics and turbulence energy transfer between the atmospheric boundary layer and large wind farms. The turbulence kinetic energy budget terms are a dominant factor in explaining the energy supply process in wind farms. Large eddy simulations are used to generate the thermally stratified wind turbine array boundary