[pt] Neste trabalho três métodos foram utilizados para estudar o escoamento turbulento em regiões de florestas. No primeiro método, a influência da vegetação no escoamento foi modelada através da adição de termos fontes nas equações de quantidade de movimento, energia cinética turbulenta e sua taxa de dissipação. No segundo, a vegetação foi considerada um meio poroso homogêneo. Finalmente, a camada do dossel foi representada por modelos 3-D de árvores, consideradas como obstáculos individuais. As equações foram resolvidas através do modelo de turbulência k −E padrão com o programa comercial FLUENT 6.2.16. As previsões dos perfis verticais da velocidade do vento médio, da intensidade da turbulência e dos tensores de Reynolds, foram comparadas com dados de experimentos de túnel de vento. Os resultados preditos dos perfis verticais da velocidade média e intensidade da turbulência, na primeira e na segunda metodologias, apresentaram boa concordância com os valores experimentais, porém, foram observadas discrepâncias nos perfis modelados do tensor de Reynolds. Entretanto, qualitativamente, a modelagem consegue capturar o comportamento físico do tensor de Reynolds no interior de florestas. Uma possível explicação para este fato, é que o modelo considera a isotropia para a viscosidade turbulenta, implicando na incapacidade de prever qualquer forte anisotropia do campo turbulento. Na terceira metodologia, as previsões dos perfis verticais de velocidade média e intensidade da turbulência apresentaram discrepâncias em relação às medições. Porém, os perfis verticais do tensor de Reynolds apresentaram boa concordância. Todos os perfis verticais da velocidade média apresentaram um ponto de inflexão na interface vegetação-atmosfera, característico de uma camada de mistura. Nas duas primeiras metodologias, este padrão foi confirmado nos perfis de tangente hiperbólica de uma camada de mistura. / [en] This work investigates different procedures in order to study the turbulent flow over the scale model of a forest region. Initially, the canopy flow was modeled by using source terms in the momentum, turbulent kinetic energy and its dissipation rate equations. After that, the forest canopy was considered a homogeneous porous medium. In the last step, the canopy boundary layer was modeled by artificial 3-D tree models. This was done by using the standard k−E turbulence model with the FLUENT commercial program. The modeled profiles of mean velocity, turbulence intensity and Reynolds stress were compared against data from wind tunnel experiments. In the two first methodologies, the model predictions of the vertical profiles of the wind speed and turbulence intensity showed good agreement with the experimental data. It was found that predictions of the Reynolds tensor were sensitive to the parameterization scheme of the standard k −E model. However, qualitatively, the model was capable of predicting the physical behavior of the Reynolds stress tensor in the canopy flow. A possible explanation for this behavior is the omission of any anisotropic eddy-viscosity effects within the k - E modelling approach. When it was considered the tree array, the model predictions for the wind speed and turbulence intensity were less satisfactory. However, it was found that the predicted results of the Reynolds stress tensor agreed well with the measured data. All the vertical profiles of the mean velocity contained an inflection point, something which is a necessary criterion for the mixing layer flow. In the tree array, the modeled results failed to the capture this behavior of the canopy flow. In the 2-D numerical simulations, it was found the characteristic hyperbolic tangent profile of a mixing layer.

Three-dimensional design and three-dimensional CFD analysis of a compressor rotor stage are performed. The design methodology followed is based on a mean line analysis and a radial equilibrium phase. The radial equilibrium is established at a selected number of radii. NACA 65 series airfoils are selected and stacked according to the experimental data available. The CFD methodology applied is based on a three-dimensional, finite difference, compressible flow Euler solver that includes the source terms belonging to rotational motion. The accuracy of the solver is shown by making use of two different test cases. The CFD solution of the designed geometry predicts the static pressure rises and flow turning angles to a good degree of accuracy.

Les machines électriques sont modélisées pour prédire leurs performances et optimiser leur rendement. Cette modélisation peut être faite par des simulations avec la méthode des éléments finis. En particulier, les machines asynchrones nécessitent des simulations 3D pour prendre en compte les courants de Foucault et les têtes de bobines. Dans le logiciel Flux®, des formulations 3D basées sur le potentiel scalaire magnétique sont utilisées avec succès depuis de nombreuses années. Néanmoins, des coupures mathématiques artificielles sont nécessaires, lorsque le domaine n'est pas simplement connexe.Afin de se libérer de ces contraintes de connexité, des formulations en potentiel vecteur magnétique ont été étudiées et développées. En 3D, leur mise en œuvre nécessite l'utilisation d'éléments finis d'arêtes afin de respecter la nature des champs. Avec les éléments d'arêtes, les formulations sont généralement résolues avec une condition de jauge pour les solveurs directs comme pour les solveurs itératifs. De nouvelles formulations en potentiel vecteur magnétique auto-jaugées ont été développées permettant la prise en compte des bobines maillées et des bobines non maillées. La prise en compte du mouvement est relativement simple à mettre en œuvre pour les formulations en potentiel scalaire magnétique avec l'interpolation nodale.Avec les éléments d'arête, l'interpolation est plus délicate. C'est pourquoi la méthode des éléments avec joints a été développée pour prendre en compte le mouvement dans un cas général. / Electric machines are modeled in order to predict their performance and to optimize their output. This modeling can be done by simulation with the finite element method. In particular, induction machines require 3D simulation to take into account eddy currents and coils overhangs. In the Flux® software, 3D formulations based on magnetic scalar potential has been used with success for many years. Nevertheless, artificial mathematical cuts are necessary, when the domain is not simply connected.In order to avoid connection constraints, magnetic vector potential formulations have been studied and developed. In 3D, their implementation requires the use of edge elements to respect the nature of fields. With edge elements, formulations are generally solved with a gauge condition for direct solvers as well as for iterative solvers. New auto-gauged magnetic vector potential formulations have been developed to take into account meshed coils and non-meshed coils. Consideration of movement is relatively simple to implement for magnetic scalar potential formulations with nodal interpolation. With edge elements, the interpolation is more delicate. For this reason, the mortar method has been developed to take into account movement in a general case.

