Cette thèse présente une extension de l'approche stochastique de l'atomisation primaire de type air assisté près d'un injecteur. Cette approche avait déjà été introduite dans les publications de Gorokhovski et al. Dans le cadre de la simulation des grandes échelles, la zone d'atomisation primaire est simulée comme un corps immergé avec une structure stochastique. Ce dernier est défini par la simulation stochastique de la position et de la courbure de l'interface entre le liquide et le gaz. La simulation de la position de l'interface est basée sur l'hypothèse de symétrie d'échelle pour la fragmentation. La normale extérieure à l'interface est modélisée en supposant une relaxation statistique vers l'isotropie. Les statistiques de la force du corps immergé servent de conditions aux limites pour le champ de vitesse issu de la LES ainsi que pour la production des gouttes de l'atomisation primaire. Celles-ci sont ensuite transportées par une approche lagrangienne. Les collisions entre les gouttes dans la zone d'atomisation primaire sont prises en compte par analogie avec l'approche standard de la théorie cinétique des gaz. Une fermeture est proposée pour la température statistique des gouttelettes. Cette approche est validée par des comparaisons avec les mesures expérimentales de la thèse de Hong. Les résultats numériques pour la vitesse et de la taille des gouttes dans le spray à différentes distances du centre du jet et de l'orifice de la buse sont relativement proches des résultats expérimentaux. Différentes conditions d'entrée pour la vitesse sont testées et comparées aux résultats expérimentaux. Par ailleurs, le rôle spécifique de la zone de recirculation devant le dard liquide est soulignée par le battement du dard liquide et la production de gouttelettes. / This thesis introduced an extension to stochastic approach for simulation of air-blast atomization closely to injector. This approach was previously proposed in publications of Gorokhovski with his PHD students. Our extension of this approach is as follows. In the framework of LES approach, the contribution of primary atomization zone is simulated as an immersed solid body with stochastic structure. The last one is defined by stochastic simulation of position-and-curvature of interface between the liquid and the gas. As it was done previously in this approach, the simulation of the interface position was based on statistical universalities of fragmentation under scaling symmetry. Additionally to this, we simulate the outwards normal to the interface, assuming its stochastic relaxation to isotropy along with propagation of spray in the down-stream direction. In this approach, the statistics of immersed body force plays role of boundary condition for LES velocity field, as well as for production of primary blobs, which are then tracked in the Lagrangian way. In this thesis, the inter-particle collisions in the primary atomisation zone are accounted also by analogy with standard kinetic approach for the ideal gas. The closure is proposed for the statistical temperature of droplets. The approach was assessed by comparison with measurements of Hong in his PHD. The results of computation showed that predicted statistics of the velocity and of the size in the spray at different distances from the center plane, at different distances from the nozzle orifice, at different inlet conditions (different gas velocity at constant gas-to-liquid momentum ratio, different gas-to-liquid momentum ratio) are relatively close to measurements. Besides, the specific role of recirculation zone in front of the liquid core was emphasized in the flapping of the liquid core and in the droplets production.

