La résolution des problèmes de contrôle optimal nécessite des temps de calcul et des capacités de stockage très élevés. Pour s’affranchir de ces contraintes, il est possible d’utiliser les méthodes de réduction de modèles comme la POD (Proper Orthogonal Decomposition). L’inconvénient de cette approche est que la base POD n’est valable que pour des paramètres situés dans un voisinage proche des paramètres pour lesquels elle a été construite. Par conséquent, en contrôle optimal, cette base peut ne pas être représentative de tous les paramètres qui seront proposés par l’algorithme de contrôle. Pour s’affranchir de cet handicap, une méthodologie de contrôle optimal utilisant des modèles réduits adaptatifs a été proposée dans ce manuscrit. Les bases réduites adaptées sont obtenues à l’aide de la méthode d’interpolation ITSGM (Interpolation on Tangent Subspace of Grassman Manifold) ou de la méthode d’enrichissement PGD (Proper Generalized Decomposition). La robustesse de cette approche en termes de précision et de temps de calcul a été démontrée pour le contrôle optimal (basé sur les équations adjointes) des équations 2D de réaction-diffusion et de Burgers. L’approche basée sur l’interpolation ITSGM a également été appliquée avec succès pour contrôler l’écoulement autour d’un cylindre 2D. Deux méthodes de réduction non intrusives, ne nécessitant pas la connaissance des équations du modèle étudié, ont également été proposées. Ces méthodes appelées NIMR (Non Intrusive Model Reduction) et HNIMR (Hyper Non Intrusive Model Reduction) ont été couplées à un algorithme génétique pour résoudre rapidement un problème de contrôle optimal. Le problème du contrôle optimal de l’écoulement autour d’un cylindre 2D a été étudié et les résultats ont montré l’efficacité de cette approche. En effet, l’algorithme génétique couplé avec la méthode HNIMR a permis d’obtenir les solutions avec une bonne précision en moins de 40 secondes. / The numerical resolution of adjoint based optimal control problems requires high computational time and storage capacities. In order to get over these high requirement, it is possible to use model reduction techniques such as POD (Proper Orthogonal Decomposition). The disadvantage of this approach is that the POD basis is valid only for parameters located in a small neighborhood to the parameters for which it was built. Therefore, this basis may not be representative for all parameters in the optimizer’s path eventually suggested by the optimal control loop. To overcome this issue, a reduced optimal control methodology using adaptive reduced order models obtained by the ITSGM (Interpolation on a Tangent Subspace of the Grassman Manifold) method or by the PGD (Proper Generalized Decomposition) method, has been proposed in this work. The robustness of this approach in terms of precision and computation time has been demonstrated for the optimal control (based on adjoint equations) of the 2D reaction-diffusion and Burgers equations. The interpolation method ITSGM has also been validated in the control of flow around a 2D cylinder. In the context of non intrusive model reduction, two non intrusive reduction methods, which do not require knowledge of the equations of the studied model, have also been proposed. These methods called NIMR (Non-Intrusive Model Reduction) and HNIMR (Hyper Non-Intrusive Model Reduction) were developed and then coupled to a genetic algorithm in order to solve an optimal control problem in quasi-real time. The problem of optimal control of the flow around a 2D cylinder has been studied and the results have shown the effectiveness of this approach. Indeed, the genetic algorithm coupled with the HNIMR method allowed to obtain the solutions with a good accuracy in less than 40 seconds.

