Le présent travail de thèse porte sur la modélisation et la simulation, en domaine temporel, de l'interaction entre un faisceau d'électrons et une onde hyperfréquence dans la structure à onde lente d'un TOP à hélice. Puisque le TOP est un instrument surdimensionné, les modèles non-stationnaires généraux utilisés dans les codes commerciaux nécessitent de trop grosses ressources de calcul pour pouvoir être utilisés en un temps raisonnable dans un but de conception. Il est donc nécessaire de faire appel à des modèles spécialisés. Durant cette thèse, nous nous sommes intéressés au "modèle discret non-stationnaire d'excitation d'un guide d'onde Périodique" de S. Kuznetsov. En 2007, N. Ryskin et al. avaient prouvé que ce modèle pouvait convenablement s'appliquer aux TOP à cavités couplées dans le cadre d'une application à une dimension du modèle. Lors de cette thèse, nous avons démontré, via le développement d'un code à une dimension (HelL-1D), que le modèle discret s'applique convenablement aux TOP à hélice. L'implémentation de ce modèle, dans un code à deux dimensions (HelL-2D) a, elle aussi, été effectuée. Enfin, nous avons développé une méthode permettant de contrôler de manière quantitative, dans le modèle discret, les phénomènes de réflexions aux extrémités de la ligne à retard, qui peuvent jouer un rôle important dans la stabilité de l'instrument. / This Ph.D. work deals with the time domain modeling and simulation of the electron beam/wave interaction in the slow-wave structure of a helix traveling-wave tube. Since a TWT is a device of which the geometry is oversized then the commercial software based on non-stationary general models needs so much computational resources that it cannot be used for design activities. During this Ph.D., we focused on the so called Kuznetsov's "discrete model of excitation of a periodic waveguide" which is a specialized model of beam/wave interaction in TWT. By 2007, N. Ryskin et al., showed that this model could conveniently apply to TWT with coupled cavities structure in one dimension. During this thesis, we first demonstrated that the discrete model could also apply to helix TWT with a sufficient (1%) accuracy via the development of a one-dimensional software called HelL-1D. We also implemented the discrete model for helix TWT in two dimensions (HelL-2D code). Finally, we developed a method, which is an extension of the discrete model, and which permits to take into account quantitatively the reflection phenomena at the terminations of a slow-wave structure in this model. This last study was very important since the stability of TWT strongly depends on this parameter.

Les tourbillons hélicoïdaux générés derrière les rotors sont étudiés. Pour les générer, une méthode basée sur le couplage entre la technique de la ligne active et un solveur des équations de Navier-Stokes (ENS), incompressibles et tridimensionnelles, a été développée. Elle consiste à modéliser la pâle par son équivalent de forces volumiques. Les équations, écrites en coordonnées cylindriques, sont résolues par un schéma de différences finies, écrit en parallèle. La méthode est d'ordre deux en temps et en espace. Le solveur des ENS a été validé par la reproduction des taux de croissance d'un écoulement de jet, instable, trouvés par la théorie d'instabilité linéaire. La comparaison avec des données expérimentales a montré que la méthode prédit bien l'aérodynamique de la pâle. Ensuite, le tourbillon de bout de pâle a été, en particulier, caractérisé. La vorticité et la vitesse azimutale ont été trouvées auto-similaire et la taille du coeur suit asymptotiquement la loi de diffusion linéaire 2D. Un modèle simple du coeur du tourbillon a été proposé. La présence d'une vitesse axiale dans le coeur du tourbillon a été montrée et a été caractérisée en fonction du rapport de vitesse au bout de la pâle. Finalement, une étude de stabilité du tourbillon a été faite en utilisant une vitesse angulaire variable pour perturber l'écoulement. Les taux de croissances des modes les plus instables sont en bon accord avec celui de l'instabilité d'appariement 2D des tourbillons. Trois types de modes ont été identifiés en fonction de la fréquence des perturbations et ont été trouvés similaires aux modes décrits par la théorie et aussi trouvés, précédemment, par l'expérience. / This present work is aimed to study helical vortices encountered in the wakes of rotating elements. For this, the generation of a helical wake of a one-bladed-rotor in a laminar velocity field, is simulated by the actuator line method. This method is a coupling of a Navier-Stokes (NS) solver with the Actuator Line Method where the blade is replaced by the body forces. This method has been implemented in a finite difference code, that we have written in parallel to solve the 3D incompressible NS equations written in cylindrical coordinates. The order of accuracy of the method is two both in time and space. The NS solver was validated comparing growth rates of an unstable jet, found numerically, and those of linear instability theory. A good agreement was found. A good agreement was also found comparing numerical results to analytical formulations and experimental data. It was shown that the method predicts well the blade aerodynamics . Then, the helical tip vortex is characterized for different Reynolds numbers and Tip Speed Ratios. The vorticity and the azimuthal velocity were found self-similar and the vortex core follows asymptotically the linear 2D diffusion law. A simple model for the helical vortex core was proposed. The presence of an axial velocity inside the vortex core was highlighted. Then, a stability study of the helical tip vortex was done using an angular velocity dependent on time to perturb the flow. The largest growth rates were found in good agreement with those of the (2D) pairing instability. Three types of modes were identified based on the perturbation frequency. The results are similar to those found in previous analytical and experimental works.

