Les principaux thèmes de cette thèse sont la théorie des trajectoires rugueuses développée par T. Lyons (1998) et ses applications, notamment à l'étude des équations différentielles stochastiques (EDS) et au calcul de sensibilités. Des applications potentielles des résultats théoriques en science du vivant et en finance y sont également développés. En premier lieu, sont étendues l'existence et l'expression des grecques Delta et Véga, sensibilités bien connues en finance, pour des EDS à coefficients bornés et dirigées par un processus gaussien multidimensionnel centré, à trajectoires continues, au-dessus duquel il existe une trajectoire géométrique naturelle. Le cas du mouvement brownien fractionnaire a particulièrement été développé afin de proposer d'une part, une application du calcul de Véga dans un modèle de marché financier à volatilité stochastique fractionnaire et d'autre part, d'effectuer des simulations. En second lieu, est étudiée une généralisation d'équation mean-reverting au cas d'un signal gaussien unidimensionnel, centré et à trajectoires continues : existence globale et unicité de la solution, intégrabilité, continuité et différentiabilité de l'application d'Itô, existence d'un schéma d'approximation convergeant dans tous les Lp avec une vitesse de convergence presque-sure, un principe de grandes déviations et, l'existence d'une densité par rapport à la mesure de Lebesgue. L'étude de cette famille d'EDS a débouché sur une application en pharmacocinétique.

Les problèmes d’optimisation discrète sont pour beaucoup difficiles à résoudre, de par leur nature combinatoire. Citons par exemple les problèmes de programmation linéaire en nombres entiers. Une approche couramment employée pour les résoudre exactement est l’approche de Séparation et Évaluation Progressive. Une approche différente appelée « Resolution Search » a été proposée par Chvátal en 1997 pour résoudre exactement des problèmes d’optimisation à variables 0-1, mais elle reste mal connue et n’a été que peu appliquée depuis. Cette thèse tente de remédier à cela, avec un succès partiel. Une première contribution consiste en la généralisation de Resolution Search à tout problème d’optimisation discrète, tout en introduisant de nouveaux concepts et définitions. Ensuite, afin de confirmer l’intérêt de cette approche, nous avons essayé de l’appliquer en pratique pour résoudre efficacement des problèmes bien connus. Bien que notre recherche n’ait pas abouti sur ce point, elle nous a amené à de nouvelles méthodes pour résoudre exactement les problèmes d’affectation généralisée et de localisation simple. Après avoir présenté ces méthodes, la thèse conclut avec un bilan et des perspectives sur l’application pratique de Resolution Search. / The combinatorial nature of discrete optimization problems often makes them diffi- cult to solve. Consider for instance integer linear programming problems, which are commonly solved using a Branch-and-Bound approach. An alternative approach, Resolution Search, was proposed by Chvátal in 1997 for solving 0-1 optimization problems, but remains little known to this day and as such has seen few practical applications. This thesis attempts to remedy this state of affairs, with partial success. Its first contribution consists in the generalization of Resolution Search to any discrete optimization problem, while introducing new definitions and concepts. Next, we tried to validate this approach by attempting to solve well-known problems efficiently with it. Although our research did not succeed in this respect, it lead us to new methods for solving the generalized assignment and uncapacitated facility location problems. After presenting these methods, this thesis concludes with a summary of our attempts at practical application of Resolution Search, along with further perspectives on this matter. / Thèse réalisée en cotutelle avec l'Université d'Avignon.