The spectral-element method (SEM) is used to study wall-bounded turbulent flowsin moderately complex geometries. The first part of the thesis is devoted to simulations of canonical flow cases, such as temporal K-type transitionand turbulent channel flow, to investigate general resolution requirements and computational efficiency of the numerical code nek5000. Large-eddy simulation (LES) is further performed of a plane asymmetric diffuser flow with an opening angle of 8.5 degrees, featuring turbulent flow separation. Good agreement with numerical studies of Herbst (2007) is obtained, and it is concluded that the use of a high-order method is advantageous for flows featuring pressure-induced separation. Moreover, it is shown, both a priori on simpler model problems and a posteriori using the full Navier--Stokes equations, that the numerical instability associated with SEM at high Reynolds numbers is cured either by employing over-integration (dealiasing) or a filter-based stabilisation, thus rendering simulations of moderate to high Reynolds number flows possible. The second part of the thesis is devoted to the first direct numerical simulation (DNS) of a truly three-dimensional, turbulent and separated diffuser flow at Re = 10 000 (based on bulk velocity and inflow-duct height), experimentally investigated by Cherry et al. (2008). The massively parallel capabilities of the spectral-element method are exploited by running the simulations on up to 32 768 processors. Very good agreement with experimental mean flow data is obtained and it is thus shown that well-resolved simulations of complex turbulent flows with high accuracy are possible at realistic Reynolds numberseven in complicated geometries. An explanation for the discovered asymmetry of the mean separated flow is provided and itis demonstrated that a large-scale quasi-periodic motion is present in the diffuser. In addition, a new diagnostic measure, based on the maximum vorticity stretching component in every spatial point, is designed and tested in a number of turbulent and transitional flows. Finally, Koopman mode decomposition is performed of a minimal channel flow and compared to classical proper orthogonal decomposition (POD). / QC 20111206

Die vorwiegend experimentelle Arbeit befasst sich mit der systematischen Untersuchung von Parametervariationen bei der aktiven Strömungskontrolle mit elektromagnetischen Kräften. An einer angestellten Platte und einem NACA0015-Profil wurde die saugseitige abgelöste Strömung durch das Einbringen einer periodischen wandparallelen Lorentzkraft an der Vorderkante beeinflusst und experimentell mittels zeitaufgelöster Particle Image Velocimetry (PIV) untersucht. Dabei wurde für verschiedene Anstellwinkel und Reynoldszahlen die Frequenz der Anregung, deren Impulseintrag und der zeitliche Kraftverlauf variiert. Strömungsmechanische Untersuchungen experimenteller und numerischer Natur wurden für eine elektrochemische Zelle und den Fall der Elektrolyse an Millieelektroden unter dem Einfluss externer Magnetfelder durchgeführt. Die Übereinstimmung der gemessenen und berechneten Geschwindigkeitsfelder war dabei sehr gut. Entgegen der Annahme, dass im Falle homogener Magnetfelder keine Strömungen induziert werden, konnte nachgewiesen werden, dass durch die lokale Krümmung der elektrischen Feldlinien in Elektrodennähe starke Lorentzkräfte generiert werden. Dies führt zu sehr komplexen Primär-und Sekundärströmungen. Die gleichen Effekte bewirken ebenfalls in der Nähe von Millieelektroden starke Lorentzkräfte in homogenen magnetischen Feldern. Die experimentellen Beobachtungen an Millieelektroden von Leventis et. al (2005), welche zum Beweis der Konzentrationsgradientenkraft herangezogen wurden, konnten alle auf das Wirken lokaler Lorentzkräfte zurückgeführt werden. Der experimentelle Nachweis der Konzentrationsgradientenkraft steht damit weiterhin aus. Zur Messung der Konzentrationen in elektrochemischen Systemen wurde erstmals das Hintergrundschlierenverfahren angewendet. Dieses Verfahren erlaubt die Bestimmung der räumlichen Konzentrationsgradienten mit erheblich weniger messtechnischen Aufwand gegenüber spektroskopischen Methoden und der Schlierentechnik.

