Dans cette thèse, on propose des méthodes analytiques ou numériques pour l'estimation de la VaR ou l'Expected Shortfall des portefeuilles linéaires, quadratiques, lorsque le vecteur des facteurs de risques suit un mélange convexe de distributions elliptiques. Aussi, on introduit pour la prémière fois la notion de "portefeuille quadratique" d'actifs de bases (ie. actions).

Nous nous intéressons à l'imagerie d'un objet enfoui dans un milieu multi-couches inhomogène, avec des données ne contenant pas la phase. Nous résolvons un problème direct modèle de propagation des ondes dans un tel milieu, à l'aide de l'analyse asymptotique et des équations intégrales. Puis nous développons des algorithmes de reconstruction à base de dérivée topologique et des techniques de l'optimisation de forme.

Cette thèse presente une étude numérique et analytique de la stabilité d'un écoulement incompressible de type jet tournant. L'entraînement du fluide externe par le jet est modélisé numériquement par l'hypothèse de frontières latérales et de sortie ouvertes, les conditions d'entrée correspondant à un profile de Grabowski. L'effet d'une faible viscosité dans le voisinage du nombre de swirl critique est étudié via une analyse asymptotique couplée à des simulations numériques axisymétriques. Un algorithme de continuation basé sur une méthode de projection récursive (RPM) a été implémenté pour capturer les états stationnaires et suivre ces branches de solutions dans l'espace de paramètres ainsi que leur stabilité. La continuation des solutions stationnaires vis-à-vis du paramètre de swirl montre l'existence d'une bifurcation pour les nombres de Reynolds assez grands. L'analyse asymptotique confirme ces résultats numériques. Le diagramme de bifurcation d'un jet tournant possédant une région de recirculation est déterminé dans le cas axisymétrique. Il est montré que l'état stationnaire subit une bifurcation de Hopf supercritique. Enfin, la stabilité globale tridimensionnelle d'un jet tournant avec une région de recirculation est étudiée numériquement par une méthode d'Arnoldi. L'état éclaté axisymétrique apparaît instable vis-à-vis de perturbations tridimensionnelles hélicoïdales. L'effet d'un gradient d'une pression extérieur sur le diagramme de bifurcation est étudié numériquement. Pour un nombre de Reynolds Re=1000, la branche colonnaire (solutions sans recirculation) existe dans le cas d'un gradient de pression favorable pour les grnads paramètres de swirl, mais disparaît quand le gradient de pression est zero. Ce résultat ouvre des perspectives pour une stratégie de contrôle pour retarder l'apparition de l'éclatement tourbillonnaire.

Le sujet de la thèse porte sur l'étude d'une nouvelle classe de modèles de transport quantique: les modèles fluides quantiques issus du principe de minimisation d'entropie. Ces modèles ont été dérivés dans deux articles publiés en 2003 et 2005 par Degond, Méhats et Ringhofer dans Journal of Statistical Physics en adaptant au cadre de la théorie quantique la méthode des moments développée par Levermore dans le cadre classique. Cette méthode consiste à prendre les moments de l'équation de Liouville quantique et à fermer ce système par un équilibre local (ou Maxwellienne quantique) défini comme minimiseur d'une certaine entropie quantique sous contrainte de conservation de certaines quantités physiques comme la masse, le courant, et l'énergie. Le principal intérêt des modèles quantiques ainsi obtenus provient du fait qu'étant macroscopiques, ils sont biens moins coûteux numériquement que des modèles microscopiques comme l'équation de Schrödinger ou l'équation de Wigner, et de plus, ils prennent en compte implicitement des effets de collision bien plus difficiles à modéliser à un niveau microscopique. Le but de cette thèse est donc de proposer des méthodes numériques pour implémenter ces modèles et de les tester sur des dispositifs physiques adéquats.<br />Nous avons donc commencé dans le chapitre I par proposer une discrétisation du plus simple de ces modèles qu'est le modèle de Dérive-Diffusion Quantique sur un domaine fermé. Puis nous avons décidé dans le chapitre II et III d'appliquer ce modèle au transport d'électrons dans les semiconducteurs en choisissant comme dispositif ouvert la diode à effet tunnel résonnant. Ensuite nous nous sommes intéressés au chapitre IV à l'étude et l'implémentation du modèle d'Euler Quantique Isotherme, avant de s'attaquer aux modèles non isothermes dans le chapitre V avec l'étude des modèles d'Hydrodynamique Quantique et de Transport d'Énergie Quantique. Enfin, le chapitre VI s'intéresse à un problème un petit peu différent en proposant un schéma asymptotiquement stable dans la limite semi-classique pour l'équation de Schrödinger écrite dans sa formulation fluide: le système de Madelung.

