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.

Tanto os motores turbo-jato quanto os turbo-fan, são os maiores responsáveis pela geração de ruído durante a decolagem, segmento de subida e de aceleração de uma aeronave. Devido a isto, o problema de ruído em jatos vem sendo intensamente investigado ao longo dos últimos anos. Já na fase do pouso, o slat é uma das fontes de ruído mais importantes. Para este caso, na maioria das aplicações práticas, existe o descolamento da camada limite no intradorso do slat a partir de onde se desenvolve uma camada de mistura. Ainda assim, existem inúmeros aspectos de tais escoamentos que precisam de investigação. Uma abordagem frequentemente feita para o estudo da instabilidade hidrodinâmica e ruído em jatos é o estudo de metade do jato. A estratégia consiste em estudar os fenômenos na camada de mistura, o que é uma aproximação razoável quando o jato tem diâmetro muito grande comparado à espessura da camada cisalhante que se desenvolve nas bordas do jato. Assim, alguns aspectos do ruído gerado pelos modos axi-simétricos de instabilidade são em grande parte reproduzidos. Um aspecto aparentemente jamais estudado antes é o efeito do emparelhamento de vórtices de diferentes geometrias na camada de mistura. Caso o efeito da modulação dos vórtices produzisse um padrão de ruído com características diferentes no emparelhamento, um controle ativo de escoamento por excitação periódica poderia ser usado para reduzir o ruído em jatos. O objetivo do presente trabalho é investigar tal efeito. A idéia é investigar este emparelhamento de vórtices na camada de mistura em desenvolvimento temporal bi-dimensional. Com isto foi possível visualizar um emparelhamento isolado de outros emparelhamentos e sem o efeito Doppler (presente na camada de mistura em desenvolvimento espacial). O método adotado foi a simulação numérica direta (DNS) das equações de Navier Stokes compressíveis na forma não-conservativa escritas na formulação característica. Os resultados mostram que a modulação dos vórtices não produz alteração significativa do ruído gerado no emparelhamento. / Turbo-fan and turbo-jet engines are the most important noise sources during the aircraft take off, climb and acceleration segments. Owing to this fact, the jet flow noise has been studied in the past years. For the landing stage, the slat is an important sound source. In this case, the slat leading edge frequently experiences a boundary layer deattachment causing the development of a mixing layer inside the slot. Nevertheless, there are many aspects of such phenomenon that have not been studied yet. Mixing layers constitutes an usual approach for jet flow instability in aeroacoustics studies. The stategy is to study the mixing layer in order to understand the jet-flow. This strategy becomes better as the ratio between the jet diameter and mixing layer thickness becomes larger. This approach is only reazonable for the jet flow axi-symetric unstable modes. The effect of vortex modulation on the vortex pairing sound production has not been found in the literature. If such effect could cause a significant change in the sound generation patterns, an active flow control system could be developed in order to enhance the jet noise performance. The purpose of the present work was to investigate such effect. It was also possible to observe a single vortex pairing inside a wide domain without the Doppler effect. The strategy was to study the vortex pairing in a bi-dimensional mixing layer under temporal development. The method used was the direct numerical simulation (DNS) of the compressible bidimensional (2D) Navier Stokes equations written in a nonconservative form of the characteristics formulation. The results showed that the vortex modulation did not produce a significant change on the vortex pairing sound.

A poluição sonora é um problema central de uma grande diversidade de aplicações industriais. Na engenharia, podemos citar diversos casos que geram ruído, como exemplos os trens, automóveis, rotores de helicópteros e o ruído aerodinâmico das aeronaves, o qual se divide em ruído gerado pelos motores a jato e a estruturas da aeronave. No presente momento o ruído dos motores aeronáuticos, principalmente os jatos, atingiu níveis de ruídos semelhantes às estruturas da aeronave, como por exemplo, eslates, flaps e trens de pouso. Desta forma, as autoridades de transporte aéreo estão exigindo também redução no ruído das estruturas. O presente trabalho apresenta a verificação das potencialidades e limitações do software CFD++, programa este adquirido pela EMBRAER e inserido como parte do projeto Aeronave Silenciosa, para assim poder compreender de uma melhor maneira o fenômeno da aeroacústica, e deste modo, poder contribuir para a redução do ruído externo das aeronaves. Para verificar as potencialidades e limitações do CFD++, foi proposto investigar o mecanismo de som do eslate. Tal fenômeno é devido ao deslocamento da camada limite no intradorso do eslate a partir de onde se desenvolve a camada de mistura, foco do presente trabalho. / Noise pollution is a central problem of a wide variety of industrial applications. In engineering, cite several cases that generate noise, as examples trains, automobiles, rotors of helicopters and the noise generated by aircraft, which is divided into noise generated by jet engines and airframe. At present the noise of aircraft engines, largely the jets reached noise levels similar structures, such as slat, flaps and landing gear. Thus, the air transport authorities are also demanding a reduction in noise of the structures of airframe. This report presents the verification of potentialities and limitations of CFD++, a program acquired by EMBRAER and inserted as part of the Silent Aircraft, so they can understand better how the phenomenon of aeroacoustics, and thus able to reduce contribute external noise from aircraft. To check the potentialities and limitations of CFD++, was proposed to investigate the mechanism of sound generated by the slat. This phenomenon is due to the displacement of the boundary layer on the lower surface of the slat from which the mixed layer develops. The mixing layer is the focus of this work.

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.

