A common alternative to full resolution of the near-wall region in computational fluid dynamics (CFD) simulations is the use of wall functions that decrease the mesh requirements in this region. This study presents two alternatives to current wall functions. The first method is based on numerically approximating a turbulent velocity profile using a one-dimensional subgrid contained within walljacent control cells. The second method is an analytical approach similar to previous wall function methods, but this method is valid both inside and outside of the fluid boundary layer. Use of both methods allows approximation of boundary layers of varying height relative to the first layer sizing. Use of these methods allows wall adjacent primary grid sizes to vary from low-Re model sizing of y+ ≈ 1 to grid sizes of y+ ~ 1000 or more without significant loss in accuracy, and with computational costs similar to currently used wall functions.

Prediction of the near wall region is crucial to the accuracy of turbulent flow computational fluid dynamics (CFD) simulation. However, sufficient near-wall resolution is often prohibitive for high Reynolds number flows with complex geometries, due to high memory and processing requirements. A common approach in these cases is to use wall functions to bridge the region from the first grid node to the wall. This thesis presents an alternative method that relaxes the near wall resolution requirement by solving one dimensional transport equations for velocity and turbulence across a locally defined subgrid contained within wall adjacent grid cells. The addition of the subgrid allows for wall adjacent primary grid sizes to vary arbitrarily from low-Re model sizing (y+ ~ 1) to wall function sizing without significant loss of accuracy or increase in computational cost.

The primary purpose of this thesis is to quantify the effects of RANS turbulence modeling on the resolution of free shear vortical flows. The simulation of aerodynamic wing-tip vortices is used as a test bed. The primary configuration is flow over an isolated finite wing with aspect ratio, , and Reynolds number, . Tip-vortex velocity profiles, vortex core and wake turbulence levels, and Reynolds stresses are compared with wind tunnel measurements. Three turbulence models for RANS closure are tested: the Lumley, Reece, and Rodi full Reynolds stress transport model and the Sparlart-Allmaras model with and without a proposed modification. The main finding is that simulations with the full Reynolds stress transport model show remarkable mean flow agreement in the vortex and wake due to the proper prediction of a laminar vortex core. Simulations with the Spalart-Allmaras model did not indicate a laminar core and predicted over-diffusion of the tip-vortex. Secondary investigations in this work include the study of wall boundary layer treatment and simulating the wake-age of an isolated rotorcraft in hover using a steady-state RANS solver. By comparing skin friction plots over the NACA 0012 airfoil, it is shown that wall functions are most effective in the trailing edge half of the airfoil, while high velocity gradient and curvature of the leading edge make them more vulnerable to discrepancies. The rotorcraft simulation uses the modified Spalart-Allmaras turbulence model and shows proper, qualitative, resolution of the interaction between the vortex sheet and the tip vortex. / Master of Science

