Cette étude s’inscrit dans le domaine de la simulation électromagnétique des problèmes de grande taille tels que la diffraction d’ondes planes par de larges plateformes et le rayonnement d’antennes aéroportées. Elle consiste à développer une méthode combinant décomposition en sous-domaines et compression hiérarchique des équations intégrales de frontière. Pour cela, nous rappelons dans un premier temps les points importants de la méthode des équations intégrales de frontière et de leur compression hiérarchique par l’algorithme ACA (Adaptive Cross Approximation). Ensuite, nous présentons la formulation IE-DDM (Integral Equations – Domain Decomposition Method) obtenue à partir d’une représentation intégrale des sous-domaines. Les matrices résultant de la discrétisation de cette formulation sont stockées au format H-matrice (matricehiérarchique). Un solveur spécialement adapté à la résolution de la formulation IE-DDM et à sa représentation hiérarchique a été conçu. Cette étude met en évidence l’efficacité de la décomposition en sous-domaines en tant que préconditionneur des équations intégrales. De plus, la méthode développée est rapide pour la résolution des problèmes à incidences multiples ainsi que la résolution des problèmes basses fréquences / This thesis is about the electromagnetic simulation of large scale problems as the wave scattering from aircrafts and the airborne antennas radiation. It consists in the development of a method combining domain decomposition and hierarchical compression of the surface integral equations. First, we remind the principles of the boundary element method and the hierarchical representation of the surface integral equations with the Adaptive Cross Approximation algorithm. Then, we present the IE-DDM formulation obtained from a sub-domain integral representation. The matrices resulting of the discretization of the formulation are stored in the H-matrix format. A solver especially fitted with the hierarchical representation of the IE-DDM formulation has been developed. This study highlights the efficiency of the sub-domain decomposition as a preconditioner of the integral equations. Moreover, the method is fast for the resolution of multiple incidences and the resolution of low frequencies problems

Cette thèse est dédiée à la résolution numérique tridimensionnelle des équations de Maxwell harmoniques, par des méthodes de décomposition de domaine couplant des résolutions par équations intégrales entre elles. Pour traiter les problèmes de diffraction d'ondes, la méthode des équations intégrales est un outil précieux. Elle consiste à paramétrer le champ électromagnétique solution par une source définie sur la surface de l'objet diffractant, solution d'une nouvelle équation linéaire (l'équation intégrale). Pour des applications à haute fréquence, le grand nombre d'inconnues (de l'ordre du million) nous oblige à utiliser un solveur itératif pour résoudre l'équation intégrale. Le problème du conditionnement des systèmes linéaires est alors crucial. De récents développements ont permis de construire une équation intégrale performante (la GCSIE) et de conditionnement stable avec la montée en fréquence. Cependant, la présence d'une cavité large et résonnante dans l'objet diffractant (telle que la cavité moteur d'un avion) dégrade le conditionnement de cette équation. Nous proposons deux méthodes de décomposition de domaine (DDM) afin de découpler le problème de la cavité du problème extérieur. La première (DDM en Y) s'exprime en fonction des opérateurs Dirichlet-to-Neumann Y, qui sont synthétisés via la résolution de problèmes métalliques par équations intégrales dans chaque sous-domaine. La seconde (DDM en S) s'exprime en fonction des opérateurs de scattering S, synthétisés par résolution de problèmes de type métal-impédant, donc bien posés à toute fréquence. La DDM en S permet ainsi de se débarrasser des phénomènes de résonance dans les cavités. Nous proposons dans un premier temps un préconditionneur analytique pour la DDM en Y, basé sur l'opérateur électromagnétique de simple couche. Nous calculons ensuite les modes guidés le long d'un cylindre infini tangent à la cavité près de l'interface, et nous diagonalisons les opérateurs Dirichlet-to-Neumann et scattering dans la base des traces de modes guidés sur l'interface. On extrait de cette étude deux préconditionneurs spectraux respectivement pour la DDM en Y et la DDM en S. Les résultats numériques confirment l'efficacité des préconditionneurs proposés / This work is dedicated to the numerical solution of the tridimensional harmonic Maxwell equations, using domain decomposition methods coupling integral equations between them. To deal with scattering problems, integral equations methods are a precious tool. They allow to look for the electromagnetic field by parameterizing it with a source only defined on the boundary of the scattering object, solution of a new linear equation (the integral equation). For applications at high frequency, the great number of unknowns forces the use of iterative methods. To accelerate the solution of integral equations, one moreover has to ensure the good condition number of the linear systems, or to propose well-suited preconditioners. An efficient method, the GCSIE, was developed in Onera. It is an intrinsically well-conditioned integral equation whose condition number remains stable whith the frequency increase. However, the existence of large and resonant cavities (such as air intakes) deteriorates the condition number. In order to circumvent this problem, we propose two domain decomposition methods (DDM) allowing to decouple the exterior problem from the problem of the cavities. The first one (Y-DDM) is based on Dirichlet-to-Neumann operators Y, which are built through the solution of metallic problems using integral equations in each subdomain. The second one (S-DDM) is based on scattering operators S, built through the solution of problems of metallic-impedant type, which are well-posed at any frequency. The S-DDM allows to avoid the resonance phenomena inside the cavities. First, we propose an analytic preconditioner for the Y-DDM, based on the electromagnetic single layer operator. We then calculate the modes guided along an artificial infinite cylinder, that is tangent to the cavity near the interface. We diagonalize the Dirichlet-to-Neumann and scattering operators in the basis of the traces of the guided modes on the interface. We deduce from this study two spectral preconditioners for the Y-DDM and the S-DDM. The numerical results confirm the efficiency of the employed preconditioners.

