scholarly journals Modelling waving crops using large-eddy simulation: comparison with experiments and a linear stability analysis

2010 ◽  
Vol 652 ◽  
pp. 5-44 ◽  
Author(s):  
S. DUPONT ◽  
F. GOSSELIN ◽  
C. PY ◽  
E. DE LANGRE ◽  
P. HEMON ◽  
...  

In order to investigate the possibility of modelling plant motion at the landscape scale, an equation for crop plant motion, forced by an instantaneous velocity field, is introduced in a large-eddy simulation (LES) airflow model, previously validated over homogeneous and heterogeneous canopies. The canopy is simply represented as a poroelastic continuous medium, which is similar in its discrete form to an infinite row of identical oscillating stems. Only one linear mode of plant vibration is considered. Two-way coupling between plant motion and the wind flow is insured through the drag force term. The coupled model is validated on the basis of a comparison with measured movements of an alfalfa crop canopy. It is also compared with the outputs of a linear stability analysis. The model is shown to reproduce the well-known phenomenon of ‘honami’ which is typical of wave-like crop motions on windy days. The wavelength of the main coherent waving patches, extracted using a bi-orthogonal decomposition (BOD) of the crop velocity fields, is in agreement with that deduced from video recordings. The main spatial and temporal characteristics of these waving patches exhibit the same variation with mean wind velocity as that observed with the measurements. However they differ from the coherent eddy structures of the wind flow at canopy top, so that coherent waving patches cannot be seen as direct signatures of coherent eddy structures. Finally, it is shown that the impact of crop motion on the wind dynamics is negligible for current wind speed values. No lock-in mechanism of coherent eddy structures on plant motion is observed, in contradiction with the linear stability analysis. This discrepancy may be attributed to the presence of a nonlinear saturation mechanism in LES.

Author(s):  
Thomas L. Kaiser ◽  
Thierry Poinsot ◽  
Kilian Oberleithner

The hydrodynamic instability in an industrial, two-staged, counter-rotative, swirled injector of highly complex geometry is under investigation. Large eddy simulations show that the complicated and strongly nonparallel flow field in the injector is superimposed by a strong precessing vortex core. Mean flow fields of large eddy simulations, validated by experimental particle image velocimetry measurements are used as input for both local and global linear stability analysis. It is shown that the origin of the instability is located at the exit plane of the primary injector. Mode shapes of both global and local linear stability analysis are compared to a dynamic mode decomposition based on large eddy simulation snapshots, showing good agreement. The estimated frequencies for the instability are in good agreement with both the experiment and the simulation. Furthermore, the adjoint mode shapes retrieved by the global approach are used to find the best location for periodic forcing in order to control the precessing vortex core.


Author(s):  
Yunfei Wang ◽  
Huanlong Chen ◽  
Huaping Liu ◽  
Yanping Song ◽  
Fu Chen

An in-house large eddy simulation (LES) code based on three-dimensional compressible N-S equations is used to research the impact of incoming wakes on unsteady evolution characteristic in a low-pressure turbine (LPT) cascade. The Mach number is 0.4 and Reynolds number is 0.6 × 105 (based on the axial chord and outlet velocity). The reduced frequency of incoming wakes is Fred = 0 (without wakes), 0.37 and 0.74. A detailed analysis of Reynolds stresses and turbulent kinetic energy inside the boundary layer has been carried out. Particular consideration is devoted to the transport process of incoming wakes and the intermittent property of the unsteady boundary layer. With the increase of reduced frequency, the inhibiting effect of wakes on boundary layer separation gradually enhances. The separation at the rear part of the suction side is weakened and the separation point moves downstream. However, incoming wakes lead to an increase in dissipation and aerodynamic losses in the main flow area. Excessive reduced frequency ( Fred = 0.74) causes the main flow area to become one of the main source areas of loss. An optimal reduced frequency exists to minimize the aerodynamic loss of the linear cascade.


2010 ◽  
Vol 10 (10) ◽  
pp. 24345-24370
Author(s):  
V. Anabor ◽  
U. Rizza ◽  
G. A. Degrazia ◽  
E. de Lima Nascimento

