Direct numerical simulation of turbulent duct flow with inclined or V-shaped ribs mounted on one wall

In this research, highly disturbed turbulent flow of distinct three-dimensional characteristics in a square duct with inclined or V-shaped ribs mounted on one wall is investigated using direct numerical simulation. The turbulence field is highly sensitive to not only the rib geometry but also the boundary layers developed over the side and top walls. In a cross-stream plane secondary flows appear as large longitudinal vortices in both inclined and V-shaped rib cases due to the confinement of four sidewalls of the square duct. However, owing to the difference in the pattern of cross-stream secondary flow motions, the flow physics is significantly different in these two ribbed duct cases. It is observed that the mean flow structures in the cross-stream directions are asymmetrical in the inclined rib case but symmetrical in the V-shaped rib case, causing substantial differences in the momentum transfer across the spanwise direction. The impacts of rib geometry on near-wall turbulence structures are investigated using vortex identifiers, joint probability density functions between the streamwise and vertical velocity fluctuations, statistical moments of different orders, spatial two-point autocorrelations and velocity spectra. It is found that near the leeward and windward rib faces, the mean inclination angle of turbulence structures in the V-shaped rib case is greater than that of the inclined rib case, which subsequently enhances momentum transport between the ribbed bottom wall and the smooth top wall.

The characteristics of weather-related wind load parameters and wind 2 power quality over the Dabanchen strong gale area

Analysis of data from cup and ultrasonic anemometers on a 100 m-tall wind mast in the Dabanchen Canyon reveals that the turbulence intensities, gust factors, and peak factors measured by cup anemometers tend to be severely undervalued compared to longitudinal values from an ultrasonic anemometer, and onsite three-dimensional measurement data are preferred for weather-related wind load calculation. This difference is related to the rotating responses of cup anemometers during wind speed acceleration or deceleration and the higher vertical fluctuation speeds driven by dynamic interactions with the sloped canyon terrain. The higher lateral turbulence is key consideration for determining wind turbine classes in Danbanchen strong wind area. The longitudinal gust and peak factors under wind speeds exceeding 25.0 m s−1 are comparable with those of typhoon-prone open, flat regions. The chances of high turbulence and unstable stratification convective processes are very high; thus, wind power collection efficiencies are not high. Except the longitudinal turbulence integral scale, China’s wind-resistance codes are not applicable due to occasional undesirable strong gales with extraordinary turbulence structures. The measured vertical power spectral density of fluctuating wind in the high-frequency domain cannot reflect the rapidly adapting features of the vast terrain under strong gales.

The characteristics of weather-related wind load parameters and wind power quality over the Dabanchen strong gale area

Analysis of data from cup and ultrasonic anemometers on a 100 m-tall wind mast in the Dabanchen Canyon reveals that the turbulence intensities, gust factors, and peak factors measured by cup anemometers tend to be severely undervalued compared to longitudinal values from an ultrasonic anemometer, and onsite three-dimensional measurement data are preferred for weather-related wind load calculation. This difference is related to the rotating responses of cup anemometers during wind speed acceleration or deceleration and the higher vertical fluctuation speeds driven by dynamic interactions with the sloped canyon terrain. The higher lateral turbulence is key consideration for determining wind turbine classes in Danbanchen strong wind area. The longitudinal gust and peak factors under wind speeds exceeding 25.0 m s−1 are comparable with those of typhoon-prone open, flat regions. The chances of high turbulence and unstable stratification convective processes are very high; thus, wind power collection efficiencies are not high. Except the longitudinal turbulence integral scale, China’s wind-resistance codes are not applicable due to occasional undesirable strong gales with extraordinary turbulence structures. The measured vertical power spectral density of fluctuating wind in the high-frequency domain cannot reflect the rapidly adapting features of the vast terrain under strong gales.

Simulation of the Scalar Transport above and within the Amazon Forest Canopy

The parallelized large-eddy simulation model (PALM) was used to understand better the turbulent exchanges of a passive scalar above and within a forested region located in the central Amazon. Weak (2 ms−1) and strong (6 ms−1) wind conditions were simulated. A passive scalar source was introduced to the forest floor for both simulations. The simulations reproduced the main characteristics of the turbulent flow and of the passive scalar transport between the forest and the atmosphere. Noteworthily, strong and weak wind conditions presented different turbulence structures that drove different patterns of scalar exchange both within and above the forest. These results show how passive scalar concentration is influenced by the wind speed at the canopy top. Additionally, higher wind speeds are related to stronger sweep and ejection regimes, generating more intense plumes that are able to reduce the passive scalar concentration inside the forest canopy. This work was the first that used PALM to investigate scalar transport between the Amazon rainforest and the atmosphere.

Turbulent flow in curved channels

