Le calcul de la distribution du champ magnétique est basé sur la formulation en potentiel scalaire magnétique dans le contexte des éléments finis nodaux. Une amélioration de la formulation a été apportée en tenant compte de la perméabilité élevée des matériaux magnétiques et du problème de connexité. La formulation en potentiel scalaire a été généralisée par l'introduction du couplage entre équations de champ et de circuits électriques. Un calcul préliminaire d'un potentiel vecteur électrique autorise la prise en compte de tout type de bobines, filaire ou massive. La résolution simultanée des équations de champ et de circuits mise en oeuvre permet traiter les problèmes à courant ou tension imposés et comportant des circuits d'alimentation avec des composants électroniques. Ces travaux sont validés sur plusieurs applications industrielles. Au cours du déroulement de ces travaux, nous avons été amené à développer une méthode de gestion robuste pour des conditions aux limites complexes et un langage dédié de haut niveau pour automatiser la mise en oeuvre de notre démarche.

Les travaux de cette thèse portent sur une nouvelle structure de machine à double excitation (MSDE) pour l'application des véhicules hybrides et électriques. Ce type de machine, ayant deux sources d'excitation, bénéficie un degré de liberté supplémentaire et un contrôle facile sur le flux. Grâce à ce degré de liberté, la machine peut être dimensionnée de manière que son meilleur rendement coïncide avec la zone de fonctionnement la plus sollicitée de la machine. Cette nouvelle structure a fait l'objet principal de ce mémoire. Le fonctionnement de la MSDE est présenté dans les deux premiers chapitres. La machine est dimensionnée suivant un cahier des charges pour véhicule hybride. La validation expérimentale a confirmé le bon fonctionnement de la structure et montré son intérêt. Une autre problématique dans le dimensionnement de la machine est l'aspect thermique car les machines sont devenues de plus en plus compactes et puissantes. Une estimation correcte des pertes est indispensable pour évaluer correctement les performances de la machine. C'est la raison pour laquelle on a décidé de consacrer une partie de cette thèse à la modélisation des pertes fer, dont l'estimation n'est pas évidente.

Le stockage holographique de données est actuellement une voie très prometteuse pour accroître de manière conséquente les capacités et débits des mémoires optiques amovibles. L'architecture de mémoire holographique visée à long terme par ces travaux consiste en une couche de données structurée en une matrice de microfibres perpendiculaires au disque. Les hologrammes sont des réseaux de Bragg multiplexés en longueur d'onde dans chaque microfibre. Le matériau d'enregistrement envisagé est un photopolymère dopé par des centres moléculaires bistables, le photopolymère servant à la structuration en microfibres et les centres à l'enregistrement des réseaux. L'impact de la répartition discrète de ces centres et de la non linéarité de leur réponse optique sur le rapport signal sur bruit de lecture a tout d'abord été évalué par de premières modélisations. Par ailleurs, une procédure d'écriture-lecture en configuration de type Lippmann a été définie afin de simplifier l'architecture de la mémoire et accroître le signal de lecture grâce à une détection homodyne des hologrammes. La mise en œuvre pratique de cette procédure a alors par la suite été étudiée au travers de deux dispositifs expérimentaux d'écriture-lecture. Le premier montage a ainsi abouti à une démonstration de la procédure dans une configuration simplifiée, limitant les difficultés relatives au matériau et à la géométrie d'enregistrement, en bon accord avec notre modèle théorique. L'acquis obtenu avec cette première démonstration a ensuite été utilisé pour mettre en place un second dispositif, afin de valider notre procédure dans une configuration plus réaliste impliquant un photopolymère comme milieu de stockage.

