Cette thèse traite avec des équations aux dérivées partielles provenant de la physique mathématique. En particulier, à partir de modèles 3D ferromagnétisme et ferroélectricité, nous obtenons des modèles 1D et 2D par l'intermédiaire de processus asymptotiques basés sur des méthodes de réduction de dimension. Le modèle 3D ferromagnétisme a été proposé par W.F. Brown depuis lesannées 40. Il est également possible d'utiliser un modèle dynamique, décrivant l'aimantation au cours du temps, en utilisant un système introduit par L.D. Landau et E.M. Lifschitz en 1935. Pour le modèle ferroélectrique, nous nous référons aux papiers de P. Chandra et P.B. Littlewood, W. Zhang et K. Bhattacharya et au livre de T. Mitsui, I. Taksuzaki et E. Nakamura.Ma thèse est constituée de trois parties :Au début, je considère l'énergie micromagnétique avec des coefficients dégénératifs dans un fil mince. Après avoir montrer l'existence de minimiseurs du problème, j'identifie l'énergie limite lorsque la section du fil tend vers zéro.Dans la deuxième partie, j'étudie le comportement asymptotique des solutions dépendant du temps des problèmes micromagnétique dans une multi-structure constituée de la jonction de deux fils minces. En supposant que les volumes des deux fils tendent vers zéro avec la même vitesse. On obtient un problème limite couplé par une condition de jonction. Le problème limite reste non-convexe, mais devient complètement local.Dans le dernier chapitre, à partir d’un modèle variationnel 3D non convexe et non-local pour la polarisation électrique dans un matériau ferroélectrique, et à l'aide d'un processus asymptotique basé sur la réduction de dimension, j'analyse des phénomènes de jonction pour deux films minces ferroélectriques joints orthogonaux. Selon la façon dont la réduction se passe, on obtienttrois modèles différents de dimension 2. On remarque qu’un effet de mémoire du processus de réduction apparaît, ce dernier dépend de la compétition entre les épaisseurs des deux films: Le paramètre de guidage est la limite du rapport des épaisseurs des deux films / This thesis deals with partial differential equations coming from mathematical physics. Particularly, starting from 3D models for ferromagnetism and ferroelectricity, we derive 1D and 2D models via asymptotic processes based on dimensional reduction methods. The 3D model for ferromagnetism was proposed by W.F. Brown in the 40s and it is based on a system introduced by L.D. Landau and E.M. Lifschitz in 1935. About the ferroelectric model, we refer tothe papers of P. Chandra and P.B. Littlewood, W. Zhang and K. Bhattacharya and to the book of T. Mitsui, I. Taksuzaki, and E. Nakamura.This thesis based on three works:At the beginning, we consider micromagnetic energy, with some degenerating coefficients, in a thin wire. After showing the existence of minimizers, we identify the limit energy as the section of the wire vanishes.In the second part, we study the asymptotic behavior of the solutions of a time dependent micromagnetic problem in a multi-structure consisting of two joined thin wires. We assume that the volumes of the two wires vanish with same rate. We obtain two 1D limit problems coupled by a junction condition on the magnetization. The limit problem remains non-convex, but now it becomes completely local.In the last chapter, starting from a non-convex and nonlocal 3D variational model for the electric polarization in a ferroelectric material, and using an asymptotic process based on dimensional reduction, we analyze junction phenomena for two orthogonal joined ferroelectric thin films. We obtain three different 2D-variational models for joined thin films, depending on how the reduction happens. We note that, a memory effect of the reduction process appears, and it depends on the competition of the relative thickness of the two films: The guide parameter is the limit of the ratio between these two small thickness