Pour les solutions transitoires répétées des structures géométriquement nonlinéaires l’effort numérique présente souvent une contrainte importante. Ainsi, l’introduction d’un modèle d’ordre réduit, qui prend en compte les effets nonlinéaires et qui accélère considérablement les calculs, s’avère souvent nécessaire.Ce travail aboutit à une méthode qui permet des solutions transitoires accélérées, fidèles et paramétrables, à travers d’un modèle réduit de la structure initiale. La structure est discrétisée et son équilibre dynamique décrit par une équation matricielle. La projection sur une base réduite est introduite afin d’obtenir un modèle réduit. Une étude numérique complète sur plusieurs bases communes démontre que la simple introduction d’une base constante ne suffit pas pour prendre en compte le comportement nonlinéaire. Trois exigences sont déduites pour une solution transitoire accélérée, fidèle et paramétrable. L’algorithme de solution doit permettre un suivi de l’évolution nonlinéaire de la solution transitoire, la solution doit être autonome des termes nonlinéaires en éléments finis et la base doit être adaptée à des paramètres externes.Trois approches sont mises en place, chacune répondant à une exigence. Ces approches sont assemblées dans la méthode intégrée. Les approches sont la mise-à-jour et augmentation de la base , la formulation polynomiale des termes nonlinéaires et l’interpolation de la base. Un algorithme de type Newmark forme le cadre de la méthode intégrée. L’application de la méthode intégrée sur des cas test en élément finis géométriquement nonlinéaires confirme qu’elle répond au but initial d’obtenir des solutions transitoires accélérées, fidèles et paramétrables. / For repeated transient solutions of geometrically nonlinear structures the numerical effort often poses a major obstacle. Thus, the introduction of a reduced order model, which takes the nonlinear effects into account and accelerates the calculations considerably, is often necessary.This work yields a method that allows for rapid, accurate and parameterisable solutions by means of a reduced model of the original structure. The structure is discretised and its dynamic equilibrium described by a matrix equation. The projection on a reduced basis is introduced to obtain the reduced model. A comprehensive numerical study on several common reduced bases shows that the simple introduction of a constant basis is not sufficient to account for the nonlinear behaviour. Three requirements for an rapid, accurate and parameterisable solution are derived. The solution algorithm has to take into account the nonlinear evolution of the solution, the solution has to be independent of the nonlinear finite element terms and the basis has to be adapted to external parameters.Three approaches are provided, each responding to one requirement. These approaches are assembled to the integrated method. The approaches are the update and augmentation of the basis, the polynomial formulation of the nonlinear terms and the interpolation of the basis. A Newmark-type time-marching algorithm provides the frame of the integrated method. The application of the integrated method on test-cases with geometrically nonlinear finite elements confirms that this method leads to the initial aim of a rapid, accurate and parameterisable transient solution.

Les sources d’énergies renouvelables représentent un gisement intéressant mais l’intermittence de leur production impose une meilleure anticipation des besoins et la mise en place d’un système de stockage d’énergie. Le stockage thermochimique par adsorption dans un système intégrant le couple zéolithe 13X/eau semble être une solution adaptée à un stockage de l’énergie à basse température pour une application aux bâtiments. Notre objectif consiste à reproduire le comportement de ce type de problèmes thermiques non-linéaires. En effet, une simulation précise et rapide du comportement du système sélectionné permettrait une régulation lors de son utilisation. Un modèle bidimensionnel de stockage d’énergie dans un adsorbeur cylindrique a été développé. La résolution numérique de ce modèle, dit d’ordre élevé, implique l’intégration d’un système de quelques centaines à quelques milliers d’équations fortement non-linéaires et couplés. Les coûts de calculs générés pouvant être prohibitifs, l’application d’une méthode de réduction a ainsi été envisagée afin de conserver les caractéristiques, le couplage des transferts de chaleur et de masse ainsi que les non-linéarités de ce modèle tout en limitant le temps de calculs. La projection de Galerkin des équations de ce dernier sur la base, obtenue grâce à une décomposition orthogonale aux valeurs propres, permet de construire un système dynamique d’ordre faible. Sa résolution est moins coûteuse que celle du modèle d’ordre élevé et reproduit correctement la dynamique de l’adsorbeur. / Renewable energy sources will play a key role in meeting future energy demand. One major criticism of those sources stands in their intermittency requiring both a more effective management of demand and efficient storage systems. We focus on thermo-chemical storage by adsorption-desorption mechanism. Eco friendly, economically viable and suitable with solar energy temperature range made the zeolite 13X - water pair ideal for buildings applications. We built an energy storage model in a cylindrical adsorber which contains the mentioned zeolite13X - water pair. Energy storage has been modeled to present coupled heat and mass transfers thanks to a bi-dimensional model. The numerical simulations lead to the time-space evolution of the heating fluid and adsorbent temperatures and pressure. These knowledge models include typically a great amount of coupled differential equations to solve and strong non linearities. The originality of this study is to build a knowledge model of coupled heat and mass transfer in an adsorber and use the Proper Orthogonal Decomposition (POD) and Galerkin projection to build a minimal model of lower dimension without significant loose of accuracy.