Ce travail de thèse a pour but de construire et d’étudier une fibration de Hitchin pour les groupes qui apparaît naturellement lorsque l’on essaie de géométriser la formule des traces. On commence par construire une telle fibration en utilisant le semi-groupe de Vinberg. Sur ce semi-groupe de Vinberg, on montre qu’il existe un certain morphisme « polynôme caractéristique » muni d’une section naturelle, de même que dans le cas des algèbres de Lie. On montre également que l’on peut construire un centralisateur régulier au-dessus de cette base des polynômes caractéristiques qui est un schéma en groupes commutatif et lisse.On s’intéresse alors à des variantes pour les groupes des fibres de Springer affines pour lesquelles on remarque que l’introduction du semi-groupe de Vinberg permet d’obtenir une condition d’intégralité analogue à celle de Kazhdan-Lusztig. Ces fibres de Springer affines sont des analogues locaux des fibres de Hitchin. On obtient alors une formule de dimension pour ces fibres.Dans un troisième temps, on s’intéresse à l’aspect global de cette fibration pour laquelle on donne une interprétation modulaire et sur laquelle on construit l’action d’un champ de Picard, issu du centralisateur régulier. L’espace total de cette fibration étant en général singulier, nous étudions son complexe d’intersection. Cet espace de Hitchin s’obtient naturellement comme l’intersection du champ de Hecke avec la diagonale du champ des G-torseurs et on démontre que sur un ouvert suffisamment gros de la base de Hitchin, le complexe d’intersection de l’espace de Hitchin s’obtient par restriction de celui du champ de Hecke corrrespondant.Enfin, dans la dernière partie de cette thèse, on établit un théorème du support dans le cas où l’espace total est singulier analogue à celui de Ngô et l’on démontre que, dans le cas de la fibration de Hitchin, les supports qui interviennent sont reliés aux strates endoscopiques. / This main goal of this work is to construct and study the properties of Hitchin fibration for groups which appears naturally when we try to geometrize the trace formula. We begin by constructing this fibration using the Vinberg’s semigroup. On this semigroup, we show that there exists a characteristic polynomial morphism equipped with a natural section, analog at the Kostant’s one in the case of Lie algebras. We also show that there exists on the base of characteristic polynomials a regular centralizer scheme, which is a smooth commutative group scheme.Then, we are interested in some variant of affine Springer fibers, for which we see that the Vinberg’s semigroup appears naturally to obtain an integrality condition analog to Kazhdan-Lusztig’s one. These affine Springer fibers are local incarnation of Hitchin fibers.In a third time, we go back to the global case and give a modular interpretation of this new Hitchin fibration on which we construct an action of a Picard stack, coming from the regular centralizer.The total space of this fibration, even on the generically regular semisimple locus will be singular and we want to understand his intersection complex. This space can be obtained as the intersection of the Hecke stack with the diagonal of the stack of G-bundles and we show that on a sufficiently big open subset of the Hitchin base, the intersection complex of the Hitchin’s space is the restriction of the corresponding intersection complex on the Hecke stack.Finally, in the last part of this work, we establish a support theorem in the case of a singular total space, generalizing Ngo’s theorem et we show that in the case of Hitchin fibration, the supports that appear are related to the endoscopic strata.