Tanto os motores turbo-jato quanto os turbo-fan, são os maiores responsáveis pela geração de ruído durante a decolagem, segmento de subida e de aceleração de uma aeronave. Devido a isto, o problema de ruído em jatos vem sendo intensamente investigado ao longo dos últimos anos. Já na fase do pouso, o slat é uma das fontes de ruído mais importantes. Para este caso, na maioria das aplicações práticas, existe o descolamento da camada limite no intradorso do slat a partir de onde se desenvolve uma camada de mistura. Ainda assim, existem inúmeros aspectos de tais escoamentos que precisam de investigação. Uma abordagem frequentemente feita para o estudo da instabilidade hidrodinâmica e ruído em jatos é o estudo de metade do jato. A estratégia consiste em estudar os fenômenos na camada de mistura, o que é uma aproximação razoável quando o jato tem diâmetro muito grande comparado à espessura da camada cisalhante que se desenvolve nas bordas do jato. Assim, alguns aspectos do ruído gerado pelos modos axi-simétricos de instabilidade são em grande parte reproduzidos. Um aspecto aparentemente jamais estudado antes é o efeito do emparelhamento de vórtices de diferentes geometrias na camada de mistura. Caso o efeito da modulação dos vórtices produzisse um padrão de ruído com características diferentes no emparelhamento, um controle ativo de escoamento por excitação periódica poderia ser usado para reduzir o ruído em jatos. O objetivo do presente trabalho é investigar tal efeito. A idéia é investigar este emparelhamento de vórtices na camada de mistura em desenvolvimento temporal bi-dimensional. Com isto foi possível visualizar um emparelhamento isolado de outros emparelhamentos e sem o efeito Doppler (presente na camada de mistura em desenvolvimento espacial). O método adotado foi a simulação numérica direta (DNS) das equações de Navier Stokes compressíveis na forma não-conservativa escritas na formulação característica. Os resultados mostram que a modulação dos vórtices não produz alteração significativa do ruído gerado no emparelhamento. / Turbo-fan and turbo-jet engines are the most important noise sources during the aircraft take off, climb and acceleration segments. Owing to this fact, the jet flow noise has been studied in the past years. For the landing stage, the slat is an important sound source. In this case, the slat leading edge frequently experiences a boundary layer deattachment causing the development of a mixing layer inside the slot. Nevertheless, there are many aspects of such phenomenon that have not been studied yet. Mixing layers constitutes an usual approach for jet flow instability in aeroacoustics studies. The stategy is to study the mixing layer in order to understand the jet-flow. This strategy becomes better as the ratio between the jet diameter and mixing layer thickness becomes larger. This approach is only reazonable for the jet flow axi-symetric unstable modes. The effect of vortex modulation on the vortex pairing sound production has not been found in the literature. If such effect could cause a significant change in the sound generation patterns, an active flow control system could be developed in order to enhance the jet noise performance. The purpose of the present work was to investigate such effect. It was also possible to observe a single vortex pairing inside a wide domain without the Doppler effect. The strategy was to study the vortex pairing in a bi-dimensional mixing layer under temporal development. The method used was the direct numerical simulation (DNS) of the compressible bidimensional (2D) Navier Stokes equations written in a nonconservative form of the characteristics formulation. The results showed that the vortex modulation did not produce a significant change on the vortex pairing sound.