Afin d’élaborer des stratégies de contrôle des écoulements en temps réel, il est nécessaire d’avoir recours à des modèles d’ordre réduit (ROMs), car la résolution des équations complètes est trop coûteuse en temps de calcul (des jours, des semaines) et en espace mémoire. Dans cette thèse, les modèles réduits ont été construits avec la méthode POD (Proper Orthogonal Decomposition). Une méthode de projection basée sur la minimisation des résidus, initiée par les travaux de Leblond et al. [134] a été proposée. Dans certaines configurations, la précision des résultats est significativement augmentée, par rapport à une projection de Galerkin classique. Dans un second temps, un algorithme d’optimisation non-linéaire, à direction de descente basée sur la méthode des équations adjointes, a été couplé avec des modèles réduits utilisant des bases POD. Deux méthodes de construction de base POD ont été employées : soit avec un paramètre (un nombre de Reynolds,. . . ), soit avec plusieurs paramètres (plusieurs nombres de Reynolds, . . . ). Les ROMs obtenus ont été utilisés pour contrôler la dispersion d’un polluant dans une cavité ventilée puis pour contrôler le champ de température dans une cavité entraînée différentiellement chauffée. Le contrôle est réalisé en temps quasi-réel et les résultats obtenus sont plutôt satisfaisants. Néanmoins, ces méthodes sont encore trop coûteuses en espace mémoire pour être aujourd’hui embarquées dans les boîtiers de contrôle utilisés dans le bâtiment. Une autre stratégie de contrôle, s’appuyant sur les contrôleurs actuels, a ainsi été développée. Celle-ci permet d’obtenir la température (ainsi que la vitesse) dans la zone d’occupation du bâtiment, en utilisant une décomposition des champs par POD et un algorithme d’optimisation de Levenberg-Marquardt. Elle a été validée sur une cavité différentiellement chauffée, puis appliquée sur une cavité ventilée 3D, proche d’un cas réel. / In order to control flows in real-time, it is necessary to resort to reduced-order models (ROMs) because the classical method of simulations is too expensive in CPU time (several days, weeks) and memory storage. In this thesis, the ROMs have been built with the POD (Proper Orthogonal Decomposition) technique. First, a projection method based on the minimization of the equations residuals and established starting from the works of Leblond et al. [134] have been developed. In some cases, the results accuracy is significantly increased. Secondly, a direct descent optimization algorithm based on adjoint-equations has been coupled with POD/ROMs. Two construction methods of POD bases has been employed: either with simulations for only one parameter (one Reynolds number, . . . ), or with simulations for several parameters (several Reynolds numbers,. . . ). The obtained ROMs have been applied in order to control the pollutant dispersion and then to control the temperature field in a lid-driven cavity heated by the left. The control is realized in quasi-real time and the results are rather satisfying. Nevertheless, these methods are still too expensive in memory storage to be embedded in the current controllers. Thus, another control strategy has been proposed, using POD and an optimization algorithm (Levenberg-Marquardt). This one enables to obtain the temperature (and the velocity) in the occupation zone of the building and has been validated on the lid-driven cavity heated by the left and applied on a 3D-ventilated cavity, similar to a real case.

