Le phénomène connu sous l'appellation "méandrement" (ou meandering) désigne les variations aléatoires de la trajectoire du sillage aérodynamique d'une éolienne. Ce phénomène est responsable de contraintes mécaniques particulières sur les éoliennes positionnées dans le sillage d'autres éoliennes et joue donc rôle dans la conception et dans la prévision de production des parcs éoliens.Ce travail propose d'étudier le méandrement par des moyens expérimentaux et numériques. La problématique est traitée en deux parties, la première portant sur la modélisation de l'écoulement de couche limite atmosphérique, avec une attention particulière portée à la modélisation des grandes échelles de la turbulence atmosphérique. La seconde partie porte sur l'étude du sillage d'un disque actuateur soumis à un écoulement atmosphérique. Chacune de ces parties comporte un volet expérimental et un volet numérique. La modélisation numérique instationnaire de l'écoulement atmosphérique fait intervenir une technique de génération stochastique de champs de vitesse turbulente avec évolution temporelle, spécialement développée au cours de la présente thèse et à laquelle un chapitre spécifique est dédié.L'un des principaux résultats est que le méandrement du sillage est fortement corrélé avec les grandes échelles de la turbulence atmosphérique. / The phenomenon known as meandering describes the unsteady trajectory variations of the wake of a wind turbine. This phenomenon is responsible for specific mechanical stresses on turbines positioned in the wake of other turbines. As such, this phenomenon must be accounted for in the design and operation of wind turbine plants.This work uses numerical fluid simulation and wind tunnel testing in order to study the meandering of the wake of a wind turbine. The subject is discussed in two parts. The first part discusses the modeling of the atmospheric boundary layer, with a focus on the large scales of the atmospheric turbulence. The second part is a study of the behavior of the wake of an actuator disc model in atmospheric wind conditions.Both parts include experimental and numerical work. The numerical simulation of the atmospheric boundary layer involves the generation of synthetic turbulent velocity time series by mean of a stochastic technique developed during this thesis, to which a chapter is dedicated.One of the main results of this work is that the meandering is highly correlated with the large scales of the atmospheric turbulence.

The spectral-element method (SEM) is used to study wall-bounded turbulent flowsin moderately complex geometries. The first part of the thesis is devoted to simulations of canonical flow cases, such as temporal K-type transitionand turbulent channel flow, to investigate general resolution requirements and computational efficiency of the numerical code nek5000. Large-eddy simulation (LES) is further performed of a plane asymmetric diffuser flow with an opening angle of 8.5 degrees, featuring turbulent flow separation. Good agreement with numerical studies of Herbst (2007) is obtained, and it is concluded that the use of a high-order method is advantageous for flows featuring pressure-induced separation. Moreover, it is shown, both a priori on simpler model problems and a posteriori using the full Navier--Stokes equations, that the numerical instability associated with SEM at high Reynolds numbers is cured either by employing over-integration (dealiasing) or a filter-based stabilisation, thus rendering simulations of moderate to high Reynolds number flows possible. The second part of the thesis is devoted to the first direct numerical simulation (DNS) of a truly three-dimensional, turbulent and separated diffuser flow at Re = 10 000 (based on bulk velocity and inflow-duct height), experimentally investigated by Cherry et al. (2008). The massively parallel capabilities of the spectral-element method are exploited by running the simulations on up to 32 768 processors. Very good agreement with experimental mean flow data is obtained and it is thus shown that well-resolved simulations of complex turbulent flows with high accuracy are possible at realistic Reynolds numberseven in complicated geometries. An explanation for the discovered asymmetry of the mean separated flow is provided and itis demonstrated that a large-scale quasi-periodic motion is present in the diffuser. In addition, a new diagnostic measure, based on the maximum vorticity stretching component in every spatial point, is designed and tested in a number of turbulent and transitional flows. Finally, Koopman mode decomposition is performed of a minimal channel flow and compared to classical proper orthogonal decomposition (POD). / QC 20111206

