Large-eddy Simulation Dataset Research Articles

A double-Gaussian wake model considering yaw misalignment

A wake steering has been known to effectively increase wind farm production by deflecting the upstream turbines’ wakes via yaw misalignment, thus minimizing their negative impacts on the downstream turbines' performances. This study presents analytical modeling of horizontal-axis wind turbine (HAWT) wake using low-cost analytical modeling as an alternative to expensive numerical and experimental trials. The existing double-Gaussian (DG) analytical wake model was modified to include the yaw misalignment effect, allowing its usability for the yawed HAWT wake modeling. The benchmark dataset produced by high-fidelity large eddy simulation (LES) of wake flowfields behind the turbine with yaw angles of 0º, 10º, 20º, and 30º were used to validate the accuracy of the DG yaw wake model. Overall, the DG yaw wake model predictions showed good agreement with the benchmark dataset under varying HAWT rotor yaw configurations. The analytical results verified by the LES dataset confirm the effectiveness of yaw misalignment in deflecting the wake trajectory, expediting the wake recovery downstream of the HAWT. In addition, a higher rotor yaw angle improves the wake recovery rate in the prevailing wind direction. Notable deviations against the benchmark dataset were found mainly within the near-wake region owing to flow acceleration arising from turbine-induced turbulence. As a result, the model’s predictions were slightly lower than the benchmark dataset, most likely due to neglecting the acceleration term in the analytical model derivation. Otherwise, the analytical model could accurately predict the mean wake velocity within the far-wake region for all evaluated cases, demonstrating its reliability in estimating wind speed potential within a practical distance for micrositing. These results were also proved quantitatively by statistical evaluations utilizing root mean square error (RMSE) and Pearson correlation coefficient R. The present study points out the importance of the upstream HAWTs’ rotor yaw controls to properly deflect their wakes away from their mainstream trajectories, thus effectively maximizing the wind speed potentials extracted by the downstream HAWTs and improving the overall wind farm production.

Improved numerical predictions of inclined negatively buoyant jet behaviour

The results of numerical simulations of inclined negatively buoyant jets are presented. These simulations address previously highlighted difficulties in capturing sufficient detail of critical flow processes to effectively predict the detailed flow behaviour. In particular, the new simulations are able to accurately capture the details of the buoyancy-induced instabilities, which are clearly evident in associated experimental investigations and that have significant impacts on the flow behaviour. This new information is captured for inclined negatively buoyant jets discharged at 45° above a horizontal reference plane. A Large Eddy Simulation (LES) approach is implemented that makes use of a Lagrangian Dynamic Sub-grid scale (SGS) model and a novel criterion for the adaptive meshing system. Comparisons with previously published simulation results and experimental data demonstrate that these new Adaptive LES simulations provide improved predictions of flow path, concentration and velocity fields, and associated mean and turbulent statistics. In addition, this study provides a set of methods for generating high-quality LES data sets for free shear flows, which are well beyond the level of detail that can be captured by current experimental systems.

Toward prediction of turbulent atmospheric flows over propagating oceanic waves via machine-learning augmented large-eddy simulation

Wind-wave interactions have important effects on the energy harvesting of offshore wind farms. High-fidelity large-eddy simulation (LES) is a powerful approach for investigating wind-wave interactions in turbulent oceanic environments. Due to the large scale of the flow domain and the high grid resolution required to resolve multi-scale flow motions, however, brute-force LES of wind-wave interactions is computationally very expensive. We propose augmenting brute-force LES via machine-learning data-driven modeling (ML-LES) to dramatically reduce the computational time required to obtain converged turbulence statistics when the brute-force approach is employed. Namely, we employ a convolutional neural network (CNN) autoencoder trained and validated with LES data sets to develop a highly efficient ML-LES approach for computing turbulence statistics from just a few snapshots of instantaneous LES flow fields. Our results demonstrate the accuracy and efficiency of ML-LES in predicting the mean velocity, velocity fluctuations, and turbulence kinetic energy in highly stretched computational grid systems required to carry out simulations in real-life oceanic environments.