Cette thèse porte sur la simulation des jets de gouttes générés par des pulvérisateurs essence haute pression, pulvérisateurs qui sont un point clef des systèmes de combustion automobile de la présente et future génération devant diminuer les émissions de CO2 et de polluants. Dans un premier temps les jets de gouttes (« sprays ») sont simulés par simulation moyennée. Les résultats de simulation d’un jet donnant des résultats en moyenne satisfaisant, l'interaction de jets en injecteurs multi-trous est alors simulée. Les résultats sont cohérents par rapport aux mesures d'entraînement d’air. La simulation permettant d'avoir accès au champ complet 3D, le mécanisme d'interaction jet à jet et de développement instationnaire du spray est décrit en détail. La formation d’un mouvement descendant au centre du spray et celle d'un point d'arrêt central sont trouvés. Finalement, Ces résultats sont étendus au cas surchauffé, cas où la pression dans la chambre est inférieure à la pression de vapeur saturante. Un modèle simple semi-empirique est proposé pour tenir compte de la modification des conditions proches de la buse d’injection. Le modèle prédit correctement les tendances des variations de paramètres et capture la forme générale du spray qui se referme sur lui-même. La seconde grande partie est consacrée au développement d’un modèle de spray par l’approche des grandes échelles (SGE), limité ici aux cas non évaporant. Il comprend la modélisation de sous-maille de la dispersion turbulente, des collisions-coalescence et des termes d’échange de quantité de mouvement de sous-maille. L'effet du choix du modèle de sous-maille pour la viscosité turbulente de sous-maille est montré, le choix retenu étant le modèle de Smagorinski dynamique. Afin d'améliorer la représentativité cruciale des conditions d’injections, un couplage faible est réalisé à partir de résultats de simulations existantes de l'écoulement interne aux buses. Les fonctions densité de probabilité simple et jointes extraits des résultats de simulations sont validés par rapport aux mesures PDA en situation pseudo-stationnaire et la pénétration liquide et la forme du spray est comparée aux visualisations par ombroscopie. Enfin, différentes zones caractéristiques sont identifiées et des longueurs sont notées pour les cas d'injection à 100 et 200bar. / Over the years numerical modelling and simulation techniques have constantly been improved with the increase in their use. While keeping the computational resources in mind, numerical simulations are usually adapted to the required degree of accuracy and quality of results. The conventional Reynolds Average Navier Stokes (RANS) is a robust, cheap but less accurate approach. Large Eddy Simulation (LES) provides very detailed and accurate results to the some of the most complex turbulence cases but at higher computational cost. On the other hand, Direct Numerical Simulation (DNS) is although the most accurate of the three approaches but at the same time it is computationally very expensive which makes it very difficult to be applied to the most of the complex industrial problems. The current work is aimed to develop a deeper understanding of multi-hole Gasoline Direct Injection (GDI) sprays which pose many complexities such as; air entrainment in the multi-hole spray cone, Jet-to-Jet interactions, and changes in the spray dynamics due to the internal flow of the injector. RANS approach is used to study multi-hole injector under cold, hot and superheated conditions. Whereas, LES is utilized to investigate the changes in the dynamics of the single spray plume due to the internal flow of the GDI injector. To reduce computational cost of the simulations, dynamic mesh refinement has been incorporated for both LES and RANS simulation. A thorough investigation of air entrainment in three and six hole GDI injectors has been carried out using RANS approach under non superheated and superheated conditions. The inter plume interactions caused by the air entrainment effects have been analysed and compared to the experimental results. Moreover, the tendencies of semi collapse and full collapse of multi-hole sprays under non superheated and superheated conditions have been investigated in detail as well. A methodology of LES has been established using different injection strategies along with various subgrid scale models for a single spray plume. In GDI multi-hole sprays, the internal flow of the injector plays a very crucial role in the outcome the spray plume. A separate already available internal flow LES simulation of the injector has been coupled with the external spray simulation in order to include the effect of nozzle geometry and the cavitation phenomenon which completely change the dynamics of the spray.

