In purely axisymmetric turbulence sustained by a linear forcing mechanism, a stable, entirely poloidal flow is observed when the toroidal component of the forcing is below a threshold. We investigate using numerical simulations whether this state persists when toroidal variations of the flow are continuously reintroduced and the forcing is purely poloidal. It is shown how the pressure–strain correlation

(2021). Special issue on rotating turbulence. Journal of Turbulence: Vol. 22, Rotating Turbulence, pp. 231-231.

Turbulence is characterised by non-linear transfer of energy and other invariants across an extended range of scales, from the injection scale to where dissipation is active. Many flow configurations in nature and in the laboratories deviate from the ideal homogeneous and isotropic turbulent cases (HIT) characterised by the presence of direct energy cascade only and strong small-scales universality

Coriolis-centrifugal convection ( C3) in a cylindrical domain constitutes an idealised model of tornadic storms, where the rotating cylinder represents the mesocyclone of a supercell thunderstorm. We present a suite of C3 direct numerical simulations, analysing the influence of centrifugal buoyancy on the formation of tornado-like vortices (TLVs). TLVs are self-consistently generated provided the flow

Non-premixed impinging jet flames with different coflow conditions are performed using PIV technology combined with numerical simulation to investigate flame instability in the vicinity of wall. Results indicate that the increase of coflow velocity results in a more chaotic flow field and higher fuel efficiency, and the increase of coflow temperature leads to ignition advance and the increase of NO

Spontaneous and stochastic reversal of large scale flow structure is an intriguing and crucial phenomenon in turbulent Rayleigh-Bénard type natural convection. This paper proposes a new control approach to eliminate the reversals through stabilising the corner flows using two small sidewall controllers. Based on a series of direct numerical simulations, it is shown that the control can successfully

Direct numerical simulations of turbulent boundary layers with roughness elements on a wall were performed to investigate the spatial characteristics of rough-wall turbulence and establish a corresponding prediction method. When the roughness height was larger than the buffer layer, the rough-wall turbulence exhibited different spatial characteristics of the turbulence structures from those pertaining

We investigate the enstrophy dynamics in relation to objective Eulerian coherent structures (OECSs) and their impact on the enstrophy and scalar transport near the turbulent/non-turbulent interface (TNTI) in flows with and without stable stratification. We confirm that vortex-stretching produces enstrophy inside the boundaries of the OECSs, while viscous diffusion transfers the enstrophy across the

In this study, a RANS model of turbulent conjugate heat transfer has been developed, which is applicable across a range of different combination of fluid and solid thermal properties. This is achieved by focusing on the transport equations for the temperature variance and its dissipation rate across the solid walls which bound the flow region. In this investigation we make use of a wider range of DNS

An experimental investigation of near-wall scale interactions in the presence of a deterministic forcing input is presented in this work. The external forcing input was generated by a wall-mounted piezoelectric (PZT) actuator, which directly introduces a dynamic perturbation into the near-wall cycle of turbulent boundary layer flow. The spectra of velocity fluctuations indicated that the fundamental

ABSTRACT This paper explores the use of a small-span direct numerical simulation for a transient, smooth-wall turbulent channel flow and then applies the small-span simulation to a transient channel flow with riblets. A flow configuration similar to that of S. He and M. Seddighi (J Fluid Mech. 2013;715:60–102) is used to study the impulse response of a half-height channel flow to an abrupt increase

ABSTRACT Hydrocyclones are widely used in industry and CFD has been used to compute them. Reynolds stress turbulence models (RSM), which are computationally costly and oftentimes hard to converge, are often recommended in these computations. The present work has selected a liquid-liquid separation hydrocyclone for which single-phase experimental tangential and axial velocity profiles are available

ABSTRACT The performance of ultra-high-lift (UHL) low-pressure turbine (LPT) is subject to complex flow phenomena (e.g. separation, transition and reattachment) which require advanced modelling for accurate numerical predictions. The feasibility and fidelity of three widely used transition-based turbulence models are evaluated in the Reynolds-Averaged Navier-Stokes (RANS) prediction of low-Reynolds

Rotating Rayleigh–Bénard convection is a simple model system used to study the interplay of buoyant forcing and rotation. Many recent studies have focused on the geostrophic regime of turbulent rotating convection where the principal balance of forces is between the Coriolis force and the pressure gradient. This regime is believed to be representative of conditions in geophysical and astrophysical