Abstract. An isolated and stationary microburst is simulated using a 3-D time-dependent, high resolution Large-Eddy Simulation (LES) model. The microburst downdraft is initiated by specifying a simplified cooling source at the top of the domain near 2 km. The modelled time scale for this damaging wind (30 m/s) is of order of few min with a spatial scale enclosing a region with 500 m radius around the impact point. These features are comparable with results obtained from full-cloud models. The simulated flow shows the principal features observed by Doppler radar and others observational full-scale downburst events. In particular are observed the expansion of the primary and secondary cores, the presence of the ring vortex at the leading edge of the cool outflow, and finally an accelerating outburst of surface winds. This result evidences the capability of LES to reproduce complexes phenomena like a Microburst and indicates the potential of LES for utilization in atmospheric phenomena situated below the storm scale and above the microscale, which generally involves high velocities in a short time scale.


2001 ◽  
Vol 436 ◽  
pp. 353-391 ◽  
Author(s):  
J. C. R. HUNT ◽  
N. D. SANDHAM ◽  
J. C. VASSILICOS ◽  
B. E. LAUNDER ◽  
P. A. MONKEWITZ ◽  
...  

Recent research is making progress in framing more precisely the basic dynamical and statistical questions about turbulence and in answering them. It is helping both to define the likely limits to current methods for modelling industrial and environmental turbulent flows, and to suggest new approaches to overcome these limitations. Our selective review is based on the themes and new results that emerged from more than 300 presentations during the Programme held in 1999 at the Isaac Newton Institute, Cambridge, UK, and on research reported elsewhere. A general conclusion is that, although turbulence is not a universal state of nature, there are certain statistical measures and kinematic features of the small-scale flow field that occur in most turbulent flows, while the large-scale eddy motions have qualitative similarities within particular types of turbulence defined by the mean flow, initial or boundary conditions, and in some cases, the range of Reynolds numbers involved. The forced transition to turbulence of laminar flows caused by strong external disturbances was shown to be highly dependent on their amplitude, location, and the type of flow. Global and elliptical instabilities explain much of the three-dimensional and sudden nature of the transition phenomena. A review of experimental results shows how the structure of turbulence, especially in shear flows, continues to change as the Reynolds number of the turbulence increases well above about 104 in ways that current numerical simulations cannot reproduce. Studies of the dynamics of small eddy structures and their mutual interactions indicate that there is a set of characteristic mechanisms in which vortices develop (vortex stretching, roll-up of instability sheets, formation of vortex tubes) and another set in which they break up (through instabilities and self- destructive interactions). Numerical simulations and theoretical arguments suggest that these often occur sequentially in randomly occurring cycles. The factors that determine the overall spectrum of turbulence were reviewed. For a narrow distribution of eddy scales, the form of the spectrum can be defined by characteristic forms of individual eddies. However, if the distribution covers a wide range of scales (as in elongated eddies in the ‘wall’ layer of turbulent boundary layers), they collectively determine the spectra (as assumed in classical theory). Mathematical analyses of the Navier–Stokes and Euler equations applied to eddy structures lead to certain limits being defined regarding the tendencies of the vorticity field to become infinitely large locally. Approximate solutions for eigen modes and Fourier components reveal striking features of the temporal, near-wall structure such as bursting, and of the very elongated, spatial spectra of sheared inhomogeneous turbulence; but other kinds of eddy concepts are needed in less structured parts of the turbulence. Renormalized perturbation methods can now calculate consistently, and in good agreement with experiment, the evolution of second- and third-order spectra of homogeneous and isotropic turbulence. The fact that these calculations do not explicitly include high-order moments and extreme events, suggests that they may play a minor role in the basic dynamics. New methods of approximate numerical simulations of the larger scales of turbulence or ‘very large eddy simulation’ (VLES) based on using statistical models for the smaller scales (as is common in meteorological modelling) enable some turbulent flows with a non-local and non-equilibrium structure, such as impinging or convective flows, to be calculated more efficiently than by using large eddy simulation (LES), and more accurately than by using ‘engineering’ models for statistics at a single point. Generally it is shown that where the turbulence in a fluid volume is changing rapidly and is very inhomogeneous there are flows where even the most complex ‘engineering’ Reynolds stress transport models are only satisfactory with some special adaptation; this may entail the use of transport equations for the third moments or non-universal modelling methods designed explicitly for particular types of flow. LES methods may also need flow-specific corrections for accurate modelling of different types of very high Reynolds number turbulent flow including those near rigid surfaces.This paper is dedicated to the memory of George Batchelor who was the inspiration of so much research in turbulence and who died on 30th March 2000. These results were presented at the last fluid mechanics seminar in DAMTP Cambridge that he attended in November 1999.


