O estudo do transporte de partículas tem uma importância fundamental em diversas áreas de pesquisas como, por exemplo, na formação de gelo em uma aeronave pois pode afetar a sua sustentação e estabilidade. Tamanha é a preocupação com a segurança de vôo que diversos estudos têm sido realizados, resultando em códigos computacionais como o LEWICE nos Estados Unidos, TRAJICE no Reino Unido, ONERA na França e CANICE no Canadá. No Brasil, um dos estudo é feito pela EMBRAER em parceria com algumas instituições. O objetivo deste trabalho é desenvolver um algoritmo que possa ser empregado na trajetória das partículas, utilizando uma formulação euleriana que elimina a dificuldade da semeadura de partículas específica da formulação lagrangiana na determinação da fração de volume da partícula. O método empregado é dos volumes finitos em malhas não-estruturadas cuja principal chave está na reconstrução mínimos quadrados de alta ordem com restrição nos contornos. O desenvolvimento do trabalho engloba 3 etapas: definição da geometria e geração das malhas; utilização de um solver para o tratamento do escoamento do ar e obtenção do campo de velocidade; implementação e utilização do esquema numérico com reconstrução mínimos quadrados de alta ordem para simular o cálculo da fração de volume com imposição de condições limites apropriadas no contorno do corpo. Os resultados dos testes realizados mostram que o esquema numérico com reconstrução mínimos quadrados pode ser empregado na resolução de equações que apresentam uma região de descontinuidade, como é o caso da região de sombra, reduzindo a largura da banda de difusão numérica e overshoots. / The particle transport study has a fundamental importance in diverse research area like in the icing accretion on an aircraft because that can affect its sustentation and stability. The concern is so big that many researches have been carried through, resulting in computational codes like the LEWICE in the United States, TRAJICE in the United Kingdom, ONERA in France and CANICE in Canada. In Brazil, one of the study has been made by the EMBRAER with some institutes. The goal of this work is to develop an algorithm that can be used in the particles trajectory study, using an Eulerian method that eliminates the difficulty particle sowing, particular of the Lagrangian method, in the determination of the droplet fraction volume. This is made by the finite volume method on unstructured meshes whose main key is the high order reconstruction with restriction on the boundary. The development of the work involves 3 stages: geometry definition and mesh generation; using code for the treatment of the air flow and obtained flow velocity; use of the high order numerical scheme least square reconstruction to simulate the droplet fraction volume result with imposition of appropriate limit conditions in the body contour. The realized simulations shown that Least Square method can be used in problem resolution that present descontinuos region like is shadow region reducing numerical diffusion and overshoots.

