Micrositing é o campo do posicionamento estratégico dos aerogeradores na área do parque eólico visando a con guração mais promissora em termos econômicos ou de produção. Dado que em terrenos complexos as características do vento variam na área do parque eólico de forma não-linear, emprega-se a modelagem numérica do vento por CFD para extrapolar os dados medidos para toda a região. O presente trabalho consiste no desenvolvimento de uma metodologia de micrositing em terrenos complexos capaz de auxiliar no projeto do layout e seleção do tipo e altura de eixo do aerogerador que maximiza o fator de capacidade (FC) utilizando softwares livres de código aberto. A metodologia consiste na simulação do vento para várias direções de incidência, assimilação das medições, convers ão de velocidade em densidade de potência, ponderação por frequência de ocorrência de cada direção, sobreposição, seleção das coordenadas dos aerogeradores e cálculo do FC para diversas con gurações de tipo e altura de eixo dos aerogeradores. Veri cação, validação e seleção das constantes do modelo de turbulência é realizada anteriormente às simulações Veri cou-se que o modelo k- produziu um escoamento horizontalmente homog êneo e que o melhor desempenho na validação foi obtido com a escolha de constantes para escoamentos atmosféricos. A metodologia foi demonstrada em uma região de terreno complexo em que o FC do parque eólico proposto apresentou caráter convergente com o re no progressivo da malha, porém oscilatório em termos do número de direções simuladas. Por m, obteve-se FC brutos superiores a 40% para as cinco melhores con gurações e de aproximadamente 52% no melhor caso, indicando bom potencial eólico. A metodologia foi capaz de preencher uma lacuna na literatura cientí ca de micrositing ao possibilitar o planejamento do layout, tipo de aerogerador e altura de eixo, bem como a estimativa da produção e FC brutos da usina em terrenos complexos. Além disso, a estrutura de trabalho com o uso de recursos computacionais livres e de código aberto reforça o caráter de desenvolvimento contínuo, compartilhamento e transparência da metodologia. / Micrositing is the eld concerned with the strategic positioning of wind turbines in the wind farm area aimed at the most promissing con guration economically- or yield-wise. Given the wind characteristics vary non-linearly across the wind farm area in complex terrain, numerical wind modeling with CFD is employed to extrapolate the measured data to the whole site. The present work consists of the development of a micrositing methodology in complex terrain capable of assisting the layout project and selection of wind turbine type and hub height that maximizes the capacity factor (CF) using free and open-source software. The methodology consists of simulating the wind for a number of incoming directions, assimilation of measurements, conversion of wind speed into power density, weighing by frequency of occurrence of each direction, overlapping, selection of wind turbine coordinates and CF calculation for a number of wind turbine types and hub heights. Veri cation, validation and selection of turbulence model constants is performed previous to the simulations It was veri ed that the k- model is able to sustain horizontally-homogeneous ow and that the classic atmospheric ow constants performed best in the validation step. The methodology was demonstrated in a complex terrain region for which the proposed wind farm CF showed converging behavior with progressive mesh re nement, however oscillating with the number of wind directions simulated. Ultimately, CF greater than 40% were obtained with the ve best performing con gurations and approximately 52% in the best case scenario, suggesting good wind potential. The methodology was capable of lling a major gap in the scienti c literature of micrositing for allowing the layout planning, selection of wind turbine type and hub height, as well as gross production estimates and CF for the wind farm in complex terrain. Additionally, the free and open-source-based framework strengthens the continuous development, sharing and transparency of the methodology.

