L'étude porte sur le contrôle actif du phénomène de tremblement transsonique sur une voilure supercritique par l'utilisation de deux stratégies de contrôle : thermique et pneumatique. Des simulations numériques de l'écoulement naturel obtenues par résolution des équations de Navier-Stokes moyennées ont permis la restitution des différents régimes (avec ou sans tremblement) en bonne cohérence avec les données expérimentales. Les observations des champs moyens, des champs turbulents et des signaux de pression permettent d'établir le principe de fonctionnement et les performances des actionneurs. Le tremblement transsonique est alors soit repoussé soit supprimé autorisant ainsi l'élargissement du domaine de vol des aéronefs. L'étude paramétrique des deux actionneurs est réalisée aussi bien en mode de fonctionnement continu que rétroactif. Une campagne d'essai consolide et approfondie les résultats obtenues avec la stratégie pneumatique.

This thesis focused on investigating flow, coherent structures, and mechanisms leading to scour around a series of three submerged groynes in a high curvature (radius of curvature (R)/channel width (B)=1.5) channel bend using a Large Eddy Simulation Numerical (LES) model. Flow was investigated during both an initial and a later stage of scour. The results showed that the groynes appeared effective in keeping the main core of high streamwise velocity away from the outer bank wall in the region where they were installed, although high potential still existed for local scour around the groynes. During the initial stage of scour, horseshoe vortices (HVs) showed the greatest propensity to induce scour immediately upstream of the groyne tips. During the later stage of scour, the HV in front of the first upstream groyne (G1) induced very high mean pressure fluctuations on the outer bank wall. Scour was also of very great concern around the tip of G1 due to severe mean bed pressure fluctuations. Downstream of the groyne field, the presence of a counter-rotating outer bank cell was capable of endangering the stability of the outer bank. The second focus of this thesis was to investigate flow in a 135° channel bend using both Reynolds Averaged Navier Stokes (RANS) and LES numerical models. The RANS study examined the effects of curvature ratio (R/B), and aspect ratio (B/H, where H is the inlet flow depth), on secondary circulation strength, and bed shear stresses. The study revealed that a decrease in R/B was associated with an increase in secondary circulation strength and peak bed shear stress. A change in B/H also substantially affected cross stream circulation strength. The LES study was conducted in a 135° (R/B = 1.5) bend flume with a fixed bed corresponding to near equilibrium scour conditions, and the results were compared to a similar high curvature 193° bend numerical study. Inner bank vortices and shear layers were present in both cases although their characteristics were substantially different. Distributions of boundary friction velocities, and turbulence were also quite different for each case.

This thesis presents the results of steady, Reynolds-averaged Navier-Stokes (RANS) computations for jet flow emanating from a generic turbofan engine exhaust. All computations were performed with commercial solver FLUENT v6.2.16. Different turbulence models were evaluated. In addition to turbulence modeling issues, a parametric study was considered. Different modeling approaches for turbulent jet flows were explained in brief, with specific attention given to the Reynolds-averaged Navier-Stokes (RANS) method used for the calculations. First, a 2D ejector problem was solved to find out the most appropriate turbulence model and solver settings for the jet flow problem under consideration. Results of one equation Spalart-Allmaras, two-equation standart k-&amp / #949 / , realizable k-&amp / #949 / , k-&amp / #969 / and SST k-&amp / #969 / turbulence models were compared with the experimental data provided and also with the results of Yoder [21]. The results of SST k-&amp / #969 / and Spalart-Allmaras turbulence models show the best agreement with the experimental data. Discrepancy with the experimental data was observed at the initial growth region of the jet, but further downstream calculated results were closer to the measurements. Comparing the flow fields for these different turbulence models, it is seen that close to the onset of mixing section, turbulence dissipation was high for models other than SST k-&amp / #969 / and Spalart-Allmaras turbulence models. Higher levels of turbulent kinetic energy were present in the SST k-&amp / #969 / and Spalart-Allmaras turbulence models which yield better results compared to other turbulence models. The results of 2D ejector problem showed that turbulence model plays an important role to define the real physics of the problem. In the second study, analyses for a generic, subsonic, axisymmetric turbofan engine exhaust were performed. A grid sensitivity study with three different grid levels was done to determine grid dimensions of which solution does not change for the parametric study. Another turbulence model sensitivity study was performed for turbofan engine exhaust analysis to have a better understanding. In order to evaluate the results of different turbulence models, both turbulent and mean flow variables were compared. Even though turbulence models produced much different results for turbulent quantities, their effects on the mean flow field were not that much significant. For the parametric study, SST k-&amp / #969 / turbulence model was used. It is seen that boundary layer thickness effect becomes important in the jet flow close to the lips of the nozzles. At far downstream regions, it does not affect the flow field. For different turbulent intensities, no significant change occurred in both mean and turbulent flow fields.

