O presente trabalho versa sobre a modelagem de uma chama turbulenta difusiva usando a técnica de simulação de grandes estruturas turbulentas (LES), juntamente com o modelo termo-químico de folha de chama (flame sheet model) e o conceito de fração de mistura como escalar conservativo. Este trabalho também é pioneiro de utilização de LES com reação química no Brasil, podendo colaborar para o desenvolvimento desta técnica na área de combustão. O trabalho consiste na construção e validação das rotinas computacionais de um código CFD, baseado em LES e com flexibilidade para uma futura utilização de cinética química detalhada de combustão (EDC/ISAT), para casos complexos onde modelos mais simples, como a fração de mistura, são falhos. O programa será validado em uma chama de difusão turbulenta não-confinada de metano (CH4), para a qual existem dados experimentais na literatura [61,62] e utilizados pela comunidade acadêmica em excelência (Stanford, TU-Darmstadt, Imperial College, Cornell University etc). As características da implementação numérica do código permitirão sua expansão futura para outras aplicações em: queima de combustíveis líquidos, combustão em câmaras fechadas e fornalhas com a inclusão de modelo de radiação. / The present work is about modeling a diffusive turbulent flame using the Large-Eddy Simulation approach (LES) and the Flame Sheet model as the chemical model with the mixture fraction concept as the conservative scalar. This work is pioneer in the sense of using LES and reactive flow in Brazil, making possible the development for LES techniques in the combustion area. The work is intended to construct and validate a CFD code based on LES and with future flexibility for a more detailed combustion chemical model (EDC/ISAT) for complex flows, where simple models are failed, like the mixture fraction. The program will be validated for a turbulent diffusion methane (CH4) flame not confined, which there are some experimental data on the specialized literature [61,62], and commonly used by the academic community (Stanford, TU-Darmstadt, Imperial College, Cornell University etc). The features of the numerical code implementation will make possible future expansion of its use in other applications: liquid fuel burning, combustion chambers and ovens with the radiation model inclusion.

Em mecânica dos fluidos, especificamente em processos turbulentos, o problema de fechamento representa um dos maiores desafios para qualquer pessoa interessada nesta área. Durante décadas, cientistas vêm usando abordagens estatísticas com o objetivo de \"fechar\" o problema ou, pelo menos, diminuir as dificuldades inerentes. Assim, o presente trabalho apresenta uma criteriosa análise com base em ferramentas estatísticas em que ondas quadradas aleatórias, aliadas a um número fixo de parâmetros, foram utilizadas para criar equações paramétricas para representar um fluxo turbulento unidimensional com uma abordagem a priori, diferenciando de outras abordagens aplicadas amplamente na área, que utilizam uma abordagem a posteriori. Em seguida, simulações foram realizadas, a fim de avaliar o comportamento do modelo. Nas simulações pôde-se reproduzir o comportamento observado na literatura e estipular a abrangência do método. Além disso, uma importante discussão acerca das condições de contorno foi desenvolvida. / In fluid mechanics, specifically in turbulent processes, the closure problem represents one of the biggest challenges for anyone interested in this area. For decades, scientists have been using statistical approaches aiming to close the problem or, at least, decrease the inherent difficulties. So, the present project presents a judicious analyze based on statistical tools in which random square waves, allied with a fixed numbers of parameters, were used to create parametric equations to represent a turbulent flow with an a priori approach, differentiating from other approaches broadly applied in the area, which use an a posteriori approach. Then simulations were done, in order to evaluate the behavior of the model. In the simulations, the behavior of some data from the literature could be followed and the scope of the method was stipulated. Besides this, an important discussion about boundary conditions was developed.