Cette étude porte sur le calcul numérique du champ acoustique rayonné par des écoulements subsoniques turbulents présentant des inhomogénéités de température. Des méthodes hybrides sont développées grâce à un développement de Janzen-Rayleigh des équations de Navier-Stokes. L'écoulement est résolu par un calcul quasi incompressible puis les perturbations acoustiques sont propagées selon deux méthodes : les équations d'Euler linéarisées (EEL) et l'approximation à faible nombre de Mach perturbée(PLMNA). Les méthodes sont validées sur des cas simples puis appliquées à une couche de mélange isotherme et anisotherme en développement spatial. / This study focuses on the numerical calculation of the acoustic field radiated by subsonic turbulent flows with temperature inhomogeneities. Hybrid methods are developed through a Rayleigh-Janzen expansion of the Navier-Stokes equations. The flow is solved in a quasi-incompressible way then the acoustic disturbances are propagated by two methods : the linearized Euler's equations (EEL) and the perturbed low Mach number approximation (PLMNA). The methods are validated on simple cases and then applied to an isothermal and non isothermal spatially evolving mixing layer.

Submitted by Pedro Silva Filho (pfsilva@ipen.br) on 2017-11-21T11:45:58Z No. of bitstreams: 0 / Made available in DSpace on 2017-11-21T11:45:58Z (GMT). No. of bitstreams: 0 / O gerenciamento de riscos em uma instalação nuclear é necessário para a segurança de trabalhadores e de populações vizinhas. Parte desse processo é a comunicação dos riscos que propicia o diálogo entre gestores da empresa e moradores das áreas de risco. A população que conhece os riscos a que está exposta, como esses riscos são gerenciados e o que deve ser feito em uma situação de emergência tende a se sentir mais segura e a confiar nas instituições responsáveis pelo plano de emergência. Sem diálogo entre empresa e público, o conhecimento dos procedimentos a serem seguidos em caso de acidente não chega à população, ou quando chega, não há confiança dessas pessoas na sua eficácia. Em Angra dos Reis, no litoral sul do Estado do Rio de Janeiro, está a Central Nuclear Almirante Álvaro Alberto. No entorno dessa Central Nuclear existe uma população que, de acordo com o Plano de Emergência Externo (PEE/RJ), deverá ser evacuada ou ficar abrigada, caso ocorra um acidente na instalação. Um trabalho de comunicação de riscos entre esses moradores é necessário para que eles conheçam o plano de emergência e os procedimentos corretos para uma situação de emergência, além de buscar esclarecer dúvidas e mitos. Esse trabalho apresenta uma análise da comunicação dos riscos feita para a população local, a percepção que ela tem dos riscos e o grau de conhecimento do plano de emergência externo por parte dessas pessoas. / Tese (Doutorado em Tecnologia Nuclear) / IPEN/T / Instituto de Pesquisas Energéticas e Nucleares - IPEN-CNEN/SP

In this work we review 12 years of developments in the field of residual based discretizations for hyperbolic problems and their application to the solution of the shallow water equations. Fundamental concepts related to the topic are recalled and he construction of second and higher order schemes for steady problems is presented. The generalization to time dependent problems by means of multi-step implicit time integration, space-time, and genuinely explicit techniques is thoroughly discussed. Finally, the issues of C-property, super consistency, and wetting/drying are analyzed in this framework showing the power of the residual based approach.

Nonlinear four-wave interactions (Snl) are critical for acquiring realistic spectra needed by operational wave models. High computational demand to calculate these interactions led to an approximation method named the Discrete Interaction Approximation (DIA) to be used broadly in the major operational wave models for a long time. However, the accuracy of the DIA has been controversial since it was first introduced and more precise approximations such as the Two Scale Approximation (TSA in short) are now available. The only issue with the initial TSA`s efficiency is performing an order of a magnitude slower than the DIA in speed. This study questions the exactness of the DIA while trying to increase the competence of the TSA by making improvements on its execution time. Particularly, in this thesis, the main effort is on the local scale term of the TSA since it is the part that consumes the most time while running the code. The findings of this work imply that the TSA can improve its operation speed significantly while maintaining its accuracy with making alterations in the code. By decreasing the number of bands in the local scale it is possible to run the TSA up to 7.5 faster than its initial version.