Flow structure over square bars at intermediate submergence : Large Eddy Simulation study of bar spacing effect

Peer reviewedPublisher PD

Impact of environmental turbulence on the performance and loadings of a tidal stream turbine

A large-eddy simulation (LES) of a laboratory-scale horizontal axis tidal stream turbine operating over an irregular bathymetry in the form of dunes is performed. The Reynolds number based on the approach velocity and the chord length of the turbine blades is approximately 60,000. The simulated turbine is a 1:30 scale model of a full-scale prototype and both turbines operate at very similar tip-speed ratio of λ ≈ 3. The simulations provide quantitative evidence of the effect of seabed-induced turbulence on the instantaneous performance and structural loadings of the turbine revealing how large-scale, energetic turbulence structures affect turbine performance and bending moments of the rotor blades. The data analysis shows that wake recovery is notably enhanced in comparison to the same turbine operating above a flat-bed and this is due to the higher turbulence levels generated by the dune. The results demonstrate the need for studying in detail the flow and turbulence characteristics at potential tidal turbine deployment sites and to incorporate observed large-scale velocity and pressure fluctuations into the structural design of the turbines

Large-eddy simulation of free-surface turbulent channel flow over square bars

The results of large-eddy simulations of free-surface turbulent channel flow over spanwise-aligned square bars are used to investigate the effects of bed roughness and water surface deformations on the root-mean-squared velocity fluctuations, dispersive shear stress, double-averaged Reynolds shear stress, wake kinetic energy and double-averaged turbulent kinetic energy. Two bar spacings, corresponding to transitional and k-type roughness, at similar Reynolds and Froude numbers are considered. The main peak of all statistical quantities occurs at the bar crest height. The effects of a standing wave at the water surface in flow over k-type roughness is marked by a local peak under the water surface for all statistical quantities considered here except wall-normal and spanwise velocity fluctuations. Quadrant analysis shows that sweeps and ejections are the strongest events contributing to both dispersive and double-averaged Reynolds shear stress but their contributions are different for the two bar spacings. Examining the budgets of dispersive and double-averaged Reynolds shear stress reveals that the dominant terms of these stresses are pressure–strain correlation and pressure transport and the contribution of wake production is similar for both of these stresses but with opposite sign. In addition to the main role of the bars in consuming or producing wake kinetic energy through production and transport and convection, the standing wave at the water surface in flow over k-type roughness induces large convection in the bulk flow too. The dominant terms in the double-averaged turbulent kinetic energy budget are similarly production, transport and convection. Large shear production renders large temporal fluctuations than spatial fluctuations of flow variables. The interaction of bars, bed and water surface is seen in the convection term in flow over k-type roughness

Reynolds and dispersive shear stress in free-surface turbulent channel flow over square bars

Reynolds and dispersive shear stresses in turbulent flow over spanwise-aligned square bars in an open channel flow are examined. Results of large-eddy simulation of flow over two different bar spacings corresponding to transitional and k -type (reattaching flow) roughness are analyzed. The Reynolds shear stress contribution to the momentum loss (or the friction factor, respectively) is greater than the dispersive shear stress contribution. By increasing the bars spacing, however, the contribution of the dispersive shear stress increases while the Reynolds shear stress contribution decreases, which is due to a standing wave at the water surface in the flow over k -type roughness which results in significant spatial variations in the time-averaged velocities. Strong sweep events take place and contribute to the friction coefficient. Investigating the dynamics of the flow reveals that there is momentum source below the crest of the bars and momentum sink above them, leading to acceleration or deceleration of flow, respectively. The contribution of dispersive shear stress is significant only in the deceleration of the flow near the crest of the bars and in the acceleration of the flow under the water surface. Quantification of the three components of total kinetic energy, i.e. mean, turbulent, and wake kinetic energy, reveals that the largest contribution is that of the mean flow in both geometries. By increasing the bar spacing, the contributions of turbulent and wake kinetic energy, which are localized at the bar height, increase, while the kinetic energy of the mean flow decreases

An immersed boundary-based large-eddy simulation approach to predict the performance of vertical axis tidal turbines

Vertical axis tidal turbines (VATTs) are perceived to be an attractive alternative to their horizontal axis counterparts in tidal streams due to their omni-directionality. The accurate prediction of VATTs demands a turbulence simulation approach that is able to predict accurately flow separation and vortex shed- ding and a numerical method that can cope with moving boundaries. Thus, in this study an immersed boundary-based large-eddy simulation (LES-IB) method is refined to allow accurate simulation of the blade vortex interaction of VATTs. The method is first introduced and validated for a VATT subjected to laminar flow. Comparisons with highly-accurate body-fitted numerical models results demonstrate the method’s ability of reproducing accurately the performance and fluid mechanics of the chosen VATT. Then, the simulation of a VATT under turbulent flow is performed and comparisons with data from exper- iments and results from RANS-based models demonstrate the accuracy of the method. The vortex-blade interaction is visualised for various tip speed ratios and together with velocity spectra detailed insights into the fluid mechanics of VATTs are provided

