On s'intéresse dans cette thèse à l'analyse haute fréquence de l'équation de Helmholtz dans un cadre dissipatif. On commence par chercher des estimations uniformes pour la résolvante de l'opérateur de Schrödinger dissipatif sur le demi-plan supérieur et près d'une énergie vérifiant une hypothèse d'amortissement sur les trajectoires classiques bornées. On généralise pour cela la méthode des commutateurs de Mourre pour des opérateurs dissipatifs. Dans une deuxième partie, on étudie les mesures semi-classiques pour la solution sortante à l'équation lorsque le terme source se concentre sur une sous-variété bornée de l'espace.

Dans la première partie de cette thèse, une simulation numérique 2D de la formation des structures chevauchantes est proposée. Elle repose sur les hypothèses des coupes équilibrées (conservation du volume dans le plan de coupe parallèle à la direction de transport) avec un plissement isopaque. Chaque coupe équilibrée est construite directement à partir de la géométrie initiale des chevauchements, du raccourcissement et de la cinématique des chevauchement. Ces paramètres initiaux du calcul sont ensuite modifiés par essai et erreur de façon à ajuster la géométrie finale calculée sur les observations géologiques. L'application informatique de ces hypothèses a été rendue possible par l'écriture d'un programme appelé CICERON. Ce modèle permet de simuler des déplacements importants le long de discontinuités cinématiques. Il a été utilisé sur deux types de structures impliquant: 1) seulement la couverture sédimentaire (Jura, Appalaches, prisme de la Barbade); 2) à la fois le socle et la couverture sédimentaire (Appenin du Sud). Dans la deuxième partie de cette thèse, l'évolution de la température pendant la formation d'une structure chevauchante est simulée numériquement grâce avec la méthode des éléments finis. A chaque incrément de déplacement sur le chevauchement, le maillage calculé par le mailleur MODULEF s'appuie sur les géométries fournies par le programme CICERON. L'étude détaillée des résultats de cette modélisation thermique montre que les transferts horizontaux conductifs de chaleur sont importants et accélèrent la mise à l'équilibre thermique de telles structures.

Un nouvel point de vue est proposé pour la simulation des plasmas utilisant le modèle cinétique qui couple les équations de Vlasov pour la distribution des particules et les équations de Maxwell pour la contribution des champs électromagnétique. On part du principe que les équations de la Physique sont des objets mathématiques qui mettent en relation des objets géométriques. Afin de conserver les propriètés géométriques des différents objets intervenant dans une équation, on utilise, pour l'étude théorique et numérique, la géométrie différentielle. Il s'avère que toutes les équations de la Physique peuvent s'écrire à l'aide des formes différentielles et que sous ce point de vue celles-ci sont indépendantes du choix des coordonnées. On propose alors une discrétisation des formes différentielles en utilisant les B-splines comme fonctions d'interpolation. Afin d'être cohérent avec la théorie, on proposera également une discrétisation des différentes opérations de la géométrie différentielle agissant sur les formes différentielles. On teste notre schéma tout d'abord sur les équations de Maxwell avec plusieurs conditions aux bords et puisque ce schéma numérique obtenu est indépendant du système de coordonnées, on le teste également lorsque l'on effectue un changement de coordonnées. Enfin, on applique la même méthode sur les équations de Vlasov-Poisson 1D et on propose plusieurs schémas numériques.

La thèse est consacrée essentiellement à l'étude de systèmes non linéaires d'évolution issus d'un modèle de Boussinesq : couplage entre les équations de Navier-stokes avec un second membre F(µ), où F est une force de gravité proportionnelle à des variations de densité qui dépendent de la température et l'équation de l'énergie.<br />Le premier chapitre nous donne un résultat d'existence d'une solution faible-renormalisée du système de Boussinesq en dimension 2, dans le cas où F est bornée.<br />Dans le chapitre 2, on aborde le cas de fonctions F plus générales : F vérifie une hypothèse de croissance. On démontre l'existence de solutions pour toutes données initiales ou pour des données initiales petites selon la croissance de F.<br />Dans le chapitre 3, nous faisons une généralisation des résultats du chapitre 2 mais sans le terme de convection.<br />Dans le chapitre 4, le manque de stabilité de l'énergie de dissipation dans L1(Q) en dimension 3, nous contraint à transformer de façon formelle le système de Boussinesq. On démontre l'existence d'une solution faible de ce nouveau système en dimension 3.