Aero-optical reduction in flat-plate turbulent boundary layer has been widely concerned due to critical application for high-speed target-seeking vehicles. From the perspective of disciplinary intersection between aerodynamic and optical engineering, the turbulence-aberrated optical behaviour, in high-Reynolds-number supersonic freestream, is numerically investigated using the in-house code High-Order

Heated horizontal plane jets find wide applications in engineering appliances such as air curtains and discharge of industrial effluents. In the present study, experimental investigations are conducted on a heated horizontal plane jet with the Reynolds numbers in the transitional regime, using a hotwire anemometer. In the far to very far-field (20 < x/d < 100) centreline velocity decay and jet spread

Asymptotic methods have been used to study a turbulent flat mixing layer comprising a viscous compressible gas flow descending from the trailing edge of a body and a gas at rest having the same chemical composition but a different temperature. In contrast to the existing methods for analysing turbulent flows using the Reynolds-averaged Navier-Stokes (RANS) equations, with some closure hypotheses, this

ABSTRACT The present numerical study aims to numerically investigate the dynamic and turbulent characteristics of a two-dimensional and turbulent offset jet. Three different parameters were investigated: The Reynolds numbers (Re) which was varied from 10000 to 30000, wall inclination angle (α) that was from −20° to +20° and finally the offset ratio (OR) which extends from 3.25–13. Ansys Fluent was

Hysteresis behaviour was reported in spanwise rotating plane Couette flow (RPCF) at Reynolds number R e w = U w h / ν = 1300 with varying rotation number R o = 2 Ω z h / U w in a recent work (Huang et al. Phys. Rev. Fluids 2019;4:052401(R)). Here, U w is half of the velocity difference between two walls, h is half of the channel width, ν is the kinematic viscosity and Ω z is the constant angular velocity

Hysteresis behaviour was reported in spanwise rotating plane Couette flow (RPCF) at Reynolds number Rew=Uwh/ν=1300 with varying rotation number Ro=2Ωzh/Uw in a recent work (Huang et al. Phys. Rev. Fluids 2019;4:052401(R)). Here, Uw is half of the velocity difference between two walls, h is half of the channel width, ν is the kinematic viscosity and Ωz is the constant angular velocity in the spanwise

The magnetohydrodynamic (MHD) control has great potential in the applications of hypersonic vehicles. Typically, the MHD flowfield in these applications is low magnetic Reynolds (Rem ) compressible turbulent flow, which is different from that without magnetic field and within the high Rem range. This paper investigated the low Rem compressible MHD isotropic turbulence in different Taylor-scale Reynolds

This paper investigates the effects of co-flow jet on the aerodynamic performance of some symmetric blade sections of wind turbines. The Reynolds number of 1.3 × 105, angles of attack between 0° and 18° and momentum coefficients of 0.03, 0.05, and 0.08 are considered for four airfoils of NACA 0012, 0015, 0018, and 0021 as the blade sections. To numerically simulate the fluid flow, the Navier-Stokes

There is a substantial need to better understand multi-physics interactions in realistic laboratory experiments, yet it is currently unfeasible to perform direct numerical simulations (DNS) of most such problems. As a result, reduced-order models such as those used to represent unclosed subgrid-scale stresses in large-eddy simulations (LES) and advanced numerical techniques, including adaptive mesh

We report of a series of laboratory experiments and numerical simulations of freely decaying rotating turbulent flows confined in domains with variable height. We show that the vertical confinement has important effects on the formation of large-scale columnar vortices, the hallmark of rotating turbulence, and in particular delays the development of the cyclone–anticyclone asymmetry. We compare the

The local flow field and wake region over a forward-facing obstacle were experimentally investigated, when the water surface above the obstacle was dipped compared to flow further upstream and further downstream. The obstacle was modelled resembling the shape of an airfoil with a flat top having stoss-side slope 35∘ and downstream lee-side slope 6∘ designed for minimal flow separation. The submerged

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

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

ABSTRACT Motivated by the needs of wall modelled Large Eddy Simulation (LES), we introduce fits to numerical solutions of the Reynolds Averaged Navier–Stokes equations in their simplest near-wall, boundary layer approximation including a mixing-length model. We formulate the problem such that independent dimensionless variables are those directly available in LES. We provide practical fits for the

(2020). Editorial. Journal of Turbulence: Vol. 21, The machine learning: Guest edited by Federico Toschi, pp. 483-483.

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

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

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