La prévention des risques industriels nécessite de simuler la dispersion turbulente de polluants. Cependant, les outils majoritairement utilisés à ce jour ne permettent pas de traiter les champs proches dans le cas de géométries complexes, et il est nécessaire d'utiliser les outils de CFD (“ Computational Fluid Dynamics ”) plus adaptés, mais plus coûteux. Afin de simuler les écoulements atmosphériques avec dispersion de polluants, les modèles CFD doivent modéliser correctement d'une part, les effets de flottabilité, et d'autre part les effets de la turbulence. Plusieurs approches existent, notamment dans la prise en compte des effets de flottabilité et la modélisation de la turbulence, et nécessitent des méthodes numériques adaptées aux spécificités mathématiques de chacune d'entre elles, ainsi que des schémas numériques précis pour ne pas polluer la modélisation. Une formulation d'ordre élevé en volumes finis, sur maillages non structurés, parallélisée, est proposée pour simuler les écoulements atmosphériques avec dispersion de polluants. L'utilisation de schémas d'ordre élevé doit permettre d'une part de réduire le nombre de cellules et diminuer les temps de simulation pour atteindre une précision donnée, et d'autre part de mieux contrôler la viscosité numérique des schémas en vue de simulations LES (Large Eddy Simulation), pour lesquelles la viscosité numérique des schémas peut masquer les effets de la modélisation. Deux schémas d'ordre élevé ont été étudiés et implémentés dans un solveur 3D Navier Stokes incompressible sur des maillages volumes finis non structurés. Nous avons développé un premier schéma d'ordre élevé, correspondant à un schéma Padé volumes finis, et nous avons étendu le schéma de reconstruction polynomiale de Carpentier (2000) aux écoulements incompressibles. Les propriétés numériques des différents schémas implémentés dans le même code de calcul sont étudiées sur différents cas tests bi-dimensionnels (calcul de flux convectifs et diffusifs sur une solution a-priori, convection d'une tâche gaussienne, décroissance d'un vortex de Taylor et cavité entraînée) et tri-dimensionnel (écoulement autour d'un obstacle cubique). Une attention particulière a été portée à l'étude de la précision et du traitement des conditions limites. L'implémentation proposée du schéma polynomial permet d'approcher, pour un maillage identique, les temps de simulation obtenus avec un schéma décentré classique d'ordre 2, mais avec une précision supérieure. Le schéma compact donne la meilleure précision. En utilisant une méthode de Jacobi sans calcul implicite de la matrice pour calculer le gradient, le temps de simulation devient intéressant uniquement lorsque la précision requise est importante. Une alternative est la résolution du système linéaire par une méthode multigrille algébrique. Cette méthode diminue considérablement le temps de calcul du gradient et le schéma Padé devient performant même pour des maillages grossiers. Enfin, pour réduire les temps de simulation, la parallélisation des schémas d'ordre élevé est réalisée par une décomposition en sous domaines. L'assemblage des flux s'effectue naturellement et différents solveurs proposés par les librairies PETSC et HYPRE (solveur multigrille algébrique et méthode de Krylov préconditionnée) permettent de résoudre les systèmes linéaires issus de notre problème. Le travail réalisé a consisté à identifier et déterminer les paramètres de résolution qui conduisent aux temps de simulation les plus faibles. Différents tests de speed-up et de scale-up ont permis de déterminer la méthode la plus efficace et ses paramètres optimaux pour la résolution en parallèle des systèmes linéaires issus de notre problème. Les résultats de ce travail ont fait l'objet d'une communication dans un congrès international “ parallel CFD juin 2008 ” et d'un article soumis à “ International Journal for Numerical Methods in Fluids ” (Analysis of high-order finite volume schemes for the incompressible Navier Stokes equations)