Cette thèse propose une méthodologie pour des études en simulation numérique en grande dimension. Elle se décline en différentes étapes : construction de plan d'expériences approprié, analyse de sensibilité et modélisation par surface de réponse. Les plans d'expériences adaptés à la simulation numérique sont les "Space Filling Designs", qui visent à répartir uniformément les points dans l'espace des variables d'entrée. Nous proposons l'algorithme WSP pour construire ces plans, rapidement, avec de bons critères d'uniformité, même en grande dimension. Ces travaux proposent la construction d'un plan polyvalent, qui sera utilisé pour les différentes étapes de l'étude : de l'analyse de sensibilité aux surfaces de réponse. L'analyse de sensibilité sera réalisée avec une approche innovante sur les points de ce plan, pour détecter le sous-ensemble de variables d'entrée réellement influentes. Basée sur le principe de la méthode de Morris, cette approche permet de hiérarchiser les variables d'entrée selon leurs effets. Le plan initial est ensuite "replié" dans le sous-espace des variables d'entrée les plus influentes, ce qui nécessite au préalable une étude pour vérifier l'uniformité de la répartition des points dans l'espace réduit et ainsi détecter d'éventuels amas et/ou lacunes. Ainsi, après réparation, ce plan est utilisé pour l'étape ultime : étude de surfaces de réponse. Nous avons alors choisi d'utiliser l'approche des Support Vector Regression, indépendante de la dimension et rapide dans sa mise en place. Obtenant des résultats comparables à l'approche classique (Krigeage), cette technique semble prometteuse pour étudier des phénomènes complexes en grande dimension. / This thesis proposes a methodology of study in numeric simulation for high dimensions. There are several steps in this methodology : setting up an experimental design, performing sensitivity analysis, then using response surface for modelling. In numeric simulation, we use a Space Filling Design that scatters the points in the entire domain. The construction of an experimental design in high dimensions must be efficient, with good uniformity properties. Moreover, this construction must be fast. We propose using the WSP algorithm to construct such an experimental design. This design is then used in all steps of the methodology, making it a versatile design, from sensitivity analysis to modelling. A sensitivity analysis allows identifying the influent factors. Adapting the Morris method principle, this approach classifies the inputs into three groups according to their effects. Then, the experimental design is folded over in the subspace of the influent inputs. This action can modify the uniformity properties of the experimental design by creating possible gaps and clusters. So, it is necessary to repair it by removing clusters and filling gaps. We propose a step-by-step approach to offer suitable repairing for each experimental design. Then, the repaired design is used for the final step: modelling from the response surface. We consider a Support Vector Machines method because dimension does not affect the construction. Easy to construct and with good results, similar to the results obtained by Kriging, the Support Vector Regression method is an alternative method for the study of complex phenomena in high dimensions.