Atmospheric models are strongly dependent on the turbulent exchange of momentum, sensible heat and moisture (latent heat) at the surface. Oceans cover about 70% of the Earth’s surface and understanding the processes that control air-sea exchange is of great importance in order to predict weather and climate. In the atmosphere, for instance, hurricane development, cyclone intensity and track depend on these processes. Ocean waves constitute an obvious example of air-sea interaction and can cause the air-flow over sea to depend on surface conditions in uniquely different ways compared to boundary layers over land. When waves are generated by wind they are called wind sea or growing sea, and when they leave their generation area or propagate faster than the generating wind they are called swell. The air-sea exchange is mediated by turbulent eddies occurring on many different scales. Field measurements and high-resolution turbulence resolving numerical simulations have here been used to study these processes. The standard method to measure turbulent fluxes is the eddy covariance method. A spatial separation is often used between instruments when measuring scalar flux; this causes an error which was investigated for the first time over sea. The error is typically smaller over ocean than over land, possibly indicating changes in turbulence structure over sea. Established and extended analysis methods to determine the dominant scales of momentum transfer was used to interpret how reduced drag and sometimes net upward momentum flux can persist in the boundary layer indirectly affected by swell. A changed turbulence structure with increased turbulence length scales and more effective mixing was found for swell. A study, using a coupled wave-atmosphere regional climate model, gave a first indication on what impact wave mixing have on atmosphere and wave parameters. Near surface wind speed and wind gradients was affected especially for shallow boundary layers, which typically increased in height from the introduced wave-mixing. A large impact may be expected in regions of the world with predominant swell. The impact of swell waves on air-sea exchange and mixing should be taken into account to develop more reliable coupled Earth system models.

Large-eddy simulation (LES) and unsteady Reynolds-averaged Navier-Stokes (URANS) calculations have been performed to investigate the effects of different mathematical models for scalar variance and its dissipation rate as applied to both a non-reacting bluff-body turbulent flow and an extension to a reacting case. In the conserved scalar formalism, the mean value of a thermo-chemical variable is obtained through the PDF-weighted integration of the local description over the conserved scalar, the mixture fraction. The scalar variance, one of the key parameters for the determination of a presumed β-function PDF, is obtained by solving its own transport equation with the unclosed scalar dissipation rate modelled using either an algebraic expression or a transport equation. The proposed approach is first applied to URANS and then extended to LES. Velocity, length and time scales associated with the URANS modelling are determined using the standard two-equation k-ε transport model. In contrast, all three scales required by the LES modelling are based on the Smagorinsky subgrid scale (SGS) algebraic model. The present study proposes a new algebraic and a new transport LES model for the scalar dissipation rate required by the transport equation for scalar variance, with a time scale consistent with the Smagorinsky SGS model.

This thesis is concerned with investigations on two related issues of turbulent jet diffusion flame, namely (a) stabilization at liftoff and (b) autoignition in a turbulent jet diffusion flame. The approach of Conditional Moment Closure (CMC) has been taken. Fully elliptic first order CMC equations are solved with detailed chemistry to simulate lifted H2/N2 flame in vitiated coflow. The same approach is further used to simulate transient autoignition process in inhomogeneous mixing layers. In Chapter 1, difficulties involved in numerical simulation of turbulent combustion problems are explained. Different numerical tools used to simulate turbulent combustion are briefly discussed. Previous experimental, theoretical and numerical studies of lifted jet diffusion flames and autoignition are reviewed. Various research issues related to objectives of the thesis are discussed. In Chapter 2, the first order CMC transport equations for the reacting flows are presented. Various closure models that are required for solving the governing equations are given. Calculation of mean reaction rate term for detailed chemistry is given with special focus on the reaction rates for pressure dependent reactions. In Chapter 3, starting with the laminar flow code, further extension is carried to include kε turbulence model and PDF model. The code is validated at each stage of inclusion of different model. In this chapter, the code is first validated for the test problem of constant density, 2D, axisymmetric turbulent jet. Further, validation of PDF model is carried out by simulating the problem of nonreacting jet of cold air issuing into a vitiated coflow. The results are compared with the published data from experiments as well as numerical simulations. It is shown that the results compare well with the data. In Chapter 4, numerical results of lifted jet diffusion flame are presented. Detailed chemistry is modelled using Mueller mechanism for H2/O2 system with 9 species and 21 reversible reactions. Simulations are carried out for different jet velocities and coflow stream temperatures. The predicted liftoff generally agrees with experimental data, as well as joint PDF results. Profiles of mean scalar fluxes in the mixture fraction space, for different coflow temperatures reveal that (1) Inside the flamezone, the chemical term balances the molecular diffusion term, and hence the structure is of a diffusion flamelet for both cases. (2) In the preflame zone, the structure depends on the coflow temperature: for low coflow temperatures, the chemical term being small, the advective term balances the axial diffusion term. However, for the high coflow temperature case, the chemical term is large and balances the advective term, the axial diffusion term being small. It is concluded that, liftoff is controlled (a) by turbulent premixed flame propagation for low cofflow temperature while (b) by autoignition for high coflow temperature. In Chapter 5, the numerical results of autoignition in inhomogeneous mixing layer are presented. The configuration consists of a fuel jet issued into hot air for which transient simulations are performed. It is found that the constants assumed in various modelling terms can severely influence the results, particularly the flame temperature. Hence, modifications to these constants are suggested to obtain improved predictions. Preliminary work is carried out to predict autoignition lengths (which may be defined by Tign × Ujet incase of jet- and coflowvelocities being equal) by varying the coflow temperature. The autoignition lengths show a reasonable agreement with the experimental data and LES results. In Chapter 6, main conclusions of this thesis are summarized. Possible future studies on this problem are suggested.