A aeroelasticidade estuda a interação mútua entre os efeitos aerodinâmicos e estruturais. É sabido que essa relação muitas vezes se comporta de maneira não linear, causando diversos problemas, tais como flutter, oscilações em ciclo limite, bifurcações e caos. Tais fenômenos são difíceis de serem diagnosticados, podendo causar problemas graves à estrutura das aeronaves e também inviabilizar as suas operações. Dentre as principais fontes de não linearidades em sistemas aeroelásticos, pode-se citar as de origem aerodinâmica e estrutural. As de origem estrutural, por sua vez, podem ter caráter distribuído ou concentrado. Sabe-se que os efeitos estruturais concentrados denominados enrijecimento e folga são os de maior impacto na aeroelasticidade não linear. Desse modo, o objetivo desse trabalho é estudar a interação não linear entre duas não linearidades estruturais, ou seja, o enrijecimento associado à rigidez em torção e a folga presente nas articulações das superfícies de controle de seções típicas aeroelásticas. Experimentos em túnel de vento são realizados utilizando um dispositivo que permite variar a intensidade do efeito de enrijecimento e do tamanho da folga na articulação da superfície de comando. O modelo numérico de seção típica aeroelástica também é utilizado e validado com dados experimentais. Análises por meio de diagramas de bifurcação de Hopf e técnicas baseadas em espectros de potência são utilizadas. Todas as respostas aeroelásticas foram caracterizadas através de ferramentas de análise nos domínios do tempo e da frequência, como técnica de reconstrução de espaço de estados e os espectros de alta ordem (HOS), os quais são importantes na identificação dos tipos de acoplamentos não lineares. Resultados indicam que a combinação dos efeitos de enrijecimento e folga são responsáveis pelo comportamento subcrítico das bifurcações de Hopf e que a intensidade do enrijecimento tem influência direta nas amplitudes de ciclo limite. / Aeroelasticity is the field of engineering that deals with the mutual interaction between the aerodynamic and structural dynamics effects. It is known that this relationship often shows nonlinear behavior, causing various problems such as flutter, limit cycle oscillations, bifurcations and chaos. Such phenomena are difficult to predict and can cause serious problems to the aircraft structure and also they can jeopardize their operations. The unsteady aerodynamic and structural dynamics provide the main sources of nonlinearities in aeroelastic systems. Structural nonlinearities can be treated as distributed or concentrated effects. It is know that the nonlinear concentrated structural effects referred as hardening and freeplay have a significant impact on nonlinear aeroelasticity. The objective of this work is to analyze an aeroelastic system under the influence of combined structural nonlinearities, i.e., the hardening nonlinearity in the pitch airfoil motion and the freeplay nonlinearity in the control surface hinge. Wind tunnel experiments are carried out using one device that allows to vary the intensity of the hardening effect and the size of the freeplay gap in the control surface hinge. The numerical model of the typical aeroelastic section is also used and validated with experimental data. All aeroelastic responses are characterized by analytical tools in time and frequency domains. It was used the state space reconstruction technique and the higher order spectral analysis (HOS) to identify types of nonlinear couplings. The results indicate that the combination of hardening and freeplay effects are responsible for inducing the subcritical behavior on the Hopf bifurcations and that the intensity of the stiffness has a direct influence on the limit cycle amplitudes.