Une courbe de fragilité sismique qui présente la probabilité de défaillance d’une structure en fonction d’une intensité sismique, est un outil performant pour l’évaluation de la vulnérabilité sismique des structures en génie nucléaire et génie civil. On se concentre dans cette thèse sur l’approche par simulations numériques pour la construction des courbes de fragilité sismique. Une étude comparative des méthodes paramétriques existantes avec l’hypothèse log-normale est d’abord réalisée. Elle permet ensuite de proposer des améliorations de la méthode du maximum de vraisemblance dans le but d’atténuer l’influence de l’excitation sismique lors de son processus de construction. Une autre amélioration est l’application de la méthode de simulations par subsets pour l’évaluation de la probabilité de défaillance faible. Enfin, en utilisant la méthode de calcul de l’évolution des fonctions de densité de probabilité qui permet d’évaluer la probabilité conjointe entre la réponse structurale et les variables aléatoires du système et de l’excitation, nous proposons également une nouvelle technique non-paramétrique de construction des courbes de fragilité sismique sans utiliser l’hypothèse de la loi log-normale. La validation des améliorations et de la nouvelle technique est réalisée sur des exemples numériques. / A seismic fragility curve that shows the failure probability of a structure in function of a seismic intensity is a powerful tool for the evaluation of the seismic vulnerability of structures in nuclear engineering and civil engineering.We focus in this thesis on the numerical simulations-based approach for the construction of seismic fragility curves. A comparative work between existent parametric methods with the lognormal assumption of the fragility curves is first performed. It then allows proposing improvements to the maximum likelihood method in order to mitigate the influence of seismic excitation during its construction process. Another improvement is the application of the subsets simulation method for the evaluation of the low probability of failure. Finally, using the Probability Density Evolution Method (PDEM) for evaluating the joint probability between a structural response and random variables of a system and/or excitations, a new technique for construction of seismic fragility curves was proposed. The seismic fragility curve can be derived without the assumption of lognormal law. The improvements and the new technique are all validated by numerical examples.

Procédé de moulage par injection de poudres est constitué de quatre étapes principales: la préparation des matières premières, moulage par injection, le déliantage et le frittage. Cette thèse présente les recherches sur deux aspects principaux: la micro-injection et frittage par micro -ondes. Les contributions principaux peuvent être conclues dans les quatre aspects suivants: Modification et complément de l'algorithme précédent pour la simulation du procédé de moulage par injection; L'évaluation et la mise en œuvre de l'effet de tension de surface en simulation pour micro-injection; Micro-ondes expériences de frittage de compacts basés sur l'acier inoxydable 17-4PH; Réalisation de la simulation de frittage à micro-ondes avec couplage de la multi-physique, y compris le chauffage à micro-ondes classique, le transfert de chaleur, et le supplément de modèle pour la densification de frittage de la poudre compacté / Powder Injection molding process consists off our main stages: feedstock preparation, injection molding, debinding and sintering. The thesis presents the research on two main aspects: micro injectionmolding and microwave sintering. The main contributions can be concluded in thefollowing four aspects: Modification and supplement of previous algorithm for the simulation ofinjection molding process; Evaluation and implementation of surface tension effect in simulation for micro injection; Microwave sintering experiments of compacts based on 17-4PH stainles ssteel; Realization of the microwave sintering simulation with the coupling of multi-physics,including the classic microwave heating, heat transfer, and the supplement of model for sintering densification of powder impacts