Investigation of flow separation control over airfoil using nanosecond pulsed plasma actuator

A detailed investigation of flow separation control over a stalled NACA 0015 airfoil model using pulsed nanosecond dielectric barrier discharge (NS-DBD) plasma actuator is conducted through combined large eddy simulation (LES), dynamic mode decomposition (DMD) and experimental particle image velocimetry (PIV) measurement possibly for the first time. Airfoil flow at chord-based Reynolds number (Re) of Re = 6.3 × 104 is considered. It is found that the interplay of the residual heat due to the nanosecond plasma discharge with external flow causes the generation of vortices through the baroclinic and vortex dilation mechanisms and finally promotes the formation of large-scale vortical structures that dominate the early-stage evolution of the flow control process. One of the most important findings of this work is that the flow control mechanisms are not universal but depend on the forcing frequencies of plasma actuation, and are categorized into three different types. Meanwhile, three-dimensional DMD analysis of LES datasets is performed and the output of modal analysis helps provide more physical insights into flow control and determine the optimum forcing frequency. Depending on the value of forcing frequency, the global or local resonance of the forced flow with the periodic plasma actuation occurs. The plasma forcing at any frequencies in the range 1≤F+≤20 produces pronounced flow control authority except for that at F+=20. Moreover, the optimum forcing frequency is identified and its relationship with the inherent characteristic frequencies of the baseline flow is established based on DMD analysis. It is found that the most favorable control effect is achieved when the plasma excitation frequency is high and very close to the characteristic frequency of the separated shear layer in the unforced baseline flow.

LiSBOA (LiDAR Statistical Barnes Objective Analysis) for optimal design of lidar scans and retrieval of wind statistics – Part 1: Theoretical framework

Abstract. A LiDAR Statistical Barnes Objective Analysis (LiSBOA) for the optimal design of lidar scans and retrieval of the velocity statistical moments is proposed. LiSBOA represents an adaptation of the classical Barnes scheme for the statistical analysis of unstructured experimental data in N-dimensional space, and it is a suitable technique for the evaluation over a structured Cartesian grid of the statistics of scalar fields sampled through scanning lidars. LiSBOA is validated and characterized via a Monte Carlo approach applied to a synthetic velocity field. This revisited theoretical framework for the Barnes objective analysis enables the formulation of guidelines for the optimal design of lidar experiments and efficient application of LiSBOA for the postprocessing of lidar measurements. The optimal design of lidar scans is formulated as a two-cost-function optimization problem, including the minimization of the percentage of the measurement volume not sampled with adequate spatial resolution and the minimization of the error on the mean of the velocity field. The optimal design of the lidar scans also guides the selection of the smoothing parameter and the total number of iterations to use for the Barnes scheme. LiSBOA is assessed against a numerical data set generated using the virtual lidar technique applied to the data obtained from a large eddy simulation (LES). The optimal sampling parameters for a scanning Doppler pulsed wind lidar are retrieved through LiSBOA, and then the estimated statistics are compared with those of the original LES data set, showing a maximum error of about 4 % for both mean velocity and turbulence intensity.

On the genesis and evolution of barchan dunes: Hydrodynamics