The present study utilizes the ultrasonic travel time technique to diagnose grid generated turbulence. Ultrasonic flow metering technology relies on the measurement and computation of small perturbations in the travel time of acoustic ultrasonic waves through the dynamic medium. The statistics of the travel time variations of ultrasonic waves that are caused by turbulence probably affect the performance of ultrasonic flow meters. Motivation for the study stems from the large travel time variations observed in typical ultrasonic flow and circulation meters. Turbulent flow data was collected downstream of a grid introduced in a uniform flow in the wind tunnel using ultrasonic techniques. Grid turbulence is well defined in literature, and is nearly homogeneous and isotropic. The experimental investigation was performed under well-controlled laboratory conditions. The grid mesh sizes varied from 0.25-0.5in, and flow velocities from 0-20m/s. The ultrasonic transducers were of 100 kHz working frequency; and all of the data was collected with them oriented perpendicular to the mean flow. Path lengths were increased from 2-10in; and the data acquisition and control system featured a very high speed data acquisition card with an analog to digital converter that enabled excellent resolution of ultrasonic signals. Experimental data was validated by comparison to other studies. The work aims to investigate the influence of the grid-generated turbulent flow on acoustic wave propagation, in terms of the variance of the travel time. The effect of turbulence on acoustic wave propagation was observed. The experimental data was used to compute average travel times, acoustic travel time variances, and standard deviation amplitude fluctuations. The data was collected in the region estimated to be homogeneous and isotropic. Average travel time data support the assumption that only the large (as compared to the wavelength ) turbulent inhomogeneities influence acoustic wave propagation. Variance data confirm the presence of a non-linear trend in the acoustic travel times with increasing path length. Amplitude fluctuations data confirm a correlation between areas of caustic formation and large amplitude fluctuations.

Questões de aeração forçada ou natural estão intimamente ligadas à capacidade de autodepuração dos corpos hídricos, ou seja, oxidar substâncias agressivas para resultar em baixo teor de toxicidade. Uma das etapas do processamento de efluentes consiste na dissolução de oxigênio em água e, para realizar essa dissolução, utilizam-se aeradores que são unidades (tanques) onde o ar é borbulhado no meio líquido, o qual se desloca em um regime contínuo de escoamento. Esta pesquisa refere-se à obtenção experimental de valores de viscosidade turbulenta para inserção em modelagem fenomenológica da transferência de oxigênio das bolhas de ar para o meio líquido. Tais modelos, se bem realísticos, podem contribuir aos estudos de gestão de recursos hídricos ou em operações nos tratamentos de efluentes líquidos. O método experimental empregado foi a velocimetria por imagem de partículas, no qual foi possível obter velocidades instantâneas do fluido (água). Estas consideram o movimento turbulento, que é o principal responsável pelo transporte de oxigênio da superfície para o seio do corpo hídrico, sendo que essa superfície pode ser livre para o ambiente, ou a superfície de uma bolha. Praticamente, o método consiste em correlacionar posições de partículas traçadoras em suspensão no fluido, as quais são assumidas ter a mesma velocidade do fluido. As posições consecutivas para fornecer a trajetória e a velocidade foram obtidas por imagens capturadas em uma freqüência definida através de uma câmera digital, onde a luz do laser contrastou as partículas em uma área desejada com uma precisão elevada. Então, nessa área (um plano), foi possível correlacionar um perfil de velocidades. Assim, os valores de viscosidade turbulenta foram obtidos para serem usados em modelagem da transferência de oxigênio, os quais poderão contribuir nos estudos de aeração em corpos hídricos. / Problems involving natural or forced aeration are intimately bind to the reaeration of water bodies. Pollutants are oxidized to yield low toxicity conditions. One of the steps of wastewater treatment consists in dissolving oxygen in water. To perform this, aeration tanks are used where bubbling air crosses the continuous liquid flow. This research focuses on the experimental determination of turbulent viscosity values to be used in modeling of oxygen transfer from air bubbles to the bulk liquid. Such models, if realistic enough, may contribute to water resources management studies or in wastewater treatment operations. Particle image velocimetry method was used, by means of which it became possible to obtain instantaneous velocities of the fluid (water). These velocities embody the turbulent flow, which is the main responsible for oxygen transport from the surface to the bulk liquid. This surface may be either facing the atmosphere or the interior of a bubble. In practice, the method consists in correlating tracking particles suspended in the liquid, which are supposed to have the same velocity of the fluid. The successive positions that give path and speed were obtained by images took in predefined intervals by a digital camera. The laser light illuminated the particles in a predefined area with high precision, making possible to determine velocity profiles. Turbulent viscosity values were so determined and may be used in the modeling of oxygen transfer, which may contribute to water body aeration studies.