The present understanding of the structure and dynamics of turbulent boundary layers on aerodynamically smooth walls has been clarified over the last decade or so. However, the dynamics of turbulent boundary layers over rough surfaces is much less well known. Nevertheless, there are many industrial and environmental flow applications that require understanding of the mean velocity and turbulence in the immediate vicinity of the roughness elements.</p> <p>This thesis reports the effects of surface roughness on the flow characteristics in a turbulent boundary layer. Both experimental and numerical investigations are used in the present study. For the experimental study, comprehensive data sets are obtained for two-dimensional zero pressure-gradient turbulent boundary layers on a smooth surface and ten different rough surfaces created from sand paper, perforated sheet, and woven wire mesh. The physical size and geometry of the roughness elements and freestream velocity were chosen to encompass both transitionally rough and fully rough flow regimes. Three different probes, namely, Pitot probe, single hot-wire, and cross hot-film, were used to measure the velocity fields in the turbulent boundary layer. A Pitot probe was used to measure the streamwise mean velocity, while the single hot-wire and cross hot-film probes were used to measure the fluctuating velocity components across the boundary layer. The flow Reynolds number based on momentum thickness, , ranged from 3730 to 13,550. The data reported include mean velocity, streamwise and wall-normal turbulence intensities, Reynolds shear stress, triple correlations, as well as skewness and flatness factors. Different scaling parameters were used to interpret and assess both the smooth- and rough-wall data at different Reynolds numbers, for approximately the same freestream velocity. The appropriateness of the logarithmic law and power law proposed by George and Castillo (1997) to describe the mean velocity in the overlap region was also investigated. The present results were interpreted within the context of the Townsends wall similarity hypothesis. </p> <p>Based on the mean velocity data, a novel correlation that relates the skin friction to the ratio of the displacement and boundary layer thicknesses, which is valid for both smooth- and rough-wall flows, was proposed. In addition, it was also found that the application of a mixed outer scale caused the velocity profile in the outer region to collapse onto the same curve, irrespective of Reynolds numbers and roughness conditions. The present results showed that there is a common region within the overlap region of the mean velocity profile where both the log law and power law are indistinguishable, irrespective of the surface conditions. For the power law formulation, functional relationships between the roughness shift, and the power law coefficient and exponent were developed for the transitionally rough flows. The present results also suggested that the effect of surface roughness on the turbulence field depends to some degree on the specific characteristics of the roughness elements and also the component of the Reynolds stress tensor being considered. </p> <p>In the case of the numerical study, a new wall function formulation based on a power law was proposed for smooth and fully rough wall turbulent pipe flow. The new formulation correctly predicted the friction factors for smooth and fully rough wall turbulent pipe flow. The existing two-layer model realistically predicted the velocity shift on a log-law plot for the fully rough turbulent boundary layer. The two-layer model results also showed the effect of roughness is to enhance the level of turbulence kinetic energy and Reynolds shear stress compared to that on a smooth wall. This enhanced level extends into the outer region of the flow, which appears to be consistent with present and recent experimental results for the boundary layer.

The present understanding of the structure and dynamics of turbulent boundary layers on aerodynamically smooth walls has been clarified over the last decade or so. However, the dynamics of turbulent boundary layers over rough surfaces is much less well known. Nevertheless, there are many industrial and environmental flow applications that require understanding of the mean velocity and turbulence in the immediate vicinity of the roughness elements.</p> <p>This thesis reports the effects of surface roughness on the flow characteristics in a turbulent boundary layer. Both experimental and numerical investigations are used in the present study. For the experimental study, comprehensive data sets are obtained for two-dimensional zero pressure-gradient turbulent boundary layers on a smooth surface and ten different rough surfaces created from sand paper, perforated sheet, and woven wire mesh. The physical size and geometry of the roughness elements and freestream velocity were chosen to encompass both transitionally rough and fully rough flow regimes. Three different probes, namely, Pitot probe, single hot-wire, and cross hot-film, were used to measure the velocity fields in the turbulent boundary layer. A Pitot probe was used to measure the streamwise mean velocity, while the single hot-wire and cross hot-film probes were used to measure the fluctuating velocity components across the boundary layer. The flow Reynolds number based on momentum thickness, , ranged from 3730 to 13,550. The data reported include mean velocity, streamwise and wall-normal turbulence intensities, Reynolds shear stress, triple correlations, as well as skewness and flatness factors. Different scaling parameters were used to interpret and assess both the smooth- and rough-wall data at different Reynolds numbers, for approximately the same freestream velocity. The appropriateness of the logarithmic law and power law proposed by George and Castillo (1997) to describe the mean velocity in the overlap region was also investigated. The present results were interpreted within the context of the Townsends wall similarity hypothesis. </p> <p>Based on the mean velocity data, a novel correlation that relates the skin friction to the ratio of the displacement and boundary layer thicknesses, which is valid for both smooth- and rough-wall flows, was proposed. In addition, it was also found that the application of a mixed outer scale caused the velocity profile in the outer region to collapse onto the same curve, irrespective of Reynolds numbers and roughness conditions. The present results showed that there is a common region within the overlap region of the mean velocity profile where both the log law and power law are indistinguishable, irrespective of the surface conditions. For the power law formulation, functional relationships between the roughness shift, and the power law coefficient and exponent were developed for the transitionally rough flows. The present results also suggested that the effect of surface roughness on the turbulence field depends to some degree on the specific characteristics of the roughness elements and also the component of the Reynolds stress tensor being considered. </p> <p>In the case of the numerical study, a new wall function formulation based on a power law was proposed for smooth and fully rough wall turbulent pipe flow. The new formulation correctly predicted the friction factors for smooth and fully rough wall turbulent pipe flow. The existing two-layer model realistically predicted the velocity shift on a log-law plot for the fully rough turbulent boundary layer. The two-layer model results also showed the effect of roughness is to enhance the level of turbulence kinetic energy and Reynolds shear stress compared to that on a smooth wall. This enhanced level extends into the outer region of the flow, which appears to be consistent with present and recent experimental results for the boundary layer.