2014 ◽  
Vol 695 ◽  
pp. 562-566
Author(s):  
Afiq Witri Muhammad Yazid ◽  
Nor Azwadi Che Sidik ◽  
Salim Mohamed Salim ◽  
Shuhaimi Mansor

This paper reports on the model sensitivity analysis of a commercial computational fluid dynamics program, ANSYS FLUENT v14. The purpose of the analysis was to determine the appropriate modeling settings for numerical model of the case study. A full scale of a simplified urban street canyon was modelled and the turbulent flow was calculated using Large Eddy Simulation (LES) techniques. The model sensitivity tests involved are mesh sensitivity, statistically steady state and sampling. Adequate numbers of cells, period time to achieve statistically steady state (SST) and sampling time to simulate wind flow and pollutant dispersion in street canyon were determined through systematic tests.


2011 ◽  
Vol 679 ◽  
pp. 156-193 ◽  
Author(s):  
EZGI S. TAŞKINOĞLU ◽  
JOSETTE BELLAN

For flows at supercritical pressure, p, the large-eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas equation of state, and the equations utilize transport properties depending on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid-scale (SGS) fluxes but also new SGS terms, each denoted as a ‘correction’. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations and represent differences, other than contributed by the convection terms, between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) which is the difference between the filtered divergence of the molecular heat flux and the divergence of the molecular heat flux computed as a function of the filtered flow field. We revisit here a previous a priori study where we only had partial success in modelling the q-correction term and show that success can be achieved using a different modelling approach. This a priori analysis, based on a temporal mixing-layer direct numerical simulation database, shows that the focus in modelling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES molecular heat-flux calculation. Furthermore, an a posteriori study is conducted for temporal mixing layers initially containing oxygen (O) in the lower stream and hydrogen (H) or helium (He) in the upper stream to examine the benefit of the new model. Results show that for any LES including SGS-flux models (constant-coefficient gradient or scale-similarity models; dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/gradient models), the inclusion of the q-correction in LES leads to the theoretical maximum reduction of the SGS molecular heat-flux difference; the remaining error in modelling this new subgrid term is thus irreducible. The impact of the q-correction model first on the molecular heat flux and then on the mean, fluctuations, second-order correlations and spatial distribution of dependent variables is also demonstrated. Discussions on the utilization of the models in general LES are presented.


Computation ◽  
2020 ◽  
Vol 8 (2) ◽  
pp. 47
Author(s):  
Gokhan Kirkil ◽  
Ching-Long Lin

A high-resolution large eddy simulation (LES) of wind flow over the Oklahoma City downtown area was performed to explain the effect of the building height on wind flow over the city. Wind flow over cities is vital for pedestrian and traffic comfort as well as urban heat effects. The average southerly wind speed of eight meters per second was used in the inflow section. It was found that heights and distribution of the buildings have the greatest impact on the wind flow patterns. The complexity of the flow field mainly depended on the location of buildings relative to each other and their heights. A strong up and downflows in the wake of tall buildings as well as large-scale coherent eddies between the low-rise buildings were observed. It was found out that high-rise buildings had the highest impact on the urban wind patterns. Other characteristics of urban canopy flows, such as wind shadows and channeling effects, are also successfully captured by the LES. The LES solver was shown to be a powerful tool for understanding urban canopy flows; therefore, it can be used in similar studies (e.g., other cities, dispersion studies, etc.) in the future.


Sign in / Sign up

Export Citation Format

Share Document