[Formulae and special characters can only be approximated here. Please see the pdf version of the Abstract for an accurate reproduction.] Since the damping functions employed by most of the low-Reynolds-number models are related to the non-dimensional distance y+[ special character ??? near-wall non-dimensional distance in y direction], which is based on local wall shear stress, these models become invalid for separated flows, because the wall shear stress is zero at the reattachment point. In addition, the pressure-velocity correlation term is neglected in most of these models, although this term is shown in this thesis to be important in the near-wall region for simple flows and large pressure gradient flows. In this thesis, two main efforts are made to improve the k ??? [special character - dissipation rate] model. First, based on Myong and Kasagi???s (1990) low-Reynolds-number model (hereafter referred to as MK model), a more general damping function [special character - turbulent viscosity damping function in LRN turbulent model] is postulated which only depends on the Reynolds numbers [formula ??? near-wall turbulence Reynolds number]. Second, a form for the pressure-velocity correlation term is postulated based on the Poisson equation for pressure fluctuations. This modified model predicts the turbulent flow over a flat plate very well. It is found that the inclusion of the pressure-velocity correlation term leads to significant improvement of the prediction of near-wall turbulence kinetic energy. When the model is applied to turbulent flow over a backward-facing step, it produces better predictions than the traditional k ??? [special character - dissipation rate] model, FLUENT???s two-layer model and the MK model. Again, the pressure-velocity correlation term improves the turbulence kinetic energy prediction in the separated region over that of other models investigated here. The studies of numerical methods concerning computational domain size and grid spacing reveal that a very large domain size is required for accurate flat plate flow computation. They also show that a fine grid distribution in the near-wall region upstream of the step is necessary for acceptable flow prediction accuracy in the downstream separated region.

Improved basic understanding, predictability, and controllability of vortex-dominated and unsteady aerodynamic flows are important in enhancement of the performance of next generation helicopters. The primary objective of this research project was improved understanding of the fundamental vorticity and turbulent flow physics for a dynamically stalling airfoil at realistic helicopter flight conditions. An experimental program was performed on a large-scale (C = 0.45 m) dynamically pitching NACA 0012 wing operating in the Texas A&M University large-scale wind tunnel. High-resolution particle image velocimetry data were acquired on the first 10-15% of the wing. Six test cases were examined including the unsteady (k>0) and steady (k=0) conditions. The relevant mechanical, shear and turbulent time-scales were all of comparable magnitude, which indicated that the flow was in a state of mechanical non-equilibrium, and the expected flow separation and reattachment hystersis was observed. Analyses of the databases provided new insights into the leading-edge Reynolds stress structure and the turbulent transport processes. Both of which were previously uncharacterized. During the upstroke motion of the wing, a bubble structure formed in the leading-edge Reynolds shear stress. The size of the bubble increased with increasing angle-of-attack before being diffused into a shear layer at full separation. The turbulent transport analyses indicated that the axial stress production was positive, where the transverse production was negative. This implied that axial turbulent stresses were being produced from the axial component of the mean flow. A significant portion of the energy was transferred to the transverse stress through the pressure-strain redistribution, and then back to the transverse mean flow through the negative transverse production. An opposite trend was observed further downstream of this region.