Le présent travail de recherche concerne l'étude de la convection de Rayleigh-Bénard-Marangoni dans un récipient cylindrique doté d'un fond en substrat solide. Le substrat solide est chauffé sur sa face inférieure par un flux de chaleur localisé. L'étude comporte deux parties : La première partie du travail consiste en une modélisation physique du problème associée à des simulations numériques. Les équations de Navier-Stokes et de l'énergie sont résolues en 3D par une méthode de volumes finis. Un transfert de chaleur conjugué solide-liquide est considéré. Des morphologies originales de cellules (type et nombre) sont observées, elles dépendent des conditions géométriques, des nombres adimensionnels qui régissent la physique de l'écoulement (nombre de Prandtl, de Rayleigh et de Marangoni ainsi que du rapport des conductivités thermiques du substrat solide et du fluide). Les transferts de chaleur sont aussi évalués pour chaque cas d'étude. Dans la deuxième partie, nous allons détaillons une étude expérimentale de la convection de Rayleigh-Bénard-Marangoni dans la même configuration que celle étudiée numériquement. Les structures convectives et leurs évolutions sont étudiées à partir d’images relevées par thermographie infra-rouge. Différents modes d'organisation des cellules convectives ont pu être mis en évidence pour ce type de chauffage à flux thermique imposé non uniforme. / The present research work concerns the study of Rayleigh-Bénard-Marangoni convection in a cylindrical container with a solid substrate base. This solid substrate is heated by a localized heat flux on its underside. The study is divided into two parts : The first part of the work consists of a physical modelling of the problem associated with numerical simulations. The Navier-Stokes and energy equations are solved by using a 3D finite volume method. A conjugate solid-liquid heat transfer is considered. Original morphology of cells (type and number) are observed, they are linked to the geometrical conditions, the dimensionless numbers which govern the physical problem (Prandtl, Rayleigh and Marangoni numbers and the ratio of solid substrate to liquid thermal conductivities). The heat transfers are also evaluated in each case. In the second part of the work, we present an experimental study of Rayleigh-Bénard-Marangoni convection in the same configuration as that studied numerically. Convective structures and their evolutions are studied from images recorded by infrared thermography. Different modes of organization of convective cells have been highlighted for this type of heating with imposed non-uniform heat flux.

Plusieurs études tant numériques qu’expérimentales font état de la présence d’instabilités thermiques dans des films liquides chauffés uniformément par le bas pour des conditions aux limites et d’écoulements particuliers. La présence de ces instabilités modifiera les transferts thermiques associés. Le sujet de ce travail de thèse consiste à étudier numériquement les instabilités thermoconvectives d’un écoulement laminaire tridimensionnel de convection mixte dans un canal horizontal à surface libre. Les variations de la tension de surface avec la température (effet Marangoni ou effet thermocapillaire) sont prises en compte. Bien que d’un intérêt certain pour de nombreuses applications industrielles, cette situation a été très peu étudiée d’un point de vue académique dans la configuration considérée ici. Dans cette de configuration plusieurs types de structures thermoconvectives sont susceptibles d’apparaître. Lorsque les forces induites par les courants de convection naturelle, forcée et thermocapillaire sont du même ordre de grandeur, les premiers résultats montrent un développement des instabilités sous forme de rouleaux convectifs longitudinaux stationnaires semblables à ceux rencontrés pour des écoulements de type Poiseuille-Rayleigh-Bénard. A notre connaissance, c’est la première fois que l’écoulement de convection de type Poiseuille-Rayleigh-Bénard associé aux effets Marangoni est étudié. Le nombre et la distribution spatiale des rouleaux convectifs le long du canal dépendent des conditions de l’écoulement. Nous proposons une étude numérique pour ces conditions particulières d’écoulement pouvant conduire à des instabilités avec une évaluation de leur effet sur les transferts de chaleur. Les équations de Navier-Stokes et de l’énergie sont résolues numériquement par la méthode de volumes finis en prenant en compte les effets thermocapillaires. Les résultats présentés concernent l’influence des paramètres contrôlant l’écoulement (nombres de Reynolds, de Rayleigh, de Biot, de Marangoni et le rapport de forme) sur les motifs de l’écoulement et les échanges thermiques. Dans une seconde partie du travail, complémentaire à la première, une analyse de stabilité linéaire de l’écoulement dans un canal ouvert à surface libre d’extension latérale infinie est réalisée en utilisant la méthode spectrale de type collocation Chebyshev pour résoudre un système aux valeurs propres. Les diagrammes de stabilité déterminant les seuils des paramètres conduisant à l’instabilité thermoconvective ont été obtenus et analysés, ainsi que les structures spatiales associées. / Several studies both numerical and experimental have reported the presence of thermal instabilities in liquid films uniformly heated from below for specific boundary conditions and flows. The presence of these instabilities modifies the associated heat transfer. The subject of this PhD thesis is to study numerically the instability of three-dimensional laminar mixed convection within a liquid flowing on a horizontal channel heated uniformly from below. The upper surface is free and assumed to be flat. The variations of the surface tension with the temperature (Marangoni effect or thermocapillary effect) are taken into account. Although of great interest for many industrial applications, this problem has received little attention from an academic point of view. In this configuration, several types of thermoconvective structures may appear. When the strength of the buoyancy, thermocapillary effects and forced convective currents are comparable, the results show the development of instabilities in the form of steady longitudinal convective rolls similar to those encountered in the Poiseuille-Rayleigh-Bénard flow. To our knowledge, this is the first time that the Poiseuille-Rayleigh-Bénard flow associated to the Marangoni effects has been investigated. The number and spatial distribution of the convective rolls along the channel depend on the flow conditions. We propose a numerical study on the flow conditions that could lead to thermal instabilities with an evaluation of their effect on the heat transfer. The coupled Navier-Stokes and energy equations are solved numerically by the finite volume method taking into account the thermocapillary effects. The results presented concern the influence of several control parameters (the Reynolds, Rayleigh, Biot and Marangoni numbers and the aspect ratio of the channel) on the flow patterns and heat transfer characteristics. In the second part of this work, complimentary to the first, a linear stability analysis of a horizontal liquid film flowing in an open channel, with infinite lateral extension and uniform heating from below, is carried out. An eigenvalue problem is obtained in the course of this analysis which is solved numerically using the Chebyshev collocation spectral method. The stability diagrams determining the threshold parameters leading to thermoconvective instabilities were obtained and analyzed as well as the associated spatial patterns.

