Étude numérique d'écoulements de fluides par une méthode vortex la marche descendante et les cavités sur les ailes de papillons

Abstract

Dans cette thèse sont considérés des écoulements isothermes de fluides visqueux (newtoniens) et incompressibles sans forces de volume. Une méthode numérique bidimensionnelle basée sur l’équation du transport de la vorticité est utilisée pour simuler les écoulements en considérant trois modèles géométriques. Plus spécifiquement, la méthode des vortex aléatoires est combinée à l’algorithme Vortex-in-cell où des transformations algébriques permettent de raffiner le maillage aux endroits critiques des écoulements. Dans cette méthode mixte eulérienne-lagrangienne, la discrétisation principale se fait au niveau du champ de vorticité, celui-ci étant représenté par un ensemble d’éléments lagrangiens (singuliers) qui transportent chacun une quantité de circulation prédéterminée. Ces éléments de vorticité sont générés sur les parois solides pour vérifier la condition d’adhérence et évoluent dans le temps selon l’équation du transport de la vorticité. Le nombre d’éléments à être générés à chaque pas de temps est déterminé par l’entremise d’un nouvel algorithme de génération qui est construit à partir d’algorithmes existants. Étant bien étudié et bien connu, l’écoulement débutant (starting flow) au-dessus d’une marche descendante s’avère un excellent modèle de validation pour des méthodes numériques bidimensionnelles à nature instationnaire : des simulations numériques ont par conséquent été effectuées dans le cadre de cette thèse en considérant un tel modèle géométrique. Les résultats de simulations sont comparés avec succès à des résultats expérimentaux où une analyse plus générale démontre les aptitudes de la méthode à simuler ce type d’écoulements. En particulier, les résultats montrent que pour des écoulements à faibles valeurs du nombre de Reynolds (Re = 97), la zone de recirculation est composée d’une seule région de vorticité pour la durée totale des simulations. Pour des écoulements à valeurs plus élevées du nombre de Reynolds (Re = 153 et 303), la zone de recirculation est respectivement composée de trois et de quatre régions de vorticité durant les stades intermédiaires de l’écoulement, tandis que pour des temps plus élevés, les structures à l’intérieur de la zone de recirculation ne sont pas très bien définies. Les résultats dévoilent aussi que, pour la gamme complète des valeurs du nombre de Reynolds étudiée, la distance qui sépare la paroi verticale et le point de recollement de la région principale de vorticité augmente quasi linéairement avec le temps durant les stades intermédiaires de l’écoulement.

Abstract: Isothermal flows of viscous (Newtonian) fluids without body forces are considered in this thesis. A two-dimensional numerical method based on the vorticity transport equation is used to simulate flows for three different model geometries. Specifically, the Random Vortex method is combined with the Vortex-In-Cell algorithm. In this mixed Eulerian-Lagrangian method, the main discretization is performed on the vorticity field, which is represented by a number of singular Lagrangian vortex elements. These vortex elements are generated on solid walls to verify the no-slip condition and are subsequently transported in time according to the vorticity transport equation. The laminar starting flow down a step is a well defined flow. Considering this flow model, instantaneous numerical simulations have been carried out and the results are compared to experimental results in the literature. Good agreement is obtained. Analysis of the numerical results and comparisons to experimental results show that for low Reynolds number flows (Re = 97), the recirculation zone is composed of one vorticity region throughout the simulations. For higher Reynolds number flows (Re = 153 and 303), the recirculation zone is respectively composed of three and four distinct vorticity regions at intermediate stages of its development, while for later times the structures inside the recirculation zone are not clearly defined. During intermediate stages of the flow development, and for the range of Reynolds numbers investigated, the distance from the step to the reattachment point of the main vorticity region increases quasi-linearly with the dimensionless time. Again, good agreement is obtained between experimental data and computed solutions. Using the same flow model, we present a numerical convergence proof by successively refining the numerical parameters. It is shown that the refinement of the vorticity generation parameters seems to primarily affect the smoothness of the solution rather than the overall structure of the flow. On the other hand, it is shown that for a given time step, higher order convection schemes improve the overall flow structure. Steady state numerical simulations were performed. Numerical streamwise velocity profiles are compared to their experimental counterparts for four Reynolds number flows (Re = 73, 125, 191 and 229). The computed length of the recirculation zone is also compared to various results. In all steady step flows, excellent agreement was obtained. Finally, both steady and unsteady flow simulations were performed on a model representing the cavities that result from the shingle-like arrangement of the scales on the upper surface of a butterfly wing. In gliding flight, the lifting force was experimentally proven to be stronger on a wing with scales. Excellent agreement is obtained in comparing results of our numerical simulations to experimental results in the literature. In particular, for low and very low Reynolds number flows (Re = 0.62, 1.00, 3.30 and 100), the recirculation zone is composed of one primary vorticity region. For higher Reynolds number flows (Re = 624), the recirculation zone area exhibits strong dynamics where coherent structures are shed at regular intervals from the vertical wall and eventually coalesce during intermediate stages of the developing flow"--Résumé abrégé par UMI