On s’intéresse ici aux erreurs de modélisation liées à l’usage de modèles de flammelette sous-maille en combustion turbulente non prémélangée. Le but de cette thèse est de développer une stratégie d’estimation d’erreur a posteriori pour déterminer le meilleur modèle parmi une hiérarchie, à un coût numérique similaire à l’utilisation de ces mêmes modèles. Dans un premier temps, une stratégie faisant appel à un estimateur basé sur les résidus pondérés est développée et testée sur un système d’équations d’advection-diffusion-réaction. Dans un deuxième temps, on teste la méthodologie d’estimation d’erreur sur un autre système d’équations, où des effets d’extinction et de réallumage sont ajoutés. Lorsqu’il n’y a pas d’advection, une analyse asymptotique rigoureuse montre l’existence de plusieurs régimes de combustion déjà observés dans les simulations numériques. Nous obtenons une approximation des paramètres de réallumage et d’extinction avec la courbe en «S», un graphe de la température maximale de la flamme en fonction du nombre de Damköhler, composée de trois branches et d’une double courbure. En ajoutant des effets advectifs, on obtient également une courbe en «S» correspondant aux régimes de combustion déjà identifiés. Nous comparons les erreurs de modélisation liées aux approximations asymptotiques dans les deux régimes stables et établissons une nouvelle hiérarchie des modèles en fonction du régime de combustion. Ces erreurs sont comparées aux estimations données par la stratégie d’estimation d’erreur. Si un seul régime stable de combustion existe, l’estimateur d’erreur l’identifie correctement ; si plus d’un régime est possible, on obtient une fac˛on systématique de choisir un régime. Pour les régimes où plus d’un modèle est approprié, la hiérarchie prédite par l’estimateur est correcte. / We are interested here in the modeling errors of subgrid flamelet models in nonpremixed turbulent combustion. The goal of this thesis is to develop an a posteriori error estimation strategy to determine the best model within a hierarchy, with a numerical cost at most that of using the models in the first place. Firstly, we develop and test a dual-weighted residual estimator strategy on a system of advection-diffusion-reaction equations. Secondly, we test that methodology on another system of equations, where quenching and ignition effects are added. In the absence of advection, a rigorous asymptotic analysis shows the existence of many combustion regimes already observed in numerical simulations. We obtain approximations of the quenching and ignition parameters, alongside the S-shaped curve, a plot of the maximal flame temperature as a function of the Damköhler number, consisting of three branches and two bends. When advection effects are added, we still obtain a S-shaped curve corresponding to the known combustion regimes. We compare the modeling errors of the asymptotic approximations in the two stable regimes and establish new model hierarchies for each combustion regime. These errors are compared with the estimations obtained by using the error estimation strategy. When only one stable combustion regime exists, the error estimator correctly identifies that regime; when two or more regimes are possible, it gives a systematic way of choosing one regime. For regimes where more than one model is appropriate, the error estimator’s predicted hierarchy is correct.

Les développements asymptotiques topologiques n'ont pas encore été étudiés pour les équations elliptiques quasilinéaires. Cette question apparaît dans la perspective d'appliquer les méthodes d'asymptotique topologique en optimisation de forme aux équations non linéaires de l'élasticité comme en imagerie pour la détection d'ensembles de codimension $\geq 2$ (points en 2D ou courbes en 3D). Dans la Partie I, notre principal résultat réside dans l'obtention du développement asymptotique topologique pour une classe d'équations elliptiques quasilinéaires, perturbées dans des sous-domaines non vides. Le gradient topologique peut être décomposé en un terme linéaire classique et en un terme nouveau, qui rend compte de la non linéarité. L'étude des difficultés spécifiques qui apparaissent avec l'équation de p-Laplace, par comparaison avec l'équation de Laplace, montre qu'un point central réside dans la possibilité de définir la variation de l'état direct à l'échelle 1 dans R^N. Nous étudions en conséquence des espaces de Sobolev à poids et quotientés, dont la semi-norme est la somme des normes L^p et L^2 du gradient dans R^N. Puis nous construisons une classe d'équations elliptiques quasilinéaires, telle que le problème définissant l'état direct à l'échelle 1 vérifie une double propriété de p- et 2- ellipticité. La méthode se poursuit par l'étude du comportement asymptotique de la solution du problème d'interface non linéaire dans R^N et par une mise en dualité appropriée des états directs et adjoints aux différentes étapes d'approximation pour les variations de l'état direct. La Partie II traite d'estimations et de développements asymptotiques de p-capacités de condensateurs, dont l'obstacle est d'intérieur vide et de codimension $\geq 2$. Après quelques résultats préliminaires, nous introduisons les condensateurs équidistants pour étudier le cas des segments. L'effet anisotrope engendré par un segment dans l'équation de p-Laplace est tel que l'inégalité de réarrangement de Pólya-Szegö pour les intégrales de type Dirichlet fournit un minorant trivial. De plus, quand p > N, on ne peut construire par extension une solution admissible pour le segment, aussi petite sa longueur soit-elle, à partir du cas du point. Nous établissons une minoration de la p-capacité N-dimensionnelle d'un segment, qui fait intervenir les p-capacités d'un point, respectivement en dimensions N et (N−1). Les cas de positivité de la p-capacité s'en déduisent. Notre méthode peut être étendue à des obstacles de dimensions supérieures et de codimension $\geq 2$. Introduisant les condensateurs elliptiques, nous montrons que le gradient topologique de la 2-capacité n'est pas un outil approprié pour distinguer les courbes et les obstacles d'intérieur non vide en 2D. Une solution pourrait être de choisir différentes valeurs de p ou bien de considérer le développement asymptotique à l'ordre 2, i.e. la hessienne topologique.

