The present thesis is an experimental study of some aspects of unsteady two dimensional boundary layers subject to adverse pressure gradient. An adverse pressure gradient usually leads to boundary layer separation or an instability which may result in transition to turbulence. Unsteady boundary layer separation is not yet fully understood and there is no specific criterion proposed in literature for its occurrence. The details of separation depend on the Reynolds number, the geometry of the body (streamlined or bluff) and the type of imposed unsteady motion (impulsive, oscillatory etc.). Similarly there are many unknowns with respect to instability and transition in unsteady boundary layers, especially those having a streamwise variation. For unsteady flows it is useful to break up the pressure gradient term in the unsteady boundary layer equation into two components:(Formula) is the velocity at the edge of the boundary layer. The first term of the right hand side of this equation may be called the temporal component (Πt) which signifies acceleration or deceleration in time of the free stream and the second term is the spatial component (Πx) which represents the spatial or convective acceleration of the free stream. Many of the studies on instability in unsteady flows found in literature are carried out in straight tubes or channels, where the Πx term is absent. However, in many cases, especially in biological systems both terms are present. An example is the unsteady flow over the moving body of a fish. To study the effects of Πt and Πx on unsteady separation and instability we have built an unsteady water tunnel where the two components can be systematically varied. The flow is created by a controlled motion of a piston. By a suitable combination of the geometry of the model and the piston motion, different types of separation bubbles may be generated. In our studies the piston motion follows a trapezoidal variation: constant acceleration from rest, followed by constant velocity and then deceleration to zero velocity. We have chosen two geometries. One is a bluff body and thus has a high value of Πx and other is a small angle diffuser with a divergence angle 6.2° and thus having a small value of Πx. Upstream and downstream of the diffuser are long lengths of constant cross section. We have performed experiments with the above mentioned geometries placed in the tunnel test section. Flow is visualized using the laser induced fluorescence technique by injecting a thin layer of fluorescein dye on the test wall. Numerical simulations have been done using the software FLUENT. Boundary layer parameters like boundary layer, displacement and momentum thicknesses are calculated from the simulations and used to analyze the experimental results. For the flow in the diffuser, quasi-steady stability analysis of the instantaneous velocity profiles gives a general idea of stability behavior of the flow. Two types of experiments have been done with the bluff body. One is the unsteady boundary layer separation and the formation of the initial vortex for a flow that is uniformly accelerated from rest. We have found some scalings for the formation time (tv) of the separation vortex. The second type of experiment was to study the vortex shedding from the separating shear layer after the boundary layer has fully separated. At high enough Reynolds number shear layer vortices are seen to shed from the separation bubble. The Strouhal number based on the momentum thickness and the velocity at the edge of the boundary layer just upstream of the separation point is found to vary between 0.004 and 0.008. This value is close to the Strouhal number value of 0.0068 found in laminar separation bubbles on a flat plate. The second part of the study concerns with the evolution of the flow in the small angle diffuser with a mild variation of the spatial component of the pressure gradient. From the experimental visualizations we have found that the ratio of Πx and Πt at the start of the deceleration phase of the piston motion is an important parameter that determines the type of instability. This value of Πx/Πt is controlled by controlling the piston deceleration: a large deceleration gives a low Πx/Πt value and a low deceleration gives a large Πx/Πt value. Three types of instabilities have been observed in our experiments. In Type I, the first vortex forms at the maximum pressure gradient point (MPGP) and which grows disproportionately with time. However, instability vortices are seen later at other locations around the MPGP. In type II an array of vortices over a certain length are observed; the vortices grow with time. In Type III, which we observe for low decelerations, we observe initial vortices only in the diffuser section in the deceleration phase of the piston motion. Type III instability is similar to the one observed in dynamic stall experiments. In all cases the instability is very localized - it occurs only over some length of the boundary layer. Transition to turbulence, which is also localized, is observed at higher Reynolds numbers. The non-dimensional time for vortex formation is not very different from that found in straight channel experiments. Quasi-steady linear stability analyses for the boundary layer at the MPGP both for the top and the bottom walls show that the flow is absolutely unstable for some cases. In summary, the thesis looks at in a unified way the separation and instability of unsteady boundary layers with reverse flow. It is hoped that the results will be useful in predicting and understanding onset of separation and instability in practically occurring unsteady flows.

