Os mecanismos de instabilidade hidrodinâmica têm um papel importante no processo da transição do escoamento de laminar para turbulento. A análise da instabilidade hidrodinâmica em uma cavidade com tampa deslizante foi realizada através da decomposição em modos globais (biglobal) para avaliar o efeito da compressibilidade neste fenômeno. O escoamento base foi obtido através de simulação numérica direta (DNS). Para tal, foi desenvolvido um código DNS compressível com discretização espacial por diferenças finitas compactas de alta resolução espectral e capacidade de processamento paralelo, com um método de decomposição de domínio que mantém a precisão das diferenças finitas compactas. O escoamento base é usado para montar o problema de autovalor oriundo das equações de Navier-Stokes linearizadas para a perturbação, discretizadas por diferenças finitas explícitas. O uso de diferenças finitas em conjunto com a implementação em matrizes esparsas reduz sensivelmente o uso de memória. Através do algoritmo de Arnoldi, a ordem do problema de autovalor é reduzida e os autovalores de interesse são recuperados. Os resultados indicam o efeito estabilizante da compressibilidade nos modos dominantes da cavidade e revelam modos inerentes ao escoamento compressível, para os quais a compressibilidade tem efeito desestabilizante. Dentre estes modos compressíveis, estão presentes modos de propagação sonora em dutos e modos relacionados à geração de som na cavidade. / Hydrodynamic instability mechanisms play an important role in laminar to turbulent transition. Hydrodynamic instability analysis of a lid-driven cavity flow was performed by global mode decomposition (biglobal) to evaluate compressibility effects on this phenomenon. The basic flow was calculated by direct numerical simulation (DNS). A compressible DNS code was developed with spectral-like compact finite difference spatial discretization. The code allows parallel processing with a domain decomposition method that preserves the compact finite difference accuracy. The basic flow is used to form the eigenvalue problem associated to the linear Navier- Stokes equations for the perturbation, which were discretized by an explicit finite difference scheme. The combination of sparse matrix techniques and finite difference discretization leads to a significant memory reduction. The order of the eigenvalue problem was reduced using the Arnoldi algorithm and the eigenvalues of interest were calculated. Results show the stabilizing effect of compressibility on the leading modes and reveal some modes intrinsic to compressible flow, for which compressibility has a destabilizing effect. Among these compressible modes, there are some related to sound propagation in ducts and to sound generation inside the cavity.

Os mecanismos de instabilidade hidrodinâmica têm um papel importante no processo da transição do escoamento de laminar para turbulento. A análise da instabilidade hidrodinâmica em uma cavidade com tampa deslizante foi realizada através da decomposição em modos globais (biglobal) para avaliar o efeito da compressibilidade neste fenômeno. O escoamento base foi obtido através de simulação numérica direta (DNS). Para tal, foi desenvolvido um código DNS compressível com discretização espacial por diferenças finitas compactas de alta resolução espectral e capacidade de processamento paralelo, com um método de decomposição de domínio que mantém a precisão das diferenças finitas compactas. O escoamento base é usado para montar o problema de autovalor oriundo das equações de Navier-Stokes linearizadas para a perturbação, discretizadas por diferenças finitas explícitas. O uso de diferenças finitas em conjunto com a implementação em matrizes esparsas reduz sensivelmente o uso de memória. Através do algoritmo de Arnoldi, a ordem do problema de autovalor é reduzida e os autovalores de interesse são recuperados. Os resultados indicam o efeito estabilizante da compressibilidade nos modos dominantes da cavidade e revelam modos inerentes ao escoamento compressível, para os quais a compressibilidade tem efeito desestabilizante. Dentre estes modos compressíveis, estão presentes modos de propagação sonora em dutos e modos relacionados à geração de som na cavidade. / Hydrodynamic instability mechanisms play an important role in laminar to turbulent transition. Hydrodynamic instability analysis of a lid-driven cavity flow was performed by global mode decomposition (biglobal) to evaluate compressibility effects on this phenomenon. The basic flow was calculated by direct numerical simulation (DNS). A compressible DNS code was developed with spectral-like compact finite difference spatial discretization. The code allows parallel processing with a domain decomposition method that preserves the compact finite difference accuracy. The basic flow is used to form the eigenvalue problem associated to the linear Navier- Stokes equations for the perturbation, which were discretized by an explicit finite difference scheme. The combination of sparse matrix techniques and finite difference discretization leads to a significant memory reduction. The order of the eigenvalue problem was reduced using the Arnoldi algorithm and the eigenvalues of interest were calculated. Results show the stabilizing effect of compressibility on the leading modes and reveal some modes intrinsic to compressible flow, for which compressibility has a destabilizing effect. Among these compressible modes, there are some related to sound propagation in ducts and to sound generation inside the cavity.