Hydrodynamic behavior in reactors used for water treatment, particularly in ozone contactors with serpentine flow, is known to strongly affect the process efficiency. However, exact flow characteristics inside these reactors are not well understood, as traditional approach either considers these reactors as black box or relies on less accurate Reynolds-Averaged Navier-Stokes (RANS) simulation. In order to provide a deep understanding of the hydrodynamics and solute transport phenomena in these reactors, high resolution numerical studies using the Large Eddy Simulation (LES) method are performed. The reactor geometries investigated in this research are Constant Baffle Spacing Multi-Chamber (CBSMC) ozone contactors and a Variable Baffle Spacing ozone contactor Model (VBSM). The LES results in two multi-chamber ozone contactors (CBSMC -Normal-Width and -Half-Width) suggest that the flow through these reactors is characterized by the presence of extensive short-circuiting and large internal recirculation. The results also suggest that the flow is highly three dimensional with a pair of symmetric counter-rotating secondary vortices. LES studies based on VBSM, the baffle spacing of which varies between 0.5 times to 5 times the size of the base chamber; suggest that the width of the recirculation zone grows at about the same rate as the baffle spacing. Instantaneous turbulent eddies are prevalent in the chamber and increase turbulent mixing. The elevated levels of turbulence are found in the short-circuiting flow path. The tracer is dispersed along the short-circuiting path and strongly into the recirculation zone due to turbulent diffusion. Baffle spacing greater than the entrance gate height, but also smaller baffle spacing, worsens the disinfection efficiency. Finally, the turbulent Schmidt number of RANS simulation was investigated by employing the previously validated LES simulation. Due to the presence of very strong turbulent diffusion in the reactors, the turbulent Schmidt number is found out to be much less than the values commonly used, and is also specific to the baffle spacing.

In this thesis we describe the development of Reynolds- averaged Navier Stokes code for the flow past two- dimensional configuration. Particularly, emphasis has been laid on the study of low Reynolds number airfoil aerodynamics. The thesis consists of five chapters covering the back ground history, problem formulation, method of solution and discussion of the results and conclusion. Chapter I deals with a detailed background history of low Reynolds number aerodynamics, problem associated with it, state of the art, its importance in practical applications in aircraft industries. Chapter II describes the mathematical model of the flow physics and various levels of approximations. Also it gives an account of complexity of the equations at low Reynolds number regarding flow separation, transition and reattachment. Chapter III describes method of solution, numerical algorithm developed, description of various upwind schemes, grid system, finite volume discrieti-zation of the governing equations described in Chapter II. Chapter IV describes the application of the newly developed Navier Stokes code for the test cases from GAMM Workshop proceedings. Also it describes validation of the code for Euler solutions, Blasius solution for the flow past flat plate and compressible Navier Stokes solution for the flow past NACA 0012 Airfoil at low Reynolds number. Chapter V describes the application of the Navier Stokes code for the more test cases of current practical interest . In this chapter laminar separation bubble characteristics are investigated in detail regarding formation, growth and shedding in an unsteady environment. Finally the conclusion is drawn regarding the robustness of the newly developed code in predicting the airfoil aerodynamic characteristics at low Reynolds number both in steady and unsteady environment. Lastly, suggestion for future work has been highlighted.