Les chocs thermiques provoquent, en général, l’endommagement et la fissuration des matériaux. Ces phénomènes sont observés, par exemple, dans le revêtement de barrière thermique pour les moteurs des turbines, le traitement des surfaces ou la soudure par laser etc. Plusieurs travaux de recherche ont été réalisés au cours des dernières décennies dans l’objectif d’améliorer les performances thermiques et/ou mécaniques des matériaux sous chargement thermique. L’étude des dommages et de la fissuration des matériaux provoqués par les chocs thermiques, tels que le décollement des interfaces et de décohésion de revêtements, a reçu également une attention considérable par les chercheurs. La majorité de ces travaux utilisent les théories classiques, tels que la loi de Fourier de conduction thermique et l'hypothèse de quasi-statique. Malheureusement ces théories ne sont pas adaptées dans le cas de charges extrêmes provoqués par le choc thermique et dans le cas des matériaux micro-fissurés. En conséquence, les théories conventionnelles doivent être enrichies.L'objectif de la thèse est de montrer le rôle crucial des termes non Fourier et les termes inertiels dans le cas de choc thermique sous conditions sévères et dans le cas où les fissures sont petites. Pour cela nous avons mené des études sur deux structures particulières soumises à des chocs thermiques. Chaque structure contient une fissure parallèle au bord libre de la structure située au voisinage de ce dernier. L’influence de la présence de fissure sur la conductivité thermique est prise en compte. Nous avons utilisé la théorie Hyperbolique de transfert de chaleur par conduction pour les champs thermique et mécanique à la place de la théorie traditionnelle classique de Fourier. Pour mener cette étude, nous avons utilisé les Transformées de Laplace et de Fourier aux équations de mouvement et à l’équation de transfert de chaleur. En s’intéressant en particulier aux champs de contrainte au voisinage de la pointe de fissure et aux facteurs d'intensité de contrainte dynamiques. Le problème se ramène à la résolution d’un système d'équations intégrales singulières dans l'espace de Laplace-Fourier. On utilise une méthode d'intégration numérique pour obtenir les différents champs. Nous résolvons ensuite un système d'équations algébriques linéaires. En effectuant des inversions numériques des transformées, nous obtenons les champs de contrainte de température et les facteurs d'intensité de contrainte dynamiques dans le domaine temporel.Les résultats numériques montrent que la conductivité thermique du milieu est affectée par l’ouverture de la fissure ce qui perturberait fortement le champ de température ainsi que l'amplitude des facteurs d'intensité de contrainte dynamiques. Les amplitudes sont supérieures à celles obtenues à partir de la théorie classique de Fourier ainsi que dans le cadre de l'hypothèse quasi-statique. On constate également qu’elles oscillent au cours du temps. La prise en compte simultanément de l’influence de la fissure sur la conductivité thermique, de l'effet non-Fourier ainsi que les effetsIVd'inertie induit un couplage entre les trois phénomènes qui rendrait le problème de choc thermique très complexe. L'effet de barrière thermique induit par la fissure affecte d’une manière significative les champs de température et des contraintes. Les effets d’inertie, et des termes non-Fourier joueraient également un rôle non négligeable lorsque la longueur de la fissure est petite. Comme dans de nombreux problèmes d'ingénierie, l'initiation et la propagation des micro-fissures sont des mécanismes dont il faut tenir compte dans les prévisions de la rupture des structures. Ces effets non conventionnels ne sont plus négligeables et doivent être inclus dans l'analyse de la fracture des structures soumises à des chocs thermiques. / Thermal shock problems occur in many engineering materials and elements, which are used in high temperature applications such as thermal barrier coatings (TBCs), solid propellant of rocket-engine, pulsed-laser processing of materials, and so on. The thermal shock resistance performances and the thermal shock damages of materials, especially the interface debonding and spallation of coatings, have received considerable attention in both analysis and design. Some conventional theories, such as the Fourier’s law of thermal conduction and the quasi-static assumption of the thermoelastic body, may no longer be appropriate because of the extreme loads provoked by the thermal shock. Therefore, these conventional theories need to be enriched or revised.The objective of this thesis is to develop the solutions of the transient temperature field and thermal stresses around a partially insulated crack in a thermoelastic strip under thermal shock loading. The crack lies parallel to the heated traction free surface. The thermal conductivity of the crack gap is taken into account. Hyperbolic heat conduction theory is used in solving the temperature field instead of the traditional Fourier thermal conduction theory. Equations of motion are applied to obtain the stress fields and the dynamic stress intensity factors of the crack. The Laplace and Fourier transforms are applied to solve the thermal-elastic governing equations such that the mixed boundary value problems are reduced to solving a singular integral equations system in Laplace-Fourier space. The numerical integration method is applied to get the temperature field and stress fields, respectively. The problems are then solved numerically by converting the singular integral equations to a linear algebraic equations system. Finally, numerical inversions of the Laplace transform are performed to obtain the temperature field and dynamic stress intensity factors in the time domain.Numerical results show that the thermal conductivity of the crack gap strongly affects the uniformity of the temperature field and consequently, the magnitude of the dynamic stress intensity factors of the crack. The stress intensity factors would have higher amplitude and oscillating feature comparing to those obtained under the conventional Fourier thermal conduction and quasi-static hypotheses. It is also observed that the interactions of the thermal conductivity of the crack gap, the non-Fourier effect and the inertia effects would make the dynamic thermal shock problem more complex. The magnitude of the thermal barrier, non-Fourier and inertia effects is estimated for some practical cases.