Les systèmes hyperboliques interviennent dans de nombreuses branches des sciences : théorie cinétique, mécanique des fluides non visqueux, magnéto hydrodynamique, dynamique des gaz non visqueux, trafic routier, flux d’une rivière ou d’un glacier, processus de sédimentation, processus d’échanges chimiques, etc. Et souvent, les systèmes qui régissent ces évènements font intervenir des petits paramètres, dont l’étude asymptotique permet d’envisager des simplifications mathématiques et/ou informatiques notoires. L’existence locale et l’existence globale de solutions, uniformément par rapport à ces paramètres, sont des questions fondamentales. Cette thèse regroupe à la fois des résultats généraux sur l’existence locale uniforme de solutions pour des systèmes hyperboliques quasi-linéaires du premier ordre ; et sur l’existence globale uniforme de solutions autour d’un équilibre constant pour ces mêmes systèmes. Le cas du système d’Euler-Maxwell ne satisfaisant pas les conditions requises pour l’existence uniforme globale, nous le traitons à part. / Hyperbolic systems arise in a large field of sciences : kinetic theory, inviscid reactive flow, magnetohydrodynamics, inviscid gas dynamics, traffic flow, river or glacier flow, sedimentation processes, chemical exchange processes, etc. In these kind of systems, small paramaters often appear, and an asymptotic study may lead to mathematical or computational simplifications. One fundamental problem that we may work on is local and global existence of solutions for these systems, uniformly with respect to these parameters. This Ph.D. thesis includes, on one hand, general results on uniform local existence of solutions for first order quasi-linear hyperbolic systems ; and on the other hand, results on uniform global existence of solutions near constant equilibriums for these same systems. In the case of Euler-Maxwell systems, required conditions are not fulfilled for uniform global existence, then we treat it separately.

Cette thèse est consacrée à l’étude mathématique des écoulements complexes dans des tubes minces. Les difficultés ne sont pas seulement liées à la rhéologie complexe, mais aussi aux conditions au bord sur la pression en entrée et en sortie (qui sont moins habituelles, mais réalistes du point de vue physique). Dans une première partie, des écoulements quasi-newtoniens stationnaires sont étudiés. D’abord, on utilise la petitesse du domaine pour montrer l’existence de la solution. Ensuite, on écrit un développement asymptotique de cette solution et on calcule formellement ses coefficients. Finalement, on justifie rigoureusement la validité de ce développement en démontrant des estimations d’erreur. Dans une deuxième partie, on considère des écoulements de fluides visco-élastiques décrits par la loi d’Oldroyd en régime stationnaire. Le modèle que nous avons choisi contient un terme diffusif en contrainte, dont l’ordre de grandeur est lié à la petitesse du domaine. Similairement à la première partie, un développement asymptotique est complètement justifié du point de vue mathématique. Dans le cas particulier de domaines axisymétriques une solution numérique est cherchée afin de la comparer à la solution obtenue via la technique asymptotique. Dans une dernière partie, on étudie les équations de Navier-Stokes non stationnaires. Un résultat d’existence des solutions fortes pour des données petites est démontré. Malheureusement, la méthode directe ne nous a pas permis pas d’avoir suffisamment de contrôle par rapport à la petitesse du domaine. Pour obtenir le résultat désiré, on utilise l’approche à la Kato, basé sur la théorie de C0 semigroupes. / This thesis is devoted to the mathematical analysis of complex flows in thin pipes. The difficulties stem not only from the complex rheology, but also from the boundary conditions used involving the pressure (which are rather atypical, but realistic from a physical point of view).In the first part, we study stationary, quasi-newtonian flows. The existence of a solution is shown using the smallness of the domain as a key ingredient. Furthermore, an asymptotic expansion of this solution is sought and its coefficients are formally computed. Lastly, the validity of this expansion is rigorously justified by proving error estimates. In the second part, we consider visco-elastic flows represented by Oldroyd’s law in stationary regime. The model which we have chosen contains a diffusive stress term, whose order of magnitude is related to the smallness of the domain. Similarly to the first part, a complete asymptotic expansion in mathematically justified. For the special case of axisymmetric domains a numerical solution is sought in order to compare it against the one obtained via the asymptotic technique. In the last part we study the non stationary Navier-Stokes equations. An existence result of strong solutions for small initial data is proven. Unfortunately, the direct method – based on energy estimates – doesn’t give us an optimal control of the smallness constant with respect to the size of the domain. To obtain the desired result, we employ the method of C 0 semigroups of linear operators.

