Le point de départ de cette thèse est un problème posé par le groupe électromagnétisme de EADS-CCR : comment résoudre plusieurs systèmes linéaires avec la même matrice mais différents seconds membres ? Pour l'application voulue, les matrices sont complexes, denses et de grande taille. Un problème standard comporte environ quelques millions d'inconnues. Comme de telles matrices ne peuvent être ni calculées, ni stockées dans un processus industriel, l'utilisation d'un produit matrice-vecteur approché est la seule alternative. En l'occurrence, le produit matrice-vecteur est effectué en utilisant la méthode multipôle rapide. Dans ce contexte, le but de cette thèse est d'adapter les méthodes itératives de type Krylov de telle sorte qu'elles traitent efficacement les nombreux seconds membres. Des travaux préliminaires avec un seul second membre ont montré que la méthode GMRES est particulièrement efficace et robuste pour cette application. En conséquence dans cette thèse nous abordons uniquement les variantes de GMRES. Les schémas d'orthogonalisation que nous avons implantés dans GMRES sont des variantes de l'algorithme de Gram-Schmidt. <br /><br />Dans une première partie, nous nous intéressons à l'influence des erreurs d'arrondi dans les algorithmes de Gram-Schmidt. Nos résultats répondent à des questions vieilles de vingt-cinq ans. Nous donnons l'explication théorique de ce qui était communément observé et accepté : <br /><br /> - l'algorithme de Gram-Schmidt modifié génère un ensemble de vecteurs bien conditionné ;<br /> - l'algorithme de Gram-Schmidt itéré deux fois fabrique un ensemble de vecteurs orthonormé.<br /><br />Ces deux propositions reposent sur l'hypothèse que la matrice de départ est "numériquement non singulière" en un sens qui est clairement défini. D'autre part, quand l'algorithme de Gram-Schmidt est itéré avec un critère de réorthogonalisation, nous proposons un nouveau critère. Nous montrons que l'algorithme obtenu est robuste alors que le critère communément utilisé est mis en défaut dans certains cas. Finalement, nous généralisons des résultats standards sur les normes en terme de valeurs singulières pour l'algorithme de Gram-Schmidt modifié. Ceci nous permet de dériver un schéma de réorthogonalisation a posteriori utilisant une matrice de rang faible. Ces résultats ont plusieurs applications directes. Nous en donnons des exemples avec les méthodes de Krylov pour résoudre des problèmes linéaires avec plusieurs seconds membres.<br /><br />Dans la deuxième partie, nous avons implémenté des variantes de la méthode GMRES pour les arithmétiques réelle et complexe, simple et double précisions. Cette implémentation convient pour des ordinateurs classiques, à mémoire partagée ou distribuée. Le code en résultant satisfait aux critères de qualité des librairies standards et son implémentation est largement détaillée. Pour des besoins de simplicité, flexibilité et efficacité, les solveurs utilisent un mécanisme de reverse communication pour les produits matrice-vecteur, les étapes de préconditionnement et les produits scalaires. Différents schémas d'orthogonalisation sont implémentés pour réduire le coût de calcul des produits scalaires, un point particulièrement important pour l'efficacité des méthodes de Krylov dans un environnement parallèle distribué. Le critère d'arrêt implémenté est basé sur l'erreur inverse normalisée. Les variantes disponibles sont GMRES-DR, seed-GMRES et block-GMRES. Ces codes s'ajoutent aux variantes déjà existantes (GMRES, flexible GMRES et SQMR). Un produit matrice-vecteur avec une décomposition LU est utilisé dans GMRES-DR de telle sorte que le stockage des approximations des vecteurs propres se fasse sur les premiers vecteurs de l'espace de Krylov. Un restart implicite et une étape de préconditionnement implicite ont été implémentés dans seed-GMRES. Nous supprimons ainsi un produit matrice-vecteur et une étape de préconditionnement par second membre et par cycle de GMRES. La version de block-GMRES permet à l'utilisateur de sélectionner différents modes de déflation. Pour terminer, des résultats reliant la norme du résidu de GMRES à la plus petite valeur singulière de l'espace construit par la méthode de Krylov ont été généralisés à la méthode block-GMRES.<br /><br />La troisième partie est consacrée à l'amélioration des techniques standards pour la résolution des systèmes linéaires dans le cadre des problèmes électromagnétiques. Après une présentation approfondie du code, nous étudions l'influence de la non-symétrie sur la convergence de l'algorithme SQMR. Nous étudions aussi le comportement de GMRES-DR sur nos problèmes. Ceci correspond à deux méthodes avec un seul second membre, le reste de cette partie concerne les cas comportant plusieurs seconds membres. Tout d'abord, nous examinons en détail les techniques qui permettent d'adapter les méthodes utilisées pour un second membre unique aux cas comportant plusieurs seconds membres. Par exemple, on peut améliorer la qualité du préconditionneur, avoir une stratégie de solution initiale, grouper les opérations de plusieurs résolutions ou encore paralléliser plusieurs résolutions. Dans le contexte du calcul de surface équivalente radar monostatique, nous avons montré que l'espace des seconds membres du problème continu était de dimension finie. La dimension donnée par notre théorie est proche de celle que nous observons en pratique. Cette propriété nous permet de réduire considérablement le nombre de systèmes linéaires à résoudre. Dans ce contexte, une version de la méthode block-GMRES est donnée. Ensuite, nous abordons certains problèmes spécifiques des méthodes seed-GMRES et block-GMRES pour lesquels nous proposons des solutions. Pour finir, des résultats plus prospectifs sont donnés. Plusieurs stratégies pour extraire et ajouter de l'information spectrale d'un cycle de GMRES à l'autre sont proposées et comparées. Puis nous utilisons le fait que la méthode multipôle rapide est un produit matrice-vecteur inexact dont la précision est réglable. Moins précis est le produit matrice-vecteur, plus rapide il est. Nous montrons comment tirer partie de cette propriété en utilisant un schéma relâché (méthode de Krylov inexacte) ou des itérations emboîtées (flexible GMRES). Enfin, le critère d'arrêt basé sur l'erreur inverse normalisée dans le cadre du calcul d'une surface équivalente radar est remis en question.