La compréhension des phénomènes thermiques et mécaniques mis en jeu lors de la mise en forme des matériaux est généralement réalisée avec l’aide de simulations numériques. Ces simulations montrent leurs limites pour les procédés qui conduisent à de très grandes déformations de la matière. Dans ce cas, de très fortes distorsions du maillage se produisent pendant le calcul, entrainant une augmentation de l’erreur, voire l’arrêt prématuré de la simulation. Cette étude porte sur le développement d’une stratégie de remaillage adaptative afin d’éviter les distorsions des éléments pendant les simulations en grandes transformations. La méthode proposée a été intégrée dans un environnement de calcul utilisant le solveur ABAQUS/Explicit, un mailleur 3D et un algorithme de transfert de champ.La méthode h-adaptative en combinaison avec un critère de contrôle basé sur l’endommagement et un estimateur d’erreur de type Zienkiewicz-Zhu Z2 (SPR-amélioré) ont été implantés. Le maillage initial est remplacé par un nouveau maillage avec le niveau de qualité désiré par l’utilisateur, tout en minimisant le nombre des degrés de liberté. Cette technique s’est montrée robuste et entièrement automatique pour déterminer la taille optimale des nouveaux éléments. Une fois le nouveau maillage généré, toutes les variables doivent être soigneusement transférées. Plusieurs techniques de transfert sont décrites et comparées. Des améliorations permettant d’augmenter leurs efficacités en termes de diffusion de l’information et de stabilité numérique ont été proposées. Une attention particulière est portée à la restauration de l'équilibre mécanique local du système. Les différentes techniques développées ont permis de modéliser différents procédés entrainant de grandes déformations élastoplastiques avec endommagement. Dans toutes les applications testées, il a été montré une amélioration de la précision et de la qualité des résultats numériques obtenus. Pour des opérations d’usinage, des mesures de champs cinématiques à travers la technique de corrélation d’images ont été réalisées afin de déterminer les champs de déformation en pointe d’outil. Ces mesures ont servi à la validation de la simulation numérique à l’échelle locale. La comparaison des champs cinématiques expérimentaux avec ceux issus du calcul éléments finis met en évidence la robustesse du processus d’adaptation du maillage proposée pour retranscrire les phénomènes locaux observés expérimentalement. En effet, la reproduction de l’écoulement de la matière sur les bords et la géométrie du copeau sont en très bonne corrélation avec les résultats expérimentaux. Ce développement a permis de proposer une description nouvelle du processus de formation des bandes de cisaillement. / In this work, a fully automated adaptive remeshing strategy, based on a tetrahedral element to simulate various 3D metal forming processes, was proposed. The aim of this work is to solve problems associated with the severe mesh distortion that occurs during the computation and which may be incompatible with the evolution of the physical behavior of the FE solution. Indeed, the quality of the mesh conditions affects the accuracy of the calculations. The proposed strategy is integrated in a computational platform which integrates a finite element solver (Abaqus/Explicit), 3D mesh generation and a field transfer algorithm.The base idea is to use the h-adaptive methodology in the combination with a damage-criterion error and Zienkiewicz-Zhu Z2 type error estimator (SPR-improved) to locally control the mesh modification-as-needed. Once a new mesh is generated, all history-dependent variables need to be carefully transferred between subsequent meshes. Therefore, different transfer techniques are described and compared. An important part of this work concerns the presentation of the proposed modification of the field transfer operator and a special attention is given to restore the local mechanical equilibrium of the system. During the large elasto-plastic deformation simulation with damage, the necessary steps for remeshing the mechanical structure are presented. The several types of applications are also given. For all studied applications, the above strategy can improve the accuracy and quality of numerical results. It also has benefits to decide how refined a mesh needs to be to reach a particular level of accuracy, or how coarse the mesh can be without unacceptably impacting solution accuracy.For the machining processes, kinematic field measurements using Digital image Correlation were performed to validate the numerical simulation at the local level. The comparison of the experimental kinematic fields and those resulting from the FE calculation highlights the robustness of the proposed mesh adaptation process which can transcribe the experimental local phenomena. Also, the reproduction of the material flow at the edges and the chip are correlated with the experimental results accurately. Finally, the physical study of the numerical results can be allowed to propose an innovative description of ASB formation.