A investigação dos mecanismos de instabilidade pode ter um papel importante no entendimento do processo laminar para turbulento de um escoamento. Análise de instabilidade de uma camada limite de uma linha de estagnação compressível foi realizada no contexto de teoria linear BiGlobal. O estudo dos mecanismos de instabilidade deste escoamento pode proporcionar uma visão útil no desenho aerodinâmico das asas. Um novo procedimento foi desenvolvido e implementado computacionalmente de maneira sequencial e paralela para o estudo de instabilidade BiGlobal. O mesmo baseia-se em formar a matriz esparsa associada ao problema discretizado por dois métodos: pontos de colocação de Chebyshev-Gauss-Lobatto e diferenças finitas, além das combinações destes métodos. Isto permitiu o uso de bibliotecas computacionais eficientes para resolver o sistema linear associado ao problema de autovalor utilizando o algoritmo de Arnoldi. O desempenho do método numérico e código computacional proposto são analisados do ponto de vista do uso de métodos de ordenação dos elementos da matriz, coeficientes de preenchimento, memória e tempo computacional a fim de determinar a solução mais eficiente para um problema físico geral com técnicas de matrizes esparsas. Um estudo paramétrico da instabilidade da camada limite de uma linha de estagnação foi realizado incluindo o estudo dos efeitos de compressibilidade. O excelente desempenho código computacional permitiu obter as curvas neutras e seus respectivos valores críticos para a faixa de número de Mach 0 \'< OU =\' Ma \'< OU =\' 1. Os resultados confirmam a teoria assintótica apresentada por (THEOFILIS; FEDOROV; COLLIS, 2004) e mostram que o incremento do número de Mach reduz o numero de Reynolds crítico e a faixa instável do número de ondas. / Investigation of linear instability mechanisms is essential for understanding the process of transition from laminar to turbulent flow. An algorithm for the numerical solution of the compressible BiGlobal eigenvalue problem is developed. This algorithm exploits the sparsity of the matrices resulting from the spatial discretization of the enigenvalue problem in order to improve the performance in terms of both memory and CPU time over previous dense algebra solutions. Spectral collocation and finite differences spatial discretization methods are implemented, and a performance study is carried out in order to determine the best practice for the efficient solution of a general physical problem with sparse matrix techniques. A combination of spectral collocation and finite differences can further improve the performance. The code developed is then applied in order to revisit and complete the parametric analyses on global instability of the compressible swept Hiemenz flow initiated in (THEOFILIS; FEDOROV; COLLIS, 2004) and obtain neutral curves of this flow as a function of the Mach number in the 0 \'< OU =\' Ma \'< OU =\' 1 range. The present numerical results fully confirm the asymptotic theory results presented in (THEOFILIS; FEDOROV; COLLIS, 2004). This work presents a complete parametric study of the instability properties of modal three dimensional disturbances in the subsonic range for the flow conguration at hand. Up to the subsonic maximum Mach number value studied, it is found that an increase in this parameter reduces the critical Reynolds number and the range of the unstable spanwise wavenumbers.