La formation du continent sud américain résulte d'une histoire complexe étalée sur plus de 3.5 Ga. Les processus qui ont façonné le continent ont impliqué l'intégralité de la lithosphère. La structure du manteau supérieur a été étudiée à l'échelle du continent en réalisant un modèle tomographique en ondes de surface, tandis qu'une étude du déphasage des ondes de cisaillement a permis d'appréhender la structure du manteau supérieur sous le sud-est du Brésil.<br /> La tomographie sismique anisotrope en ondes de surface donne une image de la structure en vitesse du manteau supérieur sous le continent et les océans environnants. Une bonne corrélation existe entre les grandes structures géologiques et les hétérogénéités de vitesse. L'anisotropie des ondes de Rayleigh suggère une absence de déformation à grande échelle au-delà de 200 km de profondeur, ce qui contraint verticalement la source de l'anisotropie mise en évidence au sud-est du Brésil. L'étude du déphasage des ondes de cisaillement y a révélé une orientation dominante du plan de polarisation de l'onde quasi-S rapide, parallèle aux chaînes péricratoniques. Un déphasage entre les ondes S rapide et lente supérieur à 2s a été mesuré à l'aplomb des décrochements majeurs. Une telle amplitude suggère que : 1) les décrochements traversent la lithosphère, 2) la lithosphère est caractérisée par une très forte anisotropie intrinsèque, 3) l'asthénosphère contribue au déphasage, par la présence de deux couches ou d'une fabrique tectonique cohérente entre lithosphère et asthénosphère, impliquant une absence de découplage depuis le Néoprotérozoïque.<br /> Une modélisation numérique tridimensionnelle de la déformation lithosphérique a permis d'étudier le couplage mécanique entre la croûte et le manteau supérieur, ainsi que la localisation de la déformation en termes de développement de zones de cisaillement.