Ce mémoire de thèse est formé de quatre chapitres. Un premier chapitre présente les différentes notions et outils mathématiques utilisés dans le corps de la thèse ainsi qu’une description des résultats principaux que nous avons obtenus. Le second chapitre présente une généralisation d’un résultat obtenu par Rousselet-Chénais en 1990 qui décrit la sensibilité des sous-espaces propres d’opérateurs compacts auto-adjoints. Rousselet-Chénais se sont limités aux sous-espaces propres de dimension 1 et nous avons étendu leur résultat aux dimensions supérieures. Nous avons appliqué nos résultats à la Décomposition par Projection Orthogonale (POD) dans le cas de variation paramétrique, temporelle ou spatiale (Gappy-POD). Le troisième chapitre traite de l’estimation du flot optique avec des énergies quadratiques ou linéaires à l’infini. On montre des résultats mathématiques de convergence de la méthode de Décomposition Progressive Généralisée (PGD) dans le cas des énergies quadratiques. Notre démonstration est basée sur la décomposition de Brézis-Lieb via la convergence presque-partout de la suite gradient PGD. Une étude numérique détaillée est faite sur différents type d’images : sur les équations de transport de scalaire passif, dont le champ de déplacement est solution des équations de Navier-Stokes. Ces équations présentent un défi pour l’estimation du flot optique à cause du faible gradient dans plusieurs régions de l’image. Nous avons appliqué notre méthode aux séquences d’images IRM pour l’estimation du mouvement des organes abdominaux. La PGD a présenté une supériorité à la fois au niveau du temps de calcul (même en 2D) et au niveau de la représentation correcte des mouvements estimés. La diffusion locale des méthodes classiques (Horn & Schunck, par exemple) ralentit leur convergence contrairement à la PGD qui est une méthode plus globale par nature. Le dernier chapitre traite de l’application de la méthode PGD dans le cas d’équations elliptiques variationnelles dont l’énergie présente tous les défis aux méthodes variationnelles classiques : manque de convexité, manque de coercivité et manque du caractère borné de l’énergie. Nous démontrons des résultats de convergence, pour la topologie faible, des suites PGD (lorsqu’elles sont bien définies) vers deux solutions extrémales sur la variété de Nehari. Plusieurs questions mathématiques concernant la PGD restent ouvertes dans ce chapitre. Ces questions font partie de nos perspectives de recherche. / This thesis is formed of four chapters. The first one presents the mathematical notions and tools used in this thesis and gives a description of the main results obtained within. The second chapter presents our generalization of a result obtained by Rousselet-Chenais in 1990 which describes the sensitivity of eigensubspaces for self-adjoint compact operators. Rousselet-Chenais were limited to sensitivity for specific subspaces of dimension 1, we have extended their result to higher dimensions. We applied our results to the Proper Orthogonal Decomposition (POD) in the case of parametric, temporal and spatial variations (Gappy- POD). The third chapter discusses the optical flow estimate with quadratic or linear energies at infinity. Mathematical results of convergence are shown for the method Progressive Generalized Decomposition (PGD) in the case of quadratic energies. Our proof is based on the decomposition of Brézis-lieb via the convergence almost everywhere of the PGD sequence gradients. A detailed numerical study is made on different types of images : on the passive scalar transport equations, whose displacement fields are solutions of the Navier-Stokes equations. These equations present a challenge for optical flow estimates because of the presence of low gradient regions in the image. We applied our method to the MRI image sequences to estimate the movement of the abdominal organs. PGD presented a superiority in both computing time level (even in 2D) and accuracy representation of the estimated motion. The local diffusion of standard methods (Horn Schunck, for example) limits the convergence rate, in contrast to the PGD which is a more global approach by construction. The last chapter deals with the application of PGD method in the case of variational elliptic equations whose energy present all challenges to classical variational methods : lack of convexity, lack of coercivity and lack of boundedness. We prove convergence results for the weak topology, the PGD sequences converge (when they are well defined) to two extremal solutions on the Nehari manifold. Several mathematical questions about PGD remain open in this chapter. These questions are part of our research perspectives.