A investigação dos mecanismos de instabilidade pode ter um papel importante no entendimento do processo laminar para turbulento de um escoamento. Análise de instabilidade de uma camada limite de uma linha de estagnação compressível foi realizada no contexto de teoria linear BiGlobal. O estudo dos mecanismos de instabilidade deste escoamento pode proporcionar uma visão útil no desenho aerodinâmico das asas. Um novo procedimento foi desenvolvido e implementado computacionalmente de maneira sequencial e paralela para o estudo de instabilidade BiGlobal. O mesmo baseia-se em formar a matriz esparsa associada ao problema discretizado por dois métodos: pontos de colocação de Chebyshev-Gauss-Lobatto e diferenças finitas, além das combinações destes métodos. Isto permitiu o uso de bibliotecas computacionais eficientes para resolver o sistema linear associado ao problema de autovalor utilizando o algoritmo de Arnoldi. O desempenho do método numérico e código computacional proposto são analisados do ponto de vista do uso de métodos de ordenação dos elementos da matriz, coeficientes de preenchimento, memória e tempo computacional a fim de determinar a solução mais eficiente para um problema físico geral com técnicas de matrizes esparsas. Um estudo paramétrico da instabilidade da camada limite de uma linha de estagnação foi realizado incluindo o estudo dos efeitos de compressibilidade. O excelente desempenho código computacional permitiu obter as curvas neutras e seus respectivos valores críticos para a faixa de número de Mach 0 \'< OU =\' Ma \'< OU =\' 1. Os resultados confirmam a teoria assintótica apresentada por (THEOFILIS; FEDOROV; COLLIS, 2004) e mostram que o incremento do número de Mach reduz o numero de Reynolds crítico e a faixa instável do número de ondas. / Investigation of linear instability mechanisms is essential for understanding the process of transition from laminar to turbulent flow. An algorithm for the numerical solution of the compressible BiGlobal eigenvalue problem is developed. This algorithm exploits the sparsity of the matrices resulting from the spatial discretization of the enigenvalue problem in order to improve the performance in terms of both memory and CPU time over previous dense algebra solutions. Spectral collocation and finite differences spatial discretization methods are implemented, and a performance study is carried out in order to determine the best practice for the efficient solution of a general physical problem with sparse matrix techniques. A combination of spectral collocation and finite differences can further improve the performance. The code developed is then applied in order to revisit and complete the parametric analyses on global instability of the compressible swept Hiemenz flow initiated in (THEOFILIS; FEDOROV; COLLIS, 2004) and obtain neutral curves of this flow as a function of the Mach number in the 0 \'< OU =\' Ma \'< OU =\' 1 range. The present numerical results fully confirm the asymptotic theory results presented in (THEOFILIS; FEDOROV; COLLIS, 2004). This work presents a complete parametric study of the instability properties of modal three dimensional disturbances in the subsonic range for the flow conguration at hand. Up to the subsonic maximum Mach number value studied, it is found that an increase in this parameter reduces the critical Reynolds number and the range of the unstable spanwise wavenumbers.

For simulations of highly complex geometries, frequently encountered in many fields of science and engineering, the process of generating a high-quality, body-fitted grid is very complicated and time-intensive. Thus, one of the principal goals of contemporary CFD is the development of numerical algorithms, which are able to deliver computationally efficient, and highly accurate solutions for a wide range of applications involving multi-physics problems, e.g. Fluid Structure Interaction (FSI). Immersed interface/boundary methods provide considerable advantages over conventional approaches, especially for flow problems containing moving boundaries.In the present work, a novel, robust, highly-accurate, Immersed Interface Method (IIM) is developed, which is based on a local Taylor-series expansion at irregular grid points enforcing numerical stability through a local stability condition. Various immersed methods have been developed in the past; however, these methods only considered the order of the local truncation error. The numerical stability of these schemes was demonstrated (in a global sense) by considering a number of different test-problems. None of these schemes used a concrete local stability condition to derive the irregular stencil coefficients. This work will demonstrate that the local stability constraint is valid as long as the DFL-number does not reach a limiting value. The IIM integrated into a newly developed Incompressible Navier-Stokes (INS) solver is used herein to simulate fully coupled FSI problems. The extension of the novel IIM to a higher-order method, the compressible Navier-Stokes equations and the Maxwell's equations demonstrate the great potential of the novel IIM.In the second part of this dissertation, the newly developed INS solver is employed to study the flow of a stalled airfoil and steady/unsteady stenotic flows. In this context, a new biglobal stability analysis approach based on solving an Initial Value Problem (IVP), instead of the traditionally used EigenValue Problem (EVP), is presented. It is demonstrated that this approach based on an IVP is computationally less expensive compared to EVP approaches while still capturing the relevant physics.