Cette thèse est dédiée à la modélisation et à l'analyse mathématique de modèles asymptotiques utilisés en océanographie décrivant la propagation des ondes internes à l'interface entre deux couches de fluides de densités différentes, soumis à la seule force de gravité.L'objectif de cette thèse est de construire et justifier de nouveaux modèles asymptotiques prenant en compte la variation de la topographie. Pour ce faire, on pose plusieurs hypothèses de petitesse sur la profondeur de l'eau et sur les déformations à l'interface et au fond. On s'intéresse plus particulièrement à deux régimes de variations topographiques, celui de moyenne amplitude et celui de lentes variations de grande amplitude.La première partie de cette thèse consiste à justifier rigoureusement et étudier mathématiquement (existence, unicité, stabilité et convergence de la solution) deux classes de modèles asymptotiques. Une classe de modèles couplés et une classe de modèles scalaires. Cette dernière classe est caractérisée par la description de la propagation unidirectionnelle des ondes internes.Dans la deuxième partie on propose un schéma numérique pour résoudre le modèle asymptotique couplé dérivé dans la première partie dans le cadre d'un font plat. Ce modèle existant dans la littérature a été reformulé d'une façon plus appropriée pour la résolution numérique en gardant le même ordre de précision que l'original et en améliorant ses propriétés de dispersion. Enfin nous présentons plusieurs simulations numériques pour valider notre schéma. / This thesis is dedicated to the modeling and the mathematical analysis of asymptotic models used in oceanography describing the propagation of internal waves at the interface between two layers of fluids of different densities, under the only influence of gravity.We aim here at constructing and justifying new asymptotic models taking into account variable topography. To this end, we assume several smallness assumptions on the depth of the water and on the deformations at the interface and at the bottom. We are interested in two topographic regimes, one for variations of medium amplitude and one for slow variations with large amplitude.In the first part of this thesis we rigorously justify and mathematically study (existence, uniqueness, stability and convergence of the solution) two classes of asymptotic models. A class of coupled models and a class of scalar models. The latter class is characterized by the description of the propagation of unidirectional internal waves. In the second part we propose a numerical resolution for the coupled asymptotic model derived in the first part restricted to the flat bottom case. This existing model in the literature has been rewritten under a new formulation more suitable for numericalresolution with the same order of precision as the standard one but with improved frequency dispersion. Finally, we present several numerical simulations to validate our scheme.