Cette thèse présente une étude expérimentale sur un voilier instrumenté, menée pour décrire le comportement aéro-élastique des voiles et du gréement pour des navigations au portant. Les formes des voiles utilisées sont des surfaces non développables avec de fortes courbures provoquant une séparation massive de l’écoulement. De plus, les spinnakers sont des voiles fines et souples rendant l’interaction fluide-structure fortement couplée. A cause du non-respect de certaines règles de similitude, le comportement dynamique d’un spinnaker se prête mal à l’étude en soufflerie et nécessite une comparaison avec des mesures in-situ. Les simulations numériques instationnaires modélisant le comportement aéro-élastique des voiles et du gréement doivent être qualifiées et demandent également des validations. C’est pourquoi un système d’instrumentation embarquée est mis en place sur un J/80, un voilier de huit mètres de long. Il s’agit de mesurer dynamiquement la forme en navigation du spinnaker, les efforts dans les gréements dormant et courant, la répartition de pression sur la voile ainsi que le vent et les attitudes du bateau. La forme du spinnaker en navigation est obtenue grâce à un système de mesure photogrammétrique développé pendant la thèse. La précision de ce système, meilleure que 1,5%, permet de mesurer la forme générale de la voile ainsi que les déformations importantes telles que celles liées au faseyement du guindant. L’effort aérodynamique produit par le spinnaker est obtenu grâce à la mesure de l’intensité des efforts et de leurs directions aux trois extrémités (drisse, amure, écoute) ainsi que par la mesure des pressions sur la voile. Le comportement général du spinnaker est analysé en fonction de l’angle du vent apparent. Une nouvelle représentation utilisant les surfaces de Bézier triangulaires est développée pour décrire la forme tridimensionnelle du spinnaker. Quelques points de contrôles suffisent pour représenter la voile et caractériser le type de voile. Un comportement dynamique propre au spinnaker est également étudié. Le réglage supposé optimal d’un spinnaker est à la limite du faseyement, en laissant le guindant se replier légèrement. Cependant ce réglage n’a jamais été scientifiquement étudié auparavant. Nous avons montré qu’il s’agit d’une forte interaction fluide-structure tridimensionnelle où une importante dépression apparaît au bord d’attaque, qui augmente temporairement les efforts, ce qui n’est pas observé avec un réglage plus bordé. / A full-scale experimental study on an instrumented sailing yacht is conducted to better assess the aero-elastic behaviour of the sails and rigging in downwind navigations. The downwind sail shape is a non-developable surface with high curvature leading to massive flow separation. In addition, spinnakers are thin and flexible sails leading to a strongly coupled Fluid-Structure Interaction. Due to the non-respect of some rules of similitude, the unsteady behaviour of downwind sails cannot be easily investigated with wind tunnel tests that would need comparison with full-scale experiments. Moreover unsteady numerical simulations modelling the aero-elastic behaviour of the sails and rigging require validations. An inboard instrumentation system has been developed on a 8 meter J/80 sailboat to simultaneously and dynamically measure the flying shape of the spinnaker, the aerodynamic loads transmitted to the rigging, the pressure distribution on the sail as well as the boat and wind data. The shape of the spinnaker while sailing is acquired by a photogrammetric system developed during this PhD. The accuracy of this new system, better than 1.5%, is used to measure the global shape and the main dynamic deformations, such as the flapping of the luff. The aerodynamic load produced by the spinnaker is assessed by the measurements of the load magnitudes and directions on the three corners of the sail (head, tack and clew), and also by the pressure distribution on the spinnaker. The global behaviour of the spinnaker is analysed according to the apparent wind angle. A new representation using Bézier triangular surfaces defines the spinnaker 3D shape. A few control points enable to represent the sail and can easily characterise the type of sail. A typical unsteady behaviour of the spinnaker is also analysed. Letting the luff of the sail flap is known by sailors as the optimal trim but has never been scientifically studied before. It is found that it is a complex three dimensional fluid-structure interaction problem where a high suction near the leading edge occurs, producing a temporary increase of the force coefficient that would not be possible otherwise.