Cette thèse est dédiée à l'étude mathématique et numérique d'une classe d'équations cinétiques collisionnelles, de type équation de Boltzmann. Nous avons porté un intérêt tout particulier à l'équation des milieux (ou gaz) granulaires, initialement introduite dans la littérature physique pour décrire le comportement hors équilibre de matériaux composés d'un grand nombre de grains, ou particules, non nécessairement microscopiques, et interagissant par des collisions dissipant l'énergie cinétique. Ces modèles se sont révélés avoir une structure mathématique très riche. Cette thèse se structure en trois partie pouvant être lues de manière indépendante, mais néanmoins en rapport avec des équations cinétiques collisionnelles en général, et l'équation des milieux granulaires en particulier. La première partie est dédiée à l'étude mathématique du comportement asymptotique de certaines équations cinétiques collisionnelles dans un cadre homogène en espace. Nous y montrons des résultats de type explosion et convergence vers la solution autosimilaire avec calcul explicite des taux, pour des opérateurs de type Boltzmann, grâce à l'utilisation (entre autre) d'une nouvelle méthode de changement de variables dépendant directement de la solution de l'équation considérée. En particulier, nous démontrons que pour un modèle de gaz granulaire - dit anormal - il est possible d'observer une explosion en temps fini. Dans la deuxième partie, orientée analyse numérique et calcul scientifique, nous nous intéressons développement et à l'étude de méthodes spectrales pour la résolution de problèmes multi-échelles, issus de la théorie des équations cinétiques collisionnelles. Les méthodes de changement de variables tiennent aussi une place importante dans cette partie, et permettent d'observer numériquement des phénomènes non triviaux qui apparaissent lors de l'étude de gaz granulaires, comme la création d'amas de matière ou la caractérisation précise du retour vers l'équilibre. La troisième et dernière partie est dédiée à l'étude spectrale de l'opérateur des milieux granulaires avec bain thermique, linéarisé au voisinage d'un équilibre homogène en espace, afin d'établir des résultats de type stabilité et convergence vers une limite hydrodynamique. Ce travail est en fait la généralisation d'un résultat célèbre dans la théorie de l'équation de Boltzmann, dû à R. Ellis et M. Pinsky, et établissant rigoureusement la première limite hydrodynamique vers les équations d'Euler compressibles linéaires puis Navier-Stokes de cette équation.

Cette thèse de doctorat est consacrée à l'étude mathématique de deux types de problèmes de jonction. Le premier modèle est obtenu comme limite variationnelle d'un problème d'élasticité avec jonction de faible épaisseur, l'adhésif occupant la jonction a une rigidité de l'ordre de l'épaisseur. La densité élastique W de l'adhésif est complétée par une densité d'énergie surfacique convexe h non nécessairement régulière. Cette densité traduit une contrainte mécanique entre adhérents et l'adhésif imparfaitement collés. On montre que le modèle limite consiste à remplacer la jonction par une contrainte qui est l'inf-convolution de h et de la densité surfacique limite de W. Dans un cadre scalaire on effectue l'analyse des concentrations de gradient à l'interface au moyen d'outils récents issus de la théorie de la mesure. Dans le second modèle, la rigidité de l'adhésif est de l'ordre inverse de l'épaisseur de le jonction, la densité élastique du matériaux adhérent a une croissance super-linéaire alors que celle de l'adhésif croit linéairement. Suivant la stratégie utilisée pour le premier problème on propose un modèle simplifié mais fiable comme limite variationnelle lorsque l'épaisseur tend vers zéro. A la limite la fine couche intermédiaire est remplacée par une interface pseudo-plastique prédisant la formation de fissures.