This thesis deals with the near-field of confluent jets, which can be of interest in many engineering applications such as design of a ventilation supply device. The physical effect of interaction between multiple closely spaced jets is studied using experimental and numerical methods. The primary aim of this study is to explore a better understanding of flow and turbulence behavior of multiple interacting jets. The main goal is to gain an insight into the confluence of jets occurring in the near-field of multiple interacting jets. The array of multiple interacting jets is studied when they are placed on a flat and a curved surface. To obtain the boundary conditions at the nozzle exits of the confluent jets on a curved surface, the results of numerical prediction of a cylindrical air supply device using two turbulence models (realizable 𝑘 − 𝜖 and Reynolds stress model) are validated with hot-wire anemometry (HWA) near different nozzles discharge in the array. A single round jet is then studied to find the appropriate turbulence models for the prediction of the three-dimensional flow field and to gain an understanding of the effect of the boundary conditions predicted at the nozzle inlet. In comparison with HWA measurements, the turbulence models with low Reynolds correction (𝑘 − 𝜖 and shear stress transport [SST] 𝑘 − 𝜔) give reasonable flow predictions for the single round jet with the prescribed inlet boundary conditions, while the transition models (𝑘 − 𝑘l − 𝜔𝜔 and transition SST 𝑘 − 𝜔) are unable to predict the flow in the turbulent region. The results of numerical prediction (low Reynolds SST 𝑘 − 𝜔 model) using the prescribed inlet boundary conditions agree well with the HWA measurement in the nearfield of confluent jets on a curved surface, except in the merging region. Instantaneous velocity measurements are performed by laser Doppler anemometry (LDA) and particle image velocimetry (PIV) in two different configurations, a single row of parallel coplanar jets and an inline array of jets on a flat surface. The results of LDA and PIV are compared, which exhibit good agreement except near the nozzle exits. The streamwise velocity profile of the jets in the initial region shows a saddle back shape with attenuated turbulence in the core region and two off-centered narrow peaks. When confluent jets issue from an array of closely spaced nozzles, they may converge, merge, and combine after a certain distance downstream of the nozzle edge. The deflection plays a salient role for the multiple interacting jets (except in the single row configuration), where all the jets are converged towards the center of the array. The jet position, such as central, side and corner jets, significantly influences the development features of the jets, such as velocity decay and lateral displacement. The flow field of confluent jets exhibits asymmetrical distributions of Reynolds stresses around the axis of the jets and highly anisotropic turbulence. The velocity decays slower in the combined regio  of confluent jets than a single jet. Using the response surface methodology, the correlations between characteristic points (merging and combined points) and the statistically significant terms of the three design factors (inlet velocity, spacing between the nozzles and diameter of the nozzles) are determined for the single row of coplanar parallel jets. The computational parametric study of the single row configuration shows that spacing has the greatest impact on the near-field characteristics.