Cette these porte sur l'utilisation du formalisme bayesien pour la repr´esentation et la manipulation des incertitudes geometriques dans les systemes de Robotique et de CAORobotique. Dans ces systemes, l'utilisation d'un modele geometrique de l'environnement est indispensable. Toutefois, la validite des calculs conduits sur ces mod`eles n´ecessite une repr´esentation des ecarts entre le modele et la realite et une prise en compte de ces ecarts lors de la resolution d'un probleme donne. L'approche proposee repr´esente une extension de la notion de specification par contraintes geometriques dans laquelle la dimension incertaine des modeles est prise en compte. Cette extension consiste a specifier les contraintes sur les positions relatives entre diff´erents corps de l'environnement non pas par de simples equations et inequations, mais par des distributions de probabilite sur les parametres de ces positions. A l'issue de cette specification, une distribution conjointe sur l'ensemble des parametres du modele est construite. Pour un probleme donne, la distribution marginale sur les parametres inconnus de ce dernier est inferee en utilisant les regles des probabilites. La resolution de ce probleme revient a optimiser cette distribution comportant, dans le cas general, une integrale portant sur un espace de grande dimension. La methode de resolution utilisee pour approcher ce double probleme d'integration/optimisation est basee sur un algorithme genetique. Cet algorithme permet en particulier de controler la precision de l'estimation numerique des integrales par une m´ethode stochastique de Monte-Carlo. L'implantation d'un systeme prototype de CAO nous a permis une experimentation assez poussee de l'approche propos´ee. La mise en oeuvre de plusieurs applications robotiques, dont les natures peuvent paraıtre tres differentes, a ete possible grace a la souplesse de la methode de specification utilisee et la robustesse de la methode de resolution implantee.

La reconstruction 3D est un outil primordial en volcanologie : études morpho-structurales, déformation du sol, modélisation analogique. Une approche originale de reconstruction 3D spécifique aux applications volcanologiques a été développée. Basée sur l'utilisation de plusieurs images numériques acquises autour d'une surface, elle se décompose en deux étapes : (1) un calibrage précis du capteur optique (2) une reconstruction 3D reposant sur la déformation itérative d'un modèle 3D initial. La méthode a été développée et testée en laboratoire sur divers modèles analogiques. Les résultats montrent qu'elle constitue un outil de quantification souple et puissant permettant de détecter des variations fines de la morphologie. Des tests sur site réel ont aussi montré le potentiel important de la méthode dans le cadre de la surveillance volcanologique. L'approche, à la frontière entre Volcanologie et Vision par Ordinateur, offre une alternative aux méthodes classiquement utilisées.

Ce travail de thèse porte sur l'analyse et la synthèse de sons impulsifs. Une des particularités de cette étude réside dans son approche pluridisciplinaire. Elle fait notamment appel à la mécanique vibratoire, aux méthodes de traitement de signaux non stationnaires et à la psychoacoustique. Dans le but d'établir des méthodes de synthèse efficaces, il est de plus en plus courant d'allier plusieurs approches qui peuvent être complémentaires. D'un point de vue physique, un son impulsif est le résultat de la mise en vibration d'un système résonant en régime d'oscillations libres. D'un point de vue perceptif, de nombreux tests psychoacoustiques montrent que l'oreille est capable d'identifier un son impulsif à partir de son amortissement (caractéristique du matériau) ainsi que de son contenu spectral qui définit la hauteur tonale (caractéristique de la géométrie de la structure). Cette étude a débouché sur la mise en oeuvre de modèles de synthèse permettant de simuler le comportement vibratoire de structures simples unidimensionnelles (du type poutre) et multidimensionnelles (du type membrane ou plaque mince) qui sont basés, lorsque cela est possible, sur la modélisation physique de la structure. Les paramètres des modèles sont définis par l'analyse de sons naturels et des grandeurs physiques caractéristiques. Par ailleurs, un modèle de synthèse a pu être mis en uvre dans le cas de structures plus complexes où les solutions analytiques ne sont pas connues. Dans ce cas, les paramètres du modèle sont directement estimés par l'analyse d'une représentation conjointe du son, du type temps-échelle, en prenant en compte des grandeurs perceptives caractéristiques. Enfin, une application concernant la synthèse du son de piano est décrite. Un des modèles de synthèse a été utilisé pour reproduire la réponse impulsionnelle de la table d'harmonie du piano, qui est une structure multidimensionnelle complexe. Les paramètres du modèle ont été estimés à partir de l'analyse de signaux expérimentaux obtenus par des mesures effectuées sur un piano placé en chambre anéchoïque. Dans le cadre de la réalité virtuelle sonore, des fonctions de transfert de la table d'harmonie et de son rayonnement acoustique ont également été estimées afin de reproduire virtuellement le champ acoustique au niveau de la tête d'un pianiste.

Dans le cadre du dimensionnement à la fatigue thermomécanique des culasses automobiles, nous proposons une méthode de simulation des structures soumises à des chargements cycliques. Cette méthode, alternative aux méthodes classiques (comme la méthode incrémentale) permet de résoudre des problèmes thermomécaniques cycliques non-linéaires. Il s'agit d'une Méthode Cyclique Directe qui consiste à rechercher directement la solution asymptotique d'une structure anélastique soumise à un chargement thermomécanique périodique, sans suivre l'histoire du chargement. Elle est fondée sur les techniques de grand incrément de temps [Ladevèze et al.], en prenant en compte la périodicité de l'état limite [Nguyen et al.] et la transformation de Fourier. Cette méthode permet des gains de temps importants pour les structures comprenant un grand nombre de degrés de liberté. Un développement de cette méthode dans une version commerciale d'Abaqus, a permis de montrer la faisabilité de cette approche dans le cas de la simulation d'une culasse.

La thèse de Monsieur Jeunechamps est intitulée "Simulation numérique, à l'aide d'algorithmes thermomécaniques implicites, de matériaux endommageables pouvant subir de grandes vitesses de déformation. Application aux structures aéronautiques soumises à impact". Elle comporte neuf chapitres et deux annexes. La bibliographie compte 285 références. Les développements informatiques ont été implémentés dans le code de calcul par éléments finis Metafor, développé au sein du département LTAS-MC&T et MN²L de l'Université de Liège. Le travail est divisé en trois parties principales. La première partie (chapitres 2 à 4) concerne la description et la modélisation thermomécanique des phénomènes à dynamique rapide sans dégradation irréversible des propriétés du matériau utilisé. La deuxième partie (chapitres 5 à 7) est consacrée à l'étude du comportement des matériaux dits endommageables éventuellement soumis à rupture, c'est-à-dire des matériaux dont les propriétés se dégradent de façon irréversible au cours de la déformation. La troisième partie (chapitre 8) est une application à l'échelle industrielle des méthodes proposées tout au long de cet ouvrage. Le chapitre 2 propose un inventaire des lois constitutives des matériaux, permettant de décrire le comportement de la structure lors de sollicitations rapides. L'accent est mis sur les principales lois d'évolution de la limite élastique implémentées dans les codes de calcul commerciaux, ainsi que sur les variantes de ces lois d'évolution. L'aspect numérique de l'intégration thermomécanique de ces lois est également abordé. Le deuxième aspect abordé dans cette première partie concerne les algorithmes d'intégration temporelle des équations de conservation de la quantité de mouvement. Le chapitre 3 décrit les algorithmes d'intégration utilisés dans ce travail, en mettant l'accent sur l'aspect dynamique et le couplage thermomécanique des phénomènes à dynamique rapide. Le chapitre 4 présente quelques applications illustrant les méthodes de calcul utilisées et permettant la validation de l'implémentation des modèles programmés dans Metafor. Les comparaisons sont effectuées dans la mesure du possible par rapport à des données expérimentales quand celles-ci sont disponibles et également par rapport à des résultats issus de codes de calcul commerciaux. Le chapitre 5 présente une modélisation de la dégradation du matériau au cours de la déformation. En effet, le matériau, sous l'effet des sollicitations et des efforts résultants, perd de ses propriétés de résistance à l'effort, et ce, de manière irréversible. Il est alors endommagé. Dans ce travail, nous avons choisi d'utiliser la théorie de l'endommagement continu pour décrire ces phénomènes. Le chapitre 5 rappelle les fondements de cette théorie ainsi que les principales lois d'endommagement continu. Une méthode générale et originale d'intégration de ces modèles d'endommagement est également proposée. Une fois que la structure est soumise à de trop fortes sollicitations, la rupture consécutive à l'endommagement du matériau apparaît. Le chapitre 6 décrit la méthode numérique développée pour modéliser le déchirement d'une structure ainsi que les différents critères de rupture utilisés. Encore une fois, nous nous limitons à une approche phénoménologique et pragmatique : il ne s'agit pas ici d'implémenter des critères complexes multi-échelles. Cependant, la structure du code de calcul est conçue pour permettre aisément de telles extensions. Le chapitre 7 présente une série d'applications permettant de valider l'implémentation des lois d'endommagement et de rupture ainsi que la formulation proposée de la théorie d'endommagement. Nous étudierons également le coût CPU engendré par la modélisation de l'endommagement et de la rupture du matériau. Enfin, le chapitre 8 présente une application industrielle proposée par la société Techspace Aero S.A. Il s'agit de l'étude du flambement d'une aube de compresseur basse pression d'un moteur d'avion lors du contact accidentel de celle-ci avec le carter du moteur. Tous les développements présentés dans les chapitres précédents sont alors utilisés pour simuler au mieux le phénomène. Les apports principaux du travail sont les suivants : utilisation d'algorithmes thermomécaniques de type étagé, dont l'intégration temporelle des équations de conservation du mouvement prend en compte les effets d'inertie ; amélioration de la technique d'intégration des lois constitutives avec endommagement, utilisant la théorie de l'endommagement continu, par l'utilisation d'un algorithme itératif ; méthode unifiée de calcul de la matrice de raideur tangente matérielle analytique pour un matériau hypoélastique avec endommagement, selon la théorie de l'endommagement continu, dans le cadre thermomécanique ; couplage thermomécanique général des modèles d'endommagement et des lois constitutives à grandes vitesses de déformation ; utilisation d'algorithmes implicites thermomécaniques dans la modélisation de la déchirure de structure par la méthode d'érosion ; établissement d'une plate-forme numérique d'accueil permettant l'implémentation future de lois matérielles, avec ou sans endommagement, ainsi que de techniques de modélisation de propagation de fissure ; processus complet de description d'un phénomène d'impact, par la modélisation du comportement thermomécanique du matériau par des lois de comportement avec endommagement adaptées au phénomène étudié et l'utilisation d'algorithmes implicites couplés à une méthode de déchirure de la structure.

Cette thèse est consacrée à certains modèles mathématiques décrivant le mouvement d'une population de cellules, qui interagissent via un signal chimique. L'accent est mis sur le modèle parabolique de Patlak-Keller-Segel, et dans une moindre mesure, sur le modèle cinétique d'Othmer-Dunbar-Alt.<br /><br />Dans une première partie nous étudions plusieurs variantes du modèle PKS classique, incluant notamment une diffusion non-linéaire des cellules, ou bien une loi de diffusion chimique à noyau de Green logarithmique. Puis nous montrons l'existence globale pour une masse sous-critique du modèle PKS classique dans tout l'espace $\mathbb{R}^2$.<br />On complexifie ensuite le modèle de base en ajoutant un intermédiaire chimique réactionnel, ce qui modifie l'homogénéité du système. Enfin les conditions d'existence globale pour le modèle cinétique ODA avec effets délocalisants sont affaiblies par rapport aux travaux précédents.<br /><br />Dans une deuxième partie nous appliquons le modèle phénoménologique de PKS, et son principe de masse critique, à un processus d'auto-organisation remarquable dans le cerveau: la sclérose concentrique de Baló. Un couplage adéquat entre un front de propagation et une instabilité de PKS décrit raisonnablement les motifs en anneaux de la maladie.<br /><br />La troisième partie adopte le point de vue du transport optimal de masse pour analyser le modèle de PKS unidimensionnel modifié auparavant (afin de partager les caractéristiques de PKS 2D). Bien que la fonctionnelle d'énergie ne soit pas convexe par déplacement, nous démontrons la convergence vers un unique état d'équilibre, lorsqu'il existe. Ces nouvelles idées sont mises en oeuvre numériquement~: un flot gradient discret pour la distance de Wasserstein est analysé, puis simulé en dimension un d'espace.<br /><br />Plusieurs annexes viennent compléter ce travail, dont une annexe qui regroupe tous les aspects numériques de la thèse.

L'énergie éolienne est en train de se développer en conséquence des problèmes environnementaux posés par les sources d'énergie traditionnelles et du progrès technologique des aérogénérateurs. La part de ce type d'énergie dans le réseau est de plus en plus importante dans les régions ventées. Par suite, l'impact sur le réseau électrique, de la qualité de la puissance produite par les aérogénérateurs augmente. D'autre part, le coût de l'énergie éolienne est encore trop élevé pour concurrencer les sources traditionnelles sur des sites moins ventés. Plus " d'intelligence " peut être introduite dans la commande des éoliennes pour résoudre ces problèmes. Cette commande doit être faite en tenant compte du comportement de l'aérogénérateur dans son ensemble. Les perturbations provenant du vent doivent également être considérées. L'objectif de cette thèse se résume ainsi: concevoir les lois de commande d'un aérogénérateur à vitesse variable et régulation pitch, en tenant compte des perturbations éoliennes et avec quatre critères d'optimisation. Ces critères sont: l'amélioration de la qualité de l'énergie électrique produite, du rendement énergétique et de la robustesse des commandes et la diminution des efforts subis par le train de puissance de l'éolienne. Des modèles de simulation permettant d'étudier les lois de commande d'aérogénérateurs existent, mais l'impact des perturbations éoliennes sur le système dans son ensemble n'est en général pas considéré avec suffisamment de précision. Ici, un modèle de simulation d'un aérogénérateur à vitesse variable et régulation pitch est présenté. Une attention particulière est portée à l'interaction entre le vent et la turbine. Ce modèle de simulation a été validé à l'aide de données réelles d'une éolienne de 180 kW et de données provenant de la littérature. Un aérogénérateur à vitesse variable et régulation pitch a différentes zones de fonctionnement. Trois zones sont étudiées ici: celle correspondante à des vitesses de vent moyens où la vitesse de rotation de l'éolienne est variable, celle liée à des vitesses plus importantes où l'éolienne tourne à la vitesse nominale et celle liée aux vitesses de vent où la puissance captée par la turbine est limitée à l'aide du pitch. En plus de la stratégie de commande générale de l'éolienne, différentes conceptions de commandes dans chacune de ces zones sont examinées. Des modèles linéaires permettant d'étudier le comportement dynamique du système éolien sont obtenus sur différents points de fonctionnement à partir du modèle de simulation non-linéaire. Certains de ces modèles sont utilisés en tant que modèle de commande. Des régulateurs numériques robustes sont synthétisés dans chacune des zones de fonctionnement. Ces lois de commande sont ensuite testées sur le modèle de simulation. Différentes commandes sont comparées en observant leur comportement face aux quatre critères liés à l'objectif de thèse. Un modèle de simulation est pratique pour tester les lois de commande, mais il convient d'essayer ces lois aussi sur un banc d'essai avant de les implémenter sur une éolienne réelle. Ici, l'aérogénérateur simulé est émulé sur un banc d'essai de 15 kW. Ensuite, quelques essais expérimentaux sont réalisés sur ce banc. Les résultats obtenus en simulations et sur le banc d'essai montrent que des lois de commande avancées permettent d'améliorer le comportement d'un aérogénérateur à vitesse variable et régulation pitch.

Ces travaux de these s'integrent dans les objectifs du projet international d'Analyse de la Mousson Africaine (AMMA). Ils abordent la caracterisation des echanges entre l'ocean et l'atmosphere dans l'Atlantique Equatorial Est (AEE) de l'echelle diurne a interannuelle a partir de donnees d'observation et de la modelisation numerique. Ces echanges sont d'une tr'es grande importance pour le declenchement de la Mousson Africaine qui determine en grande partie la saison des pluies en Afrique de l'Ouest. Nous montrons que les parametres de la couche de melange oceanique dans le GG sont principalement pilotes, a l'echelle diurne, par les flux de chaleur et les processus de subsurface (l'entraınement, le melange vertical) et nous montrons aussi qu'a partir d'une parametrisation de la turbulence oceanique (TKE) fermee a l'ordre 1.5, on parvient a restituer la turbulence oceanique observee dans cette region. Afin de quantifier le role de la tension du vent sur la variabilite des parametres de la couche de melange oceanique, nous utilisons un modele lineaire force par les anomalies de la tension du vent. Les resultats montrent que pres de l'equateur (3°S-3°N), le signal est retrouve ; par contre de part et d'autre de la bande 3°S-3°N, les amplitudes des parametres de la couche de melange sont tres largement sous estimes, ce qui montre que la tension du vent n'est pas le seul parametre a moduler la variabilite de la temperature surface 'a ces latitudes. Les bilans de chaleur integres sur la couche de melange ont permis de montrer le role des flux de chaleur et du melange vertical dans l'AEE. En effet, l'AEE est subdivise en plusieurs boıtes en tenant compte de la dynamique et de la thermodynamique et dans chaque boıte le bilan y est calcule. Les resultats montrent que les flux de chaleur et le m'elange vertical dominent le bilan aux echelles saisonniere et interannuelle. Ce melange vertical, calcule par residu, est compare avec les donnees independantes de turbulence oceanique (Dengler et al., 2010 ; Rhein et al., 2010) collectees lors des campagnes EGEE/AMMA en 2005-2007. Les resultats montrent que le residu se compare tres bien avec ces donnees independantes en terme de variabilite spatiale et temporelle. Le melange vertical est tres fort dans la region sauf dans les boıtes au sud de l'AEE et sa variabilite est largement pilotee par les flux de flottabilite et la tension du vent. En regard de ces resultats avec ceux obtenus sur le cycle diurne, on peut dire que dans l'AEE, la variabilite spatio-temporelle des parametres de la couche de melange est principalement pilotee a toutes les echelles par les flux de chaleur et les processus de subsurface.