O presente trabalho apresenta um estudo numérico dos efeitos da aplicação de geometrias geradoras de vórtices, com intuito de controlar passivamente a camada limite, em um perfil aerodinâmico que integra a asa traseira de multi elementos de um veículo de Fórmula SAE. As equações de Navier-Stokes com médias de Reynolds foram resolvidas utilizando o modelo k-ω SST (Shear Stress Transport) para o problema de fechamento da turbulência. Uma metodologia numérica padrão foi definida e utilizada nos diferentes casos analisados. Domínio de cálculo, malha, condições de contorno e critério de convergência foram escolhidos com base em norma SAE para análise numérica de escoamento externo em veículos terrestres. As camadas de volumes prismáticos próximos as superfícies com não-deslizamento foram dimensionadas de forma a resultar em um tratamento de parede adequado ao modelo de turbulência aplicado. O método GCI (Grid Convergence Index) foi utilizado para avaliar a qualidade da malha. Com o intuito de reduzir o custo computacional nos testes com diferentes configurações de geradores de vórtices, apenas parte de interesse do domínio de cálculo foi resolvido, impondo perfis de velocidade, energia cinética da turbulência e dissipação específica em sua entrada. Estas condições foram importadas da simulação com domínio completo resolvida Para verificar a correta captação dos principais efeitos físicos envolvidos, comparações com resultados experimentais foram feitas para 2 casos com escoamentos representativos: o corpo de Ahmed e um perfil aerodinâmico com geradores de vórtices. Além disso, as diferenças entre resolver o domínio completo ou parcial foram estudadas em outro comparativo com resultados experimentais. Concluiu-se que a metodologia numérica foi capaz de obter os coeficientes aerodinâmicos, e suas tendências frente a mudanças de geometria, nos casos estudados. Resolver parcialmente o domínio, impondo perfis em sua entrada, acarretou em diferença nos coeficientes obtidos na ordem de 2% para o coeficiente de sustentação e 7% para o coeficiente de arrasto. O controle passivo via geradores de vórtices foi eficaz em atrasar a separação da camada limite no flap do veículo de Fórmula SAE, as melhoras nos coeficientes de arrasto e sustentação foram da ordem de 7% e 0,3%, respectivamente. / The present work is a numerical study of the effects of the application of vortex generating geometries, in order to passively control the boundary layer, in an aerodynamic profile that integrates a multi-element rear wing of a Formula SAE vehicle. The Reynolds Averaged Navier-Stokes equations were solved using the k-ω Shear Stress Transport model for the turbulence closure problem. A standard numerical methodology was defined and used in the different cases analyzed. Computational domain, mesh, boundary conditions and convergence criteria were chosen based on SAE standard for numerical analysis of external flow in land vehicles. The layers of prismatic volumes near the non-slip surfaces were dimensioned to result in a wall treatment suitable to the applied turbulence model. The Grid Convergence Index (GCI) method was applied to evaluate the mesh quality. In order to reduce the computational cost in tests with different vortex generators configurations, only the part of interest of the calculation domain was solved, imposing velocity, turbulent kinetic energy and specific dissipation profiles on its inlet These conditions were imported from the full domain simulation already solved. To verify the correct capture of the main physical effects involved, comparisons with experimental results were made for 2 cases with representative flows: the Ahmed body and an aerodynamic profile with vortex generators. In addition, the differences between solving the complete or partial domain were studied in another comparative with experimental results. It was concluded that the numerical methodology was able to obtain the aerodynamic coefficients, and their tendencies against changes of geometry, in the cases studied. Partially solving the domain, imposing profiles at its entrance, resulted in a difference in the coefficients obtained in the order of 2% for the lift coefficient and 7% for the drag coefficient. The passive control via vortex generators was effective in delaying the separation of the boundary layer on the flap of the Formula SAE vehicle, the improvements in drag and lift coefficients were of the order of 7% and 0,3%, respectively.