Ce travail de thèse s’inscrit dans le cadre du programme R&D du CEA en support au système de conversion d’énergie à gaz du prototype industriel de Réacteur à Neutrons Rapides refroidi au Sodium (RNR-Na). Cette technologie représente une alternative aux cycles Rankine conventionnels à eau/vapeur, ayant pour avantage principal l’élimination du scenario accidentel de réaction sodium-eau. Cependant, la faible capacité de transfert de chaleur du gaz nécessite une technologie d’échangeurs compacts à plaques avec un nombre élevé de canaux à alimenter. Coté sodium, une section minimale de passage est nécessaire pour éviter le risque de bouchage par impureté. Cela induit de très faibles pertes de pression dans le faisceau qui, couplées à une condition de vitesse élevée à l’entrée, génèrent un risque réel de mauvaise distribution du débit. Les performances d’échange thermique et la tenue mécanique du composant sont alors dégradées. L’objectif principal de ce travail de thèse a été de résoudre ce problème de mauvaise distribution, en s’appuyant sur une conception innovante (BREVET FR16 57543), sur une stratégie de calcul numérique et l’établissement d’une base de données expérimentale pour la validation des travaux théoriques. Le nouveau système de distribution sodium se compose d’un collecteur d'entrée dont le design permet de guider la trajectoire du jet et d’un système de bifurcation de canaux qui augmente les pertes de pression dans le faisceau. De plus, des communications latérales entre les canaux sodium aident à homogénéiser davantage le flux. Deux installations expérimentales ont été conçues pour caractériser l'écoulement dans les canaux de bifurcation et dans le collecteur d'entrée. La conception des maquettes a permis de quantifier leur effet sur la distribution du flux entre les canaux. La base de données aérodynamiques PIV acquises a permis de valider les modèles numériques et de prouver l’efficacité du système de distribution proposé. Après avoir validé les modèles de turbulence CFD et la stratégie d'étude de la distribution dans le module SGHE, une optimisation de chaque composant du système de distribution de sodium a été réalisée. Le travail de cette thèse s’achève par la description de la conception optimale retenue pour la phase actuelle du projet ASTRID. / This PhD work was motivated by the CEA R&D program to provide solid technological basis for the use of Brayton power conversion system in Sodium-cooled Fast nuclear Reactors (SFRs). Multi-channel compact heat exchangers are necessary for the present application because of the low heat transfer capacity of the gas foreseen. In ASTRID project, a minimum size of Na channels section is required to avoid the plugging risk. However, this induces very low pressure losses in the bundle. Considering an additional inlet flow condition, a real risk of bad flow distribution remains. As a result, the thermal performance and thermal loading of the heat exchanger degrades due to it. The main goal of this work was to overcome the flow maldistribution problem by means of an innovative design of sodium distribution system (PATENT FR1657543), the development of a numerical strategy and the construction of an experimental database to validate all theoretical studies. The innovative sodium distribution system consists on an inlet header which tries to guide the evolution of the impinging jet flow while a system of bifurcating pre-distribution channels increases pressure drops in the bundle. Lateral communications between pre-distribution channels are introduced to further homogenize the flow. Two experimental facilities have been conceived to study the flow behavior in bifurcating channels and in the inlet header, respectively. At the same time, their effect on the flow distribution between channels is evaluated. The acquired PIV aerodynamic database allows to validate the numerical models and to prove the design basis for the proposed distribution system. Once having validated the CFD turbulence models and the strategy to study the flow maldistribution in the SGHE module, a decisive and trustworthy optimization of each component of the sodium distribution system has been performed. Finally, an optimal configuration has been proposed for the actual phase of ASTRID project.

Highly-loaded turbomachine blades can stall under off-design conditions. In this regime, the flow can separate close to the leading edge of the blade in a periodic manner that can lead to blade vibrations, commonly referred to as stall flutter. Prior experimental studies on stall flutter have been at large Re (Re ~ 106). In the present work, motivated by applications in Unmanned Air Vehicles (UAV) and Micro Air Vehicles (MAV), we study experimentally the forces and flow fields around an oscillating blade at low Re (Re ~ 3 x 104). At these low Re, the flow even over the stationary blade can be quite different. We experimentally study the propensity of an isolated symmetric and cambered blade (with chord c) to undergo self-excited oscillations at high angles of attack and at low Reynolds numbers (Re ~ 30, 000). We force the blade, placed at large mean angle of attack, to undergo small amplitude pitch oscillations and measure the unsteady loads on the blade. From the measured loads, the direction and magnitude of energy transfer to/from the blade is calculated. Systematic measurements have been made for varying mean blade incidence angles and for different excitation amplitudes and frequencies (f). These measurements indicate that post stall there is a possibility of excitation of the blade over a range of Strouhal Numbers (St = fc/U) with the magnitude of the exciting energy varying with amplitude, frequency and mean incidence angles. In particular, the curves for the magnitude of the exciting energy against Strouhal number (St) are found to shift to higher St values as the mean angle of attack is increased. We perform the same set of experiments on two different blade shapes, namely NACA 0012 and a compressor blade profile, SC10. Both blade profiles show qualitatively similar phenomena. The flow around both the stationary and oscillating blades is studied through Particle Image Velocimetry (PIV). PIV measurements on the stationary blade show the gradual shift of the flow separation point towards the leading edge with increasing angle of attack, which occurs at these low Re. From PIV measurements on an oscillating blade near stall, we present the flow field around the blade at different phases of the blade oscillation. These show that the boundary layer separates from the leading edge forming a shear layer, which flaps with respect to the blade. As the Strouhal number is varied, the phase between the flapping shear layer and the blade appears to change. This is likely to be the reason for the observed change in the sign of the energy transfer between the flow and the blade that is responsible for stall flutter.

