Tato dizertační práce je zaměřena na studium deformace nerovnosti uvnitř elastohydrodynamicky mazaného (EHD) kruhového kontaktu. Práce se zabývá studiem přechodu příčné nerovnosti přes kontaktní oblast, která je modelována pomocí numerických metod. Model dále uvažuje nenewtonské chování maziva. Použitý matematický model se skládá z parciální diferenciální rovnice druhého řádu pro řešení tlaku a integro-diferenciální rovnice pro řešení elastických deformací. Pro řešení tohoto modelu je použitá takzvaná multigrid (vícesíťová) metoda. Práce obsahuje popis matematického modelu EHD kontaktu a aplikované numeriké metody. Výsledky simulací jsou porovnány s experimentálně stanovenýma hodnotama tloušťky mazacího filmu. Deformace nerovnosti uvnitř kontaktní oblasti je studována pro různé provozní podmínky (střední rychlost, poměr proklzu) a různá vlastnosti maziva.

Nous nous intéressons, dans cette thèse, à l'étude des méthodes d'extrapolation polynômiales et à l'application de ces méthodes dans l'accélération de méthodes de points fixes pour des problèmes donnés. L'avantage de ces méthodes d'extrapolation est qu'elles utilisent uniquement une suite de vecteurs qui n'est pas forcément convergente, ou qui converge très lentement pour créer une nouvelle suite pouvant admettreune convergence quadratique. Le développement de méthodes cycliques permet, deplus, de limiter le coût de calculs et de stockage. Nous appliquons ces méthodes à la résolution des équations de Navier-Stokes stationnaires et incompressibles, à la résolution de la formulation Kohn-Sham de l'équation de Schrödinger et à la résolution d'équations elliptiques utilisant des méthodes multigrilles. Dans tous les cas, l'efficacité des méthodes d'extrapolation a été montrée.Nous montrons que lorsqu'elles sont appliquées à la résolution de systèmes linéaires, les méthodes d'extrapolation sont comparables aux méthodes de sous espaces de Krylov. En particulier, nous montrons l'équivalence entre la méthode MMPE et CMRH. Nous nous intéressons enfin, à la parallélisation de la méthode CMRH sur des processeurs à mémoire distribuée et à la recherche de préconditionneurs efficaces pour cette même méthode. / In this thesis, we study polynomial extrapolation methods. We discuss the design and implementation of these methods for computing solutions of fixed point methods. Extrapolation methods transform the original sequance into another sequence that converges to the same limit faster than the original one without having explicit knowledge of the sequence generator. Restarted methods permit to keep the storage requirement and the average of computational cost low. We apply these methods for computing steady state solutions of incompressible flow problems modelled by the Navier-Stokes equations, for solving the Schrödinger equation using the Kohn-Sham formulation and for solving elliptic equations using multigrid methods. In all cases, vector extrapolation methods have a useful role to play. We show that, when applied to linearly generated vector sequences, extrapolation methods are related to Krylov subspace methods. For example, we show that the MMPE approach is mathematically equivalent to CMRH method. We present an implementation of the CMRH iterative method suitable for parallel architectures with distributed memory. Finally, we present a preconditioned CMRH method.

L'objectif de ce travail est d'étudier la dé<br />formation d'un tissu posé sur un support bi- ou tri-dimensionnel et soumis à<br />son propre poids.<br />Dans la première partie, nous établissons les équations<br />d'équilibre de ce problème dans le cas général et<br />introduisons deux modèles mathématiques. Le premier est un<br />modèle membranaire non-linéaire, dont l'analyse mathématique<br />conduit au calcul de l'enveloppe quasi-convexe de la densité<br />d'énergie associée. Le deuxième modèle (modèle<br />membrane-flexion non-linéaire) est obtenu en ajoutant un terme<br />régularisant à une fonctionnelle énergie non coercive. Nous<br />prouvons l'existence d'au moins une solution de ce problème de<br />minimisation, en utilisant les techniques du Calcul des Variations. Enfin,<br />nous établissons l'existence de solutions pour le problème de<br />drapé tri-dimensionnel.<br />La seconde partie est consacrée à la résolution numérique des diffé%<br />rents modèles élaborés dans la première partie, au moyen d'une méthode ité%<br />rative de descente couplée avec une méthode multigrille, afin d'accélérer la<br />convergence de l'algorithme. Nous montrons que le problème discret admet au<br />moins une solution. Enfin, nous prouvons la convergence théorique d'une<br />sous-suite de solutions discrètes vers une solution du problème continu,<br />moyennant une hypothèse de densité.

Les polymères jouent un rôle important sur la réduction de la traînée et la modification de la structure des écoulements. Nous avons utilisé le modèle Oldroyd-B pour étudier l’effet de la viscoélasticité de solutions de polymères dilués sur des écoulements en deux dimensions autour d’obstacles dans un canal. Les obstacles sont pris en compte par la méthode de pénalisation volumique et des condition aux limites artificielles sans réflexion sont imposées à la sortie du canal. La discrétisation est effectuée par des schémas performants en différences finies et la résolution par une méthode multigrille. Les simulations numériques sont effectuées pour une large gamme de nombres de Reynolds et de nombres de Weissenberg. Les caractéristiques détaillées des écoulements viscoélastiques sont analysées et comparées entre elles et celles de l’écoulement du fluide sans polymère. En particulier les jeux de paramètres conduisant à une augmentation ou à une baisse de la traînée. Enfin les effets du polymère sur des écoulements turbulents sont aussi analysés. / Polymer plays an important role on the drag reduction and modification of the structure of flow. In this thesis, Oldroyd B model is employed to study the effectof viscoelasticity for polymer solutions diluted by two-dimensional direct numerical simulation. The obstacles are taken into account by penalization method. The artificial boundary condition is imposed without any reflection on the channel outlet. Flow pasta cylinder is investigated in detailed by present numerical methods. The numerical codes are valid for predicting the drag force and capturing the important character of viscoelastic flow by comparing with experimental and other numerical results. The drag map of the cylinder is obtained at a wide range of Reynolds number and Weissenberg number space. The detailed characteristics of viscoelastic flow are reported in the thesis. The effects of polymer on two-dimensional turbulent flow are also discussed by grid turbulent flow.