Αντικείμενο της διατριβής είναι η προσομοίωση της ροής και της μεταφοράς θερμότητας σε ηλιακά και ενεργειακά συστήματα. Η έμφαση δόθηκε στις δεξαμενές αποθήκευσης της θερμότητας που παράγεται στα συγκεκριμένα συστήματα, με στόχο τον χαρακτηρισμό των ενεργειακών απωλειών και την βελτιστοποίηση του σχεδιασμού τους. Κύριες δραστηριότητες της διατριβής θα είναι η περαιτέρω ανάπτυξη διαθέσιμων εργαλείων προσομοίωσης ροών φυσικής και μικτής συναγωγής, με διερεύνηση των νεώτερων εξελίξεων στην μοντελοποίηση με τη μέθοδο Προσομοίωσης Μεγάλων Δινών (LES). Αρχικά γίνεται εκτεταμένη επικύρωση με πειραματικά αποτελέσματα σε απλές γεωμετρικές διατάξεις (π.χ. ορθογωνικά κανάλια ή κοιλώματα με βαθμίδα θερμοκρασίας) από την βιβλιογραφία. Στη συνέχεια η μεθοδολογία εφαρμόζεται στον υπολογισμό ροών σε πιο ρεαλιστικές γεωμετρίες, επιλεγμένες από πρακτικές εφαρμογές, όπως οι δεξαμενές αποθήκευσης νερού. Αναλύονται σε βάθος οι δυναμικές διεργασίες και τα ροϊκά φαινόμενα τόσο κατά την προσαγωγή της θερμότητας στη δεξαμενή (φόρτιση) όσο και κατά την απαγωγή της (εκφόρτιση) και η επίδραση που έχουν αυτά στην αποδοτικότητα της αποθήκευσης με βάση κατάλληλους ποσοτικούς δείκτες. Από τα αποτελέσματα αναδεικνύεται η σημασία της μοντελοποίησης σε τέτοιου είδους συστήματα ως ένα σημαντικό εργαλείο στη διερεύνηση της απόδοσης τους, του ενεργειακού χαρακτηρισμού τους και ακολούθως στην προσπάθεια επίτευξης του βέλτιστου σχεδιασμού τους. / The subject of the thesis is the Simulation of Turbulent Flow and Heat Transfer in Solar and Energy Systems. Emphasis is given in the thermal storage component of these systems, with the aim of characterizing their energy losses and improve their design. Main activities of the thesis will be the further development of available computational tools for the simulation of flows in natural and mixed convection, incorporating some of the most recent developments in modeling, particularly in the Large Eddy Simulation (LES) method. Initially, an extensive validation with experimental results in simple geometric configurations is carried out (e.g. channels or differentially heated cavities). Subsequently, the methodology is applied in the calculation of flows for more realistic geometries selected from practical applications, such as various hot water storage tanks. Analysis is conducted of the dynamic processes and relevant physical phenomena during the heat supply (charging) to and removal (discharging) from the tank and their influence on the storage effectiveness using appropriate thermodynamic indices. From the simulation results, the significance of the flow and heat transfer modeling in these systems as a practical tool for studying their performance is demonstrated, by characterizing their energy content and significantly contributing to the process of optimizing their design.