Ce travail présente une contribution à la modélisation numérique des systèmes de suspensions dans le cadre de la théorie cinétique. Cette description continue des systèmes de suspensions permet de prendre en compte l'influence de la structure à l'échelle microscopique sur la cinétique de l'écoulement macroscopique. Cependant elle présente l'inconvénient majeur d'être définie sur un espace à haute dimension et rend alors difficile la résolution de ces modèles avec des approches déterministes classiques. Afin de s'affranchir, ou du moins d'alléger, le poids du caractère micro-macro des approches en théorie cinétique, plusieurs techniques de réduction dimensionnelle s'appuyant sur l'utilisation de la Décomposition Généralisée en modes Propres (PGD) sont présentées. Une étude de différents algorithmes PGD est conduite, et dont l'efficacité en termes de vitesse de convergence et d'optimalité de la solution est illustrée. La simulation de mélanges de fluides immiscibles est conduite à l'aide du Tenseur d'aire qui est un puissant outil de caractérisation du mélange. Cependant celui-ci nécessite l'introduction d'une relation de fermeture dont l'impact est évalué avec le modèle de théorie cinétique équivalent et exact. Finalement, la simulation de systèmes de suspensions colloïdales décrits par l'équation de Smoluchowski présente une approche originale de la modélisation des suspensions solides. Cette approche permet de s'affranchir avantageusement du bruit statistique inhérent aux simulations stochastiques traditionnellement mises en œuvre.

La structure des nucléons et des noyaux a été intensivement étudiée au cours duvingtième siècle au travers de la diffusion élastique d'électrons (mesure des facteurs deforme électromagnétique) et de la diffusion profondément inélastique (mesure des distributionsde partons). Le formalisme des distributions généralisées de partons (GPD)a permis d'unifier les facteurs de forme et les distributions de partons. Ce lien procureune source d'information unique sur la dynamique des partons, telle la distribution desforces nucléaires et de moment orbital au sein des hadrons. L'accès expérimental le plussimple aux GPD est la diffusion Compton profondément virtuelle (DVCS), correspondantà l'électroproduction dure d'un photon réel. Tandis que plusieurs expériences sesont déjà focalisées sur la réaction DVCS sur le nucléon, les expériences sur une ciblenucléaire s'avèrent plus rares. Cette thèse se concentre sur l'étude du canal DVCS cohérentsur l'hélium 4, avec pour objectif l'extraction des parties réelle et imaginaire dufacteur de forme Compton via l'asymétrie de spin du faisceau.

Les progrès des dispositifs d'imagerie médicale permettent l'obtention d'images volumiques, qui contiennent une grande quantité d'information. Une approche efficace de traitement de ces images consiste à utiliser la connaissance a priori de la forme des objets à analyser, et à employer des méthodes intrinsèquement tridimensionnelles. Les modèles déformables répondent à ces deux critères. Nous proposons de formaliser les modèles déformables et leur évolution dans une image dite de données, en distinguant cinq composantes : caractéristiques de liaison, représentation géométrique, déformation, déformabilité, et contrôle. Nous décrivons trois modèles déformables. Nous employons le modèle surfacique des delta-snakes pour reconstruire des objets à partir de points répartis sur leur surface. Nous approximons cette surface par une carte de distance octree-spline. Nous avons mis au point des outils interactifs pour compléter des données manquantes ou déformer directement la surface. Nous proposons ensuite pour ce même type d'application un modèle implicite à base de primitives générant un champ potentiel local. Les primitives sont placées interactivement, ou automatiquement sélectionnées dans l'axe médian discret des données. L'optimisation des paramètres des primitives mène à une représentation compacte des objets. Nous reconstruisons par ces deux modèles des objets numérisés par des capteurs de distance ou segmentés dans des images volumiques. Notre dernier modèle est volumique. Sa déformation hiérarchique par un octree-spline minimise la distance généralisée entre ses caractéristiques et celles des données, sous le contrôle de l'algorithme de Levenberg-Marquardt, et dans les limites imposées par une fonction de régularisation. Nous avons établi un algorithme de calcul de distance généralisée itérée dans un arbre k-d. Nous appliquons ce modèle à la segmentation d'images volumiques. D'autres types d'applications ont également été réalisées.