In comparison to the flow in a rigid channel, there is a multi-fold reduction in the transition Reynolds number for the flow in a micro channel when one of the walls is made sufficiently soft, due to a dynamical instability induced by the fluid-wall coupling. The flow after transition is characterized using Particle Image Velocimetry (PIV) in the x − y plane where x is the stream-wise direction and y is the cross-stream co-ordinate along the small dimension of the channel of height 0.2 − 0.3mm. For the two different soft walls of shear modulus 18 kPa and 2.19 kPaused here, the transition Reynolds number is about 250 and 330 respectively. The deformation of the microchannel due to the applied pressure gradient is measured in the experiments, and is used to predict the laminar mean velocity profiles for comparison with the experimental results. The mean velocity profiles in the microchannel are in quantitative agreement with those predicted for the laminar flow before transition, but are flatter near the centerline and have higher gradients at the wall after transition. The flow after transition is characterized by a mean velocity profile that is flatter at the center and steeper at the walls in comparison to that for a laminar flow. The root mean square of the stream-wise fluctuating velocity shows the characteristic sharp increase from the wall and a maximum close to the wall, as observed in turbulent flows in rigid-walled channels. However, the profile is asymmetric with a significantly higher maximum close to the soft wall in comparison to that close to the hard wall, and the Reynolds stress is found to be non-zero at the soft wall, indicating that there is a stress exerted by fluid velocity fluctuations on the wall. The turbulent energy production profile has a maximum at the soft wall, in contrast to the flow at a rigid surface where the turbulent energy production is zero at the wall (due to the zero Reynolds stress). The maximum of the root mean square of the velocity fluctuations and the Reynolds stress (divided by the fluid density) in the soft-walled microchannel for Reynolds numbers in the range 250-400, when scaled by suitable powers of the maximum velocity, are comparable to those in a rigid channel at Reynolds numbers in the range 5000-20000. The near-wall velocity profile shows no evidence of a viscous sub-layer for (yv∗/ν) as low as 2, but there is a logarithmic layer for (yv∗/ν) up to about 30, where the von Karman constants are very deferent from those for a rigid-walled channel. Here, v∗ is the friction velocity, ν is the kinematic viscosity and y is the distance from the soft surface. . The surface of the soft wall in contact with the fluid is marked with dye spots to monitor the deformation and motion along the fluid-wall interface. The measured displacement of the surface in the stream-wise direction, which is of the order of 5 − 12µm, is consistent with that calculated on the basis of linear elasticity. Low-frequency oscillations in the displacement of the surface are observed after transition in both the stream-wise and span-wise directions, indicating that the turbulent velocity fluctuations are dynamically coupled to motion in the solid. Modification of soft-wall turbulence in a micro channel due to the addition of small amounts of polymer The modification of soft-wall turbulence in a microchannel due to the addition of small amounts of polymer is experimentally studied using Particle Image Velocimetry (PIV) to measure the mean and the fluctuating velocities. The micro channels are of rectangular cross-section with height about 160 µm, width about 1.5 mm and length about 3 cm, with three walls made of hard Poly-dimethylsiloxane (PDMS) gel, and one wall made of soft PDMS gel with an elasticity modulus of about 18 kPa. A dynamical instabilty of the laminar flow due to the fluid-wall coupling, and a transition to turbulence, is observed at a Reynolds number of about 290 for the flow of pure water in the soft-walled microchannel (Verma and Kumaran, J. Fluid Mech., 727, 407-455, 2013). Solutions of polyacrylamide of molecular weight 5 × 106 and mass fraction up to 50 ppm, and of molecular weight 4 × 104 and mass fraction up to 1500 ppm, are used in the experiments. In all cases, the solutions are in the dilute limit be-low the critical concentration where the interactions between polymer molecules become important. The modification of the fluid viscosity due to addition of polymer molecules is small; the viscosity of the solutions with the highest polymer concentration exceed those for pure water by about 10% for the polymer with molecular weight 5 × 106, and by about 5% for the polymer with molecular weight 4 × 104. Two distinct types of flow modifications below and above a threshold mass fraction for the polymer, cTHRESHOLD , which is about 1 ppm for the polyacrylamide with molecular weight 5 × 106, and about 500 ppm for the polyacrylamide with molecular weight 4 × 104. As the polymer mass fraction increases up to the threshold value, there is no change in the transition Reynolds number, but there is significant turbulence attenuation the root mean square velocities in the stream wise and cross-stream directions decrease by a factor of 2, and the Reynolds stress decreases by a factor of 4 in comparison to that for pure water. When the polymer concentration increases beyond the threshold value, there is a decrease in the decrease in the transition Reynolds number by nearly one order of magnitude, and a further decrease in the intensity of the turbulent fluctuations. The lowest transition Reynolds number of about 35 for the solution of polyacrylamide with molecular weight 5 × 106 and mass fraction 50 ppm. For the polymer solutions with the highest concentrations, the fluctuating velocities in the stream wise and cross-stream direction are lower by a factor of 5, and the Reynolds stress is lower by a factor of 10, in comparison to pure water. Despite the significant turbulence attenuation, a sharp increase in the intensity of the fluctuating velocities is evident at transition for all polymer concentrations. Transitions to deferent kinds of turbulence in a channel with soft walls The flow in a rectangular channel with walls made of soft polyacrylamide gel is studied to examine the effect of soft walls on transition and turbulence. The width of the channel is much larger than the height, so that the flow can be considered approximately two-dimensional, the wall thickness is much larger than the channel height (smallest dimension), the bottom wall is fixed to a substrate and the top wall is unrestrained. The fluid velocity is measured using Particle Image Velocimetry, while the wall motion is studied by embedding beads in the soft wall, and measuring the time-variation of the displacement both parallel and perpendicular to the surface. As the Reynolds number increases, two different flow regimes are observed in sequence. The first is the ‘soft-wall turbulence’ resulting from a dynamical instability of the base flow due to the fluid-wall coupling. The flow in this case exhibits many of the features of the turbulent flow in a rigid channel, including the departure of the velocity profile from the parabolic profile, and the near-wall maxima in the stream-wise root mean square fluctuating velocity. However, there are also significant differences. The turbulence intensities, when scaled by suitable powers of the mean velocity, are much larger than those after the hard-wall laminar-turbulent transition at a Reynolds number of about 1000. The Reynolds stress profiles do not decrease to zero at the walls, indicating that the wall motion plays a role in the generation of turbulent fluctuations. There is no evidence of a viscous sub-layer close to the wall to within the experimental resolution. The mean velocity profile does satisfy a logarithmic law close to the surface within a region between 2-30 wall units from the surface, but the von Karman constants are very different from those for the hard-wall turbulence. The wall displacement measurements indicate that there is no observable motion perpendicular to the surface, but displacement fluctuations parallel to the surface are observed after transition, coinciding with the onset of velocity fluctuations in the fluid. The fluid velocity fluctuations are symmetric about the center line of the channel, and they show relatively little downstream variation after a flow development length of about 5 cm. As the Reynolds number is further increased, there is a second ‘wall flutter’ transition, which involves visible downstream traveling waves in the top (unrestrained) wall alone. Wall displacement fluctuations of low frequency (less than about 500 rad/s) are observed both parallel and perpendicular to the wall. The mean velocity profiles and turbulence intensities are asymmetric, with much larger turbulence intensities near the top wall. There is no evident logarithmic profile close to either the top or bottom wall. Fluctuations are initiated at the entrance of the test section, and the fluctuation intensities decrease with downstream distance, the fluctuation intensities first rapidly increase and then decrease as the Reynolds number is increased. For a channel with relatively small height (0.6 mm), the transition Reynolds number for the soft-wall instability is lower the hard-wall transition Reynolds number of about 1000, and the laminar flow becomes unstable to the soft-wall instability leading to soft-wall turbulence and then to wall flutter as the Reynolds number is increased. For a channel with relatively large height (1.8 mm), the transition Reynolds number for the soft-wall instability is higher than 1000, the flow first undergoes the hard-wall laminar-turbulent transition at a Reynolds number of about 1000, the turbulent flow undergoes the soft-wall transition leading to soft-wall turbulence, and then to wall flutter.