La demande croissante pour alléger les moteurs d’avions et diminuer les émissions polluantes de la propulsion aéronautique réclame à rendre plus compact le système de compression des moteurs, qui représente environ 40%-50% de la masse totale. Or, à taux de compression global égal, la réduction du nombre d’étage exige d’un point de vue aérodynamique une augmentation de la charge des aubes de compresseur par étage. La charge d’aube est aujourd’hui limitée car elle induit différents mécanismes de pertes tridimensionnelles très pénalisant. L’un des plus importants est le décollement de coin qui se forme à la jonction entre l’extrados de l’aube et le moyeu ou le carter. Bien que des travaux existent sur les mécanismes et paramètres intervenant dans le décollement de coin, il est encore difficile de proposer une méthode de contrôle efficace. Cela est principalement dû à deux raisons : (i) le manque de compréhension fine des mécanismes physiques, (ii) l’utilisation pour la conception de modèles de turbulence classiques de type RANS (Reynolds-averaged Navier-Stokes) qui ne sont pas capables de prédire précisément le décollement de coin, car ils ne peuvent pas décrire correctement les mécanismes de transport turbulent. Des simulations de type RANS et LES (large-eddy simulation = simulation des grandes échelles) sont présentées dans cette thèse sur une configuration de grille d’aubes de compresseur, et comparées avec les données expérimentales obtenues au LMFA (issues de travaux séparés). L’approche RANS surestime globalement le décollement de coin. Une amélioration significative est obtenue par la méthode LES, en particulier pour le coefficient de pression statique sur l’aube et les pertes de pression totale. Ces résultats montrent que la zone de décollement de coin, qui est la source principale des pertes, génère des tourbillons de grande échelle associés à de forts niveaux d’énergie. Les histogrammes bimodaux de la vitesse tangentielle qui ont été observés expérimentalement semblent confirmés par les résultats LES. En ce qui concerne les amplitudes des fluctuations de vitesse tangentielle, les résultats expérimentaux et ceux de la LES mettent en évidence deux pics sur certains profils perpendiculaires aux parois. Enfin, grâce à l’approche LES, les bilans de l’énergie cinétique turbulente sont calculés et analysés. Ils décrivent l’équilibre entre les termes de production, de dissipation et de transport. Une des perspectives de cette analyse est d’aider à améliorer la modélisation de la turbulence en approche RANS. / The increasing demand to reduce the mass of aircraft jet engines and emissions of aircraft propulsion requires to make the compression system of engines more compact, since this component accounts for about 40%-50% of the total mass. However, at a given overall pressure ratio, decreasing the number of stages will raise the compressor blade loading per stage. The blade loading is extremely restricted by different three-dimensional flow loss mechanisms. One of them is the corner separation that forms between the blade suction side and the hub or shroud. Although some works previously investigated the mechanisms and the parameters of corner separation, it is still difficult to propose an effective control method of the corner separation. That is mainly due to two reasons: (i) the lack of knowledge of the physical mechanisms, (ii) the nowadays classical RANS (Reynolds-averaged Navier-Stokes) turbulence models are not capable to accurately predict the corner separation, since they cannot correctly describe the turbulent transport mechanisms. RANS (Reynolds-averaged Navier-Stokes) and LES (large-eddy simulation) simulations are here presented on a compressor cascade configuration, in comparison with experimental data obtained at LMFA (from separate works). The RANS approach globally over-estimates the corner separation, whereas a significant improvement is achieved with the LES, especially for the blade surface static pressure coefficient and the total pressure losses. The corner separation region, which is the main source of the total pressure losses, is shown to generate large-scale energy-containing eddies. The bimodal histograms of the streamwise velocity that were observed experimentally seem to be confirmed by the LES results. Concerning the streamwise velocity fluctuations (RMS), both the experiment and the LES show some profiles with two peaks. Finally, thanks to the LES approach, the turbulent kinetic energy budget, which represents the balance between the production, dissipation and transport terms, are computed and analyzed. This may help the improvement of RANS turbulence modeling.