Čerpadla s nízkými specifickými otáčkami nacházejí uplatnění v široké škále aplikací, ale trpí nízkou účinností a rizikem nestability křivky dopravní výšky. Tato disertační práce pojednává o vylepšení hydraulických parametrů čerpadla se specifickými otáčkami ns = 32. Práce se zaměřuje na průtok v oběžném kole a ve spirále, které jsou řešeny pomocí CFD. Průzkum literatury ukázal, že přidání mezilopatek do průtokového kanálu je účinný způsob, jak zlepšit neuspokojivé proudění obecně přítomné v oběžných kolech s nízkými specifickými otáčkami. Byla zkoumána poloha mezilopatky v kanálu a vliv počtu přítomných mezilopatek. Část výsledků je ověřena experimentálně v laboratoři fakulty. Rovněž je zkoumán tvar spirály a pozornost je věnována numerickému řešení proudění u stěn a jeho vlivu na simulaci. Po numerických výpočtech na třech různých spirálách je navržen nový tvar, který zlepšuje vypočtové parametry čerpadla při nízkém i vysokém průtoku.

Despite the high cost of memory and CPU time required to resolve the boundary layer, a viscous unstructured grid solver has many advantages over a structured grid solver such as the convenience in automated grid generation and shock or vortex capturing by solution adaption. Since the geometry and flow phenomenon of a helicopter are very complex, unstructured grid-based methods are well-suited to model properly the rotor-fuselage interaction than the structured grid solver. In present study, an unstructured Cartesian grid solver is developed on the basis of the existing solver, NASCART-GT. Instead of cut-cell approach, immersed boundary approach is applied with ghost cell boundary condition, which increases the accuracy and minimizes unphysical fluctuations of the flow properties. The standard k-epsilon model by Launder and Spalding is employed for the turbulence modeling, and a new wall function approach is devised for the unstructured Cartesian grid solver. It is quite challenging and has never done before to apply wall function approach to immersed Cartesian grid. The difficulty lies in the inability to acquire smooth variation of y+ in the desired range due to the non-body-fitted cells near the solid wall. The wall function boundary condition developed in this work yields stable and reasonable solution within the accuracy of the turbulence model. The grid efficiency is also improved with respect to the conventional method. The turbulence modeling is validated and the efficiency of the developed boundary condition is tested in 2-D flow field around a flat plate, NACA0012 airfoil, axisymmetric hemispheroid, and rotorcraft applications. For rotor modeling, an actuator disk model is chosen, since it is efficient and is widely verified in the study of the rotor-fuselage interaction. This model considers the rotor as an infinitely thin disk, which carries pressure jump across the disk and allows flow to pass through it. The full three dimensional calculations of Euler and RANS equations are performed for the GT rotor model and ROBIN configuration to test implemented actuator disk model along with the developed turbulence modeling. Finally, the characteristics of the rotor-fuselage interaction are investigated by comparing the numerical solutions with the experiments.

No description available.