S požadavky materiálových věd na stále menší částice jsou potřebné i nové přístupy a metody jejich klasifikace. V disertační práci jsou zkoumány struktury turbulentního proudění a trajektorie částic uvnitř dynamického větrného třídiče. Zvyšující se výpočtový výkon a nové modely turbulence a přístupy modelování komplexních plně turbulentních problémů řešením Navier-Stokesových rovnic umožňují zkoumání stále menších lokálních proudových struktur a vlastností proudění s větší přesností. Částice menší než 10 mikronů jsou více ovlivnitelné a jejich klasifikace do hrubé nebo jemné frakce závisí na malých vírových strukturách. Práce se zaměřuje na podmínky nutné ke klasifikaci částic pod 10 mikronů, což je současná hranice možností metody větrné separace. CFD software a poslední poznatky modelování turbulence jsou použity v numerické simulaci proudových polí dynamického větrného třídiče a jsou zkoumány efekty měnících se operačních parametrů na proudová pole a klasifikaci diskrétní fáze. Experimentální verifikace numerických predikcí je realizovaná prostřednictvím částicové anemometrie na základě statistického zpracování obrazu (PIV) a proudění lopatkami rotoru je vizualizováno. Predikované trajektorie částic jsou experimentálně ověřeny třídícími testy na větrném třídiči a granulometrie je určená pomocí laserové difrakční metody. Zkoumány jsou Trompovy křivky a efektivita třídění.