Ce travail cherche a comprendre comment l'ajout d'une géométrie non euclidienne dans un problème de Schrödinger non linéaire influe sur l'existence et l'unicité des solutions explosives de masse critique. On s'inspire pour beaucoup des travaux de Merle et Raphaël sur la méthode de modulation des paramètres d'invariance géométrique pour une EDP qui possède de bonnes lois de conservations. On s'appuie ici plus particulièrement sur un article de Raphaël et Szeftel qui prouve l'existence et l'unicité d'une solution de masse critique en dimension 2 pour l'équation de Schrödinger non linéaire avec potentiel d'inhomogénéité devant la non-linéarité, et qui explose par ailleurs au maximum de l'inhomogénéité. Dans un premier temps, il s'agit de reprendre la méthode dans son ensemble afin de l'adapter à des cas où le Laplacien n'est plus plat, et est remplacé par un opérateur de type Laplace-Beltrami ou Laplacien généralisé. Ayant mis en avant le rôle de la courbure au point d'explosion, en termes de conditions sur les dérivées de termes métriques, on reprend dans un deuxième temps l'étude dans le cas plus général d'une variété riemannienne. Grâce à un ansatz sur la solution qui intègre maintenant la transformation induite par la métrique, on est capable d'énoncer un résultat d'existence et d'unicité en termes de conditions géométriques sur la variété elle même. Par soucis de simplicité, on se limite néanmoins au rôle local de la métrique, en la supposant globalement définie dans une certaine carte, et asymptotiquement équivalente a la métrique euclidienne.

La connaissance des propriétés critiques des mélanges binaires est importante pour la conception et le développement des procédés des industries chimiques, mais également, pour l'amélioration et l'extension des modèles thermodynamiques prédictifs. Cependant, des données expérimentales critiques liquide - vapeur de mélanges binaires sont relativement rare dans la littérature. C'est la raison pour laquelle, nous avons acquis un nouvel appareil permettant de mesurer les points critiques de corps purs et de mélanges. Dans une première étape, nous avons mis au point et validé le protocole expérimental, en comparant les résultats de nos mesures aux données de la littérature. Dans une deuxième étape, nous montrons l'intérêt que la mesure de points critiques de mélanges binaires peut présenter pour le développement de deux types de procédés : i) la production de biodiesel par voie supercritique et ii) l'étude de comportement de phases de nouveaux réfrigérants. Les résultats expérimentaux obtenus ont été également utilisés pour déterminer les paramètres d'interactions entre groupes du modèle PPR78 / Knowledge of mixture critical data is not only important for the design and the improvement of chemical processes but also to extend the range of applicability of predictive equations of state. However, mixture critical data comprise only a small percentage of the fluid phase equilibria data available in the open literature. This is the reason why it was decided to acquire a new synthetic-dynamic apparatus capable of determining the critical points of pure substances and multi-component mixtures. In a first step, the validation of the experimental procedure for critical point measurements was performed by comparing the experimental results obtained in this work with literature data. In a second step, the application of mixture critical data to the development of two types of process: i) biodiesel production via supercritical process and ii) phase behaviour of new refrigerants is presented. The experimental critical points measured in this work were also used for determining new group interaction parameters of the PPR78 model