Dû aux efforts apportés à la réduction des émissions de NOx dans des chambres de combustion aéronautiques il y a une tendance récente vers des systèmes à combustion pauvre. Cela résulte dans l'apparition de nouveaux types d'injecteur qui sont caractérisés par une complexité géométrique accrue et par des nouvelles stratégies pour l'injection du carburant liquide, comme des systèmes multi-point. Les deux éléments créent des exigences supplémentaires pour des outils de simulation numériques. La simulation à grandes échelles (SGE ou LES en anglais) est aujourd’hui considérée comme la méthode la plus prometteuse pour capturer les phénomènes d'écoulement complexes qui apparaissent dans une telle application. Dans le présent travail, deux sujets principaux sont abordés : Le premier est le traitement de la paroi ce qui nécessite une modélisation qui reste délicate en SGE, en particulier dans des géométries complexes. Une nouvelle méthode d'implementation pour des lois de paroi est proposée. Une étude dans une géométrie réaliste démontre que la nouvelle formulation donne de meilleurs résultats comparé à l’implémentation classique. Ensuite, la capacité d'une approche SGE typique (utilisant des lois de paroi) de prédire la perte de charge dans une géométrie représentative est analysée et des sources d'erreur sont identifiés. Le deuxième sujet est la simulation du carburant liquide dans une chambre de combustion. Avec des méthodes Eulériennes et Lagrangiennes, deux approches sont disponibles pour cette tâche. La méthode Eulérienne considère un spray de gouttelettes comme un milieu continu pour lequel on peut écrire des équations de transport. Dans la formulation Lagrangienne, des gouttes individuelles sont suivies ce qui mène à des équations simples. D’autre part, sur le plan numérique, le grand nombre de gouttes à traiter peut s’avérer délicat. La comparaison des deux méthodes sous conditions identiques (solveur gazeux, modèles physiques) est un aspect central du présent travail. Les phénomènes les plus importants dans ce contexte sont l'évaporation ainsi que le problème d'injection d'un jet liquide dans un écoulement gazeux transverse ce qui correspond à une version simplifiée d’un système multi-point. Le cas d'application final est la configuration d’un seul injecteur aéronautique, monté dans un banc d'essai expérimental. Ceci permet d'appliquer de manière simultanée tous les développements préliminaires de ce travail. L'écoulement considéré est non-réactif mais à part cela il correspond au régime ralenti d'un moteur d'avion. Dû aux conditions préchauffées, le spray issu du système d'injection multi-point s'évapore dans la chambre. Cet écoulement est simulé utilisant les approches Eulériennes et Lagrangiennes et les résultats sont comparés aux données expérimentales. / Due to efforts to reduce NOx emissions of aeronautical combustors, there is a recent trend towards lean combustion technologies. This results in novel injector designs, which are characterized by increased geometrical complexity and new injection strategies for the liquid fuel, such as multipoint systems. Both elements create additional challenges for numerical simulation tools. Large-Eddy simulation (LES) is regarded as the most promising method to capture complex flow phenomena in such an application. In the present work, two main areas of interest are considered: The first is wall modeling, which remains a challenging field in LES, in particular for complex geometries. A new implementation method for wall functions that uses a no-slip condition at the wall is proposed. It is shown that in a realistic burner geometry the new formulation yields improved results compared to a classical implementation. Furthermore, the capability of a typical LES with wall models to predict the pressure drop in a representative geometry is assessed and sources of error are identified. The second topic is the simulation of liquid fuel in a combustor. With Eulerian and Lagrangian methods, two different approaches are available for this task. The Eulerian approach considers a droplet spray as a continuum for which transport equations can be formulated. In the Lagrangian formulation, individual droplets are tracked, which leads to a simple formulation but can be challenging in terms of numerics due to the large number of particles to be treated. The comparison of these methods under identical conditions (gaseous flow solver, physical models) is a central aspect of the present work. The most important phenomena that are studied in view of the final application are evaporation and the problem of transverse liquid jets in a gaseous crossflow as a simplified representation of a multipoint system. The final application case is the configuration of a single aeronautical injector mounted in an experimental test bench. It allows to simultaneously apply all preliminary developments. The flow considered is non-reactive but otherwise corresponds to a partial load regime in an aeroengine Due to the pre-heated conditions, the spray issued by the multi-point injection undergoes evaporation. This flow is simulated using Eulerian and Lagrangian methods and the results are compared to experimental data.

No description available.