Vegetation is commonly modelled as emergent arrays of rigid, circular cylinders. However, the drag coefficient (CD) of real stems or trunks is closer to that of cylinders with a square cross-section. In this thesis, vegetation has been idealised as square cylinders in laboratory experiments with a turbulence intensity of the order of 10% which is similar to that of typical river flows. These cylinders may also represent other obstacles such as architectural structures. This research has determined CD of an isolated cylinder and cylinder pairs as a function of position as well as the average drag coefficient (CDv) of larger arrays. A strain gauge was used to measure CD whilst CDv was computed with a momentum balance which was validated by strain gauge measurements for a regularly spaced array. The velocity and turbulence intensity surrounding a pair of cylinders arranged one behind the other with respect to mean flow (in tandem) were also measured with an Acoustic Doppler Velocimeter. The isolated cylinder CD was found to be 2.11 in close agreement with other researchers. Under fixed flow conditions CD for a cylinder in a pair was found to be as low as -0.40 and as high as 3.46 depending on their relative positioning. For arrays, CDv was influenced more by the distribution of cylinders than the flow conditions over the range of conditions tested. Mean values of CDv for each array were found to be between 1.52 and 3.06. This new insight therefore suggests that CDv for vegetation in bulk may actually be much higher than the typical value of 1 which is often assumed to apply in practice. If little other information is available, a crude estimate of CDv = 2 would be reasonable for many practical applications. The validity of a 2D realizable k-epsilon turbulence model for predicting the flow around square cylinders was evaluated. The model was successful in predicting CD for an isolated cylinder. In this regard the model performed as well as Large Eddy Simulations by other authors with a significant increase in computational efficiency. However, the numerical model underestimates CD of downstream cylinders in tandem pairs and overestimates velocities in their wake. This suggests it may be necessary to expand the model to three-dimensions when attempting to simulate the flow around two or more bluff obstacles with sharp edges.