O presente trabalho apresenta um estudo numérico sobre o desempenho de trocadores de calor com microcanais. Primeiramente, para a verificação da ferramenta numérica, uma comparação com um estudo experimental conhecido é realizado. Em seguida, aperfeiçoamentos disponíveis na literatura são testados e comparados. As simulações numéricas são realizadas com o software Fluent utilizando o Método dos Volumes Finitos (MVF). O resultado numérico da resistência térmica do dissipador do estudo experimental em que a primeira parte desse trabalho foi baseada ficou em 0,097 °C/W (com uma temperatura máxima de 373 K), enquanto o resultado experimental ficou em 0,090 °C/W, o que representa uma diferença de 7,2 %. São realizados aperfeiçoamentos integrando o trocador de calor a um sistema de refrigeração. Isso reduz a temperatura de entrada e, consequentemente, a temperatura máxima encontrada no sistema (350 K para uma mesma potência dissipada) não havendo mudança na resistência térmica, porém trazendo um ganho em confiabilidade. São apresentadas mudanças na geometria do dissipador diminuindo a área de entrada, com uma sensível piora nos resultados da resistência térmica, fazendo com que a resistência aumentasse de 0,097 °C/W para o canal original para 0,272 °C/W para o canal com menor altura. Apresenta-se um trocador de calor cujos microcanais em suas dimensões longitudinais possuem ondulações, porém os efeitos esperados de melhor mistura e maior área de contato para o fluxo de calor não foram observados para as condições de contorno do experimento original; a resistência térmica calculada foi de 0,102 °C/W. A última análise apresenta uma idealização onde o fluxo de calor é distribuído uniformemente no volume de silício, encontrando-se uma resistência térmica de 0,084 °C/W. / This paper presents a numerical study on the performance of heatsinks with microchannels. For the validation of the numerical tool, a comparison with a known experimental study is carried out, then improvements available in the literature on the subject are tested and compared. The numerical simulations are performed with the Fluent software using the Finite Volumes Method (MVF). The numerical result of the thermal resistance of the heatsink from the experimental study which the first part of this paper is based was 0.097 °C/W, whereas the experimental result was 0.090 °C/W (with a maximum temperature of 373 K), representing a difference of 7.2%. An improvement is performed by integrating the heatsink to a HVAC (Heating, ventilation and air conditioning) system, which lowers the inlet temperature and consequently lowers the highest temperature found in the system (350 K for the same dissipated power) without changing the thermal resistance, but providing gains in reliability. Changes in the heatsink geometry are presented by decreasing the inlet area, with a significant worsening in the thermal resistance results, causing the resistance to vary from 0.097 °C/W for the original channel to 0.272 °C/W for the channel that has the lowest height. A heatsink whose microchannels have waves in their longitudinal dimensions is shown, but the expected effects of better mixing and greater contact area for the heat flux were not observed for the boundary conditions of the original experiment, the calculated thermal resistance was of 0.102 °C/W. The last analysis presents an idealization where the heat flux is evenly distributed in the silicon volume, with a thermal resistance of 0.084 °C/W.

Investiga-se a absorção de ondas eletromagnéticas por plasmas imersos em campos magnéticos retilíneos e inomogêneos na direção perpendicular. O formalismo empregado está dentro do limite de validade do método WKB e se baseia no conceito de tensor dielétrico efetivo, o qual é definido através de correções realizadas em um tensor dielétrico derivado a partir de uma aproximação de ondas planas, de tal forma que a equação de conservação de energia para um plasma fracamente inomogêneo seja automaticamente satisfeita. O tensor dielétrico efetivo obtido por meio deste procedimento descreve corretamente os processos de transferência de energia entre ondas e partículas em um plasma inomogêneo, o que não acontece com o tensor dielétrico obtido diretamente da aproximação de ondas planas, sem incluir as correções. Para isolar o efeito do gradiente do campo, calcula-se o tensor dielétrico em uma situação particular onde a função distribuição de equilíbrio pode ser escrita em termos das mesmas constantes de movimento existentes para um plasma homogêneo. Mostrase que o tensor efetivo obedece às relações recíprocas de Onsager e que somente as partículas que estão em ressonância com as ondas eletromagnéticas contribuem para a sua parte antiHermiteana, ligada à absorção destas ondas. Comparando com outros trabalhos que surgiram na literatura recente, mostra-se que os resultados destes trabalhos podem ser reproduzidos a partir do tensor dielétrico não corrigido, cuja parte anti-Hermiteana possui a contribuição de partículas não-ressonantes. Conclui-se, portanto, que os processos de troca de energia entre as ondas e o plasma não são adequadamente descritos nestes trabalhos. Como aplicação numérica, o tensor dielétrico efetivo é empregado na solução da relação de dispersão de modo a se obter a absorção de ondas que se propagam nos modos ordinário e extraordinário em freqüências próximas à freqüência de cíclotron eletrônica. Em particular, calcula-se a absorção local e a profundidade óptica do modo ordinário para dois conjuntos de parâmetros físicos, correspondentes a tokamaks de grande e pequena razões de aspecto. Observa-se que uma região de absorção negativa surge em freqüências próximas à freqüência de cíclotron para valores suficientemente altos do gradiente do campo. / We investigate absorption of electromagnetic waves by plasmas in ambient magnetic fields with perpendicular gradients. The basic formalism is within the range of validity of the WKB method, and is based on the concept of effective dielectric tensor, defined by the corrections dane on a dielectric tensor derived from a plane wave approximation, in such a way that the energy conservation equation for an inhomogeneous plasma is automatically satisfied. The effective dielectric tensor derived from this procedure correctly describes the energy exchange processes between waves and particles in an inhomogeneous plasma. This is not the case for the dielectric tensor derived from the plane wave approximation, without including the corrections. In arder to isolate the effect of the field gradients, we calculate the dielectric tensor for a particular situation where the equilibrium distribution function depends only on the same constants of motion that appear in the homogeneous case. We show that the effective dielectric tensor obeys the Onsager reciprocity relations and that only the particles that resonate with the waves contribute to its anti-Hermitian part, which is related to the absorption of the waves. Comparing with other works in the recent literature, we show that the results of these works can be reproduced from the non-corrected dielectric tensor, whose anti-Hermitian part has contributions from non-resonant particles. We conclude therefore that the energy exchange processes between wave and plasma are not correctly described in these works. As a numerical example, we employ the effective dielectric tensor to solve the dispersion relation in arder to calculate the absorption of ordinary and extraordinary waves in frequencies dose to the electron cyclotron frequency. Particularly, we obtain the local absorption and the optical depth of the ordinary mode for physical parameters corresponding respectively to tokamaks with large and small aspect r a tio. We find ·that a region of nega tive absorption arises near the cyclotron frequency for sufficiently high magnetic field gradients.