Porous media flows are encountered in many natural and man-made systems such as gas adsorption, filtration, heat exchangers, combustion, catalytic reactors and groundwater hydrology. This study experimentally investigates these flows as function of pore Reynolds number, Re[subscript pore]. The pore Reynolds number is based on the porous bed hydraulic diameter, D[subscript H] =φD[subscript Β]/(1−φ) where φ is bed porosity and D[subscript B] is solid phase bead diameter and average bed interstitial velocity, V[subscript int]= V[subscript Darcy]/φ, where VDarcy= Q/A[subscript bed], with Q being the volumetric flow rate and A[subscript bed] the bed cross section normal to the flow. The flow characteristics are studied through application of a particle displacement technique called particle image velocimetry, PIV. In the case of PIV, flow fields are estimated by seeding the flow with tracer particles and then evaluating their displacements. Application of quantitative imaging technique such as PIV to a complex flow domain like porous bed requires matching refractive index of liquid phase to that of the solid phase. Firstly, the effect of slight index mismatch, due to experimental uncertainties, on obtaining highly accurate PIV measurements as expressed as an experimental uncertainty was explored. Mismatch of refractive indices leads to error in estimation of particle positions and their displacements due to refraction at solid-liquid interfaces. Slight mismatch, in order of 10⁻³, in refractive indices also leads to reduction in particle density, particle signal peak intensity and degrade the particle image. These effects on velocity field estimation using PIV is studied experimentally and numerically. The numerical model, after validating against experimental results, is used to generate an expression for the error in PIV measurements as a function of refractive index mismatch for a range of bead diameters, bed widths, bed porosity, and optical magnification. After refractive index matching, planar PIV measurements were taken at discrete locations throughout a randomly packed bed with aspect ratio (bed width to bead diameter) of 4.67 for steady, low pore Reynolds number flows, Re[subscript pore] ~ 6, intermediate Re[subscript pore] of 54 and unsteady flow with high Re[subscript pore] ranging from 400-4000. Details of the measurement uncertainties as well as methods to determine local magnification and determination of the dynamic velocity range are presented. The data are analyzed using the PIV correlation averaging method for steady flows and multigrid and multipass correlation methods for unsteady turbulent flows with the largest velocity uncertainties arising from in plane image loss and out of plane motion. Results for low Re[subscript pore] flows show the correspondence of the geometric and velocity correlation functions across the bed, and that the centerline of the bed shows a random-like distribution of velocity with an integral length scale on the order of one hydraulic diameter (or 0.38 bead diameters based on the porosity for this bed). The velocity variance is shown to increase by a factor of 1.8 when comparing the center plane data versus using data across the entire bed. It is shown that the large velocity variance contributes strongly to increased dispersion estimates, and that based on the center plane data of the variance and integral length scales, the dispersion coefficient matches well with that measured in high aspect ratio beds using global data. For unsteady and turbulent flow, velocity data were used to determine the following turbulence measures: (i) turbulent kinetic energy components, (ii) turbulent shear production rate, (iii) integral Eulerian length and time scales, and (iv) energy spectra all for a range of pore Reynolds numbers, Re[subscript pore], from 418 to 3964. These measures, when scaled with the bed hydraulic diameter, DH, and average interstitial velocity, V[subscript int], all collapse for Re[subscript pore], beyond approximately 2800, except that the integral scales collapse at a lower value near 1300-1800. The results show that the pore turbulence characteristics are remarkably similar from pore to pore and that scaling based on bed averaged variables like D[subscript H] and V[subscript int] characterizes their magnitudes despite very different local mean flow conditions. In the case of high Re[subscript pore] flows, large scale structures such as stationary and convected vortices and structures resembling jets were also identified. These structures were analyzed in detail using decomposition techniques like Large Eddy Scale decomposition and critical point analysis like swirl strength analysis. Direct velocity measurements were used to estimate Lagrangian statistics through Eulerian measures and then estimate contribution of flow structures to turbulent mechanical dispersion. Results agree well with those in the literature obtained using global measurements in very high aspect ratio, long test beds. Stationary vortical or recirculation regions were seen to play a dominant role in contributing to overall dispersion in porous beds. / Graduation date: 2013

Excess global elemental sulphur (So) production has resulted in a decrease in its price. As a result, many companies, such as Syncrude Canada Ltd., have resorted to above ground storage alternatives. Geochemical reactions in these above ground blocks produce elevated concentrations of H2SO4 (acid). This acid can have potentially deleterious effects on the environment. As such, these blocks will require long-term (500 years) monitoring and maintenance. Presently the So is removed from the product stream, piped in a molten state, and poured over a low permeability liner in thin lifts. As the So cools and undergoes crystal structure change it fractures, creating preferential flow passages which are potentially highly conductive. An understanding of the liquid conductivity (Kl) of the block and knowledge regarding the spatial and temporal distribution of acid (H2SO4) within these blocks is required. In this thesis, gas pumping tests were conducted on an above ground block to determine the gas flow rates within the block and to indirectly determine the Kl of the block. Measurements of the relative humidity (RH) in the block were used to observe changes in stored acid concentrations with time and location. The results of the gas conductivity (Kg) testing showed that the block is anisotropic and is highly conductive in both the horizontal and vertical directions. Cross hole tests appeared to produce the most representative estimates of Kg due to the negation of turbulence that arises in the vicinity of the borehole. The choice of gas used in the analysis had negligible effect on the resulting Kl in contrast to choice of liquid, which resulted in larger variations in Kl. The Kl was a maximum when the liquid was pure water and decreased with increasing acid strength. The geometric mean of the resulting cross hole Kl values was 2 x 10-3 m s-1 (pure water). RH measurements were observed to fluctuate with depth and increased following precipitation. The resulting minimum pH observed within the block occurred at depths of 3 and 7 meters below the surface of the So block and increased with depth. The arithmetic mean pH value based on the daily averaged RH measurements was -1.7.