A poluição sonora é um problema central de uma grande diversidade de aplicações industriais. Na engenharia, podemos citar diversos casos que geram ruído, como exemplos os trens, automóveis, rotores de helicópteros e o ruído aerodinâmico das aeronaves, o qual se divide em ruído gerado pelos motores a jato e a estruturas da aeronave. No presente momento o ruído dos motores aeronáuticos, principalmente os jatos, atingiu níveis de ruídos semelhantes às estruturas da aeronave, como por exemplo, eslates, flaps e trens de pouso. Desta forma, as autoridades de transporte aéreo estão exigindo também redução no ruído das estruturas. O presente trabalho apresenta a verificação das potencialidades e limitações do software CFD++, programa este adquirido pela EMBRAER e inserido como parte do projeto Aeronave Silenciosa, para assim poder compreender de uma melhor maneira o fenômeno da aeroacústica, e deste modo, poder contribuir para a redução do ruído externo das aeronaves. Para verificar as potencialidades e limitações do CFD++, foi proposto investigar o mecanismo de som do eslate. Tal fenômeno é devido ao deslocamento da camada limite no intradorso do eslate a partir de onde se desenvolve a camada de mistura, foco do presente trabalho. / Noise pollution is a central problem of a wide variety of industrial applications. In engineering, cite several cases that generate noise, as examples trains, automobiles, rotors of helicopters and the noise generated by aircraft, which is divided into noise generated by jet engines and airframe. At present the noise of aircraft engines, largely the jets reached noise levels similar structures, such as slat, flaps and landing gear. Thus, the air transport authorities are also demanding a reduction in noise of the structures of airframe. This report presents the verification of potentialities and limitations of CFD++, a program acquired by EMBRAER and inserted as part of the Silent Aircraft, so they can understand better how the phenomenon of aeroacoustics, and thus able to reduce contribute external noise from aircraft. To check the potentialities and limitations of CFD++, was proposed to investigate the mechanism of sound generated by the slat. This phenomenon is due to the displacement of the boundary layer on the lower surface of the slat from which the mixed layer develops. The mixing layer is the focus of this work.

Em uma aeronave, dispositivos de hiper-sustentação à altos ângulos de ataque promovem a formação de esteiras. Tais esteiras influenciam o desempenho aerodinâmico. A esteira do eslate, por exemplo, afeta o ponto de transição do elemento principal da asa e é a fonte de ruído mais importante na estrutura do avião. Devido a geração de sustentação estas esteiras são assimétricas. Corpos imersos em escoamento cisalhante também produzem esteiras assimétricas, com importantes aplicações para a indústria petrolífera. Existem aspectos de tais escoamentos que precisam de investigação. Por exemplo, não há consenso sobre se pode ou não a assimetria suprimir desprendimento de vórtice ou como afeta o desprendimento da frequência. Na verdade, existe uma aparente discrepância entre os resultados encontrados na literatura. O objetivo do presente trabalho é contribuir para esta questão. A idéia foi investigar a influência da assimetria em um perfil de esteira bidimensional sob desenvolvimento temporal. O perfil da esteira assimétrica foi obtido por uma combinação entre um perfil de uma esteira gaussiana e um perfil tangente hiperbólico da camada de mistura. Foi desenvolvido uma análise bidimensional (2D) da teoria de estabilidade linear para o perfil. O trabalho também incluiu simulações numéricas diretas (DNS) bidimensionais da evolução da perturbação usando formulação característica das equações compressíveis de Navier-Stokes na forma não-conservativa. Os resultados mostraram que a assimetria reduz o frequência de desprendimento para a faixa investigada do parâmetro de assimetria 0 \'< OU =\' K \'< OU =\' 0,25. Para pequenos valores deste parâmetro os resultados mostraram que a assimetria promove a estabilidade. Efeitos do número de Mach na esteira assimétrica também foram investigados. Os resultados mostraram que o aumento do Ma reduziu a máxima taxa de amplificação e a banda instável. Além disso, os resultados da literatura foram reconciliados, levando-se em conta a variação do coeficiente de arrasto com a assimetria. Portanto, uma possível explicação para a controvérsia foi oferecido. / In an aircraft, high-lift devices operating at high angle of attack promote the formation wakes. Such wakes influence the aerodynamic performance. The slat wake, for example, affects the transition point of the wing main element and is the most important source of noise in the airframe. Owing to the generation of lift these wakes are asymmetries. Bodies immersed in a shear flow also produce asymmetric wakes, with important applications to the oil industry. There are aspects of such flows that need investigation. For instance, there is no consensus about whether or not the asymmetry can suppress vortex shedding or how it affects the shedding frequency. Indeed there is an apparent discrepancy between results found in literature. The aim of the present work was to contribute to this issue. The idea was to investigate the influence of asymmetry on a two-dimensional wake profile under temporal development. The asymmetric wake profile was obtained by a combination between a Gaussian wake profile and a hyperbolic tangent mixing layer profile. The bidimensional (2D) linear stability theory analysis of the profile was performed. The work also included 2D Direct Numerical Simulation (DNS) using the characteristic formulation of the compressible Navier-Stokes equations in non-conservative form. The results showed that the asymmetry reduces the shedding frequency for the range of 0 \'< OU =\' K \'< OU =\' 0,25 investigated. For small values of the asymmetric parameter the results showed that the asymmetry promotes stability. Effects of the number of Mach in asymmetric wake also were investigated. The results showed that the increase of Ma reduced the maximum rate amplification and the band unstable. Moreover, the results from literature were reconciled by taking into account the variation of the drag coefficient with the asymmetry. Therefore, a possible explanation for the controversy was offered.

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.

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. / 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.

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.