L'objectif de cette étude est de simuler numériquement le comportement des écoulements turbulents dilatables par la méthode des éléments finis. Le choix de la méthode a été adopte pour divers aspects numériques, entre autres sa capacité à décrire les géométries complexes. Après traitement statistique des équations de conservation de la masse, quantité de mouvement et l'équation d'énergie, l'analyse conduit naturellement au problème de fermeture des termes traduisant les effets de compressibilité. Tout d'abord quelques approximations de base ont été introduites afin d'analyser les termes lies aux fluctuations de masse volumique. Ensuite deux types de modélisation ont été établis pour le terme traduisant l'effet d'interaction entre les fluctuations de masse volumique et le gradient de pression moyenne. L’étude montre entre autres, le rôle joue par ce gradient pour les écoulements pesants. Enfin, l'insuffisance du modèle classique (k-) (base sur des scalaires dont aucun n'est représentatif de l'anisotropie) nous a conduit à formuler un modèle algébrique simple tenant compte de l'effet anisotrope de la gravite sur la turbulence. Les applications traitées concernent l'épanouissement d'un jet chaud dans l'atmosphère au repos et l'étude du comportement de la couche atmosphérique cinématique et thermique. Les résultats obtenus ont été confrontes a des expériences ou a des calculs tires de littérature

La méthode des éléments finis (m.e.f.) A été appliquée avec succès pour représenter l'hydrodynamique a surface libre, mais avec une modélisation de la turbulence très simplifiée. Ce travail a pour but de simuler, a l'aide de la m.e.f., des écoulements a surface libre avec des représentations de la turbulence plus complexes. Il commence par une présentation générale des modèles de turbulence. Nous avons choisi d'utiliser plus particulièrement deux types de modèle: le modèle k-, des modèles aux longueurs de mélange développés a partir de la simplification du modèle précédent. Une simplification supplémentaire est effectuée en bidimensionnalisant les modèles. Un travail de ce type a déjà été effectué par rodi et al pour le modèle k-. Nous l'avons complété pour les modèles aux longueurs de mélange. Dans le cas ou les courbures de l'écoulement horizontal sont importantes, nous proposons des variantes des deux types de modèle précédents et par ailleurs, une modification des frottements de fond intervenant dans l'équation de l'écoulement (saint-venant). Cette dernière est notamment importante pour les études sedimentologiques. La méthode de résolution des équations du modèle k- utilisée dans cette thèse est la méthode des éléments finis ; il s'agit alors de résoudre un système non linéaire complexe ; dans ce travail, nous avons fait un choix d'interpolation particulier et teste plusieurs méthodes de résolution. Le couplage avec le modèle hydrodynamique est effectue par un algorithme explicite. Sont présentés plusieurs critères pour contrôler la convergence du modèle k-. Nous terminons par une validation des modèles de turbulence utilises sur de nombreux exemples. Sur ces derniers, nous avons vérifié les méthodes employées pour résoudre les équations du modèle k-

Dans différents procédés industriels on retrouve le problème de base du mélange de deux fluides. Les caractéristiques complexes des écoulements internes dans les dispositifs de mélange utilises justifient fortement l'utilisation de la simulation numérique, qui doit allier des méthodes numériques performantes et des modèles de turbulence adaptes pour une description précise des phénomènes de transfert de masse ou de chaleur pendant le processus de mélange. L’objectif principal de ce travail est ainsi d'étudier des modèles de turbulence applicables aux situations de mélange industriel. On propose, en outre, une implantation numérique sous la forme d'un code de calcul applicable aux géométries complexes tridimensionnelles. En ce qui concerne le calcul hydrodynamique, on propose l'utilisation du modèle de turbulence a deux équations k- dans un contexte d'écoulement tridimensionnel. On évalue aussi quelques modifications de ce modèle pour le cas ou une rotation d'ensemble agit sur le champ turbulent. Concernant la modélisation du transport scalaire, deux types d'approche sont étudiés à savoir: une approche eulérienne classique avec fermeture au premier ordre, et une autre approche de type lagrangien-stochastique dont l'utilisation est un peu moins classique dans ce type de situation. Face aux géométries complexes qu'on retrouve dans les problèmes de mélange industriel, on a choisi une approche numérique s'appuyant sur la méthode des éléments finis. Un ensemble de résolution de problèmes hydrodynamiques et thermiques est présenté pour lesquels un nombre important d'informations d'origine expérimentale est disponible. On s'intéresse notamment a la résolution d'un problème complexe de mélange dans une situation de jet turbulent tridimensionnel pénétrant perpendiculairement dans une couche limite

Cette thèse a pour objectif la simulation numérique d'écoulements turbulents à l'aide d'un schéma d'éléments finis stabilisés. La méthode SUPG (Streamline Upwind Petrov Galerkin) est mise en place pour résoudre les équations de Navier-Stokes compressibles couplées au modèle turbulent de Spalart-Allmaras.

L'objet de cette étude est la recherche d'un outil efficace de prédiction des écoulements turbulents complexes, rencontrés dans les problèmes d'ingéniérie. La technique numérique choisie est la méthode des éléments finis, et seuls les écoulements incompressibles et isothermes de fluides newtoniens sont considérés. Un exposé du problème posé par la simulation numérique des écoulements turbulents ainsi que les méthodes de simulation par modélisation des structures de sous-maille, qui découlent naturellement de cette analyse, sont présentés. L'utilisation de la technique de l'homogénéisation, qui rentre dans cette catégorie, est également décrite. L'autre façon d'aborder le problème est statistique. Et la notion de description statistique de la turbulence étant introduite, une sélection des modèles statistiques en un point est entreprise. Les deux modèles retenus, modèle K- et modèle aux tensions de Reynolds algébrique (A.S.M.), font l'objet d'une analyse détaillée. L'accent est mis sur la nécéssité, pour les écoulements à géometrie complexe, d'une formulation du modèle K- non limitée aux grands nombres de Reynolds, conduisant à une utilisation plus rigoureuse de la technique des fonctions de paroi, dont une expression généralisée est proposée. Le modèle A.S.M. quant à lui, est exprimé de façon à mettre en evidence son lien avec le modèle K-, ce qui autorise par ailleurs une formulation tridimensionnelle par éléments finis commune de ces deux modèles, exposée en détail. Une discussion sur les performances de cette formulation du modèle K-, ainsi que sur les difficultés numériques liées à l'utilisation de ce modèle dans le cadre des éléments finis, est enfin conduite à partir d'exemples numériques.