This thesis is concerned with studies on spatially developing high speed mixing layers with twin objectives: (a) to provide enhanced and detailed understanding of spatial development of two-dimensional mixing layer emanating from splitter plate through large eddy simulation (LES, from now on) technique and (b) to evolve a consistent strategy for Unsteady Reynolds Averaged Navier-Stokes (URANS) approach to mixing layer calculations. The inspiration for this work arose out of the explanations that were being developed for the reduction in the mixing layer thickness with compressibility (measured by a parameter called convective Mach number, Mc). The reasons centered around increased stability, increase in compressible dissipation that was later discounted in favor of reduction in production and pressure-strain terms (with Mc, of course). These were obtained with direct numerical simulations (DNS) or LES techniques with homogeneous shear flow or temporal mixing layer. As apart, there was also a wide held view that using RANS (steady) techniques did not capture the compressibility effects when used in a way described above and so classical industrial codes for computing mixing- layer-embedded flows are unsuitable for such applications. Other important aspects that come out of the examination of literature are: the mixing layer growth is controlled in the initial stages by the double- boundary layer profile over the splitter plate and results in the mixing layer growth that is somewhat irregular due to doubling and merging of vertical structures. The view point of a smooth growth of the mixing layer is a theo- retical approximation arising out of the use of a smooth tan-hyperbolic profile that results at larger distances from the splitter plate. For all practical applications, it is inferred that the initial development is what is important because the processes of ignition and stable combustion occur close to the splitter plate. For these reasons, it was thought that understanding the development of the mixing layer is best dealt with using accurate spatial simulation with the appropriate initial profile. The LES technique used here is drawn from an OpenFOAM approach for dissimilar gases and uses one-equation Eddy Model for SGS stresses. The temporal discretization is second order accurate backward Euler and spatial discretization is fourth order least squares; the algorithm used for solving the equations is PISO and the parallelized code uses domain decomposition approach to cover large spatial domain. The calculations are performed with boundary layer profiles over the splitter plate and an initial velocity field with white noise-like fluctuations to simulate the turbulence as in the experiments. Grid independence studies are performed and several experimental cases are considered for comparison with measured data on the velocity and temperature fields as well as turbulent statistics. These comparisons are excellent for the mean field behavior and moderately acceptable for turbulent kinetic energy and shear stress. To further benefit from the LES approach, the details of the mixing layer are calculated as a function of four independent parameters on which the growth depends: convective Mach number (Mc = (U1 -U2)/ (a1 +a2)), stream speed ratio (r = U2=U1), stream density ratio (s = p2/p1) and the average velocity of the two streams ((U1+U2)=2) and examine the various terms in the equations to enable answering the questions discussed earlier. It is uncovered that r has significant influence on the attainment of self similarity (which also implies on the rate of removal of velocity defect in the double-boundary layer profile) and other parameters have a very weak influence. The minimum velocity variation with distance from the splitter plate has the 1/paxial distance behavior like in wakes; however, after a distance, departure to linear rise occurs and the distance it takes for this to appear is delayed with Mc. Other features such as the coherent structures, their merger or break up, the area of the structures, convective velocity information extraction from the coherent structures, the behavior of the pressure field in the mixing layer through the field are elucidated in detail; the behavior of the correlations between parameters (like pressure, velocity etc) at different points is used to elucidate the coherence of their fluctuating field. The effects of the parameters on the energy spectra have expected trends. An examination of the kinetic energy budget terms reveals that • the production term is the main source of the xx turbulence stress, whereas it is not significant in the yy component. • A substantial portion of this is carried by the pressure-velocity coupling from the xx direction to the yy direction, which becomes the main source term in the yy component. • Both, the production term as well as the pressure-velocity term show a clear decrease with increase in Mc. The high point of the thesis is related to using the understanding derived from an analysis of various source terms in the kinetic energy balance to evolve an unsteady Reynolds Averaged Navier Stokes (URANS) model for calculating high speed mixing layers, a subject that has eluded international research till now. It recognizes that the key feature affected by ompressibility is related to the anisotropy of the stress tensor. The relationship between stress component (_Txy) and the velocity gradient (Sxy) as obtained from LES is set out in the form of a simple relationship accounting for the effects of other parameters obtained earlier in this thesis. A minor influence due to _Tyy is extracted by describing its dependence on Sxy again as gleaned from LES studies. The needed variation of Prandtl and Schmidt numbers through the field is extracted. While the detailed variations can in fact be taken into account in URANS simulations, a simple assumption of these values being around 0.3 is chosen for the present simulations of URANS. Introduction of these features into the momentum equation gives the much expected variation of the reduction in the growth rate of the mixing layer with convective Mach number as in experiments. The relationships that can be used in high speed mixing layers are Introduction of these features into the momentum equation gives the much expected variation of the reduction in the growth rate of the mixing layer with convective Mach number as in experiments. This is then a suggested new approach to solve high speed mixing layers. While it can be thought that the principal contributions of the thesis are complete here, an additional segment is presented related to entropy view of the mixing layer. This study that considers the mixing layer with two different species expresses various terms involved in the entropy conservation equation and obtains the contribution of various terms on the entropy change for various Mc. It is first verified that the entropy derived from the conservation equation matches with those calculated from fluid properties, entropy being a state variable. It is shown that irreversible diffusion comes down the most with convective Mach number. Left: This image shows pictorially the flow of source of turbulent stress from the axial to the cross wise turbulent stress. Production (Σ) of turbulence happens mainly in the xx direction, a part of it is carried by the pressure-velocity correlation to the yy direction, which itself has a low production. With increasing Mc, both the production as well as the pressure-velocity correlation decrease. Right: This image shows the growth rate obtained from simulations scaled with the incompressible growth rate, of LES and RANS in the background of experiments (others). As is clear, the growth rate obtained is well within the band of experimental results.

This thesis focused on investigating flow, coherent structures, and mechanisms leading to scour around a series of three submerged groynes in a high curvature (radius of curvature (R)/channel width (B)=1.5) channel bend using a Large Eddy Simulation Numerical (LES) model. Flow was investigated during both an initial and a later stage of scour. The results showed that the groynes appeared effective in keeping the main core of high streamwise velocity away from the outer bank wall in the region where they were installed, although high potential still existed for local scour around the groynes. During the initial stage of scour, horseshoe vortices (HVs) showed the greatest propensity to induce scour immediately upstream of the groyne tips. During the later stage of scour, the HV in front of the first upstream groyne (G1) induced very high mean pressure fluctuations on the outer bank wall. Scour was also of very great concern around the tip of G1 due to severe mean bed pressure fluctuations. Downstream of the groyne field, the presence of a counter-rotating outer bank cell was capable of endangering the stability of the outer bank. The second focus of this thesis was to investigate flow in a 135° channel bend using both Reynolds Averaged Navier Stokes (RANS) and LES numerical models. The RANS study examined the effects of curvature ratio (R/B), and aspect ratio (B/H, where H is the inlet flow depth), on secondary circulation strength, and bed shear stresses. The study revealed that a decrease in R/B was associated with an increase in secondary circulation strength and peak bed shear stress. A change in B/H also substantially affected cross stream circulation strength. The LES study was conducted in a 135° (R/B = 1.5) bend flume with a fixed bed corresponding to near equilibrium scour conditions, and the results were compared to a similar high curvature 193° bend numerical study. Inner bank vortices and shear layers were present in both cases although their characteristics were substantially different. Distributions of boundary friction velocities, and turbulence were also quite different for each case.