Submitted by Maria Aparecida (cidazen@gmail.com) on 2017-04-17T11:57:59Z No. of bitstreams: 1 Alfredo Costa Oliveira Júnior.pdf: 789446 bytes, checksum: c5c9858983f5e6384177bda8d1ae2a0a (MD5) / Made available in DSpace on 2017-04-17T11:57:59Z (GMT). No. of bitstreams: 1 Alfredo Costa Oliveira Júnior.pdf: 789446 bytes, checksum: c5c9858983f5e6384177bda8d1ae2a0a (MD5) Previous issue date: 2017-02-16 / Atrial fibrilation(AF) is one of the most common cardiac arrhythmias worldwide. Thus, there are ample efforts to implement AF diagnosis systems. The main noninvasive way to assess cardiac health is through electrocardiogram (ECG) signal analysis, which represents the electrical activity of the cardiac muscle, and has characteristic temporal markings: P, Q, R, S and T waves. Some authors use filtering techniques, statistical analysis and even neural networks for detecting AF based on the RR interval, that is given by the temporal difference between the peaks of the R wave. However, analises of the RR interval allows for evaluating changes occurring only in the R wave of the ECG signal, not allowing to assess, for example, variations in the P wave provoked by the AF. In face of that, we propose characterize the ECG signal amplitude aiming at classifying both healthy and AF patients. The ECG signal was analyzed in the proposed methodology through the following statistics: variance, asymmetry, and kurtosis. Herein, we use the MIT-BIH Atrial Fibrillation and MIT-BIH Normal Sinus Rhythm database signals to evaluate AF and normal heartbeat intervals. Our study shown that kurtosis outperfomed variance and asymmetry with respect to sensibility (Se = 100%), specificity (Sp = 88.33%) and accuracy (Ac = 91.33%). The results were expected since kurtosis is a non-Gaussian measure and the ECG signal has sparse distribution. The proposed methodology also requires a lower number of pre-processing stages, and its simplicity allows for implementations in imbedded systems supporting the clinical diagnosis. / A Fibrilação atrial (FA) é uma das arritmias cardíacas mais comuns em todo o mundo. Por isso, amplos são os esforços para implementar sistemas que apoiem o diagnóstico de FA. A principal forma não invasiva de avaliar a saúde cardíaca, é através da análise do sinal de eletrocardiograma (ECG), o qual representa a atividade elétrica do músculo cardíaco, e possui marcações temporais características: as ondas P, Q, R, S e T. Alguns autores utilizaram técnicas de filtragem, análise estatística e até redes neurais para detectar FA com base no intervalo RR, que é dado pela diferença temporal entre os picos da onda R. Entretanto, a análise do intervalo RR permite avaliar apenas as variações que ocorrem na onda R do sinal de ECG, não permitindo avaliar, por exemplo, as alterações na onda P, provocadas pela FA. Diante disso, propõe-se caracterizar a amplitude do sinal de ECG, a fim de classificar pacientes com FA e saudáveis. Na metodologia proposta, o sinal de ECG, foi analisado por meio das seguintes estatísticas: variância, assimetria e curtose. Para avaliar o classificador proposto, usou-se sinais obtidos das bases de dados MIT-BIH Atrial Fibrillation e MIT-BIH Normal Sinus Rhythm referentes aos pacientes com FA e com ritmo cardíaco normal, respectivamente. Dentre as estatísticas analidadas, a curtose foi a que apresentou resultados superiores em termos de sensibilidade (Se = 100%), especificidade (Sp = 88, 33%) e acurácia (Ac = 91, 33%). Esses resultados são de se esperar pelo fato de que a curtose é uma medida de não-gaussianidade e que o sinal de ECG possui distribuição esparsa. A metodologia proposta também requer um número menor de etapas de pré-processamento, e sua simplicidade permite implementações em sistemas embarcados que apoiarão o diagnóstico clínico.

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.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / In this study the possible use of smart structures for noise control in a higher order acoustic duct was considered. The best option for this control was the use of axial splitters in the duct in order to prevent higher order mode propagation. It is possible to perform the active noise control in each splitter section by using a single channel control system. The use of smart structures takes advantage of the splitter plate and uses it as the control source, which substitutes the traditional loudspeakers used in active noise control systems. In order to evaluate the possibility of the noise control using smart structures, an analytical model of a thin plate with piezoelectric actuators was built then the acoustic field generated by this vibrating structure inside of the duct was obtained. However, to obtain the acoustic field inside an splitted duct, a numerical method such as the Component Mode Synthesis has to be used. Using the equation of the acoustic field generated in the duct by the plate, it was possible to obtain the acoustic field inside the splitted duct. After that, the active noise control simulations for harmonic excitations were performed and the influence of the size of the plate excited by the PZT actuators was studied. Finally the active control for random noise was simulated, in which the number of actuators in the plate was changed. In conclusion, it is possible to say that the smart structures can be used in active noise control of ducts with splitters and the advantages and disadvantages of the conveyed technique were presented. / Neste trabalho, foi estudada a proposta da utilização de estruturas inteligentes para o controle de ruído em um duto acústico com propagação de modos de alta ordem. A técnica mais adequada para este controle foi o particionamento do duto a fim de planificar as ondas que se propagam. Nesta região particionada, é possível realizar o controle ativo de ruído utilizando apenas um sensor e um atuador para cada lado da partição. A aplicação das estruturas inteligentes é proposta no sentido de aproveitar a placa particionadora para que, com a sua vibração, atuará como a fonte secundária necessária para o controle. Para a avaliação da possibilidade de controle utilizando esta técnica, primeiramente foi modelado o comportamento de uma placa instrumentada com atuadores piezoelétricos e, em seguida, obtida a modelagem analítica do campo sonoro gerado por uma estrutura vibrante no interior de um duto. Porém, a obtenção do campo acústico em um duto particionado não é facilmente obtido, sendo, então, realizada através da técnica de Síntese Modal de Componentes. Utilizando as equações do duto excitado por uma estrutura vibrante na técnica de síntese modal, foi possível obter campo acústico gerado no interior de um duto particionado. A partir disto, foram realizados simulações de controle ativo de ruído variando o trecho da placa a ser excitado para tons puros e para ruídos de banda estreita. Nesta última situação também foi avaliada a influência da quantidade de atuadores instalados. Concluiu-se deste trabalho que é possível a utilização de estruturas inteligentes no controle ativo de ruído em dutos particionados, sendo apresentadas suas vantagens e desvantagens. / Mestre em Engenharia Mecânica

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.