State of the art research in hydrodynamic stability analysis has moved from classic one-dimensional methods such as the local nonparallel approach and the parabolized stability equations to two-dimensional, biglobal, methods. The paradigm shift toward two dimensional techniques with the ability to accommodate fully three-dimensional base flows is a necessary step toward modeling complex, multidimensional flowfields in modern propulsive applications. Here, we employ a two-dimensional spatial waveform with sinusoidal temporal dependence to reduce the three-dimensional linearized Navier-Stokes equations to their biglobal form. Addressing hydrodynamic stability in this way circumvents the restrictive parallel-flow assumption and admits boundary conditions in the streamwise direction. Furthermore, the following work employs a full momentum formulation, rather than the reduced streamfunction form, accounting for a nonzero tangential mean flow velocity. This approach adds significant complexity in both formulation and implementation but renders a more general methodology applicable to a broader spectrum of mean flows. Specifically, we consider the stability of three models for bidirectional vortex flow. While a complete parametric study ensues, the stabilizing effect of the swirl velocity is evident as the injection parameter, kappa, is closely examined.

Boundary-layer transition from a laminar to a turbulent regime is a critical driver in the design of high-speed vehicles. The aerothermodynamic loads associated with transitional or fully turbulent hypersonic boundary layers are several times higher than those associated with laminar flow. The presence of isolated roughness elements on the surface of a body can accelerate the growth of incoming disturbances and introduce additional instability mechanisms in the flow field, eventually leading to a premature occurrence of transition. This dissertation studies the instabilities induced by three-dimensional discrete roughness elements located inside a high-speed boundary layer developing on a flat plate. Two-dimensional local linear stability theory (2D-LST) is employed to identify the instabilities evolving in the three-dimensional flow field that characterizes the wake induced by the roughness elements and to investigate their evolution downstream. A formulation of the disturbance energy evolution equation available for base flows depending on a single spatial direction is generalized for the first time to base flows featuring two inhomogeneous directions and perturbations depending on three spatial directions. This generalization allows to obtain a decomposition of the temporal growth rate of 2D-LST instabilities into the different contributions that lead to the production and dissipation of the total disturbance energy. This novel extension of the formulation provides an additional layer of information for understanding the energy exchange mechanisms between a three-dimensional base flow and the perturbations resulting from 2D-LST. Stability computations for a calorically perfect gas illustrate that the wake induced by the roughness elements supports the growth of different sinuous and varicose instabilities which coexist together with the Mack-mode perturbations that evolve in the flat-plate boundary layer, and which become modulated by the roughness-element wake. A single pair of sinuous and varicose disturbances is found to dominate the wake instability in the vicinity of the obstacles. The application of the newly developed decomposition of the temporal growth rate reveals that the roughness-induced wake modes extract most of their potential energy from the transport of entropy fluctuations across the base-flow temperature gradients and most of their kinetic energy from the work of the disturbance Reynolds stresses against the base-flow velocity gradients. Further downstream, the growth rate of the wake instabilities is found to be influenced by the presence of Mack-mode disturbances developing on the flat plate. Strong evidence is observed of a continuous synchronization mechanism between the wake instabilities and the Mack-mode perturbations. This phenomenon leads to an enhancement of the amplification rate of the wake modes far downstream of the roughness element, ultimately increasing the associated integrated amplification factors for some of the investigated conditions. The effects of vibrational molecular excitation and chemical non-equilibrium on the instabilities induced by a roughness element are studied for the case of a high-temperature boundary layer developing on a sharp wedge configuration. For this purpose, a 2D-LST solver for chemical non-equilibrium flows is developed for the first time, featuring a fully consistent implementation of the thermal and transport models employed for the base flow and the perturbation fields. This is achieved thanks to the automatic derivation and implementation tool (ADIT) available within the von Karman Institute extensible stability and transition analysis (VESTA) tool-kit, which enables an automatic derivation and implementation of the 2D-LST governing equations for different thermodynamic flow assumptions and models. The stability computations for this configuration show that sinuous and varicose disturbances also dominate the wake instability in the presence of vibrational molecular energy mode excitation and chemical reactions. The resulting base-flow cooling associated with the modeling of such high-temperature phenomena is found to have opposite stabilizing and destabilizing effects on the streamwise evolution of the sinuous and varicose instabilities. The modeling of vibrational excitation and chemical non-equilibrium acting exclusively on the perturbations is found to have a stabilizing influence in all cases. / Doctorat en Sciences de l'ingénieur et technologie / info:eu-repo/semantics/nonPublished