[pt] A estenose aórtica é um dos mais graves problemas decorrentes de doenças valvares. O implante da prótese valvar aórtica por cateterismo (TAVI) vem se tornando o tratamento mais indicado aos pacientes de alto risco ou inoperáveis. A estenose aórtica grave pode ser uma condição de risco à vida quando não tratada. Devido à natureza do procedimento TAVI, é esperada uma variabilidade no ângulo de inclinação da válvula implantada. O presente trabalho investigou a influência das variações de inclinação e orientação do jato transvalvar sobre o escoamento em aorta ascendente. A compreensão dos padrões hemodinâmicos do fluxo sanguíneo em aorta ascendente é importante porque eles estão intimamente relacionados ao desenvolvimento de doenças cardiovasculares. Para este fim, um modelo vascular com geometria anatômica de paciente específico foi produzido a partir de imagens de tomografia computadorizada, gerando um protótipo impresso em 3D e resina de silicone transparente. Uma configuração especial foi projetada para permitir medições tridimensionais do fluxo em diferentes seções transversais do modelo aórtico. A técnica de velocimetria por imagens estereoscópicas de partículas foi implementada para produzir informações estatísticas acerca do fluxo turbulento, tais como, campos tridimensionais de velocidade média, de energia cinética turbulenta e correlações entre os componentes de flutuação de velocidade. Os resultados obtidos indicaram que o escoamento em aorta ascendente é fortemente afetado pela direção do fluxo de entrada na aorta. / [en] Aortic stenosis is one of the most serious problems arising from valve diseases. Transcatheter Aortic Valve Implantation (TAVI) has become the preferred treatment for high-risk or inoperable patients with severe aortic stenosis that could be a lifethreatening condition when left untreated. Due to the nature of the TAVI procedure, a variability on the tilt angle of the deployed valve is expected. The present work, investigated the effects on the flow field in the ascending aorta due angle variation of the transvalvular jet. Understanding the hemodynamic patterns of blood flow in the ascending aorta is important because they are closely related to the development of cardiovascular diseases. To this end, a patient-specific vascular phantom was produced by a 3D printed model and transparent silicon resin. A special setup was designed to allow measurements of the 3D flow at different cross sections of the aorta. A stereoscopic particle image velocimetry system was implemented to yield instantaneous and averaged turbulent flow information, such as three-dimensional average velocity fields, turbulent kinetic energy, and correlations between the components of velocity fluctuation. The results obtained indicated that the velocity field in the ascending aorta is strongly affected by the inlet flow direction into the aorta.