Lorsqu’on fait face à des systèmes non linéaires, les notions classiques de stabilité ne suffisent pas à garantir un comportement approprié vis-à-vis de problématiques telles que le suivi de trajectoires, la synchronisation et la conception d’observateurs. La stabilité incrémentale a été proposée en tant qu’outil permettant de traiter de tels problèmes et de garantir que le système présente des comportements qualitatifs pertinents. Cependant, comme c’est souvent le cas avec les systèmes non linéaires, la complexité de l’analyse conduit les ingénieurs à rechercher des relaxations, ce qui introduit du conservatisme. Dans cette thèse, nous nous intéressons à la stabilité incrémentale d’une classe spécifique de systèmes, à savoir les systèmes affines par morceaux, qui pourraient fournir un outil avantageux pour aborder la stabilité incrémentale de systèmes dynamiques plus génériques.Les systèmes affines par morceaux ont un espace d’états partitionné, et sa dynamique dans chaque région est régie par une équation différentielle affine. Ils peuvent représenter des systèmes contenant des non linéarités affines par morceaux, ainsi que servir comme des approximations de systèmes non linéaires plus génériques. Ce qui est plus important, leur description est relativement proche de celle des systèmes linéaires, ce qui permet d’obtenir des conditions d’analyse exprimées comme des inégalités matricielles linéaires qui peuvent être traités numériquement de façon efficace par des solveurs existants.Dans la première partie de ce document de thèse, nous passons en revue la littérature sur l’analyse des systèmes affines par morceaux en utilisant des techniques de Lyapunov et la dissipativité. Nous proposons ensuite de nouvelles conditions pour l’analyse du gain L2 incrémental et la stabilité asymptotique incrémentale des systèmes affines par morceaux exprimés en tant qu’inégalités matricielles linéaires. Ces conditions sont montrées être moins conservatives que les résultats précédents et sont illustrées par des exemples numériques.Dans la deuxième partie, nous considérons le cas des systèmes affines par morceaux incertains représentés comme l’interconnexion entre un système nominal et un bloc d’incertitude structuré. En utilisant la théorie de la séparation des graphes, nous proposons des conditions qui étendent le cadre des contraintes quadratiques intégrales afin de considérer le cas où le système nominal est affine par morceaux, à la fois dans les cas non incrémental et incrémental. Via la théorie de la dissipativité, ces conditions sont ensuite exprimées en tant qu’inégalités matricielles linéaires.Finalement, la troisième partie de ce document de thèse est consacrée à l’analyse de systèmes non linéaires de Lur’e incertains. Nous développons une nouvelle technique d’approximation permettant de réécrire ces systèmes de façon équivalente comme des systèmes affines par morceaux incertains connectés avec l’erreur d’approximation. L’approche proposée garantit que l’erreur d’approximation est Lipschitz continue avec la garantie d’une borne supérieure prédéterminée sur la constante de Lipschitz. Cela nous permet d’utiliser les techniques susmentionnées pour analyser des classes plus génériques de systèmes non linéaires. / When dealing with nonlinear systems, regular notions of stability are not enough to ensure an appropriate behavior when dealing with problems such as tracking, synchronization and observer design. Incremental stability has been proposed as a tool to deal with such problems and ensure that the system presents relevant qualitative behavior. However, as it is often the case with nonlinear systems, the complexity of the analysis leads engineers to search for relaxations, which introduce conservatism. In this thesis, we focus on the incremental stability of a specific class of systems, namely piecewise-affine systems, which could provide a valuable tool for approaching the incremental stability of more general dynamical systems.Piecewise-affine systems have a partitioned state space, in each region of which the dynamics are governed by an affine differential equation. They can represent systems containing piecewise-affine nonlinearities, as well as serve as approximations of more general nonlinear systems. More importantly, their description is relatively close to that of linear systems, allowing us to obtain analysis conditions expressed as linear matrix inequalities that can be efficiently handled numerically by existing solvers.In the first part of this memoir, we review the literature on the analysis of piecewise-affine systems using Lyapunov and dissipativity techniques. We then propose new conditions for the analysis of incremental L2-gain and incremental asymptotic stability of piecewise-affine systems expressed as linear matrix inequalities. These conditions are shown to be less conservative than previous results and illustrated through numerical examples.In the second part, we consider the case of uncertain piecewise-affine systems represented as the interconnection between a nominal system and a structured uncertainty block. Using graph separation theory, we propose conditions that extend the framework of integral quadratic constraints to consider the case when the nominal system is piecewise affine, both in the non-incremental and incremental cases. Through dissipativity theory, these conditions are then expressed as linear matrix inequalities.Finally, the third part of this memoir is devoted to the analysis of uncertain Lur’e-type nonlinear systems. We develop a new approximation technique allowing to equivalently rewrite such systems as uncertain piecewise-affine systems connected with the approximation error. The proposed approach ensures that the approximation error is Lipschitz continuous with a guaranteed pre-specified upper bound on the Lipschitz constant. This enables us to use the aforementioned techniques to analyze more general classes of nonlinear systems.

