Neste trabalho estuda-se analiticamente um modelo de turbulência com condições de fronteira 2¼-periódicas. Para isso, considera-se uma família de modelos de Simulação em Grandes Escalas (LES) que é baseada na deconvolução aproximada de Van Cittert, onde especifica-se o filtro diferencial e suas propriedades. Para conhecer melhor os modelos de fecho, especifica-se os requisitos para um modelo de fecho ser satisfatório. Define-se uma norma e mostra-se que esta é equivalente à norma L2, sendo usada no estudo da estabilidade (solução satisfazendo uma desigualdade de energia) e da existência de soluções fracas do modelo. Para mostrar a existência de soluções, além de usar a norma definida neste trabalho, usa-se o método de Faedo-Galerkin e a propriedade antissimétrica pertinente ao termo não linear do modelo. / In this work we study analytically a model of turbulence with 2¼- periodic boundary conditions. To that end, we consider a family of Large Eddy Simulation (LES) models based on the van Cittert approximate deconvolution procedure, and we also specify a differential filter and its properties. To learn more about the closure models, we specify the requirements of a satisfactory closure model. We define a norm which is shown to be equivalent to the L2-norm and is used to study stability (solutions satisfy an energy inequality) and existence of weak solutions of the model. To show the existence of solutions, we use the Faedo-Galerkin method and the model’s nonlinear term skew-symmetry property.

Consideramos o problema de Stokes em regiões com cilindros semi-infinitos conectados por um domínio limitado qualquer. Além das condições iniciais e de contorno, é imposto uma condição sobre o comportamento da pressão quando |x| → +∞. Apresentamos um problema proposto por D. Koch envolvendo tais situações e formulamos o problema em regiões mais gerais. Técnicas usadas na resolução do conhecido problema de Leray permitem a demonstração da existência de soluções para o problema variacional. / We consider the Stokes’ problem in regions with cylindrical outlets to infinity connected by an arbitrary bounded domain. We impose an additional condition about the pressure behavior as |x| → +∞, in addition to the usual initial and boundary conditions over the velocity field. We present a problem proposed by D. Koch concerning such regions and conditions, and we formulate it in terms of more general ones. The techniques used in the well-known Leray’s problem allow the demonstration of the existence of solutions to the variational problem.

Neste trabalho são provadas algumas estimativas de erro em espaços para as aproximações de Galerkin para a solução do sistema de equações de Navier-Stokes. Mostra-se que o erro decresce em proporção inversa aos autovalores do operador de Stokes.