L’étude présentée dans cette thèse relève de la mécanique des fluides expérimentale et numérique appliquée aux écoulements pariétaux de refroidissement de chambres de combustion aéronautiques. En présence de phénomènes thermo-acoustiques, comme les instabilités de combustion, il est important d’évaluer si les capacités de l’écoulement pariétal à protéger les parois de chambre restent suffisantes. C’est ainsi que nous nous sommes intéressés aux écoulements de paroi multiperforée soumis à une excitation acoustique. Dans ce but, le banc d’essais MAVERIC a été amélioré grâce à l’installation d’un système qui permet de forcer acoustiquement l’écoulement transverse dans lequel les jets pariétaux débouchent. Nous avons pu alors mettre en évidence la forte sensibilité de ce type d’écoulements à l’excitation acoustique. Le bon accord entre les résultats expérimentaux et les simulations numériques aux grandes échelles (LES) effectuées est très encourageant dans le cas d’un forçage par onde stationnaire. Le forçage par onde progressive, étudié uniquement par simulations numériques, s’est révélé être capable de modifier significativement la topologie de l’écoulement. Enfin, à partir de l’outil numérique AVBP-AVTP qui permet le couplage de calculs fluide-solide, nous avons réalisé une étude de l’influence de la présence d’une excitation acoustique sur le comportement thermique de l’écoulement autour d’une paroi multiperforée de chambre de combustion. / This experimental and numerical study in the field of fluid mechanics deals with jets-in cross flow configurations that are relevant for the cooling of aero engine combustion chambers. Indeed, in presence of instabilities it is important to determine to which extent the film cooling is able to do its job of preserving the combustion chamber walls from the thermal load. The test facility MAVERIC has been upgraded in order to acoustically force the crossflow in which the jets are discharging. The strong sensitivity of the overall flow unsteady properties to the presence of the acoustic forcing has been clearly evidenced. The agreement between the experimental results and large-eddy simulations proved to be quite encouraging for a stationary acoustic wave whereas the case of a propagating acoustic wave investigated only numerically reveals also quite a significant change of the flow topology. In this context, the effect of the acoustic forcing on the wall thermal behavior has been analyzed thanks to the use of the fluid-solid coupled AVBP-AVTP solver.