In the last decades, the focus of internal combustion engine development has moved towards more efficient and less pollutant engines. In a Diesel engine, approximately 30-40% of the energy provided by combustion is lost through the exhaust gases. The exhaust gases are hot and therefore rich of energy. Some of this energy can be recovered by recycling the exhaust gases into turbocharger. However, the energy losses in the exhaust port are highly undesired and the mechanisms driving the total pressure losses in the exhaust manifold not fully understood. Moreover, the efficiency of the turbine is highly dependent on the upstream flow conditions. Thus, a numerical study of the flow in the exhaust port geometry of a Scania heavy-duty Diesel engine is carried out mainly by using the Large Eddy Simulation (LES) approach. The purpose is to characterize the flow in the exhaust port, analyze and identify the sources of the total pressure losses. Unsteady Reynolds Averaged Navier-Stokes (URANS) simulation results are included for comparison purposes. The calculations are performed with fixed valve and stationary boundary conditions for which experimental data are available. The simulations include a verification study of the solver using different grid resolutions and different valve lift states. The calculated numerical data are compared to existent measured pressure loss data. The results show that even global parameters like total pressure losses are predicted better by LES than by URANS. The complex three-dimensional flow structures generated in the flow field are qualitatively assessed through visualization and analyzed by statistical means. The near valve region is a major source of losses. Due to the presence of the valve, an annular, jet-like flow structure is formed where the high-velocity flow follows the valve stem into the port. Flow separation occurs immediately downstream of the valve seat on the walls of the port and also on the surface of the valve body. Strong longitudinal, non-stationary secondary flow structures (i.e. in the plane normal to the main flow direction) are observed in the exhaust manifold. Such structures can degrade the efficiency of a possible turbine of a turbocharger located downstream on the exhaust manifold. The effect of the valve and piston motion has also been studied by the Large Eddy Simulation (LES) approach. Within the exhaust process, the valves open while the piston continues moving in the combustion chamber. This process is often analyzed modeling the piston and valves at fixed locations, but conserving the total mass flow. Using advanced methods, this process can be simulated numerically in a more accurate manner. Based on LES data, the discharge coefficients are calculated following the strict definition. The results show that the discharge coefficient can be overestimated (about 20 %) when using simplified experiments, e. g. flow bench. Simple cases using fixed positions for valve and piston are contrasted with cases which consider the motion of piston and/or valves. The overall flow characteristics are compared within the cases. The comparison shows it is impossible to rebuild the dynamic flow field with the simplification with fixed valves. It is better to employ LES to simulate the dynamic flow and associated losses with valve and piston motion. / <p>QC 20131204</p>

This work deals with spatially-evolving supersonic turbulent boundary layers over adiabatic and cold walls at M∞ = 2 and up to Re0 ≈ 2600 using 3 different SGS models. The numerical methodology is based on high-order split-centered scheme to discretize the convective fluxes of the Navier-Stokes equations . For the adiabatic case, it is demonstrated that all SGS models require a comparable minimum grid-refinement in order to capture accurately the near-wall-turbulence. Overall, the models exhibit correct behavior when predictiong the dynamic properties, but show different performances for the temperature distribution in the near-wall region. For the isothermal case, it is found that the compressibility effects are not enhanced due to the wall cooling. As expected, the total temperature fluctuations are not negligible in the near-wall region. The study shows that the anti-correlation linking both velocity and temperature fields, derived from the Morkovin's hypothesis, is not satisfied.

Tidal bores are surge waves propagating upstream rivers as the tide rushes into estuaries. They induce large turbulences and mixing of the river and estuary flow of which effects remain scarcely studied. Herein, tidal bores are investigated experimentally and numerically with an idealised model of positive surges propagating upstream an initially steady flow. The experimental work estimated flow changes and typical turbulent length scale evolution induced by undular bores with and without breaking roller. The bore passage was associated with large free surface and flow velocity fluctuations, together with some variations of the integral turbulent scales. Coherent turbulent structures appeared in the wake of leading wave near the bed and moved upward into the water column during the bore propagation. The numerical simulations were based on previous experimental work on undular bores. Some test cases were realised to verify the accuracy of the numerical methods. The results gave access to the detailed flow evolution during the bore propagation. Large velocity reversals were observed close to the no-slip boundaries. In some configurations, coherent turbulent structures appeared against the walls in the wake of the bore front.

Large-eddy simulation (LES) and unsteady Reynolds-averaged Navier-Stokes (URANS) calculations have been performed to investigate the effects of different mathematical models for scalar variance and its dissipation rate as applied to both a non-reacting bluff-body turbulent flow and an extension to a reacting case. In the conserved scalar formalism, the mean value of a thermo-chemical variable is obtained through the PDF-weighted integration of the local description over the conserved scalar, the mixture fraction. The scalar variance, one of the key parameters for the determination of a presumed β-function PDF, is obtained by solving its own transport equation with the unclosed scalar dissipation rate modelled using either an algebraic expression or a transport equation. The proposed approach is first applied to URANS and then extended to LES. Velocity, length and time scales associated with the URANS modelling are determined using the standard two-equation k-ε transport model. In contrast, all three scales required by the LES modelling are based on the Smagorinsky subgrid scale (SGS) algebraic model. The present study proposes a new algebraic and a new transport LES model for the scalar dissipation rate required by the transport equation for scalar variance, with a time scale consistent with the Smagorinsky SGS model.