Les alliages à base nickel, et plus précisément les alliages 600 et 690, sont, en France, les matériaux constitutifs des tubes de générateurs de vapeur (GV) des réacteurs à eau pressurisée (REP). La corrosion généralisée résultant de l'interaction entre ces alliages et le milieu primaire des REP conduit, d'une part, à la formation d'une fine couche passive (environ 10 nm) de produits de corrosion, mais également à la dissolution de cations dans le circuit primaire, ces derniers constituant la principale source de contamination de ce circuit. L'objectif de ces travaux est d'apporter de nouveaux éléments de compréhension des phénomènes de corrosion des alliages à base nickel en milieu de type primaire, en s'attachant à décrire l'influence de paramètres métallurgiques ou physico-chimiques, sur la nature et les mécanismes de croissance de la couche d'oxyde protectrice. Dans ce contexte, le film passif formé lors de l'exposition des alliages 600, 690 et Ni-30Cr, dans des conditions simulant le milieu primaire des REP, a été analysé par plusieurs techniques de caractérisation (MEB, MET, PEC et MPEC, XPS). Le couplage de ces méthodes a permis d'obtenir une description fine, en termes de nature et de structure, de l'oxyde multicouche se formant à la surface des alliages à base nickel en milieu de type primaire. Ainsi, la part protectrice de cette couche d'oxyde est constituée d'une couche continue de chromite mixte de nickel et de fer et de nodules de Cr₂O₃ dispersés à l'interface entre le chromite mixte et l'alliage. L'étude des mécanismes de croissance de la couche protectrice notamment à partir d'expériences mettant en œuvre des marqueurs et des traceurs, a permis de conclure à une croissance de la couche de chromite mixte résultant d'une diffusion anionique, via des courts-circuits de diffusion de type joints de grains. De plus, l'impact du taux de défauts en surface de l'alliage a également été étudié, mettant en évidence un effet double de ce paramètre, qui influence à la fois la densité de courts-circuits de diffusion au sein de l'oxyde mais également la vitesse de formation des nodules de Cr₂O₃. L'ensemble de ces résultats permet de proposer un mécanisme de corrosion des alliages à base nickel en milieu de type primaire et d'aborder quelques aspects cinétiques susceptibles de constituer des éléments de réflexion en vue du développement d'une modélisation complète de ce phénomène de corrosion.

Représentations séparées pour la simulation multi-échelle du comportementmécanique et de l’endommagement des matériaux composites.Résumé: Le développement de méthodes numériques performantes pour simuler les structurescomposites est un défi en raison de la nature multi-échelle et de la complexité des mécanismed’endommagement de ce type de matériaux. Les techniques classiques de discrétisationvolumique conduisent à des coûts de calcul importants et sont restreintes en pratique à despetites structures.Dans cette thèse, un nouvelle stratégie basée sur une représentation séparée de la solution estexplorée. L’objectif est de proposer un cadre numérique efficace et fiable pour analyser les endommagementsdans les composites stratifiés sous chargements statiques et dynamiques. Ladécomposition propre généralisée (PGD) est utilisée pour construire la solution.Pour traiter l’endommagement, et plus particulière le délaminage, un modèle de zone cohésivea été implémenté dans la PGD. Une approches multi-échelle innovante est également proposéepour simuler le comportement mécanique des composites à microstructure périodique. L’idéeprincipale est de séparer deux échelles : l’échelle du motif périodique (microstructure) et l’échellemacroscopique. Les résultats de la PGD sont très proches des résultats obtenus par la méthodeéléments fini classique. Finalement, la PGD permet de réduire significativement la complexitédes modèles tout en gardant une précision satisfaisante. / Separated representations for the multiscale simulation of the mechanicalbehavior and damages of composite materials.Abstract: The development of efficient simulations for composite structures is very challengingdue to the multiscale nature and the complex damage process of this materials. When usingstandard 3D discretization techniques with advanced models for large structures, the computationalcosts are generally prohibitive.In this thesis, a new strategy based on a separated represenation of the solution is explored todevelop a computationally efficient and reliable numerical framework for the analysis of damagesin laminated composites subjected to quasi-static and dynamic loading. The PGD (Proper GeneralizedDecomposition) is used to build the solution.To treat damage, and especially delamination, a cohesive zone model has been implemented inthe PGD solver. A novel multiscale approach is also proposed to compute the mechanical behaviorof composites with periodic microstructure. The idea is to separate two scales: the scaleof periodic pattern and the macroscopic scale. The PGD results have been compared with theresults obtained with the classcial finite element method. A close agreement is found between thetwo approach and the PGD has significantly reduced the model complexity.