Suite à la découverte du boson de Higgs en Juin 2012 au Large Hadron Collider, l’accélérateur de particules situé à la frontière franco-suisse, l’étude du secteur scalaire des particules élémentaires a connu un regain d’intérêt. En particulier, le boson de Higgs étant une particule clef au sein du Modèle Standard des particules, les expérimentateurs étudient ses propriétés avec beaucoup de soin.Le Modèle Standard, dont le but est de décrire les interactions entre particules élémentaires, n’est cependant pas une théorie complète. En effet, en plus de quelques problèmes d’ordre théorique, certains phénomènes observés expérimentalement ne peuvent pas être expliqués par ce modèle. Les théoriciens en physique des particules cherchent donc à établir une nouvelle théorie venant le compléter et permettant d’expliquer pleinement les observations expérimentales.Cette thèse est axée sur l’étude du secteur scalaire de modèles au-delà du Modèle Standard des particules. J’ai plus particulièrement travaillé sur un modèle à deux doublets de Higgs – modèle purement effectif mais qui peut être inclus dans d’autres théories plus abouties – ainsi que sur un modèle construit comme une combinaison entre les théories déjà très proches de techicouleur et de Higgs composites, et ce dans le cas particulier d’une brisure de symétrie SU (4) ? Sp(4). J’ai étudié ce dernier modèle d’un point de vue effectif mais la théorie complète est capable depallier un certain nombre des limitations du Modèle Standard.Chacun de ces modèles inclut un secteur scalaire plus riche que celui du Modèle Standard et contient au moins une particule pouvant être assimilée au boson de Higgs découvert au LHC. J’ai réalisé l’étude phénoménologique de chacun de ces modèles et les ai confrontés à des contraintes tant théoriques qu’expérimentales – en particulier celles obtenues grâce aux études les plus récentes, portant sur le boson de Higgs et sur de potentielles particules scalaires additionnelles, réalisées par les équipes du LHC. Cela m’a permis de contraindre les paramètres libres des modèles et en particulier de restreindre les valeurs possibles pour la masse des autres particules scalaires, permettant de mieux cibler les zones où ces nouvelles particules, si elles existent, pourraient être détectées au LHC.Ces deux théories, bien que très contraintes par les données expérimentales, ne sont toujours pas exclues par les contraintes expérimentales les plus récentes / Following the discovery of the Higgs boson in June 2012 at the Large Hadron Collider, the particle collider located beneath the France-Switzerland border, interest in the study of the scalar sector in elementary particle physics significantly increased. In particular, as the Higgs boson plays a very special role in the Standard Model of particle physics, experimentalists study its properties with great care.The goal of the Standard Model is to describe the interactions between elementary particles. However the theory is not quite complete. Indeed, in addition to some purely theoretical problems, a number of experimental observations cannot be explained by the Standard Model. Theorists are therefore looking for a more comprehensive theory able to fully explain the observations.This thesis is based on the study of the scalar sector of two different extensions of the Standard Model of particle physics. I have worked on the Two-Higgs Doublet Model – this model is purely effective but can be included in more comprehensive theories – as well as on a model based on a combination of Technicolor and Composite Higgs theories in the framework of the SU (4) ? Sp(4) symmetry breaking pattern. I studied the latter via an effective approach but the full theory is able to get rid of some of the pitfalls of the Standard Model.These two models include a scalar sector that is richer than the one found in the Standard Model and contain at least one particle which can be assimilated to the Higgs boson discovered at the LHC.I performed a phenomenological study for these two models and tested them against both theoretical and experimental constraints. In particular I used the latest studies on the 125 GeV Higgs boson and on possible additional scalars performed by the ATLAS and CMS collaborations. The application of all these constraints drastically reduced the available parameter space of the two models. In particular it narrowed the possible mass range of the additional scalars, allowing to know more accurately where to search them experimentally in order to prove or rule out their possible existence.As of today the two theories I worked on are still not excluded by the latest experimentaldata

L’optimisation de l’efficacité énergétique des machines électriques constitue un domaine de recherche bien développé et fait partie des objectifs de plusieurs accords internationaux comme le projet Energie-Climat de la Commission Européenne visant l’amélioration de 20% d’efficacité pour 2020, encore étendu pour 2030 avec des objectifs plus importants. Ainsi, cette thèse propose un procédé d’optimisation du rendement du moteur asynchrone en agissant sur les paramètres du contrôle. Pour atteindre cet objectif, le flux dans l’entrefer est adapté selon un tableau de flux optimal calculé hors ligne pour tous les points de fonctionnement possibles. Ce flux est déterminé avec le plus haut degré de précision possible en se basant sur un modèle dynamique de la machine développé dans ces travaux. Ce dernier pallie le point faible du modèle dynamique classique, en prenant en compte l’effet des pertes fer. Le modèle des pertes fer utilisé est celui de Bertotti, qui les évalue en fonction de la fréquence et de l’amplitude du champ magnétique. Les pertes sont alors représentées par une résistance variable, continuellement évaluée selon le point de fonctionnement. Le tableau de flux optimal obtenu est fonction des conditions d’opération repérées dans le plan couple-vitesse. Ainsi l’étude montre que le flux peut être optimisé pour des valeurs de couple sensiblement inférieures à environ la moitié du couple nominal, ce seuil variant en fonction de la vitesse. La méthode d’optimisation proposée est simulée puis testée pour le contrôle scalaire et le contrôle vectoriel indirect par orientation de flux rotorique, afin de montrer la généricité de l’approche. La validation est conduite sur une maquette expérimentale d’une puissance de 5.5 kW et pour 2 machines asynchrones de générations différentes (IE2 et IE3). Les résultats obtenus montrent la réduction des pertes dans la machine et donc une amélioration du rendement global, tout en préservant un comportement dynamique satisfaisant. L’optimisation de l’efficacité énergétique est ainsi validée pour les deux structures de contrôle et pour les deux types de machine. Outre une comparaison avec la simulation, la solution proposée est comparée aux méthodes existantes afin d’en apprécier l’efficacité. / Energy efficiency optimization of electric machines is an important research field and is part of the objectives of several international projects such as the European Commission Climate and Energy package which has set itself a 20% energy savings target by 2020, and was extended for 2030 with higher targets. Therefore, this thesis proposes an efficiency optimization method of the Induction Machine (IM) through the variation of the control parameters. To achieve this goal, the flux in the airgap is modified according to an optimal flux table computed off-line for all possible operating points. The flux table is calculated with the best possible accuracy through an improved dynamic model of the IM, developed in these works. The latter avoids the main drawback of the classic dynamic model, by considering the effect of core losses. The core loss model established by Bertotti is used. It depends on the frequency and the amplitude of the magnetic field. The losses are then represented by a variable resistor, continuously evaluated according to the operating point. The established optimal flux table is a function of the operating conditions in terms of torque and speed. Indeed, the results show that the flux can be optimized for torque values less than about half the rated torque, and that this threshold is influenced by the speed. The proposed optimization method is simulated, then tested for the scalar control and the field-oriented control, in order to show the genericity of the proposed approach. The validation is carried on an experimental test bench for two 5.5 kW induction motors of different efficiency standards (IE2 and IE3). The results obtained show the reduction of the losses in the motor, thus an improvement of the overall efficiency while preserving a satisfactory dynamic behavior. Consequently, the optimization of the energy efficiency is validated for the two control structures and for the two studied motors. In addition to the validation of the simulation results, the proposed approach is compared to existing methods to assess its effectiveness