Dans cette thèse, une nouvelle loi d'endommagement en traction des roches sous chargement exploisif est établie. Basée sur l'approche micro-mécanique et énergétique, la loi d'évolution d'endommagement dynamique est conçue à l'aide de la méthode mathématique d'homogénéisation basée sur le développement asymptotique, et en tenant compte de l'effet inertiel lors de la propagation de fissure. Les simulations numériques sont présentées en vue d'illustrer la capacité du modèle à décrire les comportements connus comme les effets de taille pour la réponse structurelle, la sensibilité au taux de déformation, la transition fragile-ductile et la dispersion de l'onde.<p><p>La loi est implémentée dans le code aux éléments finis LAGAMINE pour étudier la réponse macroscopique du modèle. Plusieurs cas d'applications en dynamique ont été examinés.\\<p><p>- En 1D, les problèmes de localisation ont été étudiés dans une barre sollicitée en traction dynamique par une rampe de chargement. Selon l'amplitude de chargement, trois réponses ont été identifiées: a) comportement purement élastique pour des faibles charges, b) localisation à l'extrémité encastrée de la barre pour des chargements modérés et c) localisation à la tête de la barre pour des chargements élevés. L'influence de la taille microstrurelle sur la localisation a été examinée. Des simulations numériques de l'essai de traction dynamique par écaillage ont été éffectuées. Des essais expérimentaux de traction dynamique par écaillage ont été réalisés sur la bauxite et ont permis de valider le modèle en comparant l'analyse post mortem de l'éprouvette aux résultats des simulations numériques d'écaillage. Ces essais ont aussi permis de déterminer les caractéristiques mécaniques du matériau, la résistance mécanique en traction dynamique, l'instant de rupture et la vitesse de déformation à rupture.\\<p><p>- En 2D, des simulations numériques sont efféctuées pour reproduire le comportement d'une mine sous charge explosive. Le modèle a été capable de reproduire l'endommagement en traction sous charge explosive. L'influence des paramètres du modèle tels que: l'orientation de fissures, la taille de microstructure et la valeur initial d'endommagement sur la distribution de l'endommagement autour du trou de charge a été étudiée. Enfin, un cas d'application sur les mines de bauxite de Guinée a été étudié, incluant un calcul de l'extension de la zone endommagée ainsi qu'une prédiction numérique du niveau de vibration engendrée suite aux tirs à l'explosif. La prédiction du modèle est globalement en accord avec les résultats de la littérature. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

In this thesis, the symmetry structure of gravitational theories at null infinity is studied further, in the case of pure gravity in four dimensions and also in the case of Einstein-Yang-Mills theory in d dimensions with and without a cosmological constant.<p><p>The first part of this thesis is devoted to the presentation of asymptotic methods (symmetries, solution space and surface charges) applied to gravity in the case of the BMS gauge in three and four spacetime dimensions.<p><p>The second part of this thesis contains the original contributions.<p>Firstly, it is shown that the enhancement from Lorentz to Virasoro algebra also occurs for asymptotically flat spacetimes defined in the sense of Newman-Unti. As a first application, the transformation laws of the Newman-Penrose coefficients characterizing solution space of the Newman-Unti approach are worked out, focusing on the inhomogeneous terms that contain the information about central extensions of the theory. These transformations laws make the conformal structure particularly transparent, and constitute the main original result of the thesis.<p>Secondly, asymptotic symmetries of the Einstein-Yang-Mills system with or without cosmological constant are explicitly worked out in a unified manner in $d$ dimensions. In agreement with a recent conjecture, a Virasoro-Kac-Moody type algebra is found not only in three dimensions but also in the four dimensional asymptotically flat case.<p><p>These two parts of the thesis are supplemented by appendices. / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