Noise pollution threatens marine ecosystems, where animals rely heavily on sound for navigation and communication. The largest source of underwater noise from human activity is shipping, and propeller-induced cavitation is the dominant source of noise from ships. Mitigation strategies require accurate methods for predicting cavitation-induced noise, which remains challenging. The present thesis explores prediction and modelling strategies for cavitation-induced noise from marine propellers, and provides insight into models that can be used both during propeller design and to generate intelligent vessel control strategies. I examined three distinct approaches to predicting cavitation-induced noise, each of which is discussed in one of the three main chapters of this thesis: a high-fidelity computational fluid dynamics scheme, a parametric mapping procedure, and the use of field measurements. Each of these three chapters presents different insight into the acoustic behaviour of cavitating marine propellers, as well both real and potential strategies for mitigating this critical environmental emission. A combined experimental and numerical study of noise from a cavitating propeller, focused on both the fundamental importance of experimental findings and the effectiveness of the numerical modelling strategy used, is detailed in the first main chapter of this thesis. The experimental results highlighted that loud cavitation noise is not necessarily associated with high-power or high-speed propeller operation, affirming the need for intelligent vessel operation strategies to mitigate underwater noise pollution. Comparison of the experimental measurements and simulations revealed that the simulation strategy resulted in an over-prediction of sound levels from cavitation. Analysis of the numerical results and experiments strongly suggested that the cavitation model implemented in the simulations, a model commonly used for marine propeller simulations, was responsible for the over-prediction of sound levels. Ships are powered primarily by combustion engines, for which it is possible to generate "maps" relating the emission of pollutants to the engine’s speed and torque; the second main chapter of this thesis presents the methodology I developed for generating similar "maps" relating the level of cavitation-induced noise to the speed and torque of a ship's propeller. A proof-of-concept of the method that used the model propeller from the first main chapter is presented. To generate the maps, I used a low-order simulation technique to predict the cavitation induced by the propeller at a range of different speed and torque combinations. A pair of semi-empirical models found in the literature were combined to provide the framework for predicting noise based on cavitation patterns. The proof-of-concept map shows a clear optimal operating regime for the propeller. The final main chapter of this thesis presents an analysis of field noise measurements of coastal ferries in commercial operation, the data for which were provided by an industrial partner. The key finding was the identification of cavitation regime changes with variation in vessel speed by their acoustic signatures. The results provide a basis for remotely determining which vessels produce less noise pollution when subject to speed limits, which have been implement in critical marine habitats, and which vessels produce less noise at a specific optimum speed. / Graduate

No description available.

<p>Turbulent mixing layers have been researched for many years. Currently, research is focused on studying compressible mixing layers because of their widespread applications in high-speed flight systems. While the effect of compressibility on the shear layer growth rate is well established, there is a lack of consensus over its effect on the turbulent stresses and hence warrants additional research in this area. Computational studies on compressible shear layers could provide a deep cognizance of the dynamics of fluid structures present in these flow fields which in turn would be viable for understanding the effects of compressibility on such flows. However, performing a Direct Numerical Simulation (DNS) of a highly compressible shear layer with experimental flow conditions is extremely expensive, especially when resolving the boundary layers that lead into the mixing section. The attractive alternative is to use Large Eddy Simulation (LES), as it possesses the potential to resolve the flow physics at a reasonable computational cost. Therefore the current work deals with developing a methodology to perform LES of a compressible mixing layer with experimental flow conditions, with resolving the boundary layers that lead into the mixing section through a wall model. The wall model approach, as opposed to a wall resolved simulation, greatly reduces the computational cost associated with the boundary layer regions, especially when using an explicit time-stepping scheme. An in house LES solver which has been used previously for performing simulations of jets, has been chosen for this purpose. The solver is first verified and validated for mixing layer flows by performing simulations of laminar and incompressible turbulent mixing layer flows and comparing the results with the literature. Following this, LES of a compressible mixing layer at a convective Mach number of 0.53 is performed. The inflow profiles for the LES are taken from a precursor RANS solution based on the k-ε and RSM turbulence models. The results of the LES present good agreement with the reference experiment for the upstream boundary layer properties, the mean velocity profile of the shear layer and the shear layer growth rate. The turbulent stresses, however, have been found to be underpredicted. The anisotropy of the normal Reynolds stresses have been found to be in good agreement with the literature. Based on the present results, suggestions for future work are also discussed.</p>

This thesis discusses the application, evaluation, and extension of dimensionally adaptive meshing to the numerical solution of velocity and pressure fields inside cylindrical reactors. Due to the high length to diameter ratios of many cylindrical reactor vessels the flow field can become axisymmetric, allowing for simplification of the governing equations and significant reduction in computational time required for solution. A fully 3D solver is developed from existing computational tools at BYU and validated against theoretical velocity profiles for pipe flow at various Reynolds numbers, as well as with experimental data for an axial-fired center jet with recirculating flow. Dimensionally adaptive meshing is then incorporated into the validated 3D solver. The boundary conditions and assumptions at the dimensional boundary are discussed. The flow information is passed across the boundary through spatial mass-weighted averaging. The 3D and axisymmetric computational domains are decoupled from one another so information can only be passed from the 3D domain downstream to the axisymmetric domain. The dimensional boundary placement must meet two main requirements, the flow must be one-way and axisymmetric. It is found that the flow becomes axisymmetric early on in the reactor (~0.3-0.4 m), but recirculation exists farther downstream (until ~0.61 m) and thus governs the placement of the dimensional boundary. The resulting computational tool capable of running simulations using dimensionally adaptive meshes is called CEDAR (Computationally Efficient Dimensionally Adaptive Recirculating flow solver). Several studies are then undertaken to examine CEDAR's ability to reproduce exit velocity profiles comparable to those produced by a fully 3D mesh, including variations in pressure, firing rate, and geometry. It is found that the flow structure inside the reactor is self-similar over a wide range of operating parameters as long as the burner jets are turbulent. This observation is supported by free and confined jet theory. These theories also provide a method for placing the dimensional boundary, which is a linear function of the confining geometry diameter only (assuming that the jet diameter is less than 1/10 the diameter of the confining geometry). All exit velocity profiles produced by CEDAR are on average within 5% of the fully 3D profiles. Timing studies reveal an average 5.16 times speedup in computational time over fully 3D computations.

