111 |
Une méthode d'éléments finis adaptive pour les problèmes à surfaces libres instationnairesBenmoussa, Khalid 11 April 2018 (has links)
Le calcul des surfaces libres a fait l'objet de plusieurs ouvrages. L'objectif de ces travaux est de calculer avec précision les interfaces entre fluides immiscibles en minimisant la diffusion numérique et les pertes de masse. Plusieurs techniques existent pour effectuer de telles simulations et dans ce travail, on présente une variante améliorée de la méthode des surfaces de niveau (« level-set method »). Cette méthode requiert la résolution d'une équation de transport pour l'interface, couplée à la résolution des équations de NavierStokes. Dans la méthode proposée, l'interface est représentée par la surface de niveau 0 de la fonction distance signée. Cette interface présente une zone de transition de largeur 2e des paramètres rhéologiques des deux fluides comme la viscosité et la densité. Ces paramètres sont régularisés à l'interface pour éviter de gérer des discontinuités. La tension superficielle caractérise la résistance à la déformation de l'interface et joue un rôle très important. Il est donc primordial de bien l'imposer comme terme source dans l'équation de Navier-Stokes, ce qui est réalisé par l'introduction d'un tenseur qui agit uniquement à l'interface. Pour remédier au problème de conservation de la masse, nous avons développé une technique qui consiste à modifier la surface de niveau 0 de manière à conserver le volume. Cette correction est calculée avec la méthode de la sécante. Toutes ces stratégies sont facilitées par l'introduction d'une méthode d'adaptation de maillage instationnaire permettant de concentrer les éléments au voisinage de l'interface, là où le besoin s'en fait sentir. Cela assure une imposition très précise des forces surfaciques ainsi que de la transition des caractéristiques rhéologiques des différents fluides. Ce remaillage adaptatif est basé sur un estimateur d'erreur hiérarchique en dimension 2 et sur un estimateur basé sur une métrique en dimension 3. La résolution des équations de Navier-Stokes est faite par une méthode directe en dimension 2 et par une méthode itérative en dimension 3. L'élément de Taylor-Hood et un schéma implicite de différences arrières d'ordre 2 sont utilisés pour les discrétisations en espace et en temps. On utilise une méthode SUPG pour les équations de transport et de régularisation de l'interface. Enfin, plusieurs problèmes sont traités afin de valider la méthode et de comparer nos résultats avec ceux obtenus par d'autres chercheurs et différentes méthodes. On soulève enfin un certain nombre de questions sur les méthodes employées dans la littérature.
|
112 |
Etude expérimentale de l'interaction d'une onde de choc avec une structure mobile autour d'un axeBiamino, Laurent 30 November 2011 (has links)
Ce travail de thèse s’appuie sur une étude expérimentale en tube à choc, plus précisément, c’est une approche expérimentale de l'étude de l'interaction fluide-structure. Considérons un solide indéformable auquel on laisse un degré de liberté en rotation autour d'un axe. Cette structure ferme un espace clos. Si le contenu de l'espace clos subit le passage d'une onde de choc, ce solide va être mis mouvement et tourner autour de son axe. Concrètement, l'onde de choc va augmenter les caractéristiques physiques, en particulier sa pression, du fluide en contact avec la face impactée de cette porte. La face opposée de la porte ne subissant pas ou que très peu l'influence de l'onde de choc, une seule de ses faces est soumise à la surpression. Au moment de l'impact, le déséquilibre ainsi créé impose une action mécanique sur la porte qui va la faire accélérer et tourner autour de son axe de rotation. Jusqu'à ce stade tout est relativement simple. La difficulté intervient à l'instant où la porte commence à s'ouvrir, car les frontières du volume dans lequel le fluide évolue sont modifiées. Des fuites apparaissent et le gaz qui était maintenu dans un volume clos peut maintenant s'écouler vers un milieu libre. Une communication entre les gaz agissant de chaque coté de la porte est créée modifiant leurs propriétés et par conséquent la pression agissant sur chaque côté de la porte. Les actions mécaniques qui s'appliquent sur la porte ne sont plus les mêmes, et par conséquent l'accélération que la porte subit aussi. Au fur et à mesure que la porte change de position, le problème fluide continue d'être modifié et change en retour son action sur la porte. Cette interaction perdure soit jusqu'à ce que les limites du problème cessent d'être modifiées, la porte ne peut plus bouger, ou bien lorsque les actions mécaniques agissant sur la porte s'équilibrent, les fluides de chaque côté de la porte étant dans le même état physique. Le travail présenté ici est une étude des paramètres du fluide ou du solide en mouvement qui sont les acteurs de la loi comportementale gérant ce système complexe. Pour ce faire, nous avons réalisé une maquette expérimentale mettant en action la physique que nous venons de décrire et nous l'avons adaptée à un tube à choc. En éprouvant de nombreuses configurations expérimentales, nous avons pu déterminer comment l'écoulement interne d'un tube à choc évolue lorsqu'il est plus ou moins ouvert à son extrémité. Comment une porte fermée réagit-elle à l'impact d'une onde de choc et quelles en sont les conséquences sur l'évolution des fluides mis en jeu? Quelles sont les conséquences d'une position différente de la porte au moment de l'impact avec l'onde de choc? Ou encore, quel rôle joue l'intensité de l'onde de choc incidente ou l'inertie de la porte sur toute cette dynamique? / This thesis is based on an experimental study carried out in shock tube; in particular, this is an experimental approach to the study of fluid-structure interaction. Consider a rigid body which is allowed to rotate only around an axis and which closes a confined space. If a shock wave crosses the content of the confined space, the body will accelerate and rotate around its axis. Specifically, the shock wave will increase the physical characteristics, especially its pressure, of the fluid acting on the impacted face of the door. The opposite side of the door is not influenced by the incident shock wave, only one of its faces is subjected to overpressure. Following the first impact, the resulting imbalance imposes a mechanical action on the door that will increase its speed and make it turn around its rotation axis. The difficulty comes when the door begins to open: the volume boundaries in which the fluid is contained are modified. Leaks occur and the gas kept in this closed volume can now flow to the atmosphere. Communication between the gas acting on each side of the door is created modifying their properties and consequently the pressure acting on each side of the door.The mechanical actions that apply to the door are no more the same with time, and therefore the acceleration of the door is changing. As the door moves, the fluid problem continues to be changed and in turn it changes its action on the door. This interaction process continues until either the limits of the problem ceases to be changed, the door cannot move, or when the mechanical actions acting on the door are in equilibrium, fluids on each side of the door are in the same physical state. The presented work is a study of the parameters of the fluid or the solid motion which are main actors in the behavioral law managing this complex system. In this aim, we designed an experimental device involving the physics that we have described and we have adapted it to a shock tube. Testing many experimental configurations, we could determine how the internal flow of a shock tube evolves when the end of this shock tube is more or less open.How a closed door reacts to the impact of a shock wave and what are the implications for the evolution of the involved fluids? What are the consequences of a different position of the door at the instant of the impact with the incident shock wave? What role plays the intensity of the incident shock wave or the inertia of the door on this dynamic?
|
113 |
Unsteady dynamics of corner separation in a linear compressor cascade / Dynamiques instationnaires du décollement de coin dans une grille d’aube de compresseurZambonini, Gherardo 02 December 2016 (has links)
Le présent travail concerne l'étude du décollement de coin de compresseurs par des investigations expérimentaux sur le banc d'essai subsonique grille d'aubes situé au LMFA (Re = 3.8 x 105, M = 0,12, profil NACA 65-009). Habituellement, cette particulière séparation tridimensionnelle a lieu dans le coin entre l'aubage et la paroi du moyeu des rangées d'aubes, également dans les stators et les rotors.Les performances de l’étage sont ainsi dégradées à cause des pertes de charge et du blocage conséquent à la séparation de l’écoulement. Bien que les caractéristiques stationnaires sont bien connues par la littérature, uniquement des récentes études expérimentales avancées combinées aux améliorations de simulations numériques, telles que URANS et LES, ont permis de découvrir le comportement très instable du décollement de coin. Des études préalables sur le même banc d'essai ont rapporté un comportement intermittent du décollement, appelé comportement bimodal. Dans la présente thèse de doctorat, il est montré que le comportement bimodal correspond à deux états spécifiques du champ fluide: une séparation fermée, presque supprimée, et une séparation ouverte caractérisée par un blocage massif. Il est clair que cette commutation bimodale du décollement de coin apparaissant dans une machine réelle pourrait avoir un impact fortement déstabilisant. En utilisant la méthode de mesure PIV haute vitesse couplée à des mesures de pression résolues en temps sur la surface de l'aubage, l'écoulement d'un canal inter-aube a été étudié pour deux différentes incidences. Les mesures PIV fournissent pour la première fois des visualisations résolues en temps et étendues à toute la section de l'aubage de la variation bimodale du décollement. L'interaction des grandes structures aléatoires de la couche limite à l'amont avec le bord d'attaque de l'aubage déstabilise le front du décollement et agrandit la région de recirculation. Une séparation ainsi massive persiste jusqu'à ce que le blocage dans le canal inter-aube provoque la rupture des plus grandes structures constituant la zone de recirculation. Successivement le débit recommence à accélérer et la séparation est presque supprimée. Cette dynamique coïncide avec le régime d'écoulement hautement instationnaire et apériodique des diffuseurs, appelé large transitory stall, caractérisé par des grandes amplitudes d'énergie dans la gamme des basses fréquences du spectre. Des moyennes conditionnelles de pression et la décomposition orthogonale modale (POD) des champs de vitesse ont été appliquées pour montrer l'effet rétroactif du blocage induit par la séparation ouverte sur l'angle local au bord d'attaque. Ces résultats supportent l’hypothèse d'une instabilité auto-entretenue causée par la diffusion imposée dans le stator. En fin, des mesures de pression totale résolues en temps ont été effectuées en utilisant des capteurs à haute réponse en fréquence, situés à la même position dans les différents canaux adjacents. La cohérence et la linéarité de l'angle de phase trouvée entre les capteurs confirment que l'instabilité du décollement peut se propager d'un canal a l'autre, en analogie avec les perturbations tournantes (rotating instabilities) apparaissant dans les bancs d'essai annulaires. Ceci montre pour la première fois que, même dans une grille d'aubes linéaire de compresseur, l'instabilité intrinsèque du décollement de coin peut déclencher la propagation d'instabilités. / The present work focuses on the study of the corner separation phenomenon in compressors carried out by experimental investigations on a subsonic linear cascade test rig (Re=3.8x105, M=0.12, blade profile NACA 65-009). Usually, this particular three-dimensional separation takes place in the corner between the blade and the endwall of compressor rows, mostly at hub, both in stators and rotors.Its main features are high total pressure losses and blockage of the flow, with consequent impacts on the efficiency. Whereas time averaged characteristics are well known from the past, only recent advanced experimental studies and improvements of numerical simulations, such as URANS and LES, have permitted to uncover the highly unsteady behavior of corner separation in compressors. Precedent studies on the same test rig have reported an intermittent unsteady behavior of corner separation, called bimodal behavior. In the present thesis it is shown that the bimodal behavior corresponds to two specific states of the flow: a closed separation, almost suppressed, and an open separation characterized by massive blockage and losses. Clearly hub-separation bimodal switches appearing in a real machine could have a first order detrimental effect on the stability of the flow in the compressor. By using high speed PIV coupled with unsteady pressure measurements on the surface of the blade the flow in a single blade passage has been investigated for different incidences. The PIV measurements provide, for the first time, time-resolved flow visualizations of the size switch of the separation with an extended field of view covering the entire blade section. The interaction of random large structures of the incoming boundary layer with the blade is found to be a predominant element that destabilizes the separation boundary and enlarges the recirculation region. Such a massive separation persists until the blockage in the passage causes the breakdown of the largest structures in the aft part of the blade, reestablishing the closed separation state. Such dynamics coincide with the aperiodic intermittent flow regime of diffusers, called transitory stall regime, and the associated Fourier spectra show the largest energy amplitudes in the low frequency range. Conditional ensemble averages of pressure and proper orthogonal decomposition (POD) of velocity fields have been applied to show the feedback effect of the blockage of the separation on the flow angle around the blade leading edge. These results draw the picture of a self-sustained instability caused by the diffusion imposed by the inter-blade passage. To answer the question about the interaction between adjacent corner separations, time-resolved total pressure measurements have been carried out by using high frequency response sensors positioned in bimodal points of multiple passages. The coherent propagation velocity and the linearity of the phase angle found between the signals confirm that the unsteadiness of the separation can propagate in pitch-wise direction. It is interesting to underline that equivalent elements characterize rotating disturbances appearing in annular test rigs. This finally shows that, even in an isolated stator blade row, the intrinsic unsteadiness of corner separation can start the propagation of instabilities. It is the first time that such a propagation effect is observed in a linear compressor cascade.
|
114 |
Étude d’une méthode d’éléments finis d’ordre élevé et de son hybridation avec d’autres méthodes numériques pour la simulation électromagnétique instationnaire dans un contexte industriel / Study of a high-order finite element method and its hybridization with order numerical methods for unsteady electromagnetic simulation in an industrial contextDeymier, Nicolas 08 December 2016 (has links)
Dans cette thèse, nous nous intéressons à l’amélioration du schéma de Yee pour traiter demanière plus efficace et pertinente les problèmes industriels auxquels nous sommes confrontés à l’heureactuelle. Pour cela, nous cherchons avant tout à diminuer les erreurs numériques de dispersion et àaméliorer les modélisations des géométries courbes ainsi que des réseaux de câbles. Pour répondre àces besoins, une solution basée sur un schéma Galerkin discontinu pourrait être envisagée. Toutefois,l’utilisation d’une telle technique sur la totalité du volume de calcul est relativement coûteuse. De plus,la prise en compte de structures filaires sur un tel schéma n’est pas encore opérationnelle. C’est pourquoi,dans l’optique d’avoir un outil industriel, et après une étude bibliographique, nous nous sommes plutôtorientés sur l’étude d’un schéma éléments finis (FEM) sur maillage cartésien qui possède toutes lesbonnes propriétés du schéma de Yee. Notamment, à l’ordre d’approximation spatiale égal à 0 ce schémaFEM est exactement le schéma de Yee, et, pour des ordres supérieurs, il permet de réduire fortementl’erreur de dispersion numérique de ce dernier. Dans le travail de cette thèse, pour ce schéma, nous avons notamment donné un critère de stabilité théorique, étudié sa convergence théorique et fait une analyse de l’erreur de dispersion. Pour tenircompte des possibilités d’ordre d’approximation spatiale variable par direction, nous avons mis en placeune stratégie d’affectation des ordres suivant le maillage donné. Ceci nous a permis d’obtenir un pas detemps optimal pour une précision souhaitée tout en réduisant les coûts de calcul. Après avoir porté ceschéma sur des machines de production, différents problèmes de CEM, antennes, IEM ou foudre ont ététraités afin de montrer les avantages et le potentiel de celui-ci. En conclusion de ces expérimentationsnumériques, il s’avère que la méthode est limitée par le manque de précision pour prendre en comptedes géométries courbes. Afin d’améliorer cela, nous avons proposé une hybridation entre ce schéma et leschéma GD que l’on peut étendre aux autres schémas comme les méthodes différences finies (FDTD) etvolumes finis (FVTD). Nous avons montré que la technique d’hybridation proposée conserve l’énergie etest stable sous une condition que nous avons évaluée de manière théorique. Des exemples de validationont ensuite été montrés. Enfin, pour tenir compte des réseaux de câbles, un modèle de fils minces d’ordred’approximation spatiale élevé a été proposé. Malheureusement, celui-ci ne peut pas couvrir l’ensembledes cas industriels et pour remédier à cela, nous avons proposé une hybridation de notre approche avecune équation de ligne de transmission. L’intérêt de cette hybridation a été montré sur un certain nombred’exemples, que nous n’aurions pas pu traiter par un modèle de structure filaire simple. / In this thesis, we study the improvement of the Yee’s scheme to treat efficiently and in arelevant way the industrial issues we are facing nowadays. For that, we first of all try to reduce thenumerical errors of dispersion and then to improve the modeling of the curved surfaces and of theharness networks. To answer these needs, a solution based on a Galerkin Discontinuous (GD) methodhas been first considered. However, the use of such method on the entire modeling volume is quite costly ;moreover the wires are not taken into account in this method. That is the reason why, with the objectiveof an industrial tool and after a large bibliographic research, we headed for the study of finite elementsscheme (FEM) on a Cartesian mesh which has all the good properties of the Yee’s scheme. Especially,this scheme is exactly the Yee’s scheme when the spatial order of approximation is set to zero. Forthe higher orders, this new scheme allows to greatly reduce the numerical error of dispersion. In theframe of this thesis and for this scheme, we give a theoretical criterion of stability, study its theoreticalconvergence and we perform an analysis of the error of dispersion. To take into account the possibilityof the variable spatial orders of approximation in each direction, we put in place a strategy of orderaffectation according to the given mesh. This strategy allows to obtain an optimal time step for a givenselected precision while reducing the cost of the calculations. Once this new scheme has been adaptedto large industrial computing means, different EMC, antennas, NEMP or lightning problems are treatedto demonstrate the advantages and the potential of this scheme. As a conclusion of these numericalsimulations we demonstrate that this method is limited by a lack of precision when taking into accountcurved geometries. To improve the treatment of the curved surfaces, we propose an hybridization between this scheme andthe GD scheme. This hybridization can also be applied to other methods such as Finite Differences(FDTD) or Finite Volumes (FVTD). We demonstrate that the technique of hybridization proposed,allows to conserve the energy and is stable under a condition that we study theoretically. Some examplesare presented for validation. Finally and to take into account the cables, a thin wire model with a highorder of spatial approximation is proposed. Unfortunately, this model does not allow to cover all theindustrial cases. To solve this issue we propose an hybridization with a transmission line method. Theadvantage of this hybridization is demonstrated thanks to different cases which would not have beenfeasible with a more simple thin wire method.
|
115 |
Modélisation instationnaire de l'aérodynamique externe automobile / Unsteady computation of external aerodynamics flow in automotive industryDelassaux, François 20 December 2018 (has links)
La thèse a pour but de développer une méthodologie de calcul instationnaire permettant une étude qualitative et quantitative de l’aérodynamique externe d’une automobile. La première partie de l’étude est consacrée au développement de la méthodologie numérique sur les différents corps d’Ahmed à 25°, géométries simplifiées d’une automobile réelle, afin de valider les choix stratégiques de maillages et de méthodes numériques ainsi que s’assurer de la bonne résolution de l’écoulement. Les résultats numériques sont comparés aux données expérimentales obtenues au cours d’essais réalisés à la soufflerie La Ferté Vidame lors de travaux de thèse précédents. A l’issue de ce travail, la méthode hybride Delayed Detached Eddy Simulation Shear-Stress Transport (DDES SST) est sélectionnée pour la suite de l’étude au vu des meilleures performances obtenues (torseur aérodynamique, coefficient de pression, topologie d’écoulement) par rapport aux méthodes Reynolds Averaged Navier-Stokes (RANS), Scale Adaptive Simulation (SAS) et Stress-Blended Eddy Simulation (SBES). La seconde phase de l’étude consiste à adapter la méthodologie précédemment développée sur un véhicule réel, la Peugeot 308 SW. Au préalable, une base de données expérimentales conséquente a été réalisée au sein du Groupement d’Intérêt Economique Souffleries Aéroacoustiques Automobiles (GIE S2A) au cours de ces travaux. La géométrie est tout d’abord simplifiée afin de faciliter la mise en place de la méthodologie numérique : entrées d’air fermées, soubassement lissé, roues remplacées par des carénages. Les résultats obtenus sont encourageants et démontrent globalement la supériorité de la DDES par rapport aux méthodes RANS classiques. La topologie d’écoulement est mieux prédite (soubassement et sillage), même si la prédiction du coefficient de portance reste une difficulté majeure pour ce type de méthode hybride. / The main goal of this PhD is to develop an unsteady numerical method to study the external aerodynamic flow around real vehicles. The first part of the study focuses on the flow around simplified geometries, such as 25° Ahmed bodies (with sharp and rounded edges on the back of the body), in order to determine the optimal turbulence model, mesh setup and numerical parameters. Computational Fluid Dynamics (CFD) results are compared to experimental data reported in literature conducted in the La Ferté Vidame wind tunnel. Based on this study, Shear-Stress Transport Delayed Detached Eddy Simulation (SST DDES) demonstrates superiority over Reynolds Averaged Navier-Stokes (RANS), Scale Adaptive Simulation (SAS) and Stress Blended Eddy Simulation (SBES) turbulence models, regarding both drag and lift coefficients predictions, and flow topology.Secondly, the numerical procedure is adapted for a real vehicle, the Peugeot 308 SW estate car. A substantial experimental campaign was carried out in the Groupement d’Intérêt Economique Souffleries Aéroacoustiques Automobiles (GIE S2A) wind tunnel to provide data against which the numerical results are compared. Given the geometric complexity of a real vehicle, the car is simplified for this study as follows: the front air inlets are closed, the underbody is smoothed with additional panels and the wheels are replaced by fairings. DDES computations show encouraging results. A significant improvement of the flow topology is obtained with DDES compared to RANS models. However, the prediction of the lift coefficient remains a major difficulty with these hybrid methods.
|
116 |
Coupling of time integration schemes for compressible unsteady flows / Couplage de schémas temporels pour la simulation des écoulements compressibles instationnairesMuscat, Laurent 12 March 2019 (has links)
Dans ce travail, on s'intéresse au développement d'une méthode hybride qui couple spatialement les schémas d'intégration temporelle explicite et implicite. Afin de répondre aux contraintes du solveur industriel FLUSEPA, les schémas explicite Heun et implicite Crank-Nicolson ont été hybridés via un paramètre de transition : l'approche mise en place est appelée schéma AION. Cette dernière est étudiée en détails avec une attention particulière sur son comportement spectral et sa capacité à maintenir l'ordre de précision. On montre que le traitement hybride a d'intéressants comportements dissipatif et dispersif tout en empêchant la réflexion d'ondes parasites et en maintenant la précision attendue. De plus, l'approche hybride est validée sur plusieurs cas académiques à la fois pour les flux convectifs et pour les flux diffusifs. Et comme espéré, la méthode est plus intéressante en terme de temps de calcul que les méthodes standards d'intégration temporelle. Pour l'extension de cette approche à la méthode temporelle adaptative présente dans FLUSEPA, il a été nécessaire d'améliorer le traitement qui permet à la méthode d’être conservative tout en obtenant des propriétés spectrales acceptables. Finalement l'approche hybride a été aussi étendue pour la modélisation RANS/LES de la turbulence avec des temps de calcul intéressants tout en capturant la physique de l'écoulement / This work deals with the design of a hybrid time integrator that couples spatially explicit and implicit time integrators. In order to cope with the industrial solver of Ariane Group called FLUSEPA, the explicit scheme of Heun and the implicit scheme of Crank-Nicolson are hybridized using the transition parameter : the whole technique is called AION time integration. The latter is studied into details with special focus on spectral behaviour and on its ability to keep the accuracy. It is shown that the hybrid technique has interesting dissipation and dispersion properties while maintaining precision and avoiding spurious waves. Moreover, this hybrid approach is validated on several academic test cases for both convective and diffusive fluxes. And as expected the method is more interesting in term of computational time than standard time integrators. For the extension of this hybrid approach to the temporal adaptive method implemented in FLUSEPA, it was necessary to improve some treatments in order to maintain conservation and acceptable spectral properties. Finally the hybrid time integration was also applied to a RANS/LES turbulent test case with interesting computational time while capturing the flow physics.
|
117 |
Couplage de la vélocimétrie par images de particules en deux temps avec la décomposition en modes propres pour la caractérisation d'un écoulementFavelier, Thomas 28 February 2006 (has links) (PDF)
L'objet de cette étude est d'utiliser la décomposition en modes propres (POD) pour établir un modèle de la composante déterministe et aléatoire d'un écoulement turbulent présentant une instationnarité à grande échelle pseudo périodique.<br />Une étude expérimentale de l'écoulement bidimensionnel en moyenne en aval d'un cylindre semi-circulaire, par vélocimétrie par image de particules en deux temps (PIV2T) caractérise l'écoulement<br />Une analyse POD du champ de vitesse permet d'extraire les modes spatiaux et de définir un paramètre de phase décrivant l'instationnarité à grande échelle qui régit la partie déterministe. La modélisation de l'évolution temporelle des coefficients associés aux modes s'effectue par des fonctions soit harmoniques pour la partie déterministe, soit stochastiques pour la partie aléatoire.<br />La modélisation est en accord avec les mesures expérimentales des premiers moments statistiques en un point et des fonctions de corrélation spatio-temporelle du champ de vitesse.
|
118 |
Méthode d'éléments finis mixtes :application aux équations de la chaleur et de Stokes instationnairesKorikache, Réda 15 November 2007 (has links) (PDF)
Dans ce travail on se propose d'établir des estimations d'erreurs a priori pour les solutions approchées d'équations d'évolution obtenues par la méthode d'éléments finis mixte duale en espace et ce pour trois types de problèmes : le premier concerne le problème de Cauchy pour l'équation de diffusion de la chaleur, le second est le problème de Stokes instationnaire, et le dernier concerne le problème de Cauchy pour l'équation de diffusion de la chaleur mais avec un coefficient de diffusion aléatoire. Pour ces trois types de problèmes, il y a un certain nombre de raisons de préférer la méthode mixte duale en espace à une méthode classique en espace ; parmi elles la propriété fondamentale qu'est la conservation locale, et par suite globale, de certaines quantités physiques (la quantité de mouvement, la masse, la quantité de chaleur,...). Une autre raison bien connue pour adopter la méthode mixte duale en espace est qu'elle nous permet d'introduire des nouvelles variables : p(t) =grad u(t) le flux de chaleur à l'instant t pour l'équation de diffusion de la chaleur, p(t) = K ◊ u(t) le flux de chaleur à l'instant t pour l'équation de diffusion de la chaleur avec un coefficient de diffusion aléatoire K, ◊ dénotant le produit de Wick, σ = grad u(t) le tenseur gradient du champ des vitesses à l'instant t pour le problème de Stokes instationnaire, ces inconnues supplémentaires ayant un sens physique et une importance particulière pour plus d'une application. Il est donc important de disposer d'une méthode numérique donnant aussi de bonnes approximations de ces quantités.
|
119 |
Etude expérimentale de l'aéroélasticité d'une plaque oscillante impactée par une batterie de jets turbulentsNyirumulinga, Yohann 26 April 2011 (has links)
Les instabilités aéroélastiques de bandes d’acier constituent aujourd’hui l’un des problèmes majeurs dans les sections de refroidissement par jets impactants des lignes de recuit continues.En effet, le traitement thermique des nouveaux aciers nécessite de très fortes pentes de température impliquant constamment des augmentations de vitesse de soufflage susceptibles de mettre en jeu des instabilités aéroélastiques. Des flottements ainsi que des divergences de bande ont déjà été constatées et identifiées. Ces deux instabilités impliquent dans la plupart des cas des chocs entre la bande et les buses de soufflage ce qui engendre des défauts de surface sur la bande.Un banc d’essai a été conçu et fabriqué dans le but d’analyser ces instabilités et d’anticiper leur apparition. A partir d’observations, la dynamique structurelle de la bande a été simplifiée à un mode de rotation rigide. Le banc comporte une plaque oscillante en mouvement forcé.Celle-ci est impactée par un dispositif de plusieurs jets axisymétriques turbulents ayant une disposition identique à celle des tours industrielles. Les efforts aérodynamiques stationnaires et instationnaires agissant sur la plaque sont mesurés au moyen de capteurs de pression.L’impact de plusieurs jets en interaction crée de très importants gradients de pression sur la plaque il est donc nécessaire que la grille de prises de pression soit très fine pour que l’estimation des efforts aérodynamiques soit correcte. La plaque est donc instrumentée de 91capteurs de pression sur une surface de 18 cm². Elle peut également être translatée dans les ois directions de l’espace, ce qui permet d’obtenir la distribution des efforts instationnaires ainsi que des coefficients aéroélastiques sur une grande surface de plaque et à différentes distances d’impact.Les mesures de pression stationnaires ont permis d’établir les courbes d’évolution des efforts d’impact des jets sur la plaque en fonction de la distance jet-plaque ainsi que de la géométrie des buses. Les résultats ont permis de déterminer la stabilité statique de la plaque en mouvement de pompage. Les mesures de vitesses des jets libres ont été effectuées paranémométrie à fil chaud et ont permis de déterminer leurs propriétés statistiques.Les mesures de coefficients aéroélastiques sur la plaque en rotation ont été effectuées surune seule géométrie de soufflage, pour différentes vitesses réduites. Les résultats mettent en évidence l’importance des effets de bords sur la stabilité de plaque. Des méthodes de post traitements ont proposées afin d’extrapoler les résultats à différentes largeurs de bande. Ils sont confrontés aux travaux de Regardin et al. (réf. [1]) et mettent en évidence des désaccords avec le cas réel. Des suggestions sont apportées afin d’améliorer la représentativité du banc vis-à-vis des bandes industrielles. / Aeroelastic instabilities of steel strips impinged by arrays of cooling gas jets have becomeone of the main issues in cooling sections of continuous annealing lines. Indeed, the new steeltreatments require very high temperature variation rates which involve increases in jetvelocities that are likely to onset some aeroelastic instabilities. Strip flutter and divergencehave already been observed and identified. These two aeroelastic instabilities imply a strongrisk of contact with the blowing boxes, which can seriously blemish the strip.An experimental test rig was designed and built in order to analyze and predict of theseinstabilities. From observations, the strip’s structural dynamics was simulated by a rigidrotation mode. The rig included a forced oscillating plate which is impinged by an array ofaxisymmetric jets having the exact industrial geometry. The plate was instrumented withpressure sensors to measure the steady and unsteady surface pressures. The impingement ofinteracting jets creates very large pressure gradients on the plate and therefore a tight mesh ofpressure taps (91 over an 18cm² jet impingement surface) was necessary to allow a goodestimation of the aerodynamic loads The plate could also be moved in the three coordinatedirections as to obtain surface mappings of the unsteady jet forces and aeroelastic coefficientscan be obtained over a wide area and different jet-to-plate distances.The variation of the impinging aerodynamic forces was established as a function of the jetto-plate distance for different nozzle geometries. These results were used to determine the jetstatic stability in plunging motion. Velocity and turbulence measurements in free jets werecarried out using hot wire anemometry in order to determine their statistical properties.Aeroelastic coefficient measurements were carried out on the oscillating plate with onlyone nozzle geometry and different reduced velocities. Results show that the plate’s stability ismainly dependent on the boundary effects. Post processing methods are suggested in order toapply the results to larger plates. Results are compared to the data of Regardin et al. (réf. [1])and emphasize some discrepancies with respect to the real case. Finally some improvementsto the test-rig are suggested for it to be more representative of the industrial situation.
|
120 |
Etude expérimentale de l'aéroélasticité d'une plaque oscillante impactée par une batterie de jets turbulents / Experimental Approach to the Aeroelastic Behaviour of an Oscillating Plate Impinged by Arrays of Turbulent JetsNyirumulinga, Yohann 26 April 2011 (has links)
Les instabilités aéroélastiques de bandes d’acier constituent aujourd’hui l’un des problèmes majeurs dans les sections de refroidissement par jets impactants des lignes de recuit continues.En effet, le traitement thermique des nouveaux aciers nécessite de très fortes pentes de température impliquant constamment des augmentations de vitesse de soufflage susceptibles de mettre en jeu des instabilités aéroélastiques. Des flottements ainsi que des divergences de bande ont déjà été constatées et identifiées. Ces deux instabilités impliquent dans la plupart des cas des chocs entre la bande et les buses de soufflage ce qui engendre des défauts de surface sur la bande.Un banc d’essai a été conçu et fabriqué dans le but d’analyser ces instabilités et d’anticiper leur apparition. A partir d’observations, la dynamique structurelle de la bande a été simplifiée à un mode de rotation rigide. Le banc comporte une plaque oscillante en mouvement forcé.Celle-ci est impactée par un dispositif de plusieurs jets axisymétriques turbulents ayant une disposition identique à celle des tours industrielles. Les efforts aérodynamiques stationnaires et instationnaires agissant sur la plaque sont mesurés au moyen de capteurs de pression.L’impact de plusieurs jets en interaction crée de très importants gradients de pression sur la plaque il est donc nécessaire que la grille de prises de pression soit très fine pour que l’estimation des efforts aérodynamiques soit correcte. La plaque est donc instrumentée de 91capteurs de pression sur une surface de 18 cm². Elle peut également être translatée dans les ois directions de l’espace, ce qui permet d’obtenir la distribution des efforts instationnaires ainsi que des coefficients aéroélastiques sur une grande surface de plaque et à différentes distances d’impact.Les mesures de pression stationnaires ont permis d’établir les courbes d’évolution des efforts d’impact des jets sur la plaque en fonction de la distance jet-plaque ainsi que de la géométrie des buses. Les résultats ont permis de déterminer la stabilité statique de la plaque en mouvement de pompage. Les mesures de vitesses des jets libres ont été effectuées paranémométrie à fil chaud et ont permis de déterminer leurs propriétés statistiques.Les mesures de coefficients aéroélastiques sur la plaque en rotation ont été effectuées surune seule géométrie de soufflage, pour différentes vitesses réduites. Les résultats mettent en évidence l’importance des effets de bords sur la stabilité de plaque. Des méthodes de post traitements ont proposées afin d’extrapoler les résultats à différentes largeurs de bande. Ils sont confrontés aux travaux de Regardin et al. (réf. [1]) et mettent en évidence des désaccords avec le cas réel. Des suggestions sont apportées afin d’améliorer la représentativité du banc vis-à-vis des bandes industrielles. / Aeroelastic instabilities of steel strips impinged by arrays of cooling gas jets have becomeone of the main issues in cooling sections of continuous annealing lines. Indeed, the new steeltreatments require very high temperature variation rates which involve increases in jetvelocities that are likely to onset some aeroelastic instabilities. Strip flutter and divergencehave already been observed and identified. These two aeroelastic instabilities imply a strongrisk of contact with the blowing boxes, which can seriously blemish the strip.An experimental test rig was designed and built in order to analyze and predict of theseinstabilities. From observations, the strip’s structural dynamics was simulated by a rigidrotation mode. The rig included a forced oscillating plate which is impinged by an array ofaxisymmetric jets having the exact industrial geometry. The plate was instrumented withpressure sensors to measure the steady and unsteady surface pressures. The impingement ofinteracting jets creates very large pressure gradients on the plate and therefore a tight mesh ofpressure taps (91 over an 18cm² jet impingement surface) was necessary to allow a goodestimation of the aerodynamic loads The plate could also be moved in the three coordinatedirections as to obtain surface mappings of the unsteady jet forces and aeroelastic coefficientscan be obtained over a wide area and different jet-to-plate distances.The variation of the impinging aerodynamic forces was established as a function of the jetto-plate distance for different nozzle geometries. These results were used to determine the jetstatic stability in plunging motion. Velocity and turbulence measurements in free jets werecarried out using hot wire anemometry in order to determine their statistical properties.Aeroelastic coefficient measurements were carried out on the oscillating plate with onlyone nozzle geometry and different reduced velocities. Results show that the plate’s stability ismainly dependent on the boundary effects. Post processing methods are suggested in order toapply the results to larger plates. Results are compared to the data of Regardin et al. (réf. [1])and emphasize some discrepancies with respect to the real case. Finally some improvementsto the test-rig are suggested for it to be more representative of the industrial situation.
|
Page generated in 0.0333 seconds