Cette thèse de doctorat s’intéresse aux mécanismes d’amélioration des transferts dans des géométries de dissipateurs thermiques à plaques et ailettes. Une première partie est consacrée à l’étude d’une configuration académique à l’aide de simulations numériques visant à obtenir une amélioration du transfert de chaleur conjugué en modifiant uniquement par des découpes la forme géométrique des ailettes planes conductrices. Une analyse locale approfondie de l’écoulement et des champs thermiques a été effectuée avec notamment le principe de synergie locale, des champs de vitesse et de gradients thermiques, pour comprendre l’effet des modifications géométriques. Ce mémoire présente également le développement de dissipateurs aux performances thermo-aérauliques augmentées pour des applications de refroidissement de coffrets électronique embarqués. L’intensification des transferts thermiques est obtenue par la génération d’écoulements secondaires qui provoquent un brassage de fluide et réduisent la résistance thermique à la paroi en perturbant le développement de la couche limite thermique. Différentes configurations de dissipateurs avec deux types de générateurs d’écoulements secondaires, paires d’ailettes Delta et protrusions, ont été étudiées numériquement, en employant une modélisation de type « RANS ». Les performances thermo-aérauliques des géométries munies de générateurs de vorticité ont été comparées à celle d’un dissipateur thermique de référence « lisse ». Des prototypes ont également été fabriqués et testés sur un banc expérimental spécifiquement développé pour réaliser des mesures des performances globales en termes de puissance thermique et de pertes de charge. Les résultats expérimentaux et numériques ont été confrontés afin de qualifier les simulations réalisées. Par la suite, une étude d’optimisation employant l’analyse factorielle Taguchi a été utilisée afin d’optimiser les paramètres géométriques des dissipateurs retenus. Deux fonctions objectif ont été considérées : la maximisation du facteur de performance thermique à iso puissance de ventilation (PEC) et la réduction de la température moyenne de paroi du dissipateur par rapport au cas de référence. L’analyse des performances thermo-aérauliques globales des géométries étudiées a été complétée par une analyse qualitative locale des champs thermiques et d’écoulement notamment avec le principe de synergie. / This doctoral thesis focuses on mechanisms of heat transfer enhancement in plate and fin heat sink geometries. First part of the thesis is dedicated to study an academic configuration using numerical simulations to achieve an improvement in conjugate heat transfer by modifying only the geometrical shape (through punching) of the conductive plane fins. An in-depth local analysis of the flow and thermal fields was carried out with the local synergy principle, velocity and thermal gradients, to understand the effect of geometric modifications. This thesis also presents the development of heat sinks with increased thermo-hydraulic performance for on-board electronic box cooling applications. The intensification of the heat transfer is obtained by the generation of secondary flows which cause an intensive mixing of fluid and reduces the thermal resistance to the wall by disrupting the development of the thermal boundary layer. Different heat sink geometries with two types of secondary flow generators : delta winglet pair and protrusions were numerically studied using RANS approach. The thermo-hydraulic performances of the geometries equipped with vortex generators were compared with that of a smooth reference heat sink. The prototypes were also manufactured and tested on an experimental bench specifically designed to perform global performance measurements in terms of thermal power and pressure drops. Experimental and numerical results were compared to qualify the simulations performed. Subsequently, an optimization study using Taguchi factorial analysis was used to optimize the geometrical parameters of the chosen dissipaters. Two objective functions were considered : maximization of either iso-pumping power performance criteria (PEC) or average wall temperature of the dissipaters compared to the reference case. The global thermo-hydraulic performance analysis of the studied geometries was completed by a qualitative analysis of local flow and thermal fields, in particular with the local field synergy principle.

