Dans les centrales solaires à tour, le flux solaire est concentré vers un récepteur solaire où son énergie est transférée à un fluide caloporteur. L'écoulement au sein du récepteur solaire est turbulent, fortement anisotherme et à bas nombre de Mach. L'optimisation du récepteur solaire exige une meilleure compréhension et modélisation de l'interaction entre la température et la turbulence. Cette thèse cherche à y contribuer selon deux approches. Tout d'abord, on étudie les échanges énergétiques entre les différentes parties de l'énergie totale. On propose pour cela une nouvelle représentation des échanges énergétiques, fondée sur la moyenne de Reynolds. Cette représentation permet la caractérisation, à partir de simulations numériques directes d'un canal plan bipériodique anisotherme, de l'effet du gradient de température sur les échanges énergétiques associées à l'énergie cinétique turbulente dans les domaines spatial et spectral. Ensuite, on étudie la simulation des grandes échelles des équations de bas nombre de Mach. En utilisant les résultats de simulations numériques directes, on identifie les termes sous-mailles spécifiques à modéliser lorsque l'on utilise le filtre classique, non pondéré, et lorsque l'on utilise le filtre de Favre, pondéré par la masse volumique. Dans les deux cas, on évalue a priori la performance de différents modèles sous-mailles. La pertinence des modèles est vérifiée a posteriori par la réalisation de simulation des grandes échelles. / In solar power towers, the solar flux is concentrated towards a solar receiver, wherethrough its energy is transferred to a heat transfer fluid. The flow in the solar receiver is turbulent, strongly anisothermal and at low Mach number. The optimisation of the solar receiver requires a better understanding and modelling of the interaction between temperature and turbulence. In this thesis, this is investigated following two approaches. First, we study the energy exchanges between the different parts of total energy. To this end, a new representation of the energy exchanges, based on the Reynolds averaging, is established. The representation allows the characterisation, from direct numerical simulations of a strongly anisothermal channel flow, of the effect of the temperature gradient on the energy exchanges associated with turbulence kinetic energy in the spatial and spectral domains. Second, we study the large-eddy simulation of the low Mach number equations. Using the results of direct numerical simulations, we identify the specific subgrid terms to model when the unweighted classical filter is used and when the density-weighted Favre filter is used. In both cases, the performance of different subgrid-scale models is assessed a priori. The relevance of the subgrid-scale models is then verified a posteriori by carrying out large-eddy simulations.

Cette étude porte sur le calcul numérique du champ acoustique rayonné par des écoulements subsoniques turbulents présentant des inhomogénéités de température. Des méthodes hybrides sont développées grâce à un développement de Janzen-Rayleigh des équations de Navier-Stokes. L'écoulement est résolu par un calcul quasi incompressible puis les perturbations acoustiques sont propagées selon deux méthodes : les équations d'Euler linéarisées (EEL) et l'approximation à faible nombre de Mach perturbée(PLMNA). Les méthodes sont validées sur des cas simples puis appliquées à une couche de mélange isotherme et anisotherme en développement spatial. / This study focuses on the numerical calculation of the acoustic field radiated by subsonic turbulent flows with temperature inhomogeneities. Hybrid methods are developed through a Rayleigh-Janzen expansion of the Navier-Stokes equations. The flow is solved in a quasi-incompressible way then the acoustic disturbances are propagated by two methods : the linearized Euler's equations (EEL) and the perturbed low Mach number approximation (PLMNA). The methods are validated on simple cases and then applied to an isothermal and non isothermal spatially evolving mixing layer.

The evaluation of the noise produced by flows has reached a high level of importance in the past years. The physics surrounding flow-induced noise is quite complex and sensitive to various flow conditions like temperature, shape. Empirical models were built in the past for some special geometries but they cannot be used in a general case for a shape optimization for instance. Experimental aeroacoustic facilities represent the main tool for acoustic analyses of flow fields, but are quite expensive because extreme care must be exercised not to introduce acoustic perturbations in the flow (silent facilities). These tools allow a good analysis of the physical phenomena responsible for noise generation in the flow by a comparison of the noise sources and the flow characteristics (pressure, turbulence,). Nevertheless, the identification and location of noise sources to compare with flow structures requires quite complex methods.<p><p>The numerical approach complements the experimental one in the sense that the flow characteristics are deeply analyzed where experiments suggest noise production. For the numerical approach, the turbulence modeling is quite important. In the past, some models were appreciated for their good prediction of some aerodynamic parameters as lift and drag for instance. The challenge is now to tune these models for a correct prediction of the noise sources. In the low subsonic range, the flow field is completely decoupled from acoustics, and noise sources can be computed from a purely hydrodynamic simulation before this information is transferred to an acoustical solver which will compute the acoustic field at the listener position. This post processing of the aerodynamic results is not obvious since it can introduce non-physical noise into the solution.<p><p>This project considers the aspect of noise generation in turbulent jets and especially the noise generated by vortex pairing, as it occurs for instance in jet flows. The axisymmetric version of the flow solver SFELES has been part of this PhD research, and numerical results obtained on the jet are similar to the experimental values. Analyses performed on the numerical results are interesting to go to complete turbulence modeling for aeroacoustics since vortex pairing is one of the basic acoustical processes in vortex dynamics.<p><p>Currently, a standard static Smagorinski model is used for turbulence modeling. However, this model has well known limitations, and its influence on the noise sources extracted from the flow field is not very clear. For this reason, it is planned to adopt a dynamic procedure in which the subgrid scale model automatically adapts to the flow. We planned also to perform simulations with the variational multiscale approach to better simulate the different interactions between large and unresolved scales. The commercial software ACTRAN distributed by Free Field Technologies is used for the computation of sound propagation inside the acoustic domain. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