Ce travail présente une contribution à la modélisation numérique des systèmes de suspensions dans le cadre de la théorie cinétique. Cette description continue des systèmes de suspensions permet de prendre en compte l'influence de la structure à l'échelle microscopique sur la cinétique de l'écoulement macroscopique. Cependant elle présente l'inconvénient majeur d'être définie sur un espace à haute dimension et rend alors difficile la résolution de ces modèles avec des approches déterministes classiques. Afin de s'affranchir, ou du moins d'alléger, le poids du caractère micro-macro des approches en théorie cinétique, plusieurs techniques de réduction dimensionnelle s'appuyant sur l'utilisation de la Décomposition Généralisée en modes Propres (PGD) sont présentées. Une étude de différents algorithmes PGD est conduite, et dont l'efficacité en termes de vitesse de convergence et d'optimalité de la solution est illustrée. La simulation de mélanges de fluides immiscibles est conduite à l'aide du Tenseur d'aire qui est un puissant outil de caractérisation du mélange. Cependant celui-ci nécessite l'introduction d'une relation de fermeture dont l'impact est évalué avec le modèle de théorie cinétique équivalent et exact. Finalement, la simulation de systèmes de suspensions colloïdales décrits par l'équation de Smoluchowski présente une approche originale de la modélisation des suspensions solides. Cette approche permet de s'affranchir avantageusement du bruit statistique inhérent aux simulations stochastiques traditionnellement mises en œuvre. / This work is a contribution to the numerical modeling of suspension system in the kinetic theory framework. This continuum description of suspension system allows to account for the microstructure impact on the kinetic of the macroscopic flow. However, its main drawback is related to the high dimensional spaces in which kinetic theory models are defined and makes difficult for classical deterministic approaches to solve such systems. One possibility for circumventing, or at least alleviate, the weight of the micro-macro kinetic theory approaches lies in the use of separated representations strategies based on the Proper Generalized Decomposition (PGD). A study of different PGD algorithms is driven, illustrating the efficiency of these algorithms in terms of convergence speed and optimality of the solution obtained. The immiscible fluids blends modeling is driven using the area tensor which is a powerful numerical tool for characterizing blends. However it needs the introduction of closure relation of which impact is measured using equivalent and exact kinetic theory model. Finally, the numerical modeling of colloidal suspension system described by the Smoluchowski equation presents an original approach of the modeling of solid suspension system. This description allows to circumvent the statistical noise inherent to the stochastic approaches commonly used.

Cette thèse examine la dispersion multivariée des modelés normales stables Tweedie. Trois estimateurs de fonction variance généralisée sont discutés. Ensuite dans le cadre de la famille exponentielle naturelle deux caractérisations du modèle normal-Poisson, qui est un cas particulier de modèles normales stables Tweedie avec composante discrète, sont indiquées : d'abord par fonction variance et ensuite par fonction variance généralisée. Le dernier fournit la solution à un problème particulier d'équation de Monge-Ampère. Enfin, pour illustrer l'application de la variance généralisée des modèles Tweedie stables normales, des exemples à partir des données réelles sont fournis. / This thesis examines the multivariate dispersion of normal stable Tweedie (NST) models. Three generalize variance estimators of some NST models are discussed. Then within the framework of natural exponential family, two characterizations of normal Poisson model, which is a special case of NST models with discrete component, are shown : first by variance function and then by generalized variance function. The latter provides a solution to a particular Monge-Ampere equation problem. Finally, to illustrate the application of generalized variance of normal stable Tweedie models, examples from real data are provided.