Cette thèse analyse les évolutions liées à l’introduction des réformes néolibérales dans le secteur de l’eau à partir des histoires urbaines de trois quartiers périphériques de l’agglomération du Grand Khartoum. Articulé autour de la transition depuis des systèmes d’approvisionnement et d’accès à l’eau locaux vers les réseaux d’adduction d’eau centraux et publics, ce travail privilégie une approche scalaire et temporelle qui met en évidence la recomposition des relations de pouvoir liées aux infrastructures et au contrôle de la ressource. Cette approche met en évidence les négociations multiples qui se jouent autour de l’eau. Dans un premier temps, l’analyse d’une politique scalaire descendante de l’acteur gouvernemental, la Khartoum State Water Corporation (KSWC) contribue à la mise en place d’organisations locales de l’approvisionnement en eau par une association ou les comités populaires. Cette politique s’appuie et renforce les catégorisations de l’espace et de service discriminantes dans lesquels s’intègrent d’autres acteurs de l’approvisionnement et de l’accès à l’eau (ONG, gestionnaires locaux, population). Dans un second temps, l’intensification des réformes néolibérales entraîne la restructuration du service au niveau de l’agglomération. L’extension des réseaux et de l’administration public-es par la KSWC qui participe à la redistribution des responsabilités techniques et du pouvoir liés à la ressource en eau. Cette politique n’implique pas une désagrégation de l’autorité centrale mais une recomposition des instances étatiques au sein de la ville. Les stratégies et les pratiques scalaires des acteurs locaux et de la population lors des négociations sur le développement des réseaux unitaires témoignent des évolutions de la relation entre l’État et la société urbaine. Enfin, l’étude scalaire des relations autour de l’eau à partir des systèmes locaux met en évidence la mutation et l’intégration des formes d’exercice autoritaire de pouvoir dans la ville. / Starting from the urban history of three peripheral districts of the Greater Khartoum urban area, this thesis analyses the transformations of the water supply due to neoliberal reforms. The transition from local water supply systems to central and public networks involves embedded negotiations around water. The scalar and temporal perspective used in this work underlines important restructurings of power relations linked to infrastructures and control over the resource. First, the politic of down-scaling from the governmental actor, the Khartoum State Water Corporation (KSWC), leads to the development of local water supply systems, either by an association, or by the popular committees. This politics of down-scaling water supply reinforces discriminatory categories of space and service which are integrated by others actors like NGOs, local managers and the population. Secondly, the intensification of neoliberal reforms entails the restructuring of water service at the urban level. KSWC’s networks and administration expansion contributes to the redistribution of technical responsibilities and water powers. This turn does not lead to the disintegration of the central authority, but rather to a restructuration of the state apparatus in the urban space. When networks arrive, the scalar strategies and practices of local actors and population give evidence of changes in the State-Society relationships. Finally, the scalar approach of relations around water underlines the integration of authoritarian forms of power.