Barchans are an important signature of turbulent atmospheric and aquatic flows on parts of the Earth’s surface where the supply of sediments is scarce. Khosronejad and Sotiropoulos [“On the genesis and evolution of barchan dunes: Morphodynamics,” J. Fluid Mech. 815, 117–148 (2017)] demonstrated that high-fidelity hydro-morphodynamic large-eddy simulation (LES) can replicate sub-aqueous barchan fields exhibiting the geometric features and morphodynamic interactions observed in nature. Herein, we first validate the ability of our LES method to simulate mean flow and turbulence statistics over a frozen barchan dune created in a control laboratory experiment. We then analyze our high-fidelity LES datasets to elucidate the hydrodynamic phenomena that drive the previously identified morphodynamic interactions at various stages of subaqueous barchan dune evolution, from flatbed to quasi-equilibrium. We uncover and quantify the hydrodynamic mechanisms that cause (i) the initiation of crescent microfeatures on the initial flatbed; (ii) the development of three-dimensional barchan dunes; and (iii) several other phenomena that occur during the barchan field maturation, such as the birth of new small barchans and the exposure of the underlying bedrock in the wake of the barchans. Defining a dimensionless timescale, i.e., the saturation timescale, as a function of the saturation length scale, we attempt to categorize for the first time various temporal stages of barchan field development. We show that the barchan field has specific characteristics corresponding to each of the barchan stage developments in terms of both the flow field and the bed topology. Comparing the simulation results of transverse and barchan dunes, we illustrate that the two dune types have similar topographic and hydrodynamics characteristics until bedrock exposure occurs.

Large-eddy simulation of turbulent natural-bed flow

The turbulent flow in natural rough beds is a complex subject, still poorly understood despite the longstanding effort of several researchers. In the present work, a turbulent open-channel flow experiment, with a pebble bed at Reynolds and Froude numbers, respectively, Re = 4.65 × 104 and Fr = 0.186, has been simulated using the Large-Eddy Simulation (LES) technique, in which the wall-adapting local eddy viscosity subgrid scale closure model is used and in the presence of an air-water interface to take into account the effects of the interface deformation in the flow turbulence statistics under a low relative submergence condition. The simulations have been compared with a companion experiment, where the channel bottom is constituted by four pebble layers. For the simulations, the pebble-bed surface has been captured with a high-resolution three-dimensional laser scanner and used to morphologically characterize the numerical channel bottom. Results are presented in terms of turbulence statistics and turbulent laws, showing a good agreement with those obtained in the experiment. Since a good convergence between simulation and experimental results was obtained, the LES dataset was used to compute the Turbulent Kinetic Energy (TKE) dissipation rate across the water depth. The mesh resolution allows showing a detailed TKE dissipation rate distribution across the water depth. Moreover, the equilibrium between TKE production and dissipation was checked to verify the overlap layer existence under low relative submergence condition. Finally, a new procedure for vortex-visualization is implemented, based on the relationship between the vorticity and the TKE dissipation rate.

Two-dimensional unsteadiness map of oblique shock wave/boundary layer interaction with sidewalls

The low-frequency unsteadiness of oblique shock wave/boundary layer interactions (SBLIs) has been investigated using large-eddy simulation (LES) and high-frequency pressure measurements from experiments. Particular attention has been paid to off-centreline behaviour: the LES dataset was generated including sidewalls, and experimental pressure measurements were acquired across the entire span of the reflected shock foot. The datasets constitute the first maps of low-frequency unsteadiness in both streamwise and spanwise directions. The results reveal that significant low-frequency shock motion (with $St\approx 0.03$) occurs away from the centreline, along most of the central separation shock and in the corner regions. The most powerful low-frequency unsteadiness occurs off-centre, likely due to the separation shock being strengthened by shocks arising from the swept interactions on the sidewalls. Both simulation and experimental results exhibit asymmetry about the spanwise centre. In simulations, this may be attributed to a lack of statistical convergence; however, the fact that this is also seen in experiments is indicative that some SBLIs may exhibit some inherent asymmetry across the two spanwise halves of the separation bubble. There is also significant low-frequency power in the corner separations. The relation of the unsteadiness in the corner regions to that in the centre is investigated by means of two-point correlations: a key observation is that significant correlation does not extend across the attached flow channel between the central and corner separations.

A Three-Dimensional Scale-Adaptive Turbulent Kinetic Energy Scheme in the WRF-ARW Model