Modern gas turbine blades and vanes are operated at temperatures above their material’s melting point. Active external and internal cooling are therefore necessary to reach acceptable lifetimes. One possible internal cooling method is called matrix cooling, where a matrix of intersecting cooling air channels is integrated into a blade or vane. To further increase the efficiency of gas turbines, the amount of cooling air must be reduced. Therefore it is necessary that heat transfer inside a cooling matrix is well understood. In the first part of the thesis, a methodology for estimating heat transfer in the flow of matrix cooling channels was established using Computational Fluid Dynamics. Two four-equation RANS turbulence models based on the k-ε turbulence model showed a good correlation with experimental results, while the k-ω SST model underpredicted the heat transfer significantly. For all turbulence models, the heat transfer showed high sensitivity towards changes in the numerical setup. For the k-ω SST turbulence model, the mesh requirements were deemed too computationally expensive and it was excluded from further investigations. As the second part of the thesis, a parameter study was conducted investigating the influence of several geometric parameters on the heat transfer in a cooling matrix. The matrix was simplified as a channel flow interacting with multiple crossing flows. The highest enhancement in heat transfer was seen with changes in taper ratio, aspect ratio and matrix angle. Compared to smooth pipe flow, an increase in heat transfer of up to 60% was observed. Rounded edges of the cooling channels showed a significant influence on the heat transfer as well. In contrast, no influence of the wall thickness on the heat transfer was observed. While no direct validation is possible, the base case and the parameter sweeps show a good correlation with similar cases found in the literature.

A turbulence model was developed for computing surface velocity coefficients and discharge under steady, uniform flow conditions for rectangular and compound open-channel cross sections. Reynolds-Average Navier-Stokes (RANS) equations, Reynolds stress equations, and kinetic energy and dissipation equations were applied in the model using the finite-volume method with the SIMPLER algorithm. The models show graphical results of the velocity distributions in the longitudinal bed slope direction, secondary velocities, pressure, turbulence kinetic energy, and kinetic energy dissipation rate across the cross section. Also, the surface velocity coefficients were computed at increments of one-eighth of the base width from the vertical walls to the center of the cross section, and the submergence depth of the floating object from zero to 30 cm, with a 5-cm depth increment. Four different sets of Reynolds stress equations (one set by Boussinesq hypothesis and three sets of algebraic stress model) were used to calculate the results. Only one version of the algebraic stress model was successful in predicting the depression of the maximum streamwise velocity below the water surface. The model was calibrated and verified using laboratory data collected at Utah State University. Calculated discharges from the turbulence model had very good agreement with the laboratory data. The surface velocity coefficients from model results were generally lower than the results from the laboratory data, but higher than the values published by the United States Bureau of Reclamation. Standard cross sections of rectangular and compound cross sections were defined to simulate the model results and model sensitivity to parameter changes. The model results were summarized to show the relationship between surface velocity coefficient and channel characteristics compared with the published values by the USBR. For rectangular cross sections, the coefficients from the model are higher than the published USBR values. But the coefficients from the model and USBR are in very close agreement for the tested compound cross sections. The published coefficients by the USBR are a function of only average water depth. However, the model results show that the coefficients are also related to channel size, surface roughness height, float submergence depth, and lateral location of the float object. These factors should be included in the determination of the surface velocity coefficient to improve the discharge estimations from the application of the float method.

No description available.

No description available.

In Computational aeroacoustics, hybrid approaches first resolve the source and nearfield regions of the flow field by employing Reynolds Averaged Navier-Stokes (RANS) equations, Large Eddy Simulations (LES) or Direct Numerical Simulations (DNS).The source region data is used to form source terms, which are, in turn, applied to either empirical models or equations linearized around a mean flow. An acoustic analogy type of model is used to propagate the acoustics to the farfield regions. The aim of this research is twoold: to introduce and test a hybrid acoustic analogy, based on a coupling between the Navier-Stokes equations, applied in the source region, and the Non-linear Euler (NLE) equations applied in the acoustic propagation region; and to test and validate a recently derived generalized acoustic analogy theory in the framework of jet noise with acoustic source information obtained from RANS or LES. In the first part, the coupling between the Navier-Stokes and the NLE equations is accomplished via a buffer region, which is used to interpolate and penalize the flow variables of interest from the source region. The penalized flow variables are then applied as source terms in the NLE equations, to calculate the acoustic propagation. The non-linear Euler equations, discretized using highurate dispersion-relation preserving schemes constitute a very efficient approach for jet noise predictions in complex environments, especially for supersonic and hypersonic jets, where nonlinearities may propagate over long distances. In the second part, a RANS- or LES-informed model, which is used to provide data for Goldstein's generalized acoustic analogy, is presented. The generalized acoustic analogy of Goldstein is considered, wherein the effects of non-parallelism are taken into account and an asymptotic expansion is utilized to simplify the adjoint Green function equations. The use of the adjoint Green's function leads to a simple model for jet noise predictions for low frequencies and small observation angles, in the linear regime. Both approaches are extensively tested and validated against numerous benchmark problems and applications.