Le principe du contrôle par flux de fuite magnétique (Magnetic Flux Leakage MFL) consiste à aimanter la pièce à contrôler par un champ magnétique et à détecter à l'aide d'un capteur magnétique les fuites des lignes du champ qui résultent de la présence d'un défaut dans la pièce. Dans le but d'améliorer les performances d'un dispositif de détection, le CEA et la société Vallourec collaborent pour développer un modèle numérique dédié au contrôle virtuel des défauts longitudinaux dans les tubes ferromagnétiques. Le dispositif expérimental comprend un circuit magnétique tournant à une vitesse constante autour du tube qui défile. Dans le cadre de cette thèse, on débute le problème de la modélisation sans tenir compte des effets de la vitesse de rotation, il s'agit donc de résoudre un problème d'électromagnétisme en régime magnétostatique.Pour résoudre ce problème, on propose de comparer une approche semi-analytique basée sur le formalisme des équations intégrales (EI) et une approche purement numérique utilisant les éléments finis (EF).Dans la première partie de cette thèse, après avoir établi le formalisme théorique par EI, un premier modèle considérant des matériaux ferromagnétiques à perméabilité magnétique constante (régime linéaire) a été mis en œuvre en 2D. Ce modèle a été appliqué pour un exemple de système extrait de la littérature et validé numériquement par une comparaison des résultats EI/EF. Pour une meilleure détection, il est opportun de saturer magnétiquement la pièce. Le matériau ferromagnétique est alors caractérisé par une courbe B(H) non-linéaire. Par conséquent, la deuxième partie de la thèse a été consacrée à la mise en œuvre du modèle en régime non linéaire qui tient compte de cette caractéristique.Différentes méthodes de discrétisation ont été étudiées afin de réduire le nombre d'inconnues et le temps de calcul. L'originalité de la thèse réside dans l'utilisation des fonctions d'interpolation d'ordre élevé (polynôme de Legendre) pour une discrétisation des équations intégrales par une approche de type Galerkin. Les premiers essais de validation numérique de ce modèle ont été effectués sur un système MFL simplifié. Des premiers essais de validation expérimentale pour des données obtenues par EF ont été effectués en deux phases : La première a consisté à vérifier le distribution du champ magnétique pour un tube sain et en régime magnétostatique. La deuxième phase a consisté à calculer la réponse d'un défaut dans le tube ferromagnétique en tenant en compte les effets éventuels de la rotation du circuit magnétique par rapport au tube. / The principle of the non destructive testing by magnetic flux leakage (MFL) is to magnetize the part to be inspected by a magnetic field and to detect a flaw thanks to magnetic leakage field lines due to the strong decreasing of the magnetic permeability in the flawed region. In order to improve the performance of detection, the CEA and the Vallourec society collaborate to develop a numerical model dedicated to the virtual NDT of longitudinal defects in ferromagnetic tubes. The experimental system includes a magnetic circuit rotating at a constant speed around the tube to be inspected. The modeling task is started without considering the effects of the rotational speed, so the magnetostatic regime is considered to solve the modeling problem. In the framework of this thesis, we propose to compare a semi-analytical approach based on the formalism of integral equations method (IEM) and a purely numerical approach using finite element method (FEM).In the first part of this thesis, the theoretical formalism was established. A first simple discretization scheme is been implemented in the linear regime considering a constant magnetic permeability. This first numerical model has been validated for a simplified MFL configuration extracted and modified from the literature.For better detection, it is wishable to magnetically saturate the piece under-test. The ferromagnetic material is then characterized by a B(H) curve. Therefore, the second part of the thesis was devoted to the implementation of the model in the non-linear regime that takes into account this non-linear characteristic. Different discretization schemes have been studied in order to reduce the number of unknowns and the computational time. The originality of the thesis lies in the use of basis function of high order (Legendre polynomials) associated to a Galerkin approach for the discretization of integral equations. The first numerical result has been validated on a simplified MFL system. The first results of the experimental validation based on simulated data obtain by FEM have been performed in two steps. The first one consists to verify the distribution of the magnetic field for a ferromagnetic tube without any defect and in the magnetostatic regime. The objective of the second one was to compute the response of the flaw and to evaluate the effects of the rotational speed of the magnetic circuit around the tube.