No description available.

The hypersonic laminar to turbulent transition problem above Mach 10 is addressed experimentally in the short duration VKI Longshot gun tunnel. Reentry conditions are partially duplicated in terms of Mach and Reynolds numbers. Pure nitrogen is used as a test gas with flow enthalpies sufficiently low to avoid its dissociation, thus approaching a perfect gas behavior. The stabilizing effects of Mach number and nosetip bluntness on the development of natural boundary layer disturbances are evaluated over a 7 degrees half-angle conical geometry without angle of attack. <p><p>Emphasis is initially placed on the flow characterization of the Longshot wind tunnel where these experiments are performed. Free-stream static pressure diagnostics are implemented in order to complete existing stagnation point pressure and heat flux measurements on a hemispherical probe. An alternative method used to determine accurate free-stream flow conditions is then derived following a rigorous theoretical approach coupled to the VKI Mutation thermo-chemical library. Resulting sensitivities of free-stream quantities to the experimental inputs are determined and the corresponding uncertainties are quantified and discussed. The benefits of this different approach are underlined, revealing the severe weaknesses of traditional methods based on the measurement of reservoir conditions and the following assumptions of an isentropic and adiabatic flow through the nozzle. The operational map of the Longshot wind tunnel is redefined accordingly. The practical limits associated with the onset of nitrogen flow condensation under non-equilibrium conditions are also accounted for. <p><p>Boundary layer transition experiments are then performed in this environment with free-stream Mach numbers ranging between 10-12. Instrumentation along the 800mm long conical model includes flush-mounted thermocouples and fast-response pressure sensors. Transition locations on sharp cones compare favorably with engineering correlations. A strong stabilizing effect of nosetip bluntness is reported and no transition reversal regime is observed for Re_RN<120000. Wavelet analysis of wall pressure traces denote the presence of inviscid instabilities belonging to Mack's second mode. An excellent agreement with Linear Stability Theory results is obtained from which the N-factor of the Longshot wind tunnel in these conditions is inferred. A novel Schlieren technique using a short duration laser light source is developed, allowing for high-quality flow visualization of the boundary layer disturbances. Comparisons of these measurement techniques between each other are finally reported, providing a detailed view of the transition process above Mach 10. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Les phénomènes de transport en milieux poreux sont étudiés depuis près de deux siècles, cependant les travaux concernant les milieux fortement poreux sont encore relativement peu nombreux. Les modèles couramment utilisés pour les poreux classiques (lits de grains par exemple) sont peu applicables pour les milieux fortement poreux (les mousses par exemple), un certain nombre d’études ont été entreprises pour combler ce manque. Néanmoins, les résultats expérimentaux et numériques caractérisant les pertes de charge dans les mousses sont fortement dispersés. Du fait des progrès de l’imagerie 3D, une tendance émergente est la détermination des paramètres des lois d’écoulement à partir de simulations directes sur des géométries reconstruites. Nous présentons ici l’utilisation d’une nouvelle approche cinétique pour résoudre localement les équations de Navier-Stokes et déterminer les propriétés d’écoulement (perméabilité, dispersion, ...). / A novel kinetic scheme satisfying an entropy condition is developed, tested and implemented for the simulation of practical problems. The construction of this new entropic scheme is presented. A classical hyperbolic system is approximated by a discrete velocity vector kinetic scheme (with the simplified BGK collisional operator), but applied to an inviscid compressible gas dynamics system with a small Mach number parameter, according to the approach of Carfora and Natalini (2008). The numerical viscosity is controlled, and tends to the physical viscosity of the Navier-Stokes system. The proposed numerical scheme is analyzed and formulated as an explicit finite volume flux vector splitting (FVS) scheme that is very easy to implement. It is close in spirit to Lattice Boltzmann schemes, but it has the advantage to satisfy a discrete entropy inequality under a CFL condition and a subcharacteristic stability condition involving a cell Reynolds number. The new scheme is proved to be second-order accurate in space. We show the efficiency of the method in terms of accuracy and robustness on a variety of classical benchmark tests. Some physical problems have been studied in order to show the usefulness of both schemes. The LB code was successfully used to determine the longitudinal dispersion of metallic foams, with the use of a novel indicator. The entropic code was used to determine the permeability tensor of various porous media, from the Fontainebleau sandstone (low porosity) to a redwood tree sample (high porosity). These results are pretty accurate. Finally, the entropic framework is applied to the advection-diffusion equation as a passive scalar.