Development of modern safety standards for hydrogen storage infrastructure requires fundamental insight into the physics of buoyant gas dispersion into ambient air. Also, from a practical engineering stand-point, flow patterns and dispersion of gas originating from orifices in the side wall of circular pipe or storage tank need to be studied. In this thesis, novel configurations were considered to investigate the evolution of turbulent jets issuing from realistic pipeline geometries. First, the effect of jet densities and Reynolds numbers on vertical jets were investigated, as they emerged from the side wall of a circular pipe, through a round orifice. The resulting jet flow was thus issued through a curved surface from a source whose original velocity components were nearly perpendicular to the direction of the ensuing jets. Particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF) techniques were employed simultaneously to provide instantaneous and time-averaged flow fields of velocity and concentration. The realistic flow arrangement resulted in an asymmetric flow pattern and a significant deflection from the vertical axis of jets. The deflection was influenced by buoyancy, where heavier gases deflected more than lighter gases. These realistic jets experienced faster velocity decay, and asymmetric jet spreading compared to round jets due to significant turbulent mixing in their near field. In addition to that, horizontal multi-component jets issuing from a round orifice on the side wall of a circular tube were also investigated experimentally by the means of simultaneous velocity and concentration measurements. A range of Reynolds numbers and gas densities were considered to study the effects of buoyancy and asymmetry on the resulting flow structure. The realistic pipeline jets were always exhibited an asymmetry structure and found to deflect about the jet's streamwise axis in the near field. In the far field, the buoyancy dominated much closer to the orifice than expected in the axisymmetric round jet due to the realistic leak geometry along with the pipeline orientation considered in this study. In general, significant differences were found between the centreline trajectory, spreading rate, and velocity decay of conventional horizontal round axisymmetric jets issuing through flat plates and the pipeline leak-representative jets considered in the present study. Finally, the dispersion of turbulent multi-component jets issuing from high-aspect-ratio slots on the side wall of a circular tube were studies experimentally by employing simultaneous PIV and PLIF techniques. Two transversal & longitudinal oblong geometries in respect to the longitudinal axes of the tube , and with an aspect ratio of 10 were considered in this study. Both horizontal and vertical orientations along with broad range of Reynolds numbers and gas densities were considered to investigate the effects of buoyancy and asymmetry on the resulting flow structure. The ensuing jets were found to deflect along the jet streamwise axis, once more, due to the realistic pipeline leak-representative configuration. It was also found that increases in aspect ratio of these realistic jets caused a reduction in the angle of deflection, jet centreline decay rates and the width growth on both velocity and scalar fields compared to their round jets counterparts, most notably in the far field. These findings indicate that conventional jets (those that are issuing through flat surfaces) assumptions are inadequate to predict gas concentration, entrainment rates and, consequently, the extent of the flammability envelope of realistic gas leaks. Thus, extreme caution is required when using conventional jet assumptions to describe the physics of a buoyant jet emitted from realistic geometries. / Graduate

Les principaux objectifs de cette étude est de fournir les informations essentielles sur le comportement thermique des cellules de batterie pour une application automobile, en particulier pour les véhicules électriques et hybrides. Cette application est notre cadre de travail expérimental afin de développer un modèle électro-thermique 3D efficace pour les cellules lithium-ion et du pack batterie. L'étude expérimentale se concentre sur la distribution de température en différents points de la surface de la cellule, de l'impact de différents débits constants, et également l'importance du système de refroidissement sur le comportement en température de la batterie. Cette thèse met en évidence le comportement de température de la cellule dans des conditions de décharge agressive et de l'impact de l'empilement de plusieurs cellules à l'intérieur de la batterie. Une étude de cas sur le comportement thermique de la cellule dans une application véhicule électrique hybride série est proposée pour compléter les cycles de conduite en utilisant différentes stratégies de refroidissement. En outre, l'étude expérimentale est étendue à la caractéristique du comportement de refroidissement par flux d'air à l'intérieur de la batterie, en utilisant le système d'image de particules (PIV). Le modèle électro-thermique CFD 3D est développé sous un logiciel Open Source OpenFOAM. L'objectif principal est d'obtenir un modèle relativement simple mais précis avec un temps de calcul raisonnable. Le modèle proposé, estime la production de chaleur, à partir du courant de la batterie et la résistance interne en fonction de la température, le transfert de chaleur par conduction, convection forcée et rayonnement. / The main objectives of this study are to provide the essential information on the thermal behaviour of the battery cells for automotive purpose especially for EVs and HEVs through experimental work in order to develop an effective 3D electro-thermal model for lithium ion battery cells and pack. The experimental study is focusing on the distribution of temperature at various points of the battery cell surface, impact of different constant discharge rates, and also the importance of cooling system on the battery temperature behaviour. This thesis highlights the battery cell temperature under abuse discharge condition and the impact of stacking the battery cells inside the battery pack. Impact of different temperature and SOC on the battery cell internal resistance and a case study on the battery cell thermal behaviour used in a series HEV to complete driving cycles using different cooling strategies are also studied. Furthermore, the experimental study is extended to the characteristic of the cooling air flow behaviour inside the battery pack, using particulate image velocimetry (PIV) system. The 3D electro-thermal CFD model is implemented in a free, open source CFD software package called OpenFOAM. The target is to have a relatively simple but accurate model with reasonable computation time. This proposed model considers the heat generation from battery current and internal resistance as a function of temperature, heat transfer through conduction, forced convection and radiation.