Submitted by Renata Lopes (renatasil82@gmail.com) on 2017-06-09T17:37:21Z No. of bitstreams: 1 fabriciomatosferreira.pdf: 457154 bytes, checksum: 39a00f23ce89d7948ef54eee13b6380a (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2017-06-29T12:08:53Z (GMT) No. of bitstreams: 1 fabriciomatosferreira.pdf: 457154 bytes, checksum: 39a00f23ce89d7948ef54eee13b6380a (MD5) / Made available in DSpace on 2017-06-29T12:08:53Z (GMT). No. of bitstreams: 1 fabriciomatosferreira.pdf: 457154 bytes, checksum: 39a00f23ce89d7948ef54eee13b6380a (MD5) Previous issue date: 2015-02-27 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Teorias com derivadas de alta ordem e sua invariância sob transformações conformes são ferramentas poderosas em Teoria Quântica de Campos em espaços curvos. O principal propósito de nosso estudo é trabalhar no sentido de obter uma generalização do operador quadridimensional conforme de Paneitz para um operador análogo em seis dimensões (∆6). A consecução de um operador conformalmente invariante contendo seis derivadas justifica a escolha de se trabalhar nesta dimensão como um recurso mais simples para sua obtenção. Tal operador seria de fundamental importância na integração da anomalia conforme em teorias com derivadas de sexta ordem. A forte relação entre operadores conformes e a densidade de Euler foi investigada em D = 2 e D = 4. O termo da densidade de Euler - o qual é conhecido em 4D como termo topológico de Gauss-Bonnet - foi estendido a seis dimensões. Em dimensões pares, tal termo é representado por potências de D/2 de contrações dos tensores de Riemann e Ricci, cada um dos quais contendo D derivadas. A transformação dos escalares de terceira ordem na curvatura os quais aparecem no estudo da anomalia conforme e da densidade de Euler em seis dimensões foi obtida como um primeiro passo rumo ao nosso objetivo de obter ∆6. / Higher-derivative theories and their invariance under conformal transformati ons are powerful tools in Quantum Field Theory in curved spaces. The main purpose of our study is to work towards achievement of a generalization of the four dimensional Paneitz conformal operator to analogous six dimension one (∆6). Achieving a six deri vative conformally invariant operator justifies the choice of working in such dimension as a simpler resource for obtaining it. This operator would be of fundamental account in the integration of the conformal anomaly in sixth order derivatives theories. The strong relation between conformal operators and Euler density has been investigated in D = 2 e D = 4. Euler density term - which 4D version is known as topological Gauss-Bonnet term - was extended to the six dimension. In even dimensions, it is represented by D/2 powers of contracted Riemann and Ricci tensors, each of them carrying D derivatives. The transformation of the third order curvature scalars that arise in the study of the conformal anomaly and Euler density in six dimension was obtained as a first step to our goal of achieving ∆6.

Uma grande limitação dos métodos de diferenças finitas é que eles estão restritos a malhas e domínios retangulares. Para descrever escoamentos em domínios complexos, como, por exemplo, problemas com superfícies livres, faz-se necessário o uso de técnicas acessórias. O método de interfaces imersas é uma dessas técnicas. Nesse trabalho, primeiramente foi desenvolvido um método de projeção, totalmente livre de pressão, para as equações de Navier-Stokes com variáveis primitivas em malha deslocada. Esse método é baseado em diferenças finitas compactas, possuindo segunda ordem temporal e quarta ordem espacial. Esse método foi combinado com o método de interface imersa de Linnick e Fasel [2] para resolver numericamente as equações de Stokes com quarta ordem de precisão. A verificação do código foi feita por meio do método das soluções manufaturadas e da comparação com resultados de outros autores em problemas clássicos da literatura. / A great limitation of finite differences methods is that they are restricted to retangular meshes and domains. In order to describe flows in complex domains, e.g. free surface problems, it is necessary to use accessory techniques. The immersed interface method is one of such techniques. In the present work, firstly, a projection method was developed, which is completely pressure-free, for the Navier-Stokes equations with primitive variables in a staggered mesh. This method is based on compact finite differences, with temporal second-order precision and spatial foruth-order precision. This method was combined with the immersed interface method from Linnick e Fasel [2] in order to numerically solve the Stokes equations with fourth-order precision. The verification of the code was performed with the manufactured solutions method and by comparing results with other authors for some classical problems in the literature.

O objetivo desse trabalho consiste no estudo de métodos de identificação de parâmetros em problemas envolvendo a equação de advecção-difusão 2D. Essa equação é resolvida utilizando o método dos volumes finitos, sendo empregada métodos de reconstrução de alta ordem em malhas não-estruturadas de triângulos para calcular os fluxos nas faces dos volumes de controle. Como ferramenta de busca dos parâmetros é empregada a técnica baseadas em gradientes, sendo os mesmos calculados utilizando processos baseados em métodos adjuntos. / The aim of this work concern to study parameter identification methods on problems involving the advection-diffusion equation in two dimensions. This equation is solved employing the finite volume methods, and high-order reconstruction methods, on triangle unstructured meshes to solve the fluxes across the faces of control volumes. As parameter searching tool is employed technicals based on gradients. The gradients are solved using processes based on adjoint methods.

Em diversos escoamentos sobre superfícies há a presença de protuberâncias, como por exemplo rebites, parafusos e juntas. Estas protuberâncias podem influenciar a camada limite, acelerando a transição do escoamento do estado laminar para o estado turbulento. Em alguns casos isto pode ser indesejável, já que o escoamento turbulento implica necessariamente em uma força de atrito maior do que aquela referente ao escoamento laminar. Existem alguns aspectos neste tipo de escoamento que ainda não estão bem compreendidos. O objetivo deste trabalho é estudar a influência de uma rugosidade isolada no escoamento sobre uma superfície. Este estudo contribui para se entender o que ocorre em casos de maior complexidade. O estudo é de natureza computacional, em que se utiliza simulação numérica direta das equações de Navier-Stokes. A técnica de fronteiras imersas é utilizada para representar a rugosidade no escoamento sobre a superfície. O código numérico é verificado por meio do método de soluções manufaturadas. Comparações entre resultados experimentais, da teoria de estabilidade linear e numéricos também são utilizados para a validação do código. Resultados obtidos com diferentes alturas de rugosidade e variações no gradiente de pressão permitiram analisar a influência de elemento rugoso tridimensional em escoamentos de camada limite. / The presence of protuberances on surfaces, for example, rivets, screws and gaskets, can influence the boundary layer by accelerating the transition from laminar flow to turbulent flow. In some cases this may be undesirable, since the turbulent flow involves frictional forces greater than the ones at the laminar regime. There are some aspects of the flow in the boundary layer perturbed by a single roughness element that are not well understood. The aim of this work is to study the influence of an isolated roughness on the boundary layer. This study is a step towards to the understanding of what can happen in more complex cases. The nature of this study is computational, therefore a Direct Numerical Simulation code is used. The immersed boundary method is used to represent the roughness in the flow on the surface. The numerical code is verified via theMethod ofManufactured Solutions. Comparisons between experimental data, Linear Stability Theory and numerical results are also used for the validation of the code. Results obtained with different roughness heights and variations in the pressure gradient allowed the analysis of the influence of a three-dimensional roughness element in boundary layer flows.