A estimação de parâmetros de equações constitutiva viscoelásticas, tipicamente efetuada através de diferentes tipos de reômetros, pode ser uma tarefa árdua, especialmente no que diz respeito à estimação de parâmetros não lineares. Nas últimas décadas houve grande progresso no desenvolvimento de novas metodologias para caracterização de fluidos viscoelásticos na região de comportamento viscoelástico não linear e na simulação de escoamentos de fluidos viscoelásticos nesta região. Este trabalho tem como objetivo apresentar uma metodologia de estimação de parâmetros reológicos de equações constitutivas viscoelásticas baseada na comparação recursiva entre valores originados em experimentos com a modelagem em fluidodinâmica computacional do escoamento de referência. Esta metodologia apresenta como vantagem a possibilidade de estimar parâmetros a partir de dados originados em escoamentos mais representativos dos processos industriais, resultando em parâmetros com maior potencial de representar adequadamente o comportamento viscoelástico de um dado material nas suas condições típicas de aplicação. Dois estudos de caso foram implementados de forma a testar a metodologia: o escoamento de um fluido newtoniano entre placas paralelas e o escoamento de um fluido viscoelástico na saída do capilar. A estimação da viscosidade newtoniana foi baseada na comparação entre valores analíticos e numéricos para o gradiente de pressão ao longo do escoamento. A estimação do parâmetro não linear da equação constitutiva viscoelástica de Phan-Thien-Tanner (PTT) foi baseada na diferença entre os perfis de diâmetro numéricos com relação aos valores experimentais para o fenômeno de inchamento do filete de fluido ao emergir da saída do capilar, conhecido como inchamento do extrusado. O erro final obtido para o valor estimado da viscosidade newtoniana foi de 9 % com relação ao valor de referência. Para o parâmetro não linear do modelo PTT o erro para o valor estimado foi de 0,75 % com relação ao valor de referência obtido pela técnica de reometria rotacional. / The estimation of viscoelastic constitutive equations parameters, typically performed through different types of rheometers, can be an arduous task, especially as regards the estimation of nonlinear parameters. In recent years, there has been great progress in developing new methods for the characterization of viscoelastic fluids beyond the linear region and in the simulation of viscoelastic fluid flows in general. This work aims to present a methodology for estimation of rheological parameters of viscoelastic constitutive equations. The methodology is based on recursive comparison between experimentally measured fields with numerically computed ones by a computational fluid dynamics (CFD) model. This methodology has the advantage of estimating parameters through data coming from more representative flows of industrial processes, resulting in parameters with higher potential to represent the viscoelastic behavior of a given material under its typical conditions of application. Two case studies have been used to test this methodology: a steady state shear flow between two parallel plates of a Newtonian fluid and a viscoelastic fluid flow in a capillary. The estimation of Newtonian viscosity was based on the comparison between analytical and numerical values for the pressure gradient along the flow. The estimation of the nonlinear parameter of the viscoelastic constitutive equation Phan-Thien-Tanner was based on the difference between the experimental and numerical diameter profiles of the swelling of the fluid emerging from the capillary opening, known as extrudate swell. The final error obtained for the estimated value of Newtonian viscosity was 9 % compared to the reference value. For nonlinear parameter of PTT model, the error of the estimated value was 0.75 % compared to the reference value obtained through rotational rheometer technique.