Les systèmes cryptographiques à base de courbes elliptiques sont aujourd'hui de plus en plus employés dans les protocoles utilisant la cryptographie à clef publique. Ceci est particulièrement vrai dans le monde de l'embarqué qui est soumis à de fortes contraintes de coût, de ressources et d'efficacité, car la cryptographie à base de courbes elliptiques permet de réduire significativement la taille des clefs utilisées par rapport aux systèmes cryptographiques précédemment employés tels que RSA (Rivest-Shamir-Adleman). Les travaux qui suivent décrivent dans un premier temps l'implantation efficace et sécurisée de la cryptographie à base de courbes elliptiques sur des composants embarqués, en particulier des cartes à puce. La sécurisation de ces implantations nécessite de prendre en compte les attaques physiques dont un composant embarqué peut être la cible. Ces attaques incluent notamment les analyses par canaux auxiliaires qui consistent à observer le comportement d'un composant pendant qu'il manipule une valeur secrète pour en déduire de l'information sur celle-ci, et les analyses par faute dans lesquelles un attaquant peut perturber un composant dans le même but.Dans la seconde partie de ce mémoire de thèse, nous étudions ces attaques et leurs implications concernant l'implantation des systèmes cryptographiques à clef publique les plus répandus. De nouvelles méthodes d'analyse et de nouvelles contre-mesures sont en particulier proposées. Une étude spécifique de certaines attaques appliquées à l'algorithme de chiffrement par bloc AES est également présentée. / Elliptic curve based cryptosystems are nowadays increasingly used in protocols involving public-key cryptography. This is particularly true in the context of embedded devices which is subject to strong cost, resources, and efficiency constraints, since elliptic curve cryptography requires significantly smaller key sizes compared to other commonly used cryptosystems such as RSA.The following study focuses in a first time on secure and efficient implementation of elliptic curve cryptography in embedded devices, especially smart cards. Designing secure implementations requires to take into account physical attacks which can target embedded devices. These attacks include in particular side-channel analysis which may infer information on a secret key manipulated by a component by monitoring how it interacts with its environment, and fault analysis in which an adversary can disturb the normal functioning of a device in the same goal.In the second part of this thesis, we study these attacks and their impact on the implementation of the most used public-key cryptosystems. In particular, we propose new analysis techniques and new countermeasures for these cryptosystems, together with specific attacks on the AES block cipher.