Ce travail est consacré aux problèmes de testd’hypothèses pour les processus de Poisson nonhomogènes.L’objectif principal de ce travail est l’étude decomportement des différents tests dans le cas desmodèles statistiques singuliers. L’évolution de lasingularité de la fonction d'intensité est comme suit :régulière (l'information de Fisher finie), continue maisnon différentiable (singularité de type “cusp”),discontinue (singularité de type saut) et discontinueavec un saut de taille variable. Dans tous les cas ondécrit analytiquement les tests. Dans le cas d’un saut detaille variable, on présente également les propriétésasymptotiques des estimateurs.En particulier, on décrit les statistiques de tests, le choixdes seuils et le comportement des fonctions depuissance sous les alternatives locales. Le problèmeinitial est toujours le test d’une hypothèse simple contreune alternative unilatérale. La méthode principale est lathéorie de la convergence faible dans l’espace desfonctions discontinues. Cette théorie est appliquée àl’étude des processus de rapport de vraisemblancenormalisé dans les modèles singuliers considérés. Laconvergence faible du rapport de vraisemblance sousl’hypothèse et sous les alternatives vers les processuslimites correspondants nous permet de résoudre lesproblèmes mentionnés précédemment.Les résultats asymptotiques sont illustrés par dessimulations numériques contenant la construction destests, le choix des seuils et les fonctions de puissancessous les alternatives locales. / This work is devoted to the hypotheses testing problems for inhomogeneous Poisson processes.The main object of the work is the study of the behaviour of different tests in the case of singular statistical models. The “evolution of singularity” of the intensity function is the following: regular (finite Fisherinformation), continuous but not differentiable (“cusp”type singularity), discontinuous (jump type singularity)and discontinuous with variable jump size. In all thecases we describe analytically the tests. In the case ofvariable jump size we present as well the asymptoticproperties of the estimators.In particular we describe the test statistics, the choice ofthresholds and the form of the power functions for thelocal alternatives. The initial problem is always the testof a simple hypothesis against a one-sided alternative.The main tool is the weak convergence theory in thespace of discontinuous functions. This theory is appliedto the study of the normalized likelihood ratio processesin the considered singular models. The weakconvergence of the likelihood ratio processes underhypothesis and under alternatives to the correspondinglimit processes allows us to solve the mentioned aboveproblems.The asymptotic results are illustrated by numericalsimulations which contain the construction of the tests,the choice of the thresholds, and the power functions forlocal alternatives.

La programmation concurrente à mémoire partagée est un modèle classique de concurrence qui permet notamment de tirer parti des processeurs multicoeurs aujourd'hui très répandus dans les ordinateurs personnels. Les programmes concurrents sont sujets au problème des interblocages, notoirement difficiles à prévoir et à éliminer, en particulier dans le cas de l'utilisation du mécanisme de synchronisation très populaire que sont les mutex. Dans cette thèse nous avons travaillé à rendre plus aisée la programmation avec des mutex en étudiant des méthodes d'évitement des interblocages. Nous avons d'abord étudié une méthode utilisant une analyse statique par un système de types et d'effets, puis une variante de cette méthode dans un langage à typage dynamique. La seconde méthode est celle que nous avons le plus développée. Elle combine prévention et évitement des interblocages pour fournir une fonction de verrouillage sans interblocages expressive et utilisable. Nous l'avons implémentée sous forme d'une bibliothèque Hop (dialecte de Scheme). Ce faisant, nous avons développé un algorithme sans famine pour l'acquisition simultanée d'un nombre arbitraire de mutex, et identifié le concept d'interblocage asymptotique. Nous avons également été amenés à proposer une optimisation des exceptions (blocs finally). Nos tests de performances semblent indiquer un impact négligeable de l'utilisation de notre bibliothèque sur des applications concurrentes réelles. La majeure partie de notre recherche pourrait être appliquée à d'autres langages de programmation structurée tels que Java. / Shared-memory concurrency is a classic concurrency model which, among other things, makes it possible to take advantage of multicore processors that are now widespread in personal computers. Concurrent programs are prone to deadlocks which are notoriously hard to predict and debug. Programs using mutexes, a very popular synchronization mechanism, are no exception. In this thesis we studied deadlock avoidance methods with the aim of making programming with mutexes easier. We first studied a method that uses a static analysis by means of a type and effect system, then a variation on this method in a dynamically typed language. We developed more the second method. It mixes deadlock prevention and avoidance to provide an easy-to-use and expressive deadlock-free locking function. We implemented it as a Hop (dialect of Scheme) library. This lead us to develop a starvation-free algorithm to simultaneously acquire an arbitrary number of mutexes, and to identify the concept of asymptotic deadlock. While doing so, we also developped an optimization of exceptions(finally blocks). Our performance tests seem to show that using our library has negligible impact on theperformance of real-life applications. Most of our work could be applied to other structured programming languages such as Java.