Este trabalho apresenta um método numérico para a solução de escoamentos bi e tridimensionais de fluidos (quase) incompressíveis em torno de geometrias arredondadas. O escoamento bidimensional é analisado em torno da geometria de um cilindro (seção de um cilindro), para as equações de Euler e Navier-Stokes, e em torno da geometria aproximada de um tubarão para as equações de Euler. O escoamento tridimensional é analisado em torno de uma esfera e de um elipsóide. O método de integração empregado baseia-se no esquema explícito de Runge-Kutta de três estágios para as equações da quantidade de movimento e no de Relaxações Sucessivas para a pressão. Adota-se o esquema em diferenças finitas visando aproximações de segunda ordem no tempo e no espaço no sistema de coordenadas generalizadas. Testes numéricos são realizados para as diferentes geometrias aplicando as equações de Navier-Stokes e Euler e os resultados obtidos comparam adequadamente com dados analíticos, experimentais e/ou numéricos encontrados na literatura. / This work presents a numerical method for the solution of (almost) incompressible bi and tridimensional fl.ows for round geometries. Bidimensional fl.ows over a circular cylinder, using Euler and Navier-Stokes equations, and also for a shark approximated geometry, using Euler equations, are analyzed. Extension to tridimensional flows around a sphere and an elliptical geometry is realized. The integration method is based on the three-stage Runge-Kutta explicit scheme for momentum equations and successive under relaxation for pressure. Second order finite difference approximations for time as well as space terms in boundary fitted coordinates are employed. Numerical tests are carried out for different geometries for Euler and Navier-Stokes equations and the results showed to compare properly with analytical, numerical or experimental data found in the literature.

Neste trabalho faz-se uma análise de fluxos unidimensionais usando as equações de Burgers e de Euler. Para esta última, é obtida a solução exata e uma aproximação desta via método numérico. A obtenção da solução exata é baseada na combinação de ondas simples (uma onda de choque, uma descontinuidade de contato e uma onda de expansão) e na validade das relações de salto para as equações de Euler. Os resultados assim obtidos são utilizados para verificar (certificar) os resultados numéricos. As equações de Euler são integradas no tempo através de um esquema simplificado de Runge-Kutta; considera-se também a adição explícita de termos dissipativos ao esquema de discretização espacial. Sã.o apresentadas comparações entre as soluções exata e numérica, além de comparações da solução numérica para diferentes valores dos coeficientes de dissipação. Analisa-se também as regiões de estabilidade de métodos de Runge-Kutta para uma equação modelo, cujas propriedades são semelhantes às das equações de Burgers e de Euler. Por fim , propõe-se o estudo da convergência de um esquema semelhante ao de Runge-Kutta; faz-se uma estimativa de erro a posteriori em espaços de Banach de dimensão infinita. Além disto, são calculadas algumas estimativas a priori para a equação de Burgers que são usadas, juntamente com idéias da teoria de "shadowing" para estabelecer estimativas relativas à validade de simulações numéricas para a equação de Burgers. Além disto, são mostrados alguns estudos computacionais relevantes sobre a separação espectral, os quais poderiam ser entendidos como uma forma de estimativas a posteriori. / In this work we analyze unidimensional flows using Burgers and Euler equations. We show how to obtain both an exact and an approximate solution to the Euler equations. The exact solution is obtained by simple waves interaction (a shock wave, a contact discontinuity and an expansion fan). It is based on the jump relations for the Euler equations. \Ve use the exact solution thus obtained as a mean of comparison to the approximate (numerical) solution. A simplified Runge-Kutta method is used to integrate the Euler equations in time. Dissipation terms are added to the spatial discretized equations. vVe study the effect of these terms by plotting the exact versus the approximate solution considering different dissipation coefficients. We choose a model equation, which has similar properties to the Burgers and Euler equations, in order to study the stability regions of the Runge-Kutta scheme. Finally, we propose a convergence analysis of a Runge-Kutta-like scheme and make a posteriori error estimates in infinite dimensional Banach spaces. Furthermore, we perform some a priori estimates for the Burgers' equation and use them together with ideas of the Shadowing Theory in order to establish estimates concerning the validity of numerical simulations for Burgers' equation. As well we show a few relevant computer studies on the spectral separation which could be regarded as a form of a posteriori estimates.