Zapfluft und Wellenleistung wird den Triebwerken entnommen, um die Energie für beispielsweise die Kraftstoffpumpen, das Inflight Entertainment oder die Flügelvorderkantenenteisung zu erzeugen. Diese Energiegenerierung, hat einen Anstieg des Kraftstoffverbrauches zur Folge. Es hat sich herausgestellt, dass die Stelle der Zapfluftentnahme einen starken Einfluss auf den Gradienten des Brennstoffverbrauches hat. Das Projekt beschäftigt sich mit zwei- und dreiwelligen Turbofantriebwerken und untersucht an ihnen, die Effekte der Leistungsnahmen. Als Simulationssoftware wurde GasTurb 8.0 eingesetzt und auf die integrierten Triebwerkskonfigurationen zurückgegriffen. Ziel der Arbeit ist die Ermittlung einer mathematischen Beziehung zur Berechnung des zusätzlichen Kraftstoffmassenstromes infolge einer Zapfluft- oder Wellenleistungsentnahme. So stellt sich die Frage, welche Triebwerksparameter dafür berücksichtigt werden müssen. Eine Wellenleistungsentnahme verursacht beispielsweise einen linearen Anstieg des spezifischen Kraftstoffverbrauches. Ist diese Zunahme, identisch mit der einer Zapfluftentnahme? Am Ende der Kapitel werden die Ergebnisse mit Literaturwerten verglichen und versucht Tendenzen zu erkennen bzw. bestehende zu erhärten.

L’objectif de ce travail de thèse est de proposer la Méthodes des éléments finis de Galerkin discontinus (DGFEM) à la discrétisation des équations compressibles de Navier-Stokes. Plusieurs challenges font l’objet de ce travail. Le premier aspect a consisté à montrer l’ordre de convergence optimal de la méthode DGFEM en utilisant les polynômes d’interpolation d’ordre élevé. Le deuxième aspect concerne l’implémentation de méthodes de ‘‘shock-catpuring’’ comme les limiteurs de pentes et les méthodes de viscosité artificielle pour supprimer les oscillations numériques engendrées par l’ordre élevé (lorsque des polynômes d’interpolation de degré p>0 sont utilisés) dans les écoulements transsoniques et supersoniques. Ensuite nous avons implémenté des estimateurs d’erreur a posteriori et des procédures d ’adaptation de maillages qui permettent d’augmenter la précision de la solution et la vitesse de convergence afin d’obtenir un gain de temps considérable. Finalement, nous avons montré la capacité de la méthode DG à donner des résultats corrects à faibles nombres de Mach. Lorsque le nombre de Mach est petit pour les écoulements compressibles à la limite de l’incompressible, la solution souffre généralement de convergence et de précision. Pour pallier ce problème généralement on procède au préconditionnement qui modifie les équations d’Euler. Dans notre cas, les équations ne sont pas modifiées. Dans ce travail, nous montrons la précision et la robustesse de méthode DG proposée avec un schéma en temps implicite de second ordre et des conditions de bords adéquats. / The aim of this thesis is to deal with compressible Navier-Stokes flows discretized by Discontinuous Galerkin Finite Elements Methods. Several aspects has been considered. One is to show the optimal convergence of the DGFEM method when using high order polynomial. Second is to design shock-capturing methods such as slope limiters and artificial viscosity to suppress numerical oscillation occurring when p>0 schemes are used. Third aspect is to design an a posteriori error estimator for adaptive mesh refinement in order to optimize the mesh in the computational domain. And finally, we want to show the accuracy and the robustness of the DG method implemented when we reach very low mach numbers. Usually when simulating compressible flows at very low mach numbers at the limit of incompressible flows, there occurs many kind of problems such as accuracy and convergence of the solution. To be able to run low Mach number problems, there exists solution like preconditioning. This method usually modifies the Euler. Here the Euler equations are not modified and with a robust time scheme and good boundary conditions imposed one can have efficient and accurate results.