La thèse porte sur le développement d'algorithmes d'estimation rapides pour l'analyse modale de signaux multidimensionnels (R-D) présentant des problèmes de résolution et de complexité numérique. Un signal multidimensionnel de dimension R est la superposition de produits de R sinusoïdes. L'application visée est la spectroscopie RMN. Dans un premier temps, après un état de l'art des méthodes d'estimation dites " algébriques ", nous proposons une méthode paramétrique basée sur les tenseurs. Celle-ci utilise le treillis multidimensionnel du tenseur du signal R-D et exploite la structure des vecteurs propres du sous-espace signal obtenus en utilisant la décomposition en valeurs singulières d'ordre supérieur. Contrairement à la plupart des approches tensorielles, la méthode proposée permet d'éviter la phase d'appariement des coordonnées des modes dans chacune des dimensions ou d'une diagonalisation conjointe. Dans un deuxième temps, le problème d'estimation modale multidimensionnelle est présenté comme un problème d'approximation parcimonieuse dans lequel le dictionnaire est obtenu par la discrétisation de fonctions exponentielles complexes. Afin d'atteindre une bonne résolution spectrale, il est nécessaire de choisir une grille très fine, ce qui conduit à la manipulation d'un dictionnaire de grande taille avec tous les problèmes calculatoires sous-jacents. Nous proposons alors une méthode originale qui consiste à combiner une approximation parcimonieuse et une approche multigrille sur plusieurs niveaux de résolution. L'approche est validée au travers de plusieurs exemples 1-D et 2-D. En outre, une étude sur l'influence du choix du dictionnaire initial sur la convergence est également menée. Les méthodes développées sont ensuite appliquées à l'estimation des paramètres de signaux de spectroscopie RMN 1-D et 2-D. Afin de réduire le coût de calcul dans le cas de signaux bidimensionnels de grande taille, nous proposons également une approche exploitant la notion de parcimonie simultanée, pour estimer les coordonnées des modes sur chacune des dimensions. La procédure consiste à effectuer deux approximations parcimonieuses 1-D suivies d'une phase de reformation des paires de modes 2-D.

Cette thèse est consacrée à la définition et la mise en oeuvre d'une méthode multirésolution pour le calcul de la structure électronique d'un système composé de plusieurs noyaux et d'électrons. Dans le cadre de ce travail, nous nous intéressons à la Théorie de la Fonctionnelle de la Densité, dans laquelle le potentiel agissant sur chaque orbitale est exprimé via la densité électronique, selon un terme coulombien et un terme non linéaire. La détermination de l'état fondamental conduit au système d'équations de Kohn et Sham. La forme du potentiel de l'opérateur hamiltonien est décrite dans le chapitre deux. La résolution numérique requiert des fonctions de base possédant de bonnes propriétés physiques et algorithmiques. Les deux familles de bases couramment utilisées ne permettent pas simultanément de construire une base adaptée au système physique et de résoudre efficacement le problème auto-cohérent. C'est pourquoi l'on s'intéresse dans cette thèse à des bases d'ondelettes orthogonales et biorthogonales, et à leurs propriétés vis-à-vis du problème d'interpolation, que l'on présente dans le chapitre trois. Dans le chapitre suivant on détermine le potentiel coulombien en résolvant l'équation de Poisson, par des algorithmes itératifs utilisant le préconditionnement du laplacien en base d'ondelettes, et une méthode multigrille. On détaille ensuite la discrétisation du système d'équations par une méthode combinant formulation de Galerkin et méthode de collocation. On analyse enfin l'ordre de l'approximation pour l'oscillateur harmonique et l'hydrogène, et le comportement du système autocohérent pour différents systèmes physiques.

Le chauffage par induction électromagnétique est un procédé efficace permettant de chauffer directement une zone d'épaisseur contrôlée sous la surface de pièces métalliques en vue de les tremper. Cette thèse présente un modèle mathématique couplé électromagnétique/thermique et des approches numériques pour modéliser le procédé. Le modèle électromagnétique est basé sur une formulation en potentiel vecteur magnétique. Les courants de source sont imposés à l'aide d'une formulation en potentiel scalaire électrique permettant de modéliser des inducteurs de forme géométrique arbitraire. Le problème du transfert de chaleur est modélisé à l'aide de l'équation classique de diffusion de la chaleur. Le modèle électromagnétique est entièrement transitoire, afin de permettre l'introduction des effets non linéaires. La discrétisation spatiale est basée sur une approche éléments d'arêtes en utilisant un domaine global air/pièce/inducteur. Le système linéaire d'équations issu de la formulation implicite est creux et défini semi-positif ; il possède un noyau de taille importante. Il est démontré qu'un préconditionneur basé sur une méthode multigrille algébrique construit conjointement avec un solveur du type Krylov réduit substantiellement le temps de calcul du problème électromagnétique par rapport aux méthodes classiques de solution et peut être très efficace pour le calcul parallèle. Des exemples d'application pour le traitement thermique d'un pignon et pour un vilebrequin automobile sont présentés. Le traitement thermique des surfaces des pièces aux géométries complexes nécessite l'introduction d'un mouvement relatif de la pièce et de l'inducteur pour assurer un traitement homogène de la surface. Une nouvelle méthode est proposée, basée sur une représentation discrète d'une fonction level set du mouvement de l'inducteur qui peut être utilisée pour générer des maillages éléments finis conformes dans le cadre d'une configuration lagrangienne.

