Doctoral Thesis
DOI
https://doi.org/10.11606/T.18.2009.tde-19012011-132153
Document
Author
Full name
Ricardo Alberto Coppola Germanos
Institute/School/College
Knowledge Area
Date of Defense
Published
São Carlos, 2009
Supervisor
Committee
Medeiros, Marcello Augusto Faraco de (President)
Greco Junior, Paulo Celso
Mendonça, Márcio Teixeira de
Moser, Carlos Anissem Soares
Rodrigues, Savio Brochini
Title in Portuguese
Simulação numérica da evolução linear e não linear em uma camada de mistura compressível tridimensional
Keywords in Portuguese
Aproximações de alta ordem
Camada de mistura compressível
Instabilidade hidrodinâmica
Simulação numérica direta
Trens de onda moduladas
Abstract in Portuguese
As aplicações aeroespaciais estão frequentemente associadas a escoamentos compressíveis com altíssimos números de Reynolds. No entanto, existem no contexto aeroespacial importantes aplicações que envolvem escoamentos compressíveis a Reynolds relativamente baixos. Entre eles se destacam o escoamento em pás de turbina a gás e ao redor de dispositivos de alta sustentação como eslates e flapes em grandes ângulos de ataque. Pode-se destacar também o processo de combustão supersônica que está intimamente ligado e é fortemente beneficiado pelo presente estudo. Nas aplicações aerodinâmicas em baixos números de Reynolds frequentemente uma parcela significativa do escoamento se apresenta no regime de transição para turbulência, ou nos estágios iniciais do escoamento turbulento. O objetivo do presente projeto é a simulação numérica direta de escoamentos compressíveis transicionais com desenvolvimento de um código para simulação em três dimensões de escoamentos alto subsônicos. O escoamento a ser estudado no projeto é a evolução linear e não linear de trens de onda e pacotes de onda em uma camada de mistura compressível. A solução das equações de Navier-Stokes é obtida através do método das diferenças finitas. As derivadas espaciais são resolvidas através de um método compacto de sexta ordem, enquanto que as derivadas temporais são resolvidas através do método de Runge-Kutta de quarta ordem. Os métodos de aproximação foram modificados para trabalhar com malhas não uniformes visando refinar a malha em pontos em que o fenômeno ocorre e, consequentemente, reduzir o custo computacional. A investigação numérica inicia-se com a análise da taxa de amplificação dos trens de ondas fortemente modulados em regime linear. Os resultados obtidos foram comparados favoravelmente com a teoria linear. Os testes foram estendidos para a análise não linear, e consequentemente, foi possível reproduzir os fenômenos clássicos de instabilidade hidrodinâmica através da evolução dos trens de ondas oblíquos.
Title in English
Numerical simulation of the linear and non-linear evolution in a three-dimensional compressible mixing layer
Keywords in English
Compressible mixing layer
Direct numerical simulation
High-order schemes
Hydrodynamic instability
Wavetrains
Abstract in English
Aerospace applications are frequently associated with compressible flows at relatively high Reynolds number. Nevertheless important applications involve compressible flows at relatively low Reynolds number in the aerospace context. Among them, the flow on gas turbine blades and high lift devices such as slats and flaps at high angle of attack are particulary important. Besides, progress in aeroespace research is dependent on developing more efficient propulsion systems. In aerodynamic applications at low Reynolds number, often a substancial portion of the flow is in the transition regime, or in the initial stages of a turbulent flow. The objective of the present study is the Direct Numerical Simulation of three-dimensional transition of compressible flows in a mixing layer. Inspired on the worked devoted to modulated waves, the current work investigates the linear and nonlinear temporal evolution of wavetrains in this phenomenon. The Navier-Stokes equations were solved with a sixth-order compact finite-difference schemes. The time integration was performed by a fourth-order Runge-Kutta scheme. Moreover, the methods to solve the spatial derivatives were modified to work with non-uniform grids. This technique was implemented with the objective to improve the resolution of the grid where the phenomenon occurs and to reduce the computational cost. The numerical investigation starts with an analysis of the growth rate of the wavetrains in linear regime to verify the numerical code. The results compared favourably with linear theory. Tests were also performed in the nonlinear regime to simulate the oblique wavetrains and it was possible to reproduce the classical hydrodynamic instability phenomena.
WARNING - Viewing this document is conditioned on your acceptance of the following terms of use:
This document is only for private use for research and teaching activities. Reproduction for commercial use is forbidden. This rights cover the whole data about this document as well as its contents. Any uses or copies of this document in whole or in part must include the author's name.
Publishing Date
2011-02-07