Les structures maçonnées sont largement connues en génie civil comme constituant d’une partie des bâtiments, mais également en tant que garnissages réfractaires dans des structures utilisées à hautes températures, par exemple en sidérurgie. Malheureusement, les outils actuels ne sont pas suffisamment puissants pour prédire le comportement de ces structures avec l’apparition de fissures et pour tenir compte du comportement non linéaire d’un des deux constituants (le mortier par exemple). Ce travail de thèse contribue à la modélisation multi-échelles des maçonneries classiques et des garnissages réfractaires avec un faible coût numérique grâce à la technique d’homogénéisation périodique. Les techniques de modélisation et de simulation du comportement des maçonneries sont présentées et développées. L’influence des lois d’interface entre briques et mortier, des paramètres géométriques et matériels, ainsi que de la densité des fissures sur le comportement effectif des maçonneries est étudié. Trois approches (une extension analytique de Cecchi et Tralli, une approche numérique et un modèle micromécanique) sont proposées pour la détermination du comportement effectif d’une cellule périodique dans le cas de maçonneries avec mortier viscoélastique microfissuré et briques saines élastiques ou rigides. Les résultats des calculs sur deux exemples de maçonneries (1D et 2D) ont confirmé que l’approche multi-échelle est une solution appropriée avec une grande capacité à exprimer le comportement des maçonneries viscoélastiques microfissurées. Ce travail, limité au cas sans propagation de fissures, peut être étendu aux mortiers à comportement viscoplastique. / Masonry structures are widely used in civil engineering as part of buildings or in refractory linings of structures working at high temperatures, like in steel industry. Unfortunately, the present tools are not powerful enough to predict the behavior of these structures at their micro-cracked state and/or if one of their constituents behaves nonlinearly (e.g. the mortar). This research contributes to the multi-level modeling of classical masonries and refractory linings with a low numerical cost using basically the periodic homogenization technique. Modeling and simulation techniques of masonry behavior are presented and developped. The influence of interface law between bricks and mortar, of geometrical and material parameters, and of crack density on the effective masonry behavior is studied. Three approaches (analytical extension of Cecchi and Tralli, numerical approach and micromechanical modeling) were proposed to determine the effective behavior of a periodic masonry cell with micro-cracked viscoelastic mortar and safe elastic or rigid bricks. The results obtained on two examples of masonry (1D and 2D) confirmed that the multi-scale approach is a suitable solution with a great ability to model the effective behavior of microcracked viscoelastic masonry. This work, actually limited to the case without crack propagation, could be extended to mortars with viscoplastic behavior.

La Partie I présente l’obtention du développement asymptotique topologique pour une classe d’équations elliptiques quasilinéaires. Un point central réside dans la possibilité de définir la variation de l’état direct à l’échelle 1 dans R^N. Après avoir défini un cadre fonctionnel approprié faisant intervenir les normes L^p et L^2, et avoir justifié la classe d’équations considérée, la méthode se poursuit par l’étude du comportement asymptotique de la solution du problème d’interface non linéaire dans R^N et par une mise en dualité appropriée des états direct et adjoint aux différentes étapes d’approximation.La Partie II traite d’estimations et de développements asymptotiques de p-capacités de condensateurs, dont l’obstacle est d’intérieur vide et de codimension > ou = 2. Après les résultats préliminaires, les condensateurs équidistants permettent de donner deux illustrations de l’anisotropie engendrée par un segment dans l’équation de p-Laplace, puis d’établir une minoration de la p-capacité N-dimensionnelle d’un segment, qui fait intervenir les p-capacités d’un point, respectivement en dimensions N et (N-1). Les condensateurs elliptiques permettent d’établir que le gradient topologique de la 2-capacité n’est pas un outil approprié pour distinguer les courbes des obstacles d’intérieur non vide en 2D / Part I deals with obtaining topological asymptotic expansions for a class of quasilinear elliptic equations. A key point lies in the ability to define the variation of the direct state at scale 1 in R^N. After setting up an appropriate functional framework involving both the L^p and the L^2 norms, and then justifying the chosen class of equations, the approach goes on with the study of the asymptotic behavior of the solution of the nonlinear interface problem in R^N and by setting up an adequate duality scheme between the direct and adjoint states at each step of approximation. Part II deals with estimates and asymptotic expansions of condenser p-capacities and focuses on obstacles with empty interiors and with codimensions > ou = 2. After preliminary results, equidistant condensers are introduced to point out the anisotropy caused by a segment in the p-Laplace equation, and to provide a lower bound to the N-dimensional condenser p-capacity of a segment, by means of the N-dimensional and of the (N-1)-dimensional condenser p-capacities of apoint. Introducing elliptical condensers, it turns out that the topological gradient of the 2-capacity is not an appropriate tool to separate curves and obstacles with nonempty interior in 2D