La prévention des risques industriels nécessite de simuler la dispersion turbulente de polluants. Cependant, les outils majoritairement utilisés à ce jour ne permettent pas de traiter les champs proches dans le cas de géométries complexes, et il est nécessaire d'utiliser les outils de CFD (" Computational Fluid Dynamics ") plus adaptés, mais plus coûteux. Afin de simuler les écoulements atmosphériques avec dispersion de polluants, les modèles CFD doivent modéliser correctement d'une part, les effets de flottabilité, et d'autre part les effets de la turbulence. Plusieurs approches existent, notamment dans la prise en compte des effets de flottabilité et la modélisation de la turbulence, et nécessitent des méthodes numériques adaptées aux spécificités mathématiques de chacune d'entre elles, ainsi que des schémas numériques précis pour ne pas polluer la modélisation. Une formulation d'ordre élevé en volumes finis, sur maillages non structurés, parallélisée, est proposée pour simuler les écoulements atmosphériques avec dispersion de polluants. L'utilisation de schémas d'ordre élevé doit permettre d'une part de réduire le nombre de cellules et diminuer les temps de simulation pour atteindre une précision donnée, et d'autre part de mieux contrôler la viscosité numérique des schémas en vue de simulations LES (Large Eddy Simulation), pour lesquelles la viscosité numérique des schémas peut masquer les effets de la modélisation. Deux schémas d'ordre élevé ont été étudiés et implémentés dans un solveur 3D Navier Stokes incompressible sur des maillages volumes finis non structurés. Nous avons développé un premier schéma d'ordre élevé, correspondant à un schéma Padé volumes finis, et nous avons étendu le schéma de reconstruction polynomiale de Carpentier (2000) aux écoulements incompressibles. Les propriétés numériques des différents schémas implémentés dans le même code de calcul sont étudiées sur différents cas tests bi-dimensionnels (calcul de flux convectifs et diffusifs sur une solution a-priori, convection d'une tâche gaussienne, décroissance d'un vortex de Taylor et cavité entraînée) et tri-dimensionnel (écoulement autour d'un obstacle cubique). Une attention particulière a été portée à l'étude de la précision et du traitement des conditions limites. L'implémentation proposée du schéma polynomial permet d'approcher, pour un maillage identique, les temps de simulation obtenus avec un schéma décentré classique d'ordre 2, mais avec une précision supérieure. Le schéma compact donne la meilleure précision. En utilisant une méthode de Jacobi sans calcul implicite de la matrice pour calculer le gradient, le temps de simulation devient intéressant uniquement lorsque la précision requise est importante. Une alternative est la résolution du système linéaire par une méthode multigrille algébrique. Cette méthode diminue considérablement le temps de calcul du gradient et le schéma Padé devient performant même pour des maillages grossiers. Enfin, pour réduire les temps de simulation, la parallélisation des schémas d'ordre élevé est réalisée par une décomposition en sous domaines. L'assemblage des flux s'effectue naturellement et différents solveurs proposés par les librairies PETSC et HYPRE (solveur multigrille algébrique et méthode de Krylov préconditionnée) permettent de résoudre les systèmes linéaires issus de notre problème.

La turbulence modifie sensiblement les propriétés des écoulements évoluant en zone de proche paroi en provoquant notamment une augmentation significative du frottement pariétal. Dans le but de réduire cette traînée turbulente, on propose dans cette étude un contrôle actif de type dual dont la principale caractéristique repose sur l'interaction de deux schémas de contrôle, un schéma en boucle ouverte par soufflage localisé en espace et périodique en temps, suivi d'un schéma en boucle fermée de type suboptimal. L'idée principale est alors d'organiser la turbulence par une excitation oscillante de manière à la rendre plus prédictible et de ce fait plus contrôlable.<br />On commencera par étudier séparément l'impact de ces deux types de contrôle sur la turbulence de proche paroi. Une étude expérimentale est menée afin d'investiguer en détail l'effet de l'instationnarité induite par une excitation périodique. Les méthodes suboptimales ne pouvant être à ce jour appliquées expérimentalement, une Simulation Numérique Directe s'avère être l'outil le mieux adapté pour étudier leurs comportements. On montre finalement qu'une stratégie de contrôle duale, c'est à dire combinant les deux manipulations (Périodique + Suboptimale), peut s'avérer être une méthode de contrôle efficace afin de réduire la contrainte turbulente.