Fully developed turbulent flow in channels with mild to strong longitudinal curvature is studied by direct numerical simulations. The Reynolds based on the bulk mean velocity and channel half-width $\delta$ is fixed at $20\,000$ , resulting in a friction Reynolds number of approximately 1000. Four cases are considered with curvature varying from $\gamma = 2\delta /r_c = 0.033$ to 0.333, where $r_c$ is the curvature radius at the channel centre. Substantial differences between the mean wall shear stress on the convex and concave walls are already observed for $\gamma = 0.033$ . A log-law region is absent and a region with nearly constant mean angular momentum develops in the channel centre for strong curvatures. Spanwise and wall-normal velocity fluctuations are strongly amplified by curvature in the outer region of the concave channel side. Only near the walls, where curvature effects are relatively weak, do the mean velocity and velocity fluctuation profiles approximately collapse when scaled by wall units based on the local friction velocity. Budgets of the streamwise and wall-normal Reynolds-stress equations are presented and turbulence structures are investigated through visualizations and spectra. In the case with strongest curvature, the flow relaminarizes locally near the convex wall. On the concave channel side, large elongated streamwise vortices reminiscent of Taylor–Görtler vortices develop for all curvatures considered. The maximum in the premultiplied two-dimensional wall-normal energy spectrum and co-spectrum shifts towards larger scales with increasing curvature. The large scales substantially contribute to the wall-normal velocity fluctuations and momentum transport on the concave channel side.

Stability analysis of open-channel flows with secondary currents

This paper presents some results coming from a linear stability analysis of turbulent depth-averaged open-channel flows (OCFs) with secondary currents. The aim was to identify plausible mechanisms underpinning the formation of large-scale turbulence structures, which are commonly referred to as large-scale motions (LSMs) and very-large-scale motions (VLSMs). Results indicate that the investigated flows are subjected to a sinuous instability whose longitudinal wavelength compares very well with that pertaining to LSMs. In contrast, no unstable modes at wavelengths comparable to those associated with VLSMs could be found. This suggests that VLSMs in OCFs are triggered by nonlinear mechanisms to which the present analysis is obviously blind. We demonstrate that the existence of the sinuous instability requires two necessary conditions: (i) the circulation of the secondary currents $\omega$ must be greater than a critical value $\omega _c$ ; (ii) the presence of a dynamically responding free surface (i.e. when the free surface is modelled as a frictionless flat surface, no instabilities are detected). The present paper draws some ideas from the work by Cossu, Hwang and co-workers on other wall flows (i.e. turbulent boundary layers, pipe, channel and Couette flows) and somewhat supports their idea that LSMs and VLSMs might be governed by an outer-layer cycle also in OCFs. However, the presence of steady secondary flows makes the procedure adopted herein much simpler than that used by these authors.

The Four Stage Development of Starting Turbulent Buoyant Plumes

Turbulent heated and buoyant plumes have important applications in the atmosphere such as wildland fire plumes, volcanic plumes, and chemical plumes. The purpose of the study is to analyze the turbulence structures, and to understand the stages of the development of the starting turbulent plumes. For this purpose, data generated from an in-house Weather Research Forecast model coupled with Large-eddy simulation (WRF-bLES) with two-way feedback between the buoyant plume and the atmosphere developed has been used. The release of both dense gases (Co2, So2) and, buoyant gases (He, NH3, heated air) from a circular source at the bottom of the domain have been investigated. The simulations of the axisymmetric plume were performed at a high Reynolds number of 108. Vortex Identification methods were used to extract the Coherent structures and the large-scale features of the flow. The results have demonstrated that both the dense and the buoyant heated plumes with different initial characters exhibited universal characteristics and the development of the starting plumes occurred in four characteristic stages: Stage 1 is the plume acceleration stage, followed by stage 2 which corresponds to the formation of the head of the plume which grows spatially. Stage 3 is when the plume head is fully formed and the flow transitions to quasi-steady-state behavior. The final stage is the fully developed plume. The identification of the four-stage development of the plume in the neutral environment is the first step in studying the turbulent heated and buoyant plumes development in order to characterize realistic plumes and to quantify the extent of mixing at each of these stages. This work has important contributions to fundamental fluid dynamics of buoyant plumes with implications on forecasting the plume trajectory of smoke, wildland fire, and volcanic plumes.

Velocity Field and Turbulence Structure around Spur Dikes with Different Angles of Orientation under Ice Covered Flow Conditions

Spur dikes are well-known structures that are widely used in rivers and coastal regions. Depending on their types, sizes, and orientation angles, spur dikes can substantially change flow characteristics. Results of previous studies indicate that the presence of an ice cover in rivers can cause complicated flow structures. The present experimental study investigates velocity fields and turbulence structures in the vicinity of spur dikes under ice cover with different roughness coefficients. The spur dikes were set up at the following three angles of orientation, 90°, 60°, and 45°. Our results show that the strongest velocity fluctuation occurs immediately above the scour hole surface and very close to the dike tip. The increase in the dike angle toward upstream, the velocity component values increase, leads to a larger scour hole. Results show that an increase in dike angle of each 10° (from 45° to 90°) increases the scour depth between 5% and 10%, depending on flow conditions. Furthermore, the increase in the cover roughness coefficient and the blockage ratio of a spur dike leads to a further increase in turbulence kinetic energy and 3D velocity components values. The findings of this study imply that the appearance of an ice cover can increase turbulence intensities up to nearly 30%.