Cette thèse s'articule autour de deux thématiques : les méthodes numériques pour la propagation d'ondes acoustiques sous écoulement et les méthodes de réduction de modèles. Dans la première thématique, nous développons une méthode de couplage d'éléments finis et d'éléments de frontière pour résoudre l'équation d'Helmholtz convectée, lorsque l'écoulement est uniforme à l'extérieur d'un domaine borné. En particulier, nous proposons une formulation bien posée à toutes les fréquences de la source. Dans la deuxième thématique, nous proposons une solution au problème classique d'accumulation d'arrondis machine qui survient en calculant l'estimateur d'erreur a posteriori dans la méthode des bases réduites. Par ailleurs, nous proposons une méthode non intrusive pour calculer une approximation sous forme séparée des systèmes linéaires résultant de l'approximation en dimension finie de problèmes aux limites dépendant d'un ou plusieurs paramètres / This thesis has two topics : numerical methods for acoustic wave propagation in a flow and reduced order models. In the first topic, we develop a coupled finite element and boundary element method to solve the convected Helmholtz equation, when the flow is uniform outside a bounded domain. In particular, we propose a formulation that is well-posed at all the frequencies of the source. In the second topic, we propose a solution to the classical problem of round-off error accumulation that occurs when computing the a posteriori error bound in the reduced basis method. Furthermore, we propose a non intrusive method for the approximation, in a separated representation form, of linear systems resulting from the finite-dimensional approximation of boundary-value problems depending on one or several parameters

