Cette thèse propose la reformulation et l'implantation d'une méthode simplifiée de prévision de la transition laminaire-turbulent naturelle, la méthode des paraboles. Elle s'appuie sur une base de donnée pour fournir une expression analytique des taux d'amplification de perturbations modales de la couche limite. Ces perturbations prennent naissance dans la couche limite à travers un mécanisme de réceptivité, puis connaissent une croissance linéaire suivie d'interactions non-linéaires qui conduisent à la transition vers la turbulence.La méthode du facteur-N permet la prévision de la position de transition en supposant qu'elle a lieu quand l'amplification linéaire totale d'un mode des perturbations a atteint une amplitude critique.L'amplification totale est obtenue en intégrant les taux d'amplifications des perturbations le long de la ligne de courant extérieure à la couche limite.Cette formulation intégrale n'est pas adaptée aux codes RANS et ces travaux proposent une reformulation de la méthode du facteur-N sous forme d'équation de transport, dont le terme source dépend du taux d'amplification des perturbations donné par la méthode des paraboles.Le modèle qui en résulte (N-sigma-P ou NSP) permet de déterminer si un point donné est en amont ou en aval du point de transition et ainsi de prévoir la position de transition, pour des configurations où la transition est déclenchée par des modes longitudinaux (ondes de Tollmien-Schlichting) ou transverses (crossflow). / A simplified laminar-turbulent transition prediction method (the so-called parabolas method) is recast into a RANS compatible formulation.The parabolas method is a database approach that gives an estimation of the growth rates of modal perturbations growing in the boundary layer.These perturbations rise through a process called receptivity and undergo a linear modal growth, followed by non-linear interactions that lead to a breakdown to transition. The N-factor (or e-to-N) method predicts the transition location under the assumption that transition occurs when the total linear growth of a mode reaches a prescribed threshold. This total growth is obtained via an integration of a mode's growth rate along the boundary-layer edge streamline. This classical integral formulation is not suitable in a RANS context and is therefore recast into a set of transport equations whose source terms depend on the growth rates given by the parabolas method.The resulting model (N-sigma-P or NSP) predicts if a given point in the flow field is up- or downstream of the transition location. It has been validated for cases where transition was triggered by longitudinal (Tollmien-Schlichting) and crossflow modes.

La thèse vise à intégrer des critères de transition dans le solveur Navier-Stokes non structuré Aether utilisé chez Dassault Aviation. Une méthodologie de prévision de la transition laminaire/turbulent a été élaborée et implémentée dans le solveur RANS Aether. Deux stratégies de calcul de transition ont été testées. Soit Aether est couplé avec le code de couche limite de l’ONERA 3C3D. Soit la position de transition est calculée en utilisant directement les profils de vitesse RANS. Les deux méthodes ont été testées pour des écoulements subsoniques et transsoniques. L’influence des solveurs numériques, des critères de transition et du processus de couplage sont étudiés. L’utilisation de schémas numériques d’ordre élevé est également considérée. / This thesis present the state-of-the-art of the transition prediction numerical chain which has been developed at Dassault Aviation in the RANS solver AETHER. Two strategies for transition location estimations exist. First, AETHER is coupled with the ONERA boundary layer code 3C3D. Second, the transition location is computed by using directly the RANS velocity profiles. Both methods were preliminarily tested in subsonic and transonic. The issues of the influence of the numerical solvers, transition onset criteria and coupling process are studied. The influence of higher order numerical method, are also considered.