Dans ce travail de thèse, il s'agit essentiellement de déterminer les propriétés mécaniques et physiques linéaires effectives des composites dans lesquels l'interface entre deux phases n'est pas lisse mais très rugueuse. Une approche efficace pour surmonter les difficultés provenant de la présence de rugosités d'interface consiste d'abord à homogénéiser une zone d'interface rugueuse comme une interphase équivalente par une analyse asymptotique et ensuite à appliquer des schémas micromécaniques pour estimer les propriétés effectives en tenant en compte de la présence de l'interphase équivalente. L'objectif principal de ce travail est de développer cette approche dans un cadre général où la surface autour de laquelle l'interface oscille périodiquement et rapidement peut être courbée et les phénomènes physiques concernés peuvent être couplés. Pour atteindre cet objectif, la conduction thermique est premièrement étudiée comme un prototype des phénomènes de transport non couplés pour élaborer dans un cadre simple les éléments essentiels de notre approche. Cette étude, préliminaire mais très utile au vu de l'importance des phénomènes de transport, montre que des résultats généraux et compacts peuvent s'obtenir quand l'interface est ondulée dans une seule direction et que des méthodes numériques sont en général nécessaires dans le cas où l'interface oscille suivant deux directions. L'approche développée et les résultats obtenus pour la conduction thermique sont étendus d'abord à l'élasticité linéaire et ensuite aux phénomènes physiques linéaires couplés tels que la thermoélectricité et la piézoélectricité. Dans ces cas plus complexes, des résultats généraux sont obtenus pour les composites stratifiés avec les interfaces ondulées dans une seule direction et des méthodes numériques sont élaborées pour les composites dans lesquels les interfaces oscillent suivant deux directions / This work is essentially concerned with determining the effective linear mechanical and physical properties of composites in which the interface between two phases is not smooth but very rough. An efficient approach to overcome the difficulties arising from the presence of interfacial roughness is first to homogenize a rough interface zone as an equivalent interphase by an asymptotic analysis and then to apply micromechanical schemes to estimation of the effective properties while accounting for the equivalent interphase. The present work aims mainly to develop this approach in a general situation where the surface around which an interface oscillates periodically and quickly can be curved and the physical phenomena involved can be coupled. To achieve this goal, thermal conduction is first studied as a prototype of transport phenomena so as to elaborate key elements of our approach in a simple situation. This study,even preliminary but very useful in view of the importance of transport phenomena, shows that general and compact results can be obtained when the interface is corrugated in only one direction and that numerical methods are generally required when an interface is curved along two directions. The approach developed and the results obtained for thermal conduction are extended first to linear elasticity and then to linear coupled physical phenomena such as thermoelectricity and piezoelectricity. In these more complex cases, general results are obtained for composite laminates with interfaces oscillating in only one direction, and numerical methods are elaborated for composites in which the interfaces oscillate in two directions

Le procédé de soudage par friction et malaxage est un procédé récent qui a été développé au sein de l'institut de soudure britannique "The Welding Institute" au début des années 90. Ce procédé, utilisé généralement en aéronautique, est sans apport de matière et permet de souder principalement des alliages d'aluminium difficilement soudables par les procédés classiques de soudage. Il consiste à malaxer le matériau de base à l'aide d'un outil constitué d'un pion et d'un épaulement frottant sur les faces supérieures des tôles à souder. La modélisation de ce procédé est très complexe puisque ce dernier implique des couplages entre des phénomènes mécaniques, thermiques et métallurgiques. Le malaxage dans le procédé de soudage FSW est difficile à simuler à l'aide de la méthode des éléments finis (en lagrangien) puisque la zone proche de l'outil de soudage est le siège de grandes déformations. Donc le remaillage est nécessaire. Cependant, le remaillage est cher et très difficile pour les problèmes tridimensionnels. Par ailleurs, après un remaillage, il est nécessaire d'interpoler les champs (vitesses, contraintes,...) correspondant à la solution courante, ce qui peut introduire des erreurs supplémentaires dans le calcul (on parle de diffusion numérique). Nous proposons dans ce travail des modèles basés sur la méthode sans maillage dite "Smoothed Particle Hydrodynamics SPH" et la méthode des moindres carrés mobiles (Moving Least Square MLS) pour la simulation de ce procédé. Ces modèles sont formulés dans le cadre lagrangien et utilisent la forme forte des équations aux dérivées partielles. Le premier modèle basé sur SPH considère la zone de soudure comme un fluide non newtonien faiblement compressible et dont la viscosité dépend de la température. Ce modèle est proposé pour la simulation numérique du comportement thermomécanique d'un matériau soudé par le procédé FSW. Dans le deuxième modèle, un algorithme itératif implicite de premier ordre a été proposé, pour simuler le malaxage de la matière dans le cas d'un matériau viscoplastique, en utilisant la méthode MLS et la technique de collocation. Le troisième modèle est un algorithme implicite d'ordre élevé basée sur le couplage de la méthode MLS et la Méthode Asymptotique Numérique MAN. Cet algorithme permet de réduire le temps de calcul par rapport à l'algorithme itératif implicite de premier ordre. La validation de ces trois modèles proposés a été faite par le code industriel Fluent / Friction stir welding is a recent process that has been developed by the British Welding Institute TWI "The Welding Institute" since 1990s. This process, generally used in aerospace, does not need additional material and allows mainly joining plates of aluminum alloys which are difficult to weld by the classical welding processes. It consists in mixing the base material using a tool comprising a pin and a shoulder which heats the plates to be welded by friction. The modeling of this process is very complex since it involves the coupling between mechanical, thermal and metallurgical phenomena. The mixing in welding process FSW is difficult to simulate using finite element method in lagrangian framework since the area near the welding tool is submitted to large deformations. So remeshing procedure is often required. However, remeshing can be very expensive and difficult to perform for three-dimensional problems. Moreover, after remeshing step, it is necessary to interpolate the fields (velocities, constraints ...) corresponding to the current solution, which may lead to additional errors in the calculation (called numerical diffusion). We propose in this work models based on meshless methods called "Smoothed Particle Hydrodynamics SPH" and Moving Least Square method for the simulation of this welding process. These models are formulated in lagrangian framework and use the strong form of partial differential equations. The first model based on SPH considers the welding zone as a weakly compressible non-newtonian fluid and whose viscosity depends on the temperature. This model is proposed for the numerical simulation of thermo-mechanical behavior of a welded material by FSW process. The second model is a first order implicit iterative algorithm proposed to simulate material mixing in the case a visco-plastic behavior using the MLS method and the collocation technique. The third model is a high order implicit algorithm based on the coupling of MLS method and Asymptotic Numerical Method (ANM). This algorithm allows reducing the computation time by comparison with the first order implicit iterative algorithm. The validation of these three proposed models was done by the industrial code Fluent