De nos jours, il n'existe encore aucune théorie capable de proposer une description précise et quantitative du transfert de proton en solution. En effet, ce problème est complexe du fait de la grande diversité des interactions existant dans l'eau liquide, à savoir: des interactions non liantes de type Van der Waals, des liaisons faiblement covalentes et des liaisons hydrogènes remarquablement fortes. Ces dernières sont d'ailleurs à l'origine des nombreuses propriétés fascinantes de l'eau à l'échelle macroscopique. À cela s'ajoutent les effets quantiques nucléaires dus à la faible masse de l'hydrogène, qui modifient profondément la nature de la surface d'énergie potentielle décrivant le transfert de proton le long de sa coordonnée de réaction. Nous proposons dans cette thèse une approche tout quantique basée sur une description quasi exacte de la fonction d'onde du système par l'utilisation de méthodes stochastiques de type Monte Carlo Quantique. Cette technique, combinée avec le formalisme des équations de Langevin et des intégrales de chemin de Feynman, permet de simuler à un niveau de précision inédit, n'importe quel système chimique en phase gaz ou en solution. Nous appliquons cette méthodologie à des agrégats d'eau neutres ou protonés pour apporter de nouveaux éclaircissements sur les phénomènes microscopiques régissant la diffusion du proton hydraté dans de tels systèmes. Il est mis en évidence que la mobilité du proton est optimale pour des températures proches des conditions ambiantes, du fait de la compétition subtile entre les effets thermiques et quantiques nucléaires. / There is no theory up to now able to provide an accurate and quantitative description of the proton transfer (PT) yet. Indeed, the complexity of the problem stems from the large diversity of the existing interactions in liquid water, namely: non bonding Van der Waals interactions, weakly covalent bonds and remarkably strong H-bonds. The latter ones are at the origin of the numerous fascinating properties of water at the macroscopic scale. In addition to such interactions, the nuclear quantum effects arising from the hydrogen light mass deeply modify the potential energy surface, and must be taken into account. In this thesis, we propose a fully quantum approach based on an almost exact description of the electronic wave function by means of Quantum Monte Carlo (QMC) methods. Our novel technique combines QMC with a Langevin-based Molecular Dynamics and the Feynman's path integral formalism. This allows one to perform fully quantum simulations of systems in gas or condensed phase, at an unprecedented level of accuracy,. We apply our approach to neutral or charged protonated water clusters to shed light on the microscopic phenomena driving the proton diffusion in such systems. We discovered that the proton hopping is optimal for temperatures close to ambient conditions, due to the subtle competition between thermal and nuclear quantum effects. This is highly suggestive of the importance of quantum nuclear effects to make PT processes - relevant for life - most efficient at room temperature.

Ce travail porte sur l'enseignement et l'apprentissage des intégrales multiples pour les calculs d'aires et de volumes, dans des universités brésiliennes et les classes préparatoires technologiques françaises. L'objectif est de mieux comprendre les difficultés rencontrées par les étudiants, avec ces calculs et d'étudier dans quelle mesure l'utilisation d'un logiciel comme Maple peut les aider à surmonter ces difficultés, en favorisant les interactions entre représentation graphique et représentation analytique des objets concernés. <br />En nous appuyant sur l'approche instrumentale de l'intégration des outils informatiques à l'enseignement, nous partons de l'hypothèse que la réussite de ce type de tâche mathématique avec Maple nécessite, au-delà de la connaissance des commandes pertinentes et de leur syntaxe, l'apprentissage de techniques instrumentées spécifiques, un apprentissage qui doit être pris en charge par l'institution. Dans la thèse, une technique de représentation dite de crible géométrique occupe de ce point de vue une place centrale, en étant vue comme un moyen de forcer l'interaction entre représentations graphique et analytique.<br /><br />Le travail comprend une analyse institutionnelle (étude écologique de programmes et de manuels) du point de vue de l'enseignement des intégrales multiples et de celui de l'intégration de l'outil Maple, et deux expérimentations sans et avec Maple.

On s'intéresse à la géométrie des variétés de courbure de Ricci presque supérieure à $k$ (i.e. telle qu'une norme $L^p$---locale ou globale---de la fonction $(\underline(\rm Ric)-k)^-$ soit petite, où $\underline(\rm Ric)(x)$ est la plus petite valeur propre de la courbure de Ricci en $x$). On démontre sous cette hypothèse les équivalents des inégalités géométriques classiques de Myers, de Bishop-Gromov, de Lichnerowicz,... puis on caractérise les variétés qui réalisent presque les cas d'égalité (généralisant des travaux de T.~Colding et de P.~Petersen). Sur une variété compacte $M^n$ de courbure presque positive, le laplacien sur les 1-formes a au plus $n$ petites valeurs propres. S'il a exactement $n$ petites valeurs propres ($n-1$ suffisent si $M$ est orientable) alors $M$ est difféomorphe à une Nilvariété et la métrique est presque invariante à gauche. Ces résultats découlent d'estimées analytiques établies dans la première partie de la thèse.