Cette thèse a pour objectif la construction de schémas d’ordre élevé et peu diffusifs pour le transport d’un scalaire dans les écoulements à surface libre, en deux ou trois dimensions. On souhaite en particulier obtenir des schémas robustes, qui gardent au niveau discret les propriétés mathématiques de l’équation de transport avec une faible diffusion numérique, et les utiliser sur des cas industriels. Dans ce travail deux méthodes numériques sont envisagées : une méthode aux volumes finis (VF) et une méthode aux résidus distribués (RD). Dans les deux cas, l’équation de transport est résolue avec une approche découplée, qui est la solution la plus avantageuse en termes de précision et de coûts de calcul. Pour ce qui concerne la méthode aux volumes finis, les équations de Saint-Venant couplées à l’équation du transport sont d’abord résolues avec un schéma dit vertex-centred où le flux numérique est approximé avec un solveur de Riemann appelé Harten-Lax-Van Leer-Contact [135]. A partir de cette approche, une formulation découplée est proposée. Cette dernière permet de résoudre l’équation du transport avec un pas de temps plus grand que celui de la formulation couplée. Cette idée a été d’abord proposée pour d’autres schémas dans [13]. Pour augmenter l’ordre de précision en espace, la technique MUSCL [89] est utilisée en combinaison avec l’approche découplée. Finalement, la problématique des zones sèches est abordée. Dans le cas de la méthode aux résidus distribués, les équations de Saint-Venant sont résolues avec une méthode éléments finis, et la méthode RD est utilisée seulement pour discrétiser l’équation du transport, en focalisant l’attention sur les problèmes non stationnaires. L’équation de continuité du fluide discrétisée est employée pour garantir la conservation de la masse et le principe du maximum. Pour obtenir des schémas d’ordre deux dans les problèmes non stationnaires, un schéma prédicteur-correcteur [112] est utilisé, en l’adaptant au cas de concentration moyennée sur la verticale. Une version d’ordre 1 mais peu diffusive, est aussi présentée dans ce travail. De plus, un schéma localement implicite, complètement nouveau, est aussi formulé pour pouvoir traiter le problème des bancs découvrant. Les deux techniques sont validées d’abord sur des cas simples, pour évaluer l’ordre de précision des schémas et ensuite sur des cas plus complexes pour vérifier aussi les autres propriétés numériques. Les résultats montrent que les nouveaux schémas sont à la fois précis et conservatifs, tout en gardant la monotonie comme le prévoient les démonstrations. Un cas d’application industriel est aussi présenté en conclusion. Le schéma prédicteur-correcteur RD est adapté aussi au cas 3D, sans aucun problème théorique nouveau, par rapport au cas 2D. Les propriétés de base des schémas sont validées sur des cas test préliminaires / The purpose of this thesis is to build higher order and less diffusive schemes for pollutant transport in shallow water flows or 3D free surface flows. We want robust schemes which respect the main mathematical properties of the advection equation with relatively low numerical diffusion and apply them to environmental industrial applications. Two techniques are tested in this work: a classical finite volume method and a residual distribution technique combined with a finite element method. For both methods we propose a decoupled approach since it is the most advantageous in terms of accuracy and CPU time. Concerning the first technique, a vertex-centred finite volume method is used to solve the augmented shallow water system where the numerical flux is computed through an Harten-Lax-Van Leer-Contact Riemannsolver [135]. Starting from this solution, a decoupled approach is formulated and is preferred since it allows to compute with a larger time step the advection of a tracer. This idea was inspired by [13]. The Monotonic Upwind Scheme for Conservation Law [89], combined with the decoupled approach, is then used for the second order extension in space. The wetting and drying problem is also analysed and a possible solution is presented. In the second case, the shallow water system is entirely solved using the finite element technique and the residual distribution method is applied to the solution of the tracer equation, focusing on the case of time-dependent problems. However, for consistency reasons the resolution of the continuity equation must be considered in the numerical discretization of the tracer. In order to get second order schemes for unsteady cases a predictor-corrector scheme [112] is used in this work. A first order but less diffusive version of the predictor-corrector scheme is also introduced. Moreover, we also present a new locally semi-implicit version of the residual distribution method which, in addition to good properties in terms of accuracy and stability, has the advantage to cope with dry zones. The two methods are first validated on academical test cases with analytical solution in order to assess the order of the schemes. Then more complex cases are addressed to test the robustness of the schemes and their performance under different flow conditions. Finally a real test case for which real data are available is carried out. An extension of the predictor-corrector residual distribution schemes to the 3D case is presented as final contribution. Even in this case the RD technique is completely compatible with the finite element framework used for the Navier-Stokes equations, thus its extension to the 3D case does not present any extra theoretical problem. The method is tested on preliminary cases