Influence of bubble size, diffuser width and flow rate on the integral behaviour of bubble plumes

A large-eddy simulation based Eulerian-Lagrangian model is employed to quantify the impact of bubble size, diffuser diameter, and gas flow rate on integral properties of bubble plumes, such as the plume's width, centerline velocity, and mass flux. Calculated quantities are compared with experimental data and integral model predictions. Furthermore, the LES data were used to assess the behavior of the entrainment coefficient, the momentum amplification factor, and the bubble-to-momentum spread ratio. It is found that bubble plumes with constant bubble size and smaller diameter behave in accordance with integral plume models. Plumes comprising larger and non-uniform bubble sizes appear to deviate from past observations and model predictions. In multi-diameter bubble plumes, a bubble self-organisation takes place, i.e., small bubbles cluster in the center of the plume whilst large bubbles are found at the periphery of the plume. Multi-diameter bubble plumes also feature a greater entrainment rate than single-size bubble plumes, as well as a higher spread ratio and lower turbulent momentum rate. Once the plume is fully established, the size of the diffuser does not appear to affect integral properties of bubble plumes. However, plume development is affected by the diffuser width, as larger release areas lead to a delayed asymptotic behavior of the plume and consequently to a lower entrainment and higher spread ratio. Finally, the effect of the gas flow rate on the integral plume is studied and is deemed very relevant with regards to most integral plume properties and coefficients. This effect is already fairly well described by integral plume models

Effect of blade cambering on dynamic stall in view of designing vertical axis turbines

This paper presents large-eddy simulations of symmetric and asymmetric (cambered) airfoils forced to undergo deep dynamic stall due to a prescribed pitching motion. Experimental data in terms of lift, drag, and moment coefficients are available for the symmetric NACA 0012 airfoil and these are used to validate the large-eddy simulations. Good agreement between computed and experimentally observed coefficients is found confirming the accuracy of the method. The influence of foil asymmetry on the aerodynamic coefficients is analysed by subjecting a NACA 4412 airfoil to the same flow and pitching motion conditions. Flow visualisations and analysis of aerodynamic forces allow an understanding and quantification of dynamic stall on both straight and cambered foils. The results confirm that cambered airfoils provide an increased lift-to-drag ratio and a decreased force hysteresis cycle in comparison to their symmetric counterpart. This may translate into increased performance and lower fatigue loads when using cambered airfoils in the design of vertical axis turbines operating at low tip-speed ratios

Physically realistic roughness closure scheme to simulate turbulent channel flow over rough beds within the framework of LES

A physically realistic roughness closure method for the simulation of turbulent open-channel flow over natural beds within the framework of large-eddy simulation (LES) is proposed. The description of bed roughness in LES is accomplished through a roughness geometry function together with forcing terms in the momentum equations. The major benefit of this method is that the roughness is generated from one physically measurable parameter, i.e., the mean grain diameter of the bed material. A series of flows over rough beds, for which mean flow and turbulence statistics are available from experiments, is simulated. Measured and computed values are compared to validate the proposed roughness closure approach. It is found that predicted streamwise velocity profiles, turbulence intensities, and turbulent shear stress profiles match the measured values fairly well. Furthermore, the effect of roughness on the overall flow resistance is predicted in reasonable agreement with experimental values

Effect of three-dimensional mixing conditions on water treatment reaction process

The performance of water disinfection facilities traditionally relies on Hydraulic Efficiency Indicators (HEIs), extracted from experimentally derived Residence Time Distribution (RTD) curves. This approach has often been undertaken numerically through computational fluid dynamics (CFD) models, which can be calibrated to predict accurately RTDs, enabling the assessment of disinfection facilities prior to the construction of disinfection tanks. However, a significant drawback of the conventional efficiency methodology prescribed for disinfection tanks is associated with the respective indicators, as they are predominantly linked to the internal flow characteristics developed in the reactor, rather than the disinfection chemistry which should be optimized. In this study three-dimensional numerical models were refined to simulate the processes of chlorine decay, pathogen inactivation and the by-product formation in disinfection contact tanks (CTs). The main objective of this study was to examine the effect of three-dimensional mixing on the reaction processes which were modelled through finite-rate kinetic models. Comparisons have been made between pathogen inactivation and disinfection by-product accumulation results produced by a RANS approach against the findings of a Segregated Flow Analysis (SFA) of conservative tracer transport. CFD Results confirm that three-dimensional mixing does have an effect on the reaction processes, which, however, is not apparent through the SFA approach