Cette thèse se consacre à l'étude du modèle de dénaturation de l'ADN introduit par Poland et Scheraga dans les années soixante. Les modèles de dépiégeage en milieu aléatoire, avec lesquels la correspondance a été établie, sont également traités. Dans le cas où les interactions entre le système et l’environnement sont homogènes, le problème a été résolu : selon la valeur d'un paramètre géométrique, une transition de phase d'ordre un ou deux se produit. En revanche, lorsque les interactions sont prises aléatoires (on parle d'un système en présence de désordre), nous ne connaissons ni le point critique, ni l'ordre de la transition en régime defort désordre. Pour simplifier le problème, de nombreux auteurs font usage d'une représentation hiérarchique grâce à laquelle une renormalisation exacte de la fonction de partition peut être écrite. Mais à nouveau, la question du point critique et de l'ordre de la transition n'a pas été résolue. Nous avons introduit un nouveau système (Toymodel) plus simple que la version hiérarchique en changeant la forme de la renormalisation. Le problème, ainsi posé, permet de mettre en évidence une famille de distributions qui ne varient presque pas lors d'une renormalisation, avec lesquelles nous avons pu dériver des équations du type Berezinskii-Kosterlitz- Thouless. Aussi, en présence de désordre, la transition de phase n'admet pas de point fixe critique. Ces deux éléments, en accord avec nos résultats numériques, nous poussent à croire que nous sommes en présence d'une transition de phase d'ordre infini. La seconde partie de la thèse rapporte un travail sur le processus simple d'exclusion symétrique, qui est l'un des modèles les plus simples de physique statistique hors d'équilibre pour lequel un état stationnaire est connu. La fonction de grandes déviations a été calculée dans le passé par les approches microscopiques et macroscopiques et ici, nous en avons calculé la première correction de taille finie. Le résultat a ensuite été comparé aux corrections similaires pour des systèmes à l'équilibre. / This thesis is a study of a DNA denaturation model, introduced by Poland and Scheraga during the 1960s. The depinning models with random environment, with which the similarity has been made, are also concerned. If the interactions between the system and the environment are homogeneous, the problem has been solved: depending on the value of a geometrical parameter, a first or a second order phase transition happens. On the other hand, when the interactions are random, we know neither the critical point nor the phase transition order in the case of strong disorder. In order to simplify the problem, some authors have used a hierarchical representation through which an exact renormalization can be written. Despite this simplification, the critical point and the transition order have not been found. By changing the renormalization relation, we introduced a Toy-model which is simpler than the hierarchical version. The new problem leaded us to a family of distributions, which stay almost the same under renormalization, and allow us to derive the Berezinskii-Kosterlitz- Thouless equations. Also, with strong disorder, the phase transition does not have a critical fixed point. These two elements, according to our numerical results, predict that the order transition is infinite. The second part of this thesis reports on a work about the simple symmetric exclusion process, which is one of the simplest out of equilibrium models for which a stationary state is known. The large deviation function has been calculated in the past through microscopic and macroscopic approaches. Here, we calculated the leading finite-size correction. Then the result has been compared to similar corrections for equilibrium systems.

Ce travail concerne le développement d’une fibre optique à structure originale permettant la montée en puissance moyenne et crête dans les systèmes lasers à 2μm. La gestion des effets thermiques est devenue aujourd’hui un enjeu primordial notamment pour essayer de repousser le seuil d’apparition des instabilités modales transverses. Le principe mis en avant ici est basé sur la rupture de symétrie de la gaine microstructurée de la fibre afin d’améliorer la délocalisation des modes d’ordre supérieur vers l’extérieur du milieu à gain. Ainsi une propagation quasi-monomode est obtenue dans une fibre apériodique passive avec un coeur de 140μm à 2μm. Un travail de modélisation a été mené de manière à proposer d’autres structures basées sur une modulation contrôlée de l’indice de réfraction dans le milieu à gain afin de repousser encore plus le seuil du régime multimode. En outre, la fabrication de la première fibre dopée thulium à large aire modale (Dcoeur = 18μm) par la méthode REPUSIL a montré une efficacité de 50%. Ensuite, la fabrication de la première fibre rigide complètement apériodique à gaine réduite dopée thulium a été réalisée. Une fibre avec un coeur de 29μm et un diamètre extérieur de 769μm et une longueur de 86cm a été caractérisée en configuration laser et a permis d’obtenir une émission laser à 2μm de 3,8W (puissance limitée par la puissance de la diode de pompe disponible) avec un rendement de 20% et une qualité de faisceau quasi-monomode. De plus, une fibre passive complètement apériodique à polarisation unique avec un coeur de 140μm a également été réalisée et a permis d’obtenir un ratio d’extinction de polarisation de 16,5dB à 2μm. / This work deals with the development of an original leaky structure of optical fiber aiming at generating a high power laser radiation at 2μm in CW and pulsed regimes. The management of thermal effects in high power/energy regime became a major issue, in particular to push further the transverse modal instabilities threshold. Our approach consists in breaking the symmetry of the microstructured fiber cladding for to improve the delocalization of the high order modes outside of the gain medium. Thus, an effective single-mode propagation at a 2μm operating wavelength was first demonstrated into a passive aperiodic fibers whose the core diameter reaches up to 140 μm. After implementing some Stress Applying Parts over our aperiodic design, a PER of 16.5dB was achieved at 2μm using a single polarization passive FA-LPF with a core of 140 μm. Thence, a numerical study on the tailoring the active core refractive index has been carried out so as to fend off the threshold of multimodedness. Additionally, the first LMA Tm-doped fiber (Dcoeur = 18μm) fabricated by the REPUSIL method and showing an efficiency of 50% is demonstrated. Thereafter, the fabrication of the first rod-type Tm-doped FA-LPF with reduced cladding is shown. A 29 μm core FA-LPF was characterized in laser configuration, leading to an effective single-mode emission of 3.8W of average power at 2μm strictly restricted by the available pump power and an efficiency of 20%.

Le transport anormal de l’energie observé en régime turbulent joue un rôle majeur dans les propriétés de stabilite des plasmas de fusion par confinement magnétique, dans des machines comme ITER. En effet, la turbulence plasma est intimement corrélée au temps de confinement de l’energie, un point clé des recherches en fusion thermonucléaire.<p>Du point de vue théorique, la turbulence plasma est décrite par les équations gyrocinétiques, un ensemble d équations aux dérivées partielles non linéaires couplées. Par suite des très différentes échelles spatiales mises en jeu dans des conditions expérimentales réelles, une simulation numérique directe et complète (DNS) de la turbulence gyrocinétique est totalement hors de portée des plus puissants calculateurs actuels, de sorte que démontrer la faisabilité d’une alternative permettant de réduire l’effort numérique est primordiale. En particulier, les simulations de grandes échelles (”Large-Eddy Simulations” - LES) constituent un candidat pertinent pour permettre une telle r éduction. Les techniques LES ont initialement été développées pour les simulations de fluides turbulents à haut nombre de Reynolds. Dans ces simulations, les plus grandes échelles sont explicitement simulées numériquement, alors que l’influence des plus petites est prise en compte via un modèle implémenté dans le code.<p>Cette thèse présente les premiers développements de techniques LES dans le cadre des équations gyrocinétiques (GyroLES). La modélisation des plus petites échelles est basée sur des bilans d’énergie libre. En effet, l’energie libre joue un rôle important dans la théorie gyrocinétique car elle en est un invariant non lin éaire bien connu. Il est démontré que sa dynamique partage de nombreuses propriétés avec le transfert d’energie dans la turbulence fluide. En particulier, il est montré l’existence d’une cascade d énergie libre, fortement locale et dirigée des grandes échelles vers les petites, dans le plan perpendiculaire â celui du champ magnétique ambiant.<p>La technique GyroLES est aujourd’hui implantée dans le code GENE et a été testée avec succès pour les instabilités de gradient de température ionique (ITG), connues pour jouer un rôle crucial dans la micro-turbulence gyrocinétique. A l’aide des GyroLES, le spectre du flux de chaleur obtenu dans des simulations à très hautes résolutions est correctement reproduit, et ce avec un gain d’un facteur 20 en termes de coût numérique. Pour ces raisons, les simulations gyrocinétiques GyroLES sont potentiellement un excellent candidat pour réduire l’effort numérique des codes gyrocinétiques actuels. <p>/ Anomalous transport due to plasma micro-turbulence is known to play an important role in confinement properties of magnetically confined fusion plasma devices such as ITER. Indeed, plasma turbulence is strongly connected to the energy confinement time, a key issue in thermonuclear fusion research. Plasma turbulence is described by the gyrokinetic equations, a set of nonlinear partial differential equations. Due to the various scales characterizing the turbulent fluctuations in realistic experimental conditions, Direct Numerical Simulations (DNS) of gyrokinetic turbulence remain close to the computational limit of current supercomputers, so that any alternative is welcome to decrease the numerical effort. In particular, Large-Eddy Simulations (LES) are a good candidate for such a decrease. LES techniques have been devised for simulating turbulent fluids at high Reynolds number. In these simulations, the large scales are computed explicitly while the influence of the smallest scales is modeled.<p>In this thesis, we present for the first time the development of the LES for gyrokinetics (GyroLES). The modeling of the smallest scales is based on free energy diagnostics. Indeed, free energy plays an important role in gyrokinetic theory, since it is known to be a nonlinear invariant. It is shown that its dynamics share many properties with the energy transfer in fluid turbulence. In particular, one finds a (strongly) local, forward (from large to small scales) cascade of free energy in the plane perpendicular to the background magnetic field.<p>The GyroLES technique is implemented in the gyrokinetic code Gene and successfully tested for the ion temperature gradient instability (ITG), since ITG is suspected to play a crucial role in gyrokinetic micro-turbulence. Employing GyroLES, the heat flux spectra obtained from highly resolved direct numerical simulations are recovered. It is shown that the gain of GyroLES runs is 20 in terms of computational time. For this reason, Gyrokinetic Large Eddy Simulations can be considered a serious candidate to reduce the numerical cost of gyrokinetic simulations. / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

Le présent travail de thèse s’inscrit dans le cadre d’un projet industriel proposé par une entreprise luxembourgeoise et en collaboration avec le Centre de Recherche Public Henry Tudor du Luxembourg (Laboratoire des Technologies Industriels (LTI)). L’objectif consiste à mettre au point une méthode numérique de détermination de la forme des outils d’emboutissage et du flan de pièces de précision en tôlerie fine pour que ce dernier, une fois déformé en une ou plusieurs opérations, correspond à la définition tridimensionnelle du cahier des charges. La méthode a pour objectif de remplacer une démarche expérimentale coûteuse par essais–erreur. Deux démarches numériques sont proposées, la première est relative à la détermination de la forme du flan. Elle consiste à estimer la forme du flan par Approche Inverse en partant de la forme 3D demandée. Puis, un logiciel de simulation incrémental par éléments finis en 3D est utilisé dans une procédure d’optimisation heuristique pour déterminer la forme du flan. Dans la deuxième démarche, il s’agit de déterminer la forme des outils d’emboutissage en utilisant le logiciel de simulation incrémental couplé avec une méthode de compensation du retour élastique en 2D. La démarche numérique est validée expérimentalement dans le cas d’un emboutissage réalisé en une ou plusieurs passes, à l’aide d’une presse manuelle, sans serre flan et avec des outils de forme très complexe. / The present study deals with an industrial project proposed by a luxembourgian enterprise and in collaboration with the luxembourgian research centre Henry Tudor (Laboratory of Industrial Technologies (LTI)). The main objective is to build a numerical approaches for the determination of the initial blank shape contour and tools shape for 3D thin metallic precision parts obtained by stamping, knowing the 3D CAD geometry of the final part. The purpose of the present procedure is to replace the expensive and time consuming experimental trial and error optimization method. Two numerical approaches have been proposed, the first is regarding the determination of the blank shape. An estimation of the blank shape can be given using the Inverse Approach. Update of the blank shape will then be continued by iterations combining heuristic optimization algorithms and incremental stamping codes. The second approach is based on precise finite element models and on spring-back compensation algorithm. The numerical approaches are tested in the case of a special stamping process where the parts are pressed in one or more steps using a manual press, without blank holder and by the mean of complex shape tools.

La simulation de la combustion turbulente connait un nouvel essor avec l'introduction de la Simulation aux Grandes Échelles (SGE) qui permet de prédire l'évolution in stationnaire de l'écoulement réactif turbulent. Dans ce contexte la prise en compte du rayonnement soulève des questions d'ordre a la fois fondamental et pratique. En effet les processus physiques du rayonnement et de la combustion sont de nature radicalement différente : la combustion est contrôlée par des échanges locaux sur une durée finie, alors que le rayonnement est instantané et fait intervenir des échanges a distance. En premier lieu il convient de s'interroger sur l'impact de la modélisation SGE de la combustion turbulente sur le rayonnement. Cette question est traitée dans le cadre plus général de l'interaction rayonnement-turbulence. A partir d'études théoriques et numériques, il est montre que cette interaction est faible et qu'une solution SGE peut être directement utilisée pour un calcul radiatif, sans modélisation supplémentaire. Il s'agit ensuite de mettre en place de façon pratique le couplage in stationnaire rayonnement-combustion turbulente. Un point clé est la réduction du temps de calcul pour le rayonnement, et diverses stratégies sont proposées. En particulier un nouveau modèle spectral est introduit, utilisant une technique de tabulation et garantissant un niveau de précision suffisant. Le temps de calcul radiatif a ainsi été réduit de deux ordres de grandeur, permettant la réalisation d'un calcul couple sur une configuration de flamme pré-melangée turbulente. / Simulation of turbulent combustion has gained high potential with the Large Eddy Simulation (LES) approach, allowing to predict unsteady turbulent reactive flows. In this context, taking into account radiation rises new fundamental and practical questions. Indeed the physics involved in radiation and in combustion are completely different : combustion is controlled by local exchanges and finite times whereas radiation is instantaneous and is based on non-local exchanges. In a first step, the impact of LES modelling of turbulent combustion on radiation is regarded. This question is treated in the more general frame of the turbulence-radiation interaction. From theoretical and numerical studies, it is shown that this interaction is weak in the LES context so that LES solutions can be directly coupled to radiative calculations, without further modelling. Then the unsteady coupling of radiation and turbulent combustion is realised. A key point is the reduction of calculation time of radiation, and several strategies are proposed. In particular a new global spectral model is introduced, using a tabulation technique and ensuring a sufficient level of accuracy. The radiative time calculation is finally decreased by two orders of magnitude, enabling the realization of a coupled calculation of a turbulent premixed flame

Ce travail a pour cadre une collaboration entre le LMA et Thales Alenia Space. Nous nous intéressons au comportement des structures flexibles et plus particulièrement des mètres rubans qui ont la particularité de pouvoir, grâce à l’aplatissement de la section, s’enrouler ou développer des pliages localisés. Une première thèse a permis d’une part la mise au point d’un nouveau type de mètre ruban au déroulement maîtrisable thermiquement et d’autre part le développement d’un modèle plan de poutre à section flexible. Dans le travail de thèse présenté ici, nous proposons une version étendue de ce modèle adaptée à la simulation du comportement dynamique tridimensionnel des mètres rubans en grands déplacements et en grandes rotations. Ce modèle est dérivé de la théorie des coques et repose sur l’introduction d’hypothèses cinématiques et sthéniques adaptées. La déformation de la section est caractérisée par celle de sa ligne moyenne qui peut se déformer dans son plan par flexion et torsion mais non par extension, ainsi que hors de son plan par gauchissement de torsion. Les fortes variations de forme de la section dans son plan peuvent alors être décrites par une cinématique de type Elastica, tandis qu’une cinématique de type Vlassov est utilisée pour définir le gauchissement dans le repère local attaché à la section. Le modèle unidimensionnel est obtenu par intégration sur la section des expressions de la théorie des coques, une approche énergétique permet ensuite de formuler le problème associé qui est résolu grâce au logiciel de modélisation par éléments finis COMSOL. / This work was carried out within the framework of a collaboration between the LMA and Thales Alenia Space. We focus on the behaviour of flexible structures and more specifically of tape springs, whose particularity lies in their capacity to coil up or to form localized folds through the flattening of their cross-section. A first thesis led to the development of a new type of tape spring whose uncoiling is controlled thermically on one hand and of a planar rod model with a flexible thin-walled cross-section on the other hand. In this thesis, we offer an extended version of this model dedicated to the simulation of three-dimensional dynamic behavior of tape springs in large displacements and large rotations. This model is derived from shell theory and is based on the introduction of adapted kinematic and sthenic hypotheses. The deformation of the cross-section is characterized by that of its average line which can deform in its own plane by flexion and twisting but not by extension, as well as out of its plane through torsional warping. The large changes of the cross-section shape in its plane can then be described by an Elastica kinematics, whereas a Vlassov kinematics is used to define the warping in the local frame attached to the section. The unidimensionnal model is obtained by integration over the cross-section of the expressions of the shell theory, an energetic approach then allows to express the associated problem which is solved thanks to the finite element modeling software COMSOL.

Ce mémoire de thèse traite de l’analyse, l’intensification et l’optimisation du transfert thermique convectif dans les échangeurs à plaques et ailettes à pas décalés utilisés notamment dans le domaine automobile comme refroidisseurs d’air de suralimentation. Deux approches complémentaires sont abordées dans cette étude : des simulations numériques visant à l’analyse fine locale des caractéristiques de l’écoulement et des mécanismes de transfert, et une modélisation de type nodale permettant une caractérisation globale des performances thermo-aérauliques. Sur la gamme de Reynolds considérée différentes modélisations de la turbulence sont mises en œuvre et comparées. Ainsi des simulations aux grandes échelles (LES) permettent de qualifier des simulations de type RANS classiquement utilisées jusque-là : de fortes différences tant au niveau structuration de l’écoulement qu’au niveau performances globales sont ainsi mises en évidence selon le régime d’écoulement considéré. La mise au point d’un modèle nodal est ensuite abordée dans le but de mener des optimisations de géométries d’échangeurs non-conventionnels à pas décalés. Les différents scénarii d’optimisation considérés montrent l’intérêt de cette approche autorisant l’évaluation d’un nombre élevé de configurations géométriques. Dans une dernière partie une nouvelle géométrie innovante permettant de générer des tourbillons longitudinaux sur ce type d’ailettes est proposée et étudiée. / This thesis deals with the analysis, intensification and optimization of convective heat transfer in offset strip fins (OSF) heat exchangers used, for example, in the automotive field as water-cooled charge air coolers. Two complementary approaches are carried out in this study: CFD simulations to perform local fine analysis of the flow characteristics and transfer mechanisms, and a nodal type modeling allowing calculation of global aerothermal performance. Over the range of Reynolds numbers considered, different turbulence modeling approaches are implemented and compared: Large Eddy Simulations (LES) and RANS simulations which are usually used. The qualification of the RANS models shows that strong differences, both in the flow structure and at the overall performance evaluation level, may beobserved, depending on the flow regime considered. Then the development of a nodal model is presented. It aims at carrying out rapid optimization of geometries of unconventional OSF heat exchangers. The various optimization scenarios considered show the interest of this approach allowing the evaluation of a large number of geometric configurations. In a last part, an innovative new geometry that generates longitudinal vortices on this type of fins is proposed and studied.