O estudo do comportamento reológico de fluidos não Newtonianos tem grande importância em diversas áreas da engenharia. O aumento na demanda destes fluidos por exemplo, no uso doméstico, pessoal e em processos químicos acarreta em dificuldades que vão desde o processo de sua mistura ao seu manuseio. Dentre os fluidos não Newtonianos estão os viscoelásticos, os quais exibem deformação aparente quando os níveis de tensões são inferiores ao limite de escoamento do material e, dentro desta classe, alguns ainda podem apresentam comportamento elástico quando submetidos à baixa taxa de cisalhamento. Juntamente com os efeitos elasto-viscoplásticos, os materiais podem apresentar comportamento tixotrópico, onde devido as tensões sua reestruturação não é instantânea. Na presente Tese, fez-se um estudo numérico a fim de analisar o problema especifico de escoamentos de fluidos elasto-viscoplásticos com comportamento tixotrópico em uma expansão planar abrupta na razão de aspecto 4 : 1. O modelo mecânico aplicado consiste de uma equação viscoelástica para o campo de tensões, uma evolutiva para o parâmetro de estrutura do material, bem como as equações de conservação de massa e momento. O modelo mecânico aplicado mostrou-se capaz de prever o comportamento tixotrópico. A aproximação numérica do modelo aplicado foi feita através do método estabilizado de elementos finitos, especificamente o método Galerkin Mínimos-Quadrados (GLS), o qual foi implementado no código de elementos finitos para fluidos não Newtonianos em desenvolvimento no Laboratório de Mecânica dos Fluidos Aplicada e Computacional (LAMAC) da Universidade Federal do Rio Grande do Sul. Os fenômenos reológicos presentes no problema foram analisados a partir da influência da cinemática do escoamento, da elasticidade e da tixotropia, no nível de estruturação do material, na posição e tamanho das regiões não escoadas e na deformação elástica do material. Os resultados mostraram-se satisfatórios, pois condizem com os apresentados na literatura, apontando a boa predição do modelo mecânico aplicado bem como a robustez de sua implementação computacional. / The rheological behavior of non-Newtonian uids study is of great then importance in many areas of engineering. The increase in demand of these uids - for example, domestic use, personal and processes chemical - causes di culties ranging from the process of mixing it to handling. Among the non-Newtonian uids are viscoelastic, which exhibit apparent deformation when stress levels are lower than the yield limit of the material and, within this class, some still have elastic behavior when subjected to low shear rates. Together with the elastic-viscoplastic e ects, materials may exhibit thixotropic behavior, ie, due the restructuring strain is not instantaneous. In this thesis was made a numerical study to simulate the speci c problem of ow of elastic-viscoplastic uids with thixotropic behavior in abrupt planar expansion { common geometry in industrial systems associated with elastic-viscoplastic uids, whose ratio the aspect is 4:1. The mechanical behavior of most of these structured materials, are highly non-Newtonian, with this, there is need to make them more pseudoplastic, causing undesirable behaviors such as thixotropic which is a phenomenon somewhat characterized and modeled in the literature. The mechanical model applied is able to predict thixotropic behavior and is composed of a viscoelastic equation for the stress eld and an evolutionary to the material structure parameter in addition to the mass and momentum conservation equations. This mechanical model is approximated by a stabilized nite element model, called the Galerkin method of least squares (GLS). In order to study the rheological phenomena present in the problem is analyzed the in uence of the ow kinematics, elasticity and thixotropy in the level of structure of the material, in position and size of unyielded regions and the elastic deformation of the material. The results were satisfactory, since the the study of intensity the U the results agree with those reported in the literature pointing to good prediction of the mechanical model applied well as the hardiness of their computational implementation.The results of the elasticity showed quite spectacular behavior of unyielded regions of the material since for high values of relaxation time, the unyielded region of channel larger has the form of " ngers"because of the high exibility of the material along the line of symmetry of the channel. The results of the analysis of the elastic deformation show, that model correctly dosing the elasticity in the unyielded regions. The thixotropic e ects reported a slower structuring of the material with increasing characteristic time, in response to strain change caused by the expansion, results indicate that for high values of relaxation time, higher the distance between the unyielded regions of smaller channel and larger.

Analisam-se as instabilidades eletrostáticas ocasiona das pela injeção de um feixe de íons num plasma magnetizado, in finito, homogé- neo e sem colisões, como um possível mecanismo de aquecimento do plasma. A direção do feixe é suposta formar um ângulo arbitrário com relação ao campo magnético. As razões de cres cimento das instabilidades são calculadas para diversas freqbencias naturais de oscilação do plasma, discutindo-se as restrições que essas instabilidades impõem ao ingulo de injeção e i densidade do feixe. Para freqbencias próximas i freqbencia hibri da inferior, analisam-se os modos íon-íon transversal ao campo e a instabilidade das duas correntes modificada. Para freqüencias na região da fregüéncia acústica iOnica, consideram-se quatro ca sos limites do modo íon-acústico, com distintas características. Estuda-se também um modo íon-acústico modificado, cujas ondas, para parãmetros adequadamente escolhidos, possuem velocidade de fase comparável a das ondas acústicas lentas. As equações dos momenta da teoria quasilinear são usa das para o cálculo da variação temporal do momentum e da energia térmica de cada componente do plasma, na presença das instabili dades. Os Tons do plasma são aquecidos pela maior parte das on das, excessões feitas is instabilidades no modo das duas correntes modificado e is instabilidades em um dos modos íon-acústicos. / Electrostatic instabilities excited by an energetic ion beam injected at an arbitrary angle with respect to the magnetic field in a infinite,homogeneous, collisionless plasma are examined as a possible plasma heating mechanism. The instability growth rates are calculated for several natural oscillation frequencies of the plasma, and threshold conditions on the angle of injection and on the beam density are derived for the instabilities. For frequencies near the lower-hybrid frequency, the cross-field ion-ion mode and the modified two-stream instability are analized. When the frequency is Glose to the ion-acoustic frequency, four limiting cases of the ion-acoustic mode, with distinct characteristics, are considered. A modified ion-acoustic mode, which propagates with a phase velocity comparable to the phase velocity of the slow acoustic waves, for appropriately chosen parameters, is also derived. The quasilinear moment equations are used to calculate the time evolution of the momentum and the thermal energy of each plasma component in the presence of the instabilities. The background ions are heated by most of the instabilities, exceptions being the modified two-stream mode and one of the ion-acoustic modes.

O presente trabalho apresenta um estudo sobre ondas de gravidade. São descritas as equações que governam o movimento de ondas lineares, fracamente não-lineares e completamente não-lineares. Foi implementado um método espectral que resolve, empregando o método de Newton, as equações completamente não-lineares que descrevem o movimento das ondas estacionárias em água. Com isso, é possível simular o empinamento das ondas, ou seja, o comportamento destas quando a profundidade da água diminui. Calcula-se a altura das ondas em função da profundidade, o coeficiente de empinamento e o ângulo de quebra das ondas. Assim, pôde-se observar as relações destas propriedades com a quebra das ondas e associar os resultados numéricos com aqueles previstos pela teoria. / A study on the propagation of gravity waves is presented. The governing equations for linear, weakly nonlinear and fully nonlinear water waves are described. A spectral method has been implemented which solves, employing the Newton method, the fully nonlinear equations describing the motion of steady waves. It is then possible to simulate the shoaling of waves, that is, the behaviour of them when the water depth decreases. The wave height as a function of depth, the shoaling coefficient and the angle of wave breaking are calculated. It was possible to observe the relationship of these properties with the breaking of the waves and associate the numerical results with those predicted by the theory.

Tese (doutorado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Elétrica, 2014. / Submitted by Ana Cristina Barbosa da Silva (annabds@hotmail.com) on 2015-02-20T17:57:17Z No. of bitstreams: 1 2014_ViniciusdeCarvalhoRispoli.pdf: 16939736 bytes, checksum: 8dfa17c58667403d827192cadef858c7 (MD5) / Approved for entry into archive by Raquel Viana(raquelviana@bce.unb.br) on 2015-04-29T21:10:38Z (GMT) No. of bitstreams: 1 2014_ViniciusdeCarvalhoRispoli.pdf: 16939736 bytes, checksum: 8dfa17c58667403d827192cadef858c7 (MD5) / Made available in DSpace on 2015-04-29T21:10:38Z (GMT). No. of bitstreams: 1 2014_ViniciusdeCarvalhoRispoli.pdf: 16939736 bytes, checksum: 8dfa17c58667403d827192cadef858c7 (MD5) / As doenças cardiovasculares representam hoje uma das maiores causas de morte no mundo. Desta forma, o conhecimento dos padrões de escoamento sanguíneo no corpo humano se faz importante no diagnóstico e pesquisa de algumas doenças. Utilizando a ressonância magnética (RM), o escoamento sanguíneo pode ser mensurado in vivo diretamente usando ressonância magnética (RM) com contraste de fase (CF) ou através da codificação de velocidade em Fourier (CVF). Por outro lado, o escoamento sanguíneo pode também ser simulado utilizando uma abordagem baseada em modelos da dinâmica dos fluidos computacional (DFC). O CF tem como desvantagens os longos tempos de aquisição, a resolução espaço-temporal limitada, efeitos de volume parcial e baixa relação sinal-ruído (SNR). Por outro lado, a codificação de velocidade em Fourier é um método de imageamento por ressonância magnética promissor para realizar medidas do escoamento cardiovascular. A CVF oferece SNR consideravelmente maior que o CF e também é robusto aos efeitos de volume parcial. Os conjuntos de dados com CVF são adquiridos com baixa resolução espacial devido aos longos tempos de aquisição relacionados a sua alta dimensionalidade, (x,y,v,t), mas proporcionam a distribuição de velocidade associada a cada um dos seus grandes voxels. Entretanto, não um mapa de velocidades do escoamento. Finalmente, a DFC proporciona alta resolução temporal e espacial, além de tempos de aquisição reduzidos. Porém a sua acurácia é dependente das hipóteses do modelo. O objetivo deste trabalho é apresentar um método capaz de integrar medidas diretas de RM (CF ou CVF) a um algoritmo para determinar solução das equações da dinâmica dos fluidos, tendo como principal interesse a redução no tempo de exame num ambiente clínico. Desta forma, este trabalho apresenta duas contribuições originais: (i) um método para extrair mapas de velocidade de alta resolução espacial a partir de dados com CVF de baixa resolução espacial; e (ii) um método para reconstruir campos de velocidade que são influenciados ao mesmo tempo pelas medidas diretas por RM (CF ou CVF) e pelos modelos físicos da dinâmica dos fluidos. _________________________________________________________________________________ ABSTRACT / Cardiovascular diseases represents today one of the leading causes of death worldwide. Thus, knowledge of the blood flow patterns in the human body is important for diagnosis and research of certain diseases. Using magnetic resonance (MRI), in vivo 3D blood flow patterns can be either measured directly using phase–contrast (PC) MRI or Fourier velocity encoding (FVE) MRI. On the other hand, blood flow patterns can be obtained from model-based computational fluid dynamics (CFD) calculations. PC–MRI suffers from long scan times, limited spatio–temporal resolution, partial–volume effects and low signal–to–noise ratio (SNR). Fourier velocity encoding (FVE) is a promising magnetic resonance imaging method for measurement of cardiovascular blood flow. FVE provides considerably higher SNR than phase contrast imaging, and is robust to partial–volume effects. FVE data is usually acquired with low spatial resolution, due to scantime restrictions associated with its higher dimensionality, (x; y; v; t), and provides the velocity distribution associated with a large voxel, but does not directly provides a velocity map. Since the acquisition time is a disadvantage of FVE, then, preferably, it should be acquired with rapid spiral trajectories. Finally, CFD provides arbitrarily high spatial and temporal resolution and reduced scan times, but its accuracy hinges on the model assumptions. The objective of this work is to present a method capable of integrate direct MRI (PC or FVE) measures inside a CFD solver on the way to reduce scan time in a clinical environment. This work presents two main original contributions: (i) a method to derive velocity maps with high spatial resolution from low spatial resolution FVE data; and (ii) a numerical framework for constructing a flow field that is influenced by both direct measurements (PC–MRI or FVE) and a fluid physics model.