L’enjeu de la simulation de la dynamique atmosphérique et océanographique a pris ces dernières années une importance accrue avec la question du réchauffement climatique. Le modèle à simuler est complexe. Il combine les équations de la mécanique des fluides avec celles de la thermodynamique. Au 19ème siècle, le mathématicien Adhémar Barré de Saint-Venant formule un système d’équations aux dérivées partielles décrivant les mouvements d’un fluide soumis à la gravité et de faible épaisseur. Il s’agit des équations Shallow Water. L’objectif de cette thèse est de développer et d’analyser un algorithme de résolution des équations Shallow Water sur une sphère en rotation. Dans un premier temps, j’étudie différents aspects mathématiques des opérateurs aux différences finis utilisés par la suite en géométrie sphérique. Les schémas aux différences obtenus sont utilisés pour résoudre l’équation de transport, l’équation des ondes et l’équation de Burgers. Les propriétés de stabilité précision et conservation sont analysées. Dans un second temps, la grille Cubed-Sphere est introduite et analysée. La structure de ce maillage est analogue à celle d’un cube. L'interprétation de la Cubed-Sphere à l’aide de grands cercles permet de construire des opérateurs sphériques discrets gradient, divergence et vorticité d'ordre au moins égal à 3 (en pratique d'ordre 4). La troisième partie de la thèse est dédiée à différents tests pour le système d’équations Shallow Water ainsi que pour l’équation d’advection. Les résultats démontrent une précision proche de celle obtenue par les algorithmes conservatifs d'ordre 4 les plus récents / The problem to obtain accurate simulations of the atmospheric and oceanic equations has become essential in recent years for a proper understanding of the climate change. The full mathematical model to simulate is rather complex. It consists of the coupling of several equations involving fluid dynamics and thermodynamics. In the 19th century, Adhémar Barré de Saint-Venant first formulated the equations describing the dynamic of a fluid subject to gravity and bottom topography. This system is Shallow Water equations. The goal of this thesis is to develop and analyze a numerical scheme to solve the shallow water equation on a rotating sphere. First, a mathematical analsysis of finite difference operators that will be used on the sphere is presented. These schemes are then used to solve various equations in a spehreical setting, in particular the advection equation, the wave equation and the Burgers equation. Stability, accuracy and conservation properties are studied. In a second part, I consider in detail the Cubed-Sphere grid. This particular spherical grid has the mesh topology of a cube. Another interpretation makes use of great circles, this allows to obtain spherical discret operators gradient, divergence and curl of a preved third order. These operators are numercially of 4th order. Numerial results are show in particular for the SW equations an acurracy similar to the one of conservative schemes of 4th order published recently

Une couche limite turbulente se développant le long d’une paroi présente des fluctuationsde vitesses et de pression importantes. Si la paroi du profil est suffisammentsouple, les fluctuations de pression pariétale peuvent la faire rentrer en vibration cequi induit un rayonnement acoustique de chaque côté de la paroi. Ce scenario estl’un des mécanismes de génération de bruit interne dans les aéronefs. Le but de cettethèse est de proposer un modèle de reconstruction des fluctuations de pression pariétaleafin de prévoir in fine le bruit rayonné. Plutôt que de reposer sur une approchesemi-empirique, les modèles développés dans cette thèse se basent sur la résolutionanalytique de l’équation de Poisson liant les fluctuations de pression aux fluctuationsde vitesses. Ces dernières sont modélisées par exemple à l’aide des profils moyensde la couche limite obtenus grâce à un calcul RANS. La résolution de l’équation dePoisson dans ce contexte a déjà été entreprise en particulier par Lysak et Aupoixet leurs travaux sont le point de départ de cette thèse. Cependant, leur modèle nedonne qu’une description temporelle des fluctuations de pression pariétale alors queles aspects spatiaux sont nécessaires pour une application vibro-acoustique. L’apportde cette thèse consiste donc en une modification de leur modèle afin de palliercette difficulté. En parallèle de ces travaux de modélisation, une expérience de validationen soufflerie a été élaborée et mise en place. Les fluctuations de vitesses ontété mesurées par vélocimétrie laser tandis que les fluctuations de pression pariétaleont été mesurées à l’aide de micro-tiges mobiles. Le modèle initialement développéà été affiné à l’aide de ces mesures. En particulier, une description anisotrope desfluctuations de vitesses a été développée, ce qui est plus cohérent pour un écoulementcisaillé que la description homogène isotrope utilisée jusqu’alors. Les modèlesdéveloppés ont un large recoupement avec le modèle semi-empirique de Corcos quiest la référence utilisée pour les applications en vibro-acoustique. Cependant, desdifférences comportementales importantes aux hautes et basses fréquences ont étémises en évidence. Le modèle de Corcos peut donc être remis en question pour cesplages fréquentielles. Ces résultats théoriques doivent néanmoins être confortés pardes mesures. / Large pressure and velocity fluctuations are present inside a turbulent boundarylayer developing along a wall. A non-rigid wall can be excited by the wall pressurefluctuations and thus acoustic radiations will be emitted above and bellow the wall.This scenario is one the mechanism of intern noise generation inside aircraft. Theaim of this thesis is to elaborate a turbulent wall pressure fluctuations model inorder to compute the noise radiated by the vibrating wall at the end. Modelsdeveloped during this thesis do not rely on an semi-empirical approach as they arebased on the analytical resolution of Poisson’s equation relating pressure and velocityfluctuations. This kind of approach has already been used by Lysak and Aupoix andtheir work was the starting point of this thesis. However, their model only gives atemporal description of the wall pressure fluctuations while a temporal and spatialdescription is needed for vibro-acoustics application. The major contribution of thisthesis is to modify their model in order to overcome this incapacity. In parallelto this theoretical modelling, a wind-tunnel experimental campaign dedicated tovalidation was designed and implemented. Velocity fluctuations were measured usingLaser Doppler Velocimetry while wall pressure fluctuations were measured usingmobile wall-mounted microphones. The initial model was improved using thesemeasurements. In particular, an anisotropic description of the velocity fluctuationswas developed which is more consistent for a sheared flow than the homogeneous andisotropic description used by Lysak and Aupoix. For a large frequency range, thefinal model behaves similarly to Corcos’ model which is the most used reference forvibro-acoustics applications. However, large differences were highlighted for low andhigh frequencies. Therefore the validity of Corcos’ model is questionable for thesefrequency range. These theoretical results must still be confirmed by experimentaldata.

La propagation d'une fissure chargée par un écoulement de fluide visqueux est un phénomène complexe où la compréhension des phénomènes mécaniques mis en jeu en pointe de fissures reste encore partielle. C'est le cas de la zone de décollement entre le solide et le fluide qui apparaît pour un certain choix de débit d'injection, de viscosité du fluide et de ténacité du matériau. Cette thèse propose une modélisation simplifiée de ce problème d'interaction fortement couplé. Dans un premier chapitre, on étudie un modèle simplifié unidimensionnel de film élastique collé sur un substrat rigide et on considère une injection de fluide visqueux entre le film et le substrat. On suppose que la propagation de la fissure est régie par la loi de Griffith. On néglige l'existence du retard possible entre le fluide et le solide et on choisit une loi de comportement non-linéaire pour le fluide visqueux. A partir d'une analyse asymptotique pour une faible viscosité, on établit une solution approchée du problème. On montre que le champ pression est singulier en pointe de fissure et on montre l'influence du débit d'injection sur la cinétique du trajet de fissuration. Dans le deuxième chapitre on propose de prendre en compte l'existence de la zone de décollement en modifiant la formulation du modèle et en le réécrivant sous la forme d'un problème d'optimisation en temps discret où les zones de décollement font partie des inconnues du problème. On valide la formulation proposée sur l'exemple analytique de l'écrasement d'une goutte par une barre rigide. On montre ensuite que cette formulation et l'algorithme lié à son implémentation sont capables de gérer l'évolution de l'écrasement de plusieurs gouttes de forme quelconque en capturant correctement les phase d'étalement des gouttes ainsi que de leur coalescence. On étend ensuite cette formulation au cas de l'écrasement d'une goutte par un film élastique. Dans le dernier chapitre, on examine la validité de l'hypothèse de lubrification utilisée en fracturation hydraulique. A l'aide de la méthode de développement asymptotique, on construit une équation de Reynolds régularisée avec des termes de gradient supérieur tenant compte de la variation spatiale de la hauteur des parois. On compare alors le comportement des champs de pression donnés par les équations de Reynolds classique et régularisée sur des exemples d'écoulement entre des conduits de formes multiples. / The crack evolution under a viscous fluid action is a complex phenomenon where the understanding of the mechanical phenomena near the crack tip is still largely limited. This is the case for the lag between the solid and the fluid front propagation which appears for some configurations of injection rate, fluid viscosity and material toughness. This thesis proposes a simplified model for this strongly coupled interaction problem.The first chapter studies a simplified one-dimensional model of a elastic film bonded to a rigid substrate. We consider a viscous fluid injection between the film and the substrate. The crack propagation is assumed to follow the Griffith's law. The existence of the lag is neglected and a non-linear behavior law is chosen for the viscous fluid. Using an asymptotic analysis, an approximate solution is established for the low viscosity case. It is shown that the pressure field diverges at the crack tip and that the kinetics of the crack is influenced by the injection rate. The second chapter proposes to take into account the existence of the lag by modifying the model formulation and rewriting it as a discrete time optimisation problem where the delamination zones are part of the unknowns of the problem. This formulation is validated for the analytical example of a drop crushed by a rigid bar. It is shown that this formulation and its implementation can manage the evolution of several drops of any shape and correctly captures the drops spreading and coalescence. This formulation is then extended to the case of a drop crushed by an elastic film. In the last chapter, the validity of the lubrication hypothesis is examinated. Using an asymptotic analysis, a regularized Reynolds equation is constructed with higher gradient terms taking into account the spatial variation of the walls height. A comparison between the pressure fields behaviour given by the classical and the regularized Reynolds equation is shown for different conducts.

La mise en œuvre de la méthode des éléments finis de frontière nécessite l'évaluation d'intégrales comportant un intégrand singulier. Un calcul fiable et précis de ces intégrales peut dans certains cas se révéler à la fois crucial et difficile. La méthode que nous proposons consiste en une réduction récursive de la dimension du domaine d'intégration et aboutit à une représentation de l'intégrale sous la forme d'une combinaison linéaire d'intégrales mono-dimensionnelles dont l'intégrand est régulier et qui peuvent s'évaluer numériquement mais aussi explicitement. L'équation de Helmholtz 3-D sert d'équation modèle mais ces résultats peuvent être utilisés pour les équations de Laplace et de Maxwell 3-D. L'intégrand est décomposé en une partie homogène et une partie régulière ; cette dernière peut être traitée par les méthodes usuelles d'intégration numérique. Pour la discrétisation du domaine, des triangles plans sont utilisés ; par conséquent, nous évaluons des intégrales sur le produit de deux triangles. La technique que nous avons développée nécessite de distinguer entre diverses configurations géométriques ; c'est pourquoi nous traitons séparément le cas de triangles coplanaires, dans des plans sécants ou parallèles. Divers prolongements significatifs de la méthode sont présentés : son extension à l'électromagnétisme, l'évaluation de l'intégrale du noyau de Green complet pour les coefficients d'auto-influence, et le calcul de la partie finie d'intégrales hypersingulières. / The implementation of the boundary element method requires the evaluation of integrals with a singular integrand. A reliable and accurate calculation of these integrals can in some cases be crucial and difficult. The proposed method is a recursive reduction of the dimension of the integration domain and leads to a representation of the integral as a linear combination of one-dimensional integrals whose integrand is regular and that can be evaluated numerically and even explicitly. The 3-D Helmholtz equation is used as a model equation, but these results can be used for the Laplace and the Maxwell equations in 3-D. The integrand is decomposed into a homogeneous part and a regular part, the latter can be treated by conventional numerical integration methods. For the discretization of the domain we use planar triangles, so we evaluate integrals over the product of two triangles. The technique we have developped requires to distinguish between several geometric configurations, that's why we treat separately the case of triangles in the same plane, in secant planes and in parallel planes.