This document is a report summering the master thesis work dealing with the Computational Fluid Dynamic (CFD) study of the flow around a high-speed train. The model is a scaled 1:50 generic train with two cars, one inter-car gap and simplified bogies. A platform is set on the side of the train since one of the aim of the study is to look at the consequences of the phenomena in the wake on people or objects standing on the platform. The slipstream is one of this phenomena, it is due to the fact that the viscous air is dragged when the train is passing. If too strong, it can move or destabilize people or objects on the platform. In addition of the slipstream study, a velocity profile study, a drag and lift coefficients analyze as well as a Q-factor study and a frequency study have been realized. Some results of these different studies are compared with the ones obtained on the same model with a Delayed Detached Eddy Simulation (DDES). Since the flow is turbulent, for those different studies, the flow has been simulated with a Reynolds Averaged Navier-Stokes equation model (RANS) which is the k-ω SST model for the turbulence. The study of the slipstream allowed to calculate the Technical Specification for Interoperability (TSI) which must not be higher that the European Union requirement set at 15.5 m/s, the result obtained is 8.1 m/s which is then lower than the limit. The velocity profile shows similarities with the DDES results even though it is less detailed. The same conclusion is done for the Q-plot where is clearly visible the two counter-rotating vortices in the wake. Finally, a Fast Fourier Transform algorithm has been applied to instantaneous velocity results in the wake of the train in order to get the frequency of the aerodynamic phenomena in that wake. The main frequency is 25 Hz and corresponds to a Strouhal number of 0.1, quite closed to the results obtained with DDES which is 0.085. The results of the RANS and DDES are reasonably similar and by regarding at the large difference between the cell numbers (respectively 8 500 000 and 20 000 000) it can be conclude that in some ways the RANS model can be preferred at the DDES to save time for the computation but it does not contain the small scales resolved by the DDES.

The aim of this master thesis is to characterize the fluid forces applied to a fuel assembly inthe core of a nuclear power plant in case of seism. The forces are studied with a simplifiedtwo-dimensional model constituted of an array of 3 by 3 infinite cylinders oscillating in aclosed box. The axial flow of water, which convects the heat in the core of a nuclear powerplant, is also taken into account. The velocity of the axial flow reaches 4m/s in the middle ofthe assembly and modifies the forces features when the cylinders move laterally.The seism is modeled as a lateral displacement with high amplitude (several cylinderdiameters) and low frequencies (below 20 Hz). In order to study the effects of the amplitudeand of the frequency of the displacement, the displacement taken is a sine function withboth controlled amplitude and frequency. Four degrees of freedom of the system will bestudied: the amplitude of the displacement, its frequency, the axial velocity amplitude andthe confinement (due to the closed box).The fluid forces exerted on the cylinders can be seen as a combination of three terms: anadded mass, related to the acceleration of cylinders, a drift force, related to the damping ofthe fluid and a force due to the interaction of the cylinder with residual vortices. The firsttwo components will be characterized through the Morison expansion, and their evolutionwith the variation of the degree of freedom of the system will be quantified. The effect ofthe interaction with the residual vortices will be observed in the plots of the forces vs. timebut also in the velocity and vorticity map of the fluid.The fluid forces are calculated with the CFD code Code_Saturne, which uses a second orderaccurate finite volume method. Unsteady Reynolds Averaged Navier Stokes simulations arerealized with a k-epsilon turbulence model. The Arbitrary Lagrange Euler model is used todescribe the structure displacement. The domain is meshed with hexahedra with thesoftware gmsh [1] and the flow is visualized with Paraview [2]. The modeling techniquesused for the simulations are described in the first part of this master thesis.