A new three-dimensional (3D) turbulent kinetic energy (TKE) subgrid mixing scheme is developed using the Advanced Research version of the Weather Research and Forecasting (WRF) Model (WRF-ARW) to address the gray-zone problem in the parameterization of subgrid turbulent mixing. The new scheme combines the horizontal and vertical subgrid turbulent mixing into a single energetically consistent framework, in contrast to the conventionally separate treatment of the vertical and horizontal mixing. The new scheme is self-adaptive to the grid-size change between the large-eddy simulation (LES) and mesoscale limits. A series of dry convective boundary layer (CBL) idealized simulations are carried out to compare the performance of the new scheme and the conventional treatment of subgrid mixing to the WRF-ARW LES dataset. The importance of including the nonlocal component in the vertical buoyancy specification in the newly developed general TKE-based scheme is illustrated in the comparison. The improvements of the new scheme with the conventional treatment of subgrid mixing across the gray-zone model resolutions are demonstrated through the partitioning of the total vertical flux profiles. Results from real-case simulations show the feasibility of using the new scheme in the WRF Model in lieu of the conventional treatment of subgrid mixing.

Editorial: IFP Energies nouvelles International Conference LES4ICE 2014 – Large-Eddy Simulation for Internal Combustion Engine Flows

Further improving the environmental performances of Internal Combustion Engines (ICE) increasingly requires moving beyond traditional multidimensional simulation tools based on a cycle averaged approach (Reynolds-Averaged Navier–Stokes, RANS), and to reliably predict and control individual engine cycles under realistic operating conditions. Large-Eddy Simulation (LES) offers this unique potential by predicting spatially filtered flow realizations, thus opening up new perspectives for extending the scope of application of Computational Fluid Dynamics (CFD) for ICE. Since its 1 edition in 2008, the LES4ICE conference provides a forum for exchange concerning research and development of LES and related experimental techniques for their application to ICE flows. It brings together researchers and engineers working in the field of piston engine combustion to debate the state of the art in LES applied to ICE and examine advanced experimental techniques capable of supporting and validating its development. The present special issue proposes an extract from the contributions to the 2014 edition of LES4ICE, following a selection by the conference’s Scientific Committee. The selected four articles provide a view on some recent advances in the domain of LES and related experimental techniques, and on their application to ICE flows. The study of Cyclic Combustion Variability (CCV) in Spark-Ignition Engines (SIE) has become a prominent topic for LES research in the past decade. It has been demonstrated how combining advanced optical diagnostics and state-of-the-art LES can allow identifying the causes of CCV, and that the acquired understanding could efficiently be capitalized in the form of reduced models for 1D CFD simulations of complete engine systems. This on-going research opened exciting perspectives for the application of LES to other yet poorly understood and mastered non-cyclic phenomena in SIE, as in particular the link between CCV and engine knock. A longstanding issue when studying the causes of CCV in SI engines is the relative importance of the variability of large-scale coherent structures, and of small scale turbulent fluctuations. A method aiming at addressing this question was proposed by S. Buhl, F. Hartmann and C. Hasse, in their contribution “Identification of Large-Scale Structure Fluctuations in IC Engines using POD-Based Conditional Averaging” [1]. They propose a method combining POD and conditional averaging, and discuss its application to the analysis of the origins of CCV for two different LES datasets. Oil & Gas Science and Technology – Rev. IFP Energies nouvelles (2016) 71, E1 C. Angelberger and C. Mounaim-Rousselle, published by IFP Energies nouvelles, 2016 DOI: 10.2516/ogst/2015048 IFP Energies nouvelles International Conference Rencontres Scientifiques d'IFP Energies nouvelles

A Year-Long Large-Eddy Simulation of the Weather over Cabauw: An Overview