Sous l'hypothèse gaussienne, la relation entre indépendance conditionnelle et parcimonie permet de justifier la construction d'estimateurs de l'inverse de la matrice de covariance -- également appelée matrice de précision -- à partir d'approches régularisées. Cette thèse, motivée à l'origine par la problématique de classification d'images, vise à développer une méthode d'estimation de la matrice de précision en grande dimension, lorsque le nombre $n$ d'observations est petit devant la dimension $p$ du modèle. Notre approche repose essentiellement sur les liens qu'entretiennent la matrice de précision et le modèle de régression linéaire. Elle consiste à estimer la matrice de précision en deux temps. Les éléments non diagonaux sont tout d'abord estimés en considérant $p$ problèmes de minimisation du type racine carrée des moindres carrés pénalisés par la norme $ell_1$.Les éléments diagonaux sont ensuite obtenus à partir du résultat de l'étape précédente, par analyse résiduelle ou maximum de vraisemblance. Nous comparons ces différents estimateurs des termes diagonaux en fonction de leur risque d'estimation. De plus, nous proposons un nouvel estimateur, conçu de sorte à tenir compte de la possible contamination des données par des {em outliers}, grâce à l'ajout d'un terme de régularisation en norme mixte $ell_2/ell_1$. L'analyse non-asymptotique de la convergence de notre estimateur souligne la pertinence de notre méthode / Under the Gaussian assumption, the relationship between conditional independence and sparsity allows to justify the construction of estimators of the inverse of the covariance matrix -- also called precision matrix -- from regularized approaches. This thesis, originally motivated by the problem of image classification, aims at developing a method to estimate the precision matrix in high dimension, that is when the sample size $n$ is small compared to the dimension $p$ of the model. Our approach relies basically on the connection of the precision matrix to the linear regression model. It consists of estimating the precision matrix in two steps. The off-diagonal elements are first estimated by solving $p$ minimization problems of the type $ell_1$-penalized square-root of least-squares. The diagonal entries are then obtained from the result of the previous step, by residual analysis of likelihood maximization. This various estimators of the diagonal entries are compared in terms of estimation risk. Moreover, we propose a new estimator, designed to consider the possible contamination of data by outliers, thanks to the addition of a $ell_2/ell_1$ mixed norm regularization term. The nonasymptotic analysis of the consistency of our estimator points out the relevance of our method