La pollution atmosphérique et ses conséquences sur la santé et l'environnement constituent un domaine d'étude complexe à cause du nombre de phénomènes physiques mis en jeu. L'objectif de ce travail est d'étudier les principales caractéristiques de l'écoulement et de la dispersion d'un scalaire passif au sein de la canopée urbaine. Pour ce faire, un code numérique a été développé. Il résout les équations de Navier-Stokes dans le cadre d'un écoulement incompressible pour une atmosphère neutre en faisant usage de la méthode de la compressibilité artificielle selon la méthode des volumes finis. Le modèle de Spalart-Allmaras a été utilisé pour la modélisation de la turbulence. La canopée urbaine est explicitement prise en compte et est modélisée par un groupe d'obstacles de forme cubique. Le code fut d'abord testé pour des configurations bidimensionnelles avec un seul et 4 obstacles en configuration alignée selon deux approches : une simulation stationnaire RANS et instationnaire URANS qui reproduit le décrochement tourbillonnaire. La prise en compte du décrochement tourbillonnaire se traduit par une diffusion dans le sillage turbulent du groupe d'obstacles. Les résultats ont été comparés à des mesures expérimentales et d'autres résultats numériques de référence dans la bibliographie et montrent l'amélioration du champ de vitesse moyen par l'approche code fut ensuite testé sur un cas tridimensionnel avec un groupe d'obstacles organisés selon 2 configurations géométriques: alignée et en quinconce. Afin d'éliminer les effets des conditions aux limites, l'écoulement fut calculé sur un volume élémentaire de calcul en utilisant des conditions aux limites périodiques. Deux types de simulations furent réalisés: l'approche RANS classique et la version DDES du modèle de Spalart-Allmaras. L'écoulement obtenu par la DDES améliore de façon significative les résultats par rapport au RANS en comparaison de mesures expérimentales de simulation directe et montrent la bonne potentialité du modèle. La dispersion d'un scalaire passif émis au sein de la canopée fut obtenue sur un domaine plus important comprenant 16 volumes élémentaires par le biais des conditions aux limites périodiques utilisées. Une analyse du champ de concentration a ensuite été réalisée et des comparaisons effectuées en fonction du type de calcul et de la configuration géométrique. / Atmospheric pollution and its impact on health and the environment depend on many physical phenomena, and this makes it a difficult subject to study. The main objective of this work is to investigate the main characteristics of the flow and dispersion of a passive scalar in the urban canopy. Specifically, the urban canopy is simulated by a group of cubical obstacles in a neutrally-buoyant atmospheric boundary layer. A numerical code bas been developed as a tool to aid in this study; flow is computed by solving the Navier-Stokes equations for an incompressible flow, using a finite volume approach, and the method of artificial compressibility. The turbulence is modeled using the method proposed by Spalart and Allmaras. The code was tested first in a 2-D configuration, for flow over a single obstacle, and over a group of 4 obstacles; in both cases two types of simulation were studied: a stationary RANS simulation, and an unsteady RANS (URANS), which reproduced vortex shedding from the obstacles. The explicit inclusion of vortex shedding in the URANS simulation leads to diffusion in the obstacle wakes, and the results compare better with experimental measurements and other published numerical simulations than do those for the RANS simulations. The code was then tested for some 3-D cases consisting of a group of obstacles arrangcd either in aligned or staggered configurations. In order to avoid the influence of boundary conditions, the flow field was simulated using periodic boundary conditions and a small sub-unit from the group of obstacles. Two types of simulation were performed: a classical RANS type calculation and the DDES proposed by Spalart and Allmaras. The results obtained using the ODES agree much more closely with experimental measurements and the results of other numerical simulations than do those obtained using RANS, and indicate the potential of this approach. The dispersion of a passive scalar in the urban canopy was simulated on a much larger domain consisting of 16 of the sub-units used to compute the flow field. The concentration fields were analyzed to show the influence of the geometrical configuration and the type of model.