Abstract Results are presented of two large-eddy simulation (LES) runs of the entire year 2012 centered at the Cabauw observational supersite in the Netherlands. The LES is coupled to a regional weather model that provides the large-scale information. The simulations provide three-dimensional continuous time series of LES-generated turbulence and clouds, which can be compared in detail to the extensive observational dataset of Cabauw. The LES dataset is available from the authors on request. This type of LES setup has a number of advantages. First, it can provide a more statistical approach to the study of turbulent atmospheric flow than the more common case studies, since a diverse but representative set of conditions is covered, including numerous transitions. This has advantages in the design and evaluation of parameterizations. Second, the setup can provide valuable information on the quality of the LES model when applied to such a wide range of conditions. Last, it also provides the possibility to emulate observation techniques. This might help detect limitations and potential problems of a variety of measurement techniques. The LES runs are validated through a comparison with observations from the observational supersite and with results from the “parent” large-scale model. The long time series that are generated, in combination with information on the spatial structure, provide a novel opportunity to study time scales ranging from seconds to seasons. This facilitates a study of the power spectrum of horizontal and vertical wind speed variance to identify the dominant variance-containing time scales.

Evaluating Large-Eddy Simulation (LES) and High-Speed Particle Image Velocimetry (PIV) with Phase-Invariant Proper Orthogonal Decomposition (POD)

This study is part of a program to understand the stochastic variations in IC engine flows; in particular, it is a comparison of measured (PIV) and computed (LES) velocity from multiple cycles of the same motored engine. Comparison procedures included traditional RANS (Reynolds Averaged Navier-Stokes) decomposition (ensemble-averaged and RMS (Root Mean Square) velocity), phase-dependent, and phase- invariant POD. Phase-dependent POD was performed on the PIV and LES samples separately and on the combined samples, thus creating separate or a single POD mode sets, respectively. The phase- invariant POD was performed both on normalized snapshots and on snapshots where the original energy was conserved. Initial comparisons of the mass-specific kinetic energies of the ensemble average and RMS velocities revealed that the PIV and LES data sets differed significantly during most of the intake stroke. This discrepancy was quantifiedfirst by comparing the relevance indices calculated between ensemble average velocity fields and, second, using phase-dependent POD, which quantified cycle-to-cycle flow variations of the RANS average and turbulence. Phase-dependent POD was applied separately to the PIV and LES data sets during the intake stroke (76° ATDCE, After Top Dead Center Exhaust), where the intake-valve jet is strong and the PIV and LES data were earlier found to be significantly different. The cyclic variability of the LES ensemble average was estimated to be significantly higher than that of the PIV data set. POD was also applied to the combined sample of LES and PIV snapshots for quantitative comparison creating a single set of modes, so that comparison could be made with POD coefficients alone. Example comparisons were made at again for data at 76° and also 330° ATDCE, which is a viable spark timing in a fired engine. The results at 76° ATDCE were similar to those obtained with the POD analysis of the separate samples. At 330° ATDCE, the PIV ensemble average and RMS velocities showed somewhat more cyclic variability. Phase-invariant POD was applied to the combined PIV and LES velocity data from all crank angles to study the flow evolution over the crank angle range. Conclusions derived from the two different energy transformations are contrasted. When the energy of each velocity field is normalized, the phase-invariant POD results focus on differences in flow structures and their evolution. On the other hand, when the energy of each velocity field is conserved, the phase-invariant POD results also take into account differences in energy between the PIV and LES data.

A Large Eddy Simulation Intercomparison Study of Shallow Cumulus Convection

This paper reports an intercomparison study on undisturbed trade wind cumulus convection under steadystate conditions as observed during the Barbados Oceanographic and Meteorological Experiment (BOMEX) with 10 large eddy simulation (LES) models. A main objective of this study is to obtain a quantitative assessment of the quality of the turbulent dynamics for this type of boundary layer clouds as produced by the different LES codes. A 6-h simulation shows excellent model-to-model agreement of the observed vertical thermodynamical structure, reasonable agreement of variances and turbulent fluxes, and good agreement of quantities conditionally sampled within the model clouds, such as cloud cover, liquid water, and cloud updraft strength. In the second part of this paper the LES dataset is used to evaluate simple models that are used in parameterizations of current general circulation models (GCMs). Finally, the relation of this work to subsequent LES studies of more complicated regimes is discussed, and guidance is given for the design of future observational studies of shallow cumulus boundary layers.