In modern aircraft, one of the most difficult issues has been how to provide avionics with adequate cooling. Future versions of the fighter aircraft JAS 39 Gripen is equipped with new applications that have increased heat loads. In previous versions of the JAS 39 Gripen avionics was cooled by zero degree air and fuel, but in the next version a liquid loop will be installed to cool the new radar.The fluid in the liquid loop is cooled to proper temperature by pressurized bleed air from the engine which is cooled by ram air. The air to cool the avionics is produced the same way and this is a very expensive process for the airplane which lowers its performance. It is important to minimize the production of cooling air and therefore three new adjustable valves that provide various components of cooling air are installed in the next version of the JAS 39 Gripen. The cooled and pressure controlled air from the engine is distributed between different avionic shelves, each containing a set of components. Depending on the type of tasks performed and current flight mode of the aircraft the requirement of functions which should be active varies and therefore also the cooling demand to avionics. The first part of this thesis studies the overall priority of how the engine bleed shall be used. This part of the thesis results in a decision basis for the distribution of cooling air to be regulated in the absence of full cooling capacity. The amount of cooling which must be distributed to the radar is proportional to its developed power which varies widely depending on the radar’s operational mode. Since the pump which determines the liquid flow velocity operates at a constant speed is the regulation of cooling to the radar is controlled by varying the bleed air flow into the heat exchanger which cools the fluid and thus the temperature of the fluid has when it reaches the radar. This part of the thesis creates a control algorithm for controlling the airflow into the heat exchanger. The regulation keeps the fluid inlet temperature to the radar within the range of +25 ± 5 ˚ C and the gradient of the temperature less than 0.5° C per second. The PI-controller with the feed-forward filter succeeded in controlling the temperature of the liquid as it reached the radar within +25 ± 1° C, the temperature gradient requirement, 0.5° C per second, was also passed in all flight cases which were used to evaluate the controller. The PI-controller with feed-forward has a low convergence time and no static error. It also performs well when the measurement signals contain a lot of noise because of the controllers integrated low pass filter.   The three new adjustable valves saves 12 to 97 g/s of cooling air for the different valve positions studied in this thesis, this corresponds to 9 - 70% of the total amount of controllable air to the avionics. Since the production of cooling air is a costly process for the aircraft, the use of all 3 valves is recommended.

Improved basic understanding, predictability, and controllability of vortex-dominated and unsteady aerodynamic flows are important in enhancement of the performance of next generation helicopters. The primary objective of this research project was improved understanding of the fundamental vorticity and turbulent flow physics for a dynamically stalling airfoil at realistic helicopter flight conditions. An experimental program was performed on a large-scale (C = 0.45 m) dynamically pitching NACA 0012 wing operating in the Texas A&M University large-scale wind tunnel. High-resolution particle image velocimetry data were acquired on the first 10-15% of the wing. Six test cases were examined including the unsteady (k>0) and steady (k=0) conditions. The relevant mechanical, shear and turbulent time-scales were all of comparable magnitude, which indicated that the flow was in a state of mechanical non-equilibrium, and the expected flow separation and reattachment hystersis was observed. Analyses of the databases provided new insights into the leading-edge Reynolds stress structure and the turbulent transport processes. Both of which were previously uncharacterized. During the upstroke motion of the wing, a bubble structure formed in the leading-edge Reynolds shear stress. The size of the bubble increased with increasing angle-of-attack before being diffused into a shear layer at full separation. The turbulent transport analyses indicated that the axial stress production was positive, where the transverse production was negative. This implied that axial turbulent stresses were being produced from the axial component of the mean flow. A significant portion of the energy was transferred to the transverse stress through the pressure-strain redistribution, and then back to the transverse mean flow through the negative transverse production. An opposite trend was observed further downstream of this region.