This thesis deals with study of dynamics of the swirling flow. The swirling flow occurs frequently in hydraulic machinery (e.g., vortex rope in draft tube of the hydraulic turbine) and often influences operation of these machines. For this reason, sufficient knowledge regarding this characteristic flow is necessary for subsequent improvement of hydraulic machines. The theoretical part of this thesis contains description of flow instabilities and their manifestations, notably Kármán vortex street and vortex rope. In the next part, two methods are applied to these two transient flows in order to identify their specific dynamic properties. The first method, Fourier transform, enables to find frequencies of transient flow. The second method, proper orthogonal decomposition (POD), enables to identify planar or spatial eigenmodes of a specific swirling flow. Proper orthogonal decomposition is used in this thesis to identify planar eigenmodes of Kármán vortex street and spatial eigenmodes of vortex rope.

La machine à réluctance variable peut être utilisée dans les véhicules électriques où pour des considérations d’autonomie, le rendement est crucial. En raison du fort champ de fuite dans la région de l’entrefer de la machine à réluctance variable due à sa géométrie particulière à pôles saillants, les pertes « cuivre » peuvent devenir conséquentes. Il est alors recommandé de ne pas placer les conducteurs au voisinage de l’entrefer. Cependant, des instructions concrètes pour la conception d’un enroulement optimal sont manquantes. Généralement, les pertes « cuivre » dans les machines électriques sont la somme des pertes Ohm DC classiques et des pertes additionnelles dites par courants de Foucault. Les pertes DC étant constantes à un point de fonctionnement donné, l’optimisation est axée alors sur la réduction des pertes par courants de Foucault en jouant sur la configuration géométrique de l’enroulement. Dans le cas de calculs répétitifs fastidieux, rencontrés par exemple lors des processus de conception et d’optimisation du bobinage des machines électriques, il y a un intérêt significatif à réduire le temps de calcul. Dans ce travail, on présente trois techniques de réduction du modèle et leurs applications dans les analyses par la méthode des éléments finis. Outre l’influence de la fréquence d’alimentation et de la section du conducteur, plusieurs facteurs liés à la configuration de l’enroulement influent sur les pertes additionnelles par courants de Foucault :i) la position du conducteur dans l’encoche au voisinage de la dent du stator ou de la zone de l’entrefer .ii) la disposition des conducteurs envers les lignes du champ magnétique bidimensionnelles de l’encoche .iii) l’utilisation d’un conducteur massif ou multi filamentaire; les filaments sont connectés en parallèle et peuvent permuter leurs positions périodiquement au sein du conducteur tout au long du bobinage. Dans cette thèse, on étudie principalement l’influence de la disposition géométrique des spires dans l’encoche et du type du conducteur utilisé s’il s’agit d’un conducteur massif, en fils de Litz ou en fils torsadés. Les pertes par courants de Foucault sont la conséquence d’un couplage fort électrique-magnétique entre la densité du courant et la variation en fonction du temps du champ magnétique. En utilisant le modèle de Maxwell, ce couplage est décrit par une équation différentielle à dérivée partielle qui ne peut être résolue simplement. La résolution de cette équation utilisant l’approche analytique n’est possible que sous certaines hypothèses simplificatrices qui peuvent dégrader la fiabilité de la solution. La modélisation par la méthode des éléments finis permet quant à elle de prendre en compte le mouvement du rotor et la non-linéarité du circuit magnétique garantissant ainsi une meilleure précision. Néanmoins, cela conduit à une large capacité de stockage et à un temps de calcul substantiel qui peut entraver tout processus de conception ou d’optimisation. Pour surmonter ce problème, on propose dans ce manuscrit trois techniques de réduction du modèle. Ces techniques assurent une réduction efficace de la taille du système matriciel associé à la modélisation par la méthode des éléments finis et diminuent par conséquent le temps de calcul : i) une réduction spatiale qui évite une modélisation en 3D des conducteurs complexes en fils torsadés et en fils de Litz et propose une modélisation 2D satisfaisante .ii) la technique de la perturbation. iii) la réduction de l’ordre du modèle utilisant la méthode de la décomposition orthogonale aux valeurs propres combinée à la méthode d’interpolation empirique discrète. La comparaison du modèle réduit à un modèle complet de référence montre l’efficacité de la réduction du modèle à réduire le temps de calcul tout en restant en deçà d’une erreur de précision acceptable. / The switched reluctance machine can be used in hybrid or electric vehicle where, for autonomy considerations, energy efficiency is crucial. Because of the strong stray field in the air-gap region of the switched reluctance machine due to its salient pole geometry, the copper losses can become substantial. It is firmly recommended therefore not to place the coil conductors near the air-gap region. Nevertheless, concrete instructions for optimal winding design are missing. The copper losses in electrical machines are subdivided into classical DC ohmic losses and additional eddy current losses occurring due to the time varying magnetic fields penetrating the copper conductors. Based on the fact that the DC losses are constant at a given operating point, the optimization is focused on reducing the eddy current losses by modifying the winding geometry configuration. In the case of tedious repetitive calculations, met for example during design and optimization processes of electrical machine windings, there is a significant interest in reducing the computation time. This work suggests three model reduction techniques and their applications in the finite element analyses.Besides the frequency of the excitation current and the cross section of the coil conductors, several factors related to the winding configuration can affect the addition al eddy current losses:i) the coil conductor position in the winding slot especially near the stator pole or close to the air gapii) the disposition of the coil conductor against the two-dimensional flux lines in the slot windingiii) the subdivision of the solid conductor into multiple parallel strands swapping their positions periodically in the conductor cross section throughout the length of the machine winding.This thesis mainly studies the influence of the geometric coils disposition in the slot windings and the type of the conductor used whether it is solid or stranded, with Litz or twisted wires.The eddy current losses exit through the strong electro-magnetic coupling between the electric current density and the time dependent magnetic flux lines penetrating the conductors; it is described mathematically by a partial differential equation that cannot be solved easily. The analytical approach, which is used practically for a quick resolution of the strong electro-magnetic coupling equation, is only possible under certain simplifying assumptions that deteriorate brutally the reliability of the copper losses calculation. The finite element modeling as for it, allows taking into account the rotor motion and the non-linear behavior of the magnetic circuit, thus ensuring a higher accuracy. However, it leads under these conditions to a substantial calculation time and requires large storage capacity. These constraints are critical and may hinder therefore any process of conception or optimization. In this thesis, we suggest three different model reduction techniques that can be effective in reducing the size of large scale complete finite element models and enable therefore to shorten the computational time:i) the spatial reduction avoiding the 3D modeling which seems required in the case of twisted and Litz wires and suggesting an alternative satisfactory 2D modeling.ii) the perturbation technique.iii) the model order reduction using the proper orthogonal decomposition combined with the discrete empirical interpolation method.The comparison between the reduced model solutions to that of the complete finite element model has proved the effectiveness of the proposed model reduction techniques; they allow shrinking the required computational time while staying below an acceptable error of accuracy.