Nessa dissertação apresenta-se o desenvolvimento de um modelo espectral para o decaimento da turbulência na Camada Limite Convectiva (CLC). A equação dinâmica para o espectro de energia é obtida a partir das equações de Navier-Stokes [Rinze, 1975J, na qual o termo de produção de energia por efeito mecânico foi desprezado. Durante o período de transição dia-noite foi considerado que o fluxo de calor da superfície é interrompido instântaneamente, o que permite desconsiderar o termo de produção ou perda de energia por efeito térmico. Para parametrizar o termo de transferência de energia cinética, para um fluxo turbulento e isotrópico, foi utilizado o modelo de Heisenberg, que descreve o decaimento da turbulência como um processo em que os turbilhões de um certo tamanho transferem energia por efeito de uma viscosidade cinemática turbulenta, aos turbilhões menores através de interações em cascata. Para resolver a equação diferencial não linear, que descreve o decaimento da turbulência na Camada Limite Convectiva (CLC), fixa-se o parâmetro "z" e considera-se que o espectro de energia pode ser dividido como um produto de duas funções. Neste trabalho, foi utilizado o modelo de Kristensen [1989] para obter o espectro inicial tridimensional.Foi calculada a componente vertical do coeficiente difusão, considerando o espectro vertical, que foi obtido usando-se uma função peso, a qual informa a contribuição de cada componente na formação do espectro tridimensional. Finalmente, os resultados obtidos no trabalho foram confrontados com os dados de LES - Large Eddy Simulation [Nieuwstadt e Brost, 1986] existentes na literatura. / This work presents the development of a spectral model for the decay of turbulence in the Convective Boundary Layer (CBL). The dynamic equation for the energy spectrum is obtained from the Navier-Stokes equations [Hinze, 1975J, in which the term of energy production by mechanical effect was disregarded. During the period of transition from day to night it was considered that the flow of heat from the surface is instantaneously interrupted, which aIlowed the dismissal of the term of production or loss of energy instantaneously. For the purpose of parameters for the kinetic energy term of transference in a turbulent isotropic flow, the Heisenberg model was used, which describes the decay of turbulence as a process in which certain sized eddy transfer energy to smaIler eddy by means of turbulent kinetic viscosity through cascading interactions. For the purpose of solving the non-linear differential equation describing the decay of turbulence in the Convective Boundary Layer (CBL), we established the "z" parameter and considered that the energy spectrum can be divided as a product of two functions. This paper uses the Kristensen model [1989] for obtaining the initial tridimensional spectrum. The vertical component of the diffusion coefficient was calculated considering the vertical spectrum, which in turn was obtained by means of a weight function, which informs the contribution of each component in the whole of the tridimensional spectrum.FinaIly, the results obtained were confronted with the LES - Large Eddy Simulation [Nieuwstadt e Brost, 1986] existing in literature, for validation.

Este trabalho consiste de duas partes. Na primeira, estendemos o resultado de Braz e Silva e Zingano [2], [3] sobre soluções u(•; t) ε C°([0; T[;Lp(Rn)) de equações de advecção-difusão em meios heterogêneos para classes mais gerais de equações parabólicas, aplicando os resultados nas equações de Navier-Stokes incompressíveis no plano formuladas em termos da vorticidade do escoamento. Em particular, estabelecemos estimativas mostrando o decaimento em certas normas do campo de velocidade u(•; t) em caso de escoamentos de energia infinita. Na segunda parte, consideramos as equações de Navier-Stokes em dimensão n = 2; 3 examinando soluções u(•; t) de energia finita. Inicialmente, obtemos uma nova derivação, mais simples, do resultado obtido originalmente por Kato [20] estabelecendo o decaimento assintótico (t → ∞) de ||u(•; t)||L²(Rn), para estados iniciais u0 ε H¹(Rn) (com divergente nulo) arbitrários. Na linha deste argumento obtemos uma formula»c~ao mais forte dos resultados fundamentais de Wiegner [36] relacionando u(•; t) com soluções evΔtu0 da equação do calor, adaptando o método recentemente introduzido em [22], [23] para a derivação destes resultados. O método de [22], [23] também é utilizado para estabelecermos (dimensão n=3) que, ocorrendo "blow- up"de u(•; t) em tempo finito t*, necessariamente t* < 0:159||u0||4Lp(Rn)º-5, sendo ν a viscosidade dinãmica do escoamento. / In the first part of the this work, we extend results by Braz e Silva e Zingano [2], [3] concerning Lp solutions u(•; t) ε C°([0; T[;Lp(Rn)) of advection-dicusion equations in heterogeneous media to broader classes of quasilinear parabolic equations, applying the results to incompressible Navier-Stokes flows in the plane by way of the vorticity formulation. In particular, we obtain some decay rates (as t → ∞) for certain norms of the velocity field u(•; t) in case of flow with infinity energy. In the second part, we consider the Navier-Stokes equations in dimension n = 2; 3 and examine solutions u(•; t) with finite energy. First, we give a new (and simpler) derivation of the time asymptotic result originally obtained by Kato [20] and Masuda [28] showing the decay of the L2 norm of divergence-free, finite- energy solutions. Following these footsteps, we give a stronger formulation of the fundamental results obtained by Wiegner [36] relating the velocity field u(•; t) to solutions evΔtu0 of the heat equation, adapting the approach introduced in [22], [23] for the derivation of Wiergner's results. The analysis in [22], [23] is also used to obtain an interesting bound for the blow-up time t* in 3-D flows, in case solutions cease to be smooth: one must have t* < 0:159||u0||4Lp(Rn)º-5, where v is the dynamic viscosity.

Neste trabalho, consideramos um problema de valor inicial e de fronteira, semelhante ao sistema de equações de Navier-Stokes, que descreve o movimento de um meio granular com densidade constante. Utilizando o método de Galerkin: provamos a existência de soluções fracas em espaços de Sobolev. / In this work, we consider an initial-boundary value problem similar to the Navier-Stokes system of equations describing the motion of a granulated medium with constant density. We prove, using Galerkin method, the existence of weak solutions in Sobolev spaces.

Neste trabalho são provadas algumas estimativas de erro em espaços para as aproximações de Galerkin para a solução do sistema de equações de Navier-Stokes. Mostra-se que o erro decresce em proporção inversa aos autovalores do operador de Stokes.

O estudo das equações de Navier-Stokes desperta interesse dos estudiosos da área da análise numérica, visto que a partir destas pode-se determinar os campos de velocidade e pressão de um escoamento. Com estas equações também pode-se aproximar coeficientes aerodinâmicos, fato de grande interesse nas indústria automobilística e aeronáutica. Por isso, propõe-se estudar a aproximação das equações de Navier Stokes via o método de elementos finitos, que tem se mostrado um bom método de resolução de problemas envolvendo fiuidos. Estudam-se métodos de linearização do termo convectivo. Apresentam-se três propostas de métodos de discretização temporal para as equações dadas. Através da análise de erro e taxas de convergência para problemas com solução exata conhecida, comparam-se diferentes espaços de discretização espacial e diferentes métodos de discretização temporal. Introduz-se um modelo de regularização e através do cálculo dos coeficientes de arrasto e sustentação comprova-se a sua efetividade (no sentido de que ele permite trabalhar com malhas mais grossas, mas ainda obter soluções comparáveis àquelas obtidas por DNS). / The study of the Navier-Stokes equations arouses interest of researchers in the area of numerical analysis, since from these one can determine the velocity and pressure fields of a fiow. These equations approach aerodynamic coefficients, a fact of great interest in the aeronautical and automotive industries. Therefore, we propose to study an approximation of the Navier Stokes equations through the finite element method, which has shown to be a good method for solving problems involving fiuids. Ve study methods of linearizing the convective term, and present three methods for time discretization. Through error analysis and convergence rates for problems with known exact solution, we compare different spatial and time discretization methods. By calculating the drag and lift coefficients around a cilinder we confirm the effectiveness of the Leray-deconvolution model in comparison to direct simulation (DNS) for solving problems with coarser meshes.

Made available in DSpace on 2016-08-29T15:33:15Z (GMT). No. of bitstreams: 1 tese_4176_.pdf: 6013163 bytes, checksum: 2e413f36a79ecc7f2ea7e35d0744e463 (MD5) Previous issue date: 2011-08-25 / As unidades de processamento gráfico (Graphics Processing Unit -- GPU) surgiram como um poderoso dispositivo computacional e a plataforma Compute Unified Device Architecture (CUDA) como um ambiente adequado para a implementação de um código na GPU. Especializada inicialmente em processamento gráfico, a GPU vem sendo designada à otimização de cálculos lógicos e aritméticos beneficiando diversas áreas de pesquisa com a redução do tempo de computação. O objetivo deste trabalho é mostrar como aplicações em mecânica dos fluidos, discretizadas pelo método das diferenças finitas, podem lucrar bastante com esta tecnologia. Implementações paralelas na GPU em C+CUDA das equações de Navier-Stokes e de transporte são comparadas com uma versão sequencial implementada na CPU em C. É utilizada uma formulação em diferenças finitas implícita-explícita, sendo o algoritmo caracterizado como sendo explícito nas velocidades e temperatura e implícito na pressão. A resolução dos sistemas lineares resultantes é feita utilizando um esquema de coloração Red-Black para as células internas da malha e o método iterativo successive-over-relaxation (SOR), denominado Red-Black-SOR. É discutido neste trabalho os impactos do uso de tipos de dados double e float e também a utilização de memórias shared e global existentes na GPU. O algoritmo C+CUDA é verificado para o seguinte conjunto de problemas conhecidos da literatura: cavidade com cobertura deslizante, escoamento sobre um degrau, escoamento laminar com um obstáculo cilíndrico, convecção natural e convecção de Rayleigh-Bénard, considerando casos bidimensionais e tridimensionais. O tempo de processamento é comparado com o mesmo algoritmo implementado em C. Os resultados numéricos mostraram que é possível alcançar speedups da ordem de 85 vezes para dados float e 61 vezes para dados double utlizando C+CUDA.

Dissertação (mestrado) - Universidade Federal de Santa Catarina, Centro de Ciências Físicas e Matemáticas. Programa de Pós-Graduação em Matemática e Computação Científica / Made available in DSpace on 2012-10-21T15:12:44Z (GMT). No. of bitstreams: 1 221522.pdf: 893831 bytes, checksum: e0f7a25edf5c463738b416dc7b1239f0 (MD5) / O trabalho aqui presente destina-se a fazer uma análise comparativa, no contexo do método de Newton inexato, os desempenhos das metodologias iterativas baseadas em subespaços de Krylov: GMRES (Generalized Minimum Residual Method) e Bi-CGStab (Biconjugate Gradient Stabilized) e um método direto (LU esparso). As características das desempenhos (número de iterações e tempo computacional) das metodologias investigadas são acessadas com o uso de alguns testes padrão largamente utilizados como "benchmark" em mecânica dos fluidos computacional. O método de Newton inexato baseado em GMRES e Bi-CGStab é aplicado no sistema não linear gerado pelo método de elementos finitos (MEF) sobre o problema de valor de contorno composto pelas equações de Navier-Stokes. Uma importante observação diz respeito a condição necessária e suficiente de Brezzi-Babuka (ou condições inf-sup), a qual é satisfeita com o uso de parâmetros de estabilização.