L’apparition accidentelle d’un arc électrique dans le système de distribution électrique d’un aéronef peut compromettre la sécurité du vol. Il existe peu de travaux liés à cette problématique.Le but de ce travail est donc d’étudier le comportement d’un arc électrique, en conditions aéronautiques,par des approches théorique, numérique, et expérimentale. Dans ce travail, un modèle MHD de la colonne d’arc à l’ETL a été utilisé, et résolu à l’aide du logiciel commercial comsolMultiphysics. Afin de décrire l’interaction plasma-électrodes, le modèle a dû étendu pour inclure les écarts à l’équilibre près des électrodes. Ces zones ont été prises en compte en considérant la conservation du courant et de l’énergie dans la zone hors-équilibre. L’approche choisie et le développement du modèle ont été détaillés. La validation du modèle dans le cas d’un arc libre a montré un excellent accord avec les résultats numériques et expérimentaux de la littérature.Ce modèle d’arc libre a été étendu au cas de l’arc se propageant entre des électrodes en configuration rails et en géométrie 3D. Une description auto-cohérente du déplacement de l’arc entre les électrodes a été réalisée. La simulation numérique a été faite pour des arcs en régimes DC, pulsé et AC à des pressions atmosphériques et inférieures. Les principales caractéristiques de l’arc ont été analysées et discutées. Les résultats obtenus ont été comparés avec les résultats expérimentaux et ont montré un bon accord.Ce modèle d’arc électrique est capable de prédire le comportement d’un arc de défaut dans des conditions aéronautiques. Des améliorations du modèle sont discutées comme perspectives de ce travail. / The ignition of an electric arc in the electric distribution system of an aircraft can be a serious problem for flight safety. The amount of information on this topic is limited, however. Therefore,the aim of this work is to investigate the electric arc behavior by means of experiment and numerical simulations.The MHD model of the LTE arc column was used and resolved numerically using the commercial software comsol Multiphysics. In order to describe plasma-electride interaction, the model had to be extended to include non-equilibrium effects near the electrodes. These zones were taken into account by means of current and energy conservation in the non-equilibrium layer. The correct matching conditions were developed and are described in the work. Validation of the model in the case of a free burning arc showed excellent agreement between comprehensive models and the experiment.This model was then extended to the case of the electric arc between rail electrodes in a 3D geometry. Due to electromagnetic forces the electric arc displaces along the electrodes. A self-consistent description of this phenomenon was established. The calculation was performed for DC, pulsed and AC current conditions at atmospheric and lower pressures. The main characteristics of the arc were analyzed and discussed. The results obtained were compared with the experimental measurements and showed good agreement.The model of electric arcs between busbar electrodes is able to predict the behavior of a fault arc in aeronautical conditions. Further improvements of the model are discussed as an outlook of the research.