Ce travail de recherche est une contribution au développement de la méthode Décomposition Propre Généralisée (PGD) à la résolution de problèmes de diffusion-réaction raides dédiés à la mécanique des matériaux. Ce type d’équations est notamment rencontré lors de l’oxydation des matériaux polymères et il est donc nécessaire de mettre en place un outil pour simuler ce phénomène afin de prédire numériquement le vieillissement de certains matériaux composites à matrice organique utilisés dans l’aéronautique. La méthode PGD a été choisie dans cette thèse car elle permet un gain en temps de calcul notable par rapport à la méthode des éléments finis. Néanmoins cette famille d’équations n’a jamais été traitée avec cette méthode. Cette dernière se résume à la recherche de solutions d’Équations aux Dérivées Partielles sous forme séparée. Dans le cas d’un problème 1D transitoire, cela revient à chercher la solution sous la forme d’une représentation séparée espace-temps. Dans le cadre de cette thèse, un outil numérique a été mis en place permettant une flexibilité telle que différents algorithmes peuvent être testés. La diffusion Fickienne 1D est tout d’abord évaluée avec en particulier une discussion sur l’utilisation d’un schéma de type Euler ou Runge-Kutta à pas adaptatif pour la détermination des fonctions temporelles. Le schéma de Runge-Kutta permet de réduire notablement le temps de calcul des simulations.Ensuite, la mise en place de l’outil pour les systèmes d’équation de type diffusion-réaction nécessite des algorithmes de résolution de systèmes non linéaires, couplés et raides. Pour cela, différents algorithmes ont été implémentés et discutés.Dans le cas d’un système non linéaire, l’utilisation de la méthode de Newton-Raphson dans les itérations pour la recherche du nouveau mode permet de réduire le temps de calcul en limitant le nombre de modes à considérer pour une erreur donnée. En ce qui concerne les couplages, deux stratégies de résolution ont été évaluées. Le couplage fort mène aux mêmes conclusions que dans le cas non linéaire. Les systèmes raides mais linéaires ont ensuite été traités en implémentant l’algorithme de Rosenbrock pour la détermination des fonctions temporelles. Cet algorithme permet contrairement à Euler et à Runge-Kutta de construire une solution avec un temps de calcul raisonnable liée à l’adaptation du maillage temporel sous-jacent à l’utilisation de cette méthode. La résolution d’un système d’équations de diffusion-réaction raides non linéaires utilisée pour la prédiction de l’oxydation d’un composite issu de la littérature a été testée en utilisant les différents algorithmes mis en place. Néanmoins, les non linéarités et la raideur du système génèrent des équations différentielles intermédiaires à coefficients variables pour lesquelles la méthode de Rosenbrock montre ses limites. Il sera donc nécessaire de tester ou développer d’autres algorithmes pour lever ce verrou.Mots / This work presents the development of the Proper Generalized Decomposition (PGD) method for solving stiff reaction-diffusion equations in the framework of mechanics of materials. These equations are particularly encountered in the oxidation of polymers and it is therefore necessary to develop a tool to simulate this phenomenon for example for the ageing of organic matrix composites in aircraft application. The PGD method has been chosen in this work since it allows a large time saving compared to the finite element method. However this family of equations has never been dealt with this method. The PGD method consists in approximating a solution of a Partial Differential Equation with a separated representation. The solution is sought under a space-time separated representation for a 1D transient equation.In this work, a numerical tool has been developed allowing a flexibility to test different algorithms. The 1D Fickian diffusion is first evaluated and two numerical schemes, Euler and Runge-Kutta adaptive methods, are discussed for the determination of the time modes. The Runge-Kutta method allows a large time saving. The implementation of the numerical tool for reaction-diffusion equations requires the use of specific algorithms dedicated to nonlinearity, couplingand stiffness. For this reason, different algorithms have been implemented and discussed. For nonlinear systems, the use of the Newton-Raphson algorithm at the level of the iterations to compute the new mode allows time saving by decreasing the number of modes required for a given precision. Concerning the couplings, two strategies have been evaluated. The strong coupling leads to the same conclusions as the nonlinear case. The linear stiff systems are then studied by considering a dedicated method, the Rosenbrock method, for the determination of the time modes. This algorithm allows time saving compared to the Runge-Kutta method. The solution of a realistic nonlinear stiff reaction-diffusionsystem used for the prediction of the oxidation of a composite obtained from the literature has been tested by using the various implemented algorithms. However, the nonlinearities and the stiffness of the system generate differential equations with variable coefficients for which the Rosenbrock method is limited. It will be necessary to test or develop other algorithms to overcome this barrier.