La prévention des risques industriels nécessite de simuler la dispersion turbulente de polluants. Cependant, les outils majoritairement utilisés à ce jour ne permettent pas de traiter les champs proches dans le cas de géométries complexes, et il est nécessaire d'utiliser les outils de CFD (“ Computational Fluid Dynamics ”) plus adaptés, mais plus coûteux. Afin de simuler les écoulements atmosphériques avec dispersion de polluants, les modèles CFD doivent modéliser correctement d'une part, les effets de flottabilité, et d'autre part les effets de la turbulence. Plusieurs approches existent, notamment dans la prise en compte des effets de flottabilité et la modélisation de la turbulence, et nécessitent des méthodes numériques adaptées aux spécificités mathématiques de chacune d'entre elles, ainsi que des schémas numériques précis pour ne pas polluer la modélisation. Une formulation d'ordre élevé en volumes finis, sur maillages non structurés, parallélisée, est proposée pour simuler les écoulements atmosphériques avec dispersion de polluants. L'utilisation de schémas d'ordre élevé doit permettre d'une part de réduire le nombre de cellules et diminuer les temps de simulation pour atteindre une précision donnée, et d'autre part de mieux contrôler la viscosité numérique des schémas en vue de simulations LES (Large Eddy Simulation), pour lesquelles la viscosité numérique des schémas peut masquer les effets de la modélisation. Deux schémas d'ordre élevé ont été étudiés et implémentés dans un solveur 3D Navier Stokes incompressible sur des maillages volumes finis non structurés. Nous avons développé un premier schéma d'ordre élevé, correspondant à un schéma Padé volumes finis, et nous avons étendu le schéma de reconstruction polynomiale de Carpentier (2000) aux écoulements incompressibles. Les propriétés numériques des différents schémas implémentés dans le même code de calcul sont étudiées sur différents cas tests bi-dimensionnels (calcul de flux convectifs et diffusifs sur une solution a-priori, convection d'une tâche gaussienne, décroissance d'un vortex de Taylor et cavité entraînée) et tri-dimensionnel (écoulement autour d'un obstacle cubique). Une attention particulière a été portée à l'étude de la précision et du traitement des conditions limites. L'implémentation proposée du schéma polynomial permet d'approcher, pour un maillage identique, les temps de simulation obtenus avec un schéma décentré classique d'ordre 2, mais avec une précision supérieure. Le schéma compact donne la meilleure précision. En utilisant une méthode de Jacobi sans calcul implicite de la matrice pour calculer le gradient, le temps de simulation devient intéressant uniquement lorsque la précision requise est importante. Une alternative est la résolution du système linéaire par une méthode multigrille algébrique. Cette méthode diminue considérablement le temps de calcul du gradient et le schéma Padé devient performant même pour des maillages grossiers. Enfin, pour réduire les temps de simulation, la parallélisation des schémas d'ordre élevé est réalisée par une décomposition en sous domaines. L'assemblage des flux s'effectue naturellement et différents solveurs proposés par les librairies PETSC et HYPRE (solveur multigrille algébrique et méthode de Krylov préconditionnée) permettent de résoudre les systèmes linéaires issus de notre problème. / The prevention of industrial risks requires simulating turbulent dispersion of pollutants. However, the tools mostly used so far do not allow near fields treated in the case of complex geometries, and it is necessary to utilize the tools of CFD (Computational Fluid Dynamics ") more suitable but more expensive. To simulate atmospheric flows with dispersion of pollutants, the CFD models must correctly model the one hand, the effects of buoyancy, and secondly the effects of turbulence. Several approaches exist, including taking into account the effects of buoyancy and turbulence modeling, and require numerical methods adapted to the specific mathematics of each, and accurate numerical schemes to avoid pollution modeling. A formulation of high order finite volume on unstructured meshes, parallelized, is proposed to simulate the atmospheric flows with dispersion of pollutants. The use of high order schemes allow one hand to reduce the number of cells and decrease the simulation time to achieve a given accuracy, and secondly to better control the viscosity numerical schemes for simulation LES (Large Eddy Simulation), for which the numerical viscosity patterns may mask the effects of modeling. Two high-order schemes have been studied and implemented in a 3D Navier Stokes solver on unstructured mesh finite volume. We developed the first high-order scheme, corresponding to a Padé finite volume scheme, and we have extended the scheme of reconstruction polynomial Carpentier (2000) for incompressible flows. The numerical properties of the various schemes implemented in the same computer code are studied different two-dimensional test cases (calculation of diffusive and convective flow on a solution a priori, a task Gaussian convection, decay of a vortex of Taylor and driven cavity) and tri-dimensional (flow past an obstacle cubic). Particular attention has been paid to the study of the accuracy and treatment of boundary conditions. The implementation of the polynomial allows to obtain quasi identical simulation time compared to a classical upwind scheme of order 2, but with higher accuracy. The compact layout gives the best accuracy. Using a Jacobi method without calculation implied matrix to calculate the gradient, the simulation time becomes interesting only when the required accuracy is important. An alternative is the resolution of linear system by an algebraic multigrid method. This method significantly reduces the computation time of the gradient and the Padé scheme is effective even for coarse meshes. Finally, to reduce simulation time, the parallelization schemes of high order is achieved by a decomposition into subdomains. The assembly flow occurs naturally and different solvers provided by PETSc libraries and HYORE (algebraic multigrid solver and preconditioned Krylov method) used to solve linear systems from our problem. The work was to identify and determine the parameters that lead to lowest time resolution simulation. Various tests of speed-up and scale-up were used to determine the most effective and optimal parameters for solving linear systems in parallel from our problem. The results of this work have been the subject of a communication in an international conference "Parallel CFD 2008" and an article submitted to "International Journal for Numerical Methods in Fluids" (Analysis of high-order finite volume schemes for the incompressible Navier Stokes equations)

Le chauffage par induction électromagnétique est un procédé efficace permettant de chauffer directement une zone d'épaisseur contrôlée sous la surface de pièces métalliques en vue de les tremper. Cette thèse présente un modèle mathématique couplé électromagnétique/thermique et des approches numériques pour modéliser le procédé. Le modèle électromagnétique est basé sur une formulation en potentiel vecteur magnétique. Les courants de source sont imposés à l'aide d'une formulation en potentiel scalaire électrique permettant de modéliser des inducteurs de forme géométrique arbitraire. Le problème du transfert de chaleur est modélisé à l'aide de l'équation classique de diffusion de la chaleur. Le modèle électromagnétique est entièrement transitoire, afin de permettre l'introduction des effets non linéaires. La discrétisation spatiale est basée sur une approche éléments d'arêtes en utilisant un domaine global air/pièce/inducteur. Le système linéaire d'équations issu de la formulation implicite est creux et défini semi-positif ; il possède un noyau de taille importante. Il est démontré qu'un préconditionneur basé sur une méthode multigrille algébrique construit conjointement avec un solveur du type Krylov réduit substantiellement le temps de calcul du problème électromagnétique par rapport aux méthodes classiques de solution et peut être très efficace pour le calcul parallèle. Des exemples d'application pour le traitement thermique d'un pignon et pour un vilebrequin automobile sont présentés. Le traitement thermique des surfaces des pièces aux géométries complexes nécessite l'introduction d'un mouvement relatif de la pièce et de l'inducteur pour assurer un traitement homogène de la surface. Une nouvelle méthode est proposée, basée sur une représentation discrète d'une fonction level set du mouvement de l'inducteur qui peut être utilisée pour générer des maillages éléments finis conformes dans le cadre d'une configuration lagrangienne. / Electromagnetic induction heating is an efficient process allowing to directly heat up a prescribed area beneath the surface of metallic workpieces to enable quenching. This work presents a mathematical model for the coupled electromagnetic/heat transfer process as well as numerical solution methods. The electromagnetic model is based on a magnetic vector potential formulation. The source currents are prescribed using a voltage potential formulation enabling the modelling of arbitrary inductor geometries. The heat transfer problem is modelled using the classical heat diffusion equation. The electromagnetic model is fully transient, in order to allow the introduction of non-linear effects. The space discretisation is based on an edge finite element approach using a global domain including air, workpiece and inductor. The resulting linear system of equations of the implicit formulation is sparse and semi-definite, including a large kernel. It is demonstrated that a preconditioner based on the auxiliary space algebraic multigrid method in connection with a Krylov solver substantially reduces the solution time of the electromagnetic problem in comparison to classical solution methods and can be effectively applied in parallel. Applications for the heat treatment of a gearwheel and for an automotive crankshaft are presented. The surface heat treatment of complex geometrical parts requires the introduction of a relative movement of workpiece and inductor to ensure a homogeneous surface treatment. A novel method is proposed, which is based on a discrete level set representation of the inductor motion that can be used to generate conforming finite element meshes in a Lagrangian setting.

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.