<p dir="ltr">The Boeing-AFOSR Mach-6 Quiet Tunnel (BAM6QT) was used for a set of experiments aiming to characterize the stability regimes of stagnation point injection in noisy and quiet flow across an array of different injected gases. Four gases were used in this experiment: air, helium, carbon dioxide, and argon. These gases were injected at varying thrust coefficients, ranging from 0.0516 to 0.5666, using a 7 degree half-angle cone with a 19 mm radius spherical nose and a single 1.93 mm-radius sonic jet in the center of the model. The primary data collected consists of schlieren images gathered at a sample rate of 76 kHz. These data were then analyzed using a shock tracking software to measure the physical locations of flow features as well as through spectral proper orthogonal decomposition (SPOD) to analyze specific modes in the flow.</p><p dir="ltr">Through this analysis, it was observed that three principle modes exist in stagnation point injection regardless of the injecting gas: a high frequency vortex-coupled mode, a low frequency Mach-shock-rigid mode, and a hybrid mode residing between these two modes. The first two modes were observed in all stability regimes, whereas the hybrid mode was only observed in the bifurcated regime. Furthermore, the unsteady regime was observed to be mostly characterized by this first, vortex-coupled mode. Conversely, the steady regime was observed to be driven by the Mach-shock-rigid mode instead. This transition was measured to occur as the thrust coefficient was increased.</p><p dir="ltr">This research also found that freestream noise resulted in an amplified and widened frequency range within the Mach-shock-rigid mode. This same freestream noise did not appear to have an impact on the other two principle modes; however, in some cases the noise produced in the Mach-shock-rigid mode due to this freestream noise did in fact mask the other principle modes.</p><p dir="ltr">Lastly, it was observed that the thrust coefficient, in and of itself, is not the sole indicator of stability in stagnation point injection. Across the different injected gases in this research, transition between the stability regimes did not in fact occur at a constant thrust coefficient value. Additionally, even within the same injected gas, this transition did not occur at the same thrust coefficient value between noisy and quiet runs—indicating an effect of freestream noise on stability.</p>

This thesis presents the implementation and validation of a new methodology developed by Glegg et al. (2004) for solving the trailing edge noise problem. This method is based on the premises that the noise produced by a surface can be computed by the integral of the cross product between the velocity and vorticity fields, of the boundary layer and shed vorticity (Howe (1978)). To extract the source terms, proper orthogonal decomposition is applied to the velocity cross spectrum to extract modes of the unsteady velocity and vorticity. The new formulation of the trailing edge noise problem by Glegg et al. (2004) is attractive because it applies to the high frequencies of interest but does not require an excessive computational effort. Also, the nature of the formulation permits the identification of the modes producing the noise and their associated velocity fluctuations as well as the regions of the boundary layer responsible for the noise production. The source terms were obtained using the direct numerical simulation of a turbulent channel flow by Moser et al. (1998). Two-point velocity and vorticity statistics of this data set were obtained by averaging 41 instantaneous fields. For comparisons purposes, experimental boundary layer data by Adrian et al. (2000) was chosen. Statistical reduction of 50 velocity fields obtained by particle image velocimetry was performed and analysis of the two-point correlation function showed features similar to the DNS data case. Also, proper orthogonal decomposition revealed identical dominant modes and eddy structures in the flow, therefore justifying considering the channel flow as an external boundary layer for noise calculations. Comparison of noise predictions with experimental data from Brooks et al. (1989) showed realistic results with the largest discrepancies, on the order of 5 dB, occurring at the lowest frequencies. The DNS results are least applicable at these frequencies, since these correspond to the longest streamwise lengthscales, which are the most affected by the periodicity conditions used in the DNS and also are the least representative of the turbulence in an external boundary layer flow. Most of the noise was shown to be produced by low-frequency streamwise velocity modes in the bottom 10% of the boundary layer and locations closest to the wall. Only 6 modes were required to obtain noise levels within 1 dB of the total noise. Finally, the method for predicting spatial velocity correlation from Reynolds stress data in wake flows, originally developed by Devenport et al. (1999, 2001) and Devenport and Glegg (2001), was adapted to boundary-layer type flows. This method, using Reynolds stresses and the prescription of a lengthscale to extrapolate the full two-point correlation, was shown to produce best results for a lengthscale prescribed as proportional to the turbulent macroscale. Noise predictions using modeled two-point statistics showed good agreement with the DNS inferred data in all but frequency magnitude, a probable consequence of the modeling of the correlation function in the streamwise direction. Other quantities associated to noise were seen to be similar to the ones obtained using the DNS. / Master of Science

The aim of this work was to understand the physical requirements as well as to develop methodology required to employ Time Resolved Digital Particle Image Velocimetry (TRDPIV) for measuring high speed, high magnification, near wall flow fields. Previous attempts to perform measurements such as this have been unsuccessful because of both limitations in equipment as well as proper methodology for processing of the data. This work addresses those issues and successfully demonstrates a test inside of a transonic turbine cascade as well as a high speed high magnification wall jet. From previous studies it was established that flow tracer delivery is not a trivial task in a high speed high back pressure environment. Any TRDPIV measurement requires uniform spatial seeding density, but time-resolved measurements require uniform temporal seeding density as well. To this end, a high pressure particle generator was developed. This advancement enhanced current capability beyond what was previously attainable. Unfortunately, this was not sufficient to resolve the issue of seeding all together, and an advanced data reconstruction methodology was developed to reconstruct areas of the flow field that where lost do to inhomogeneous seeding. This reconstruction methodology, based on Proper Orthogonal Decomposition (POD), has been shown to produce errors in corrected velocities below tradition spatial techniques alone. The combination of both particle generator and reconstruction methodology was instrumental for successfully acquiring TRDPIV measurements in a high speed high pressure environment such as a transonic wind tunnel facility. This work also investigates the development of a turbulent wall jet. This experiment helped in demonstrating the capability of taking high speed high magnification TRDPIV measurements. This experiment was very unique in that it is one of only a few experiments that studied the developing region of these jets. The Reynolds number ranged for this experiment from 150 – 10,000 which corresponded to velocities of 1 - 80 m/s. The results from this experiment showed good agreement with currently published time averaged data. Using scaling laws for fully developed jets a new scaling law was found for the developing region of the jet that could be applied to all Reynolds numbers in this study. A temporal investigation was also carried out using the temporal coefficients from POD. A vortex identification scheme was also applied to all of the Reynolds numbers showing clear trends as Reynolds number increased. / Master of Science

In this study the complete four dimensional space time correlation function was measured in the wake of an untripped circular cylinder at a Reynolds number of 60 000. This correlation serves as the complete inflow boundary condition for an open rotor ingesting inhomogeneous turbulence. An important aspect of the turbulence ingestion problem is understanding how different inflow boundary conditions effect the sound produced by a rotor. In the present study the turbulence structure of two plane wakes were compared. Measurements completed by a previous study in the wake of a NACA 0012 airfoil were compared with the measurements completed by the present study in the wake of a cylinder. The mean flows of both plane wakes were found to be very similar, however the Reynolds stress profiles show that the cylinder wake is substantially more turbulent. The structures of the two-point correlation function in each wake are also similar, although the cylinder wake had greater maximum correlation values and was correlated at greater separations. The two-point correlation was used along with proper orthogonal decomposition to compute the average instantaneous velocity fields of both wake flows. These velocity fields represent the average eddy structures present in each wake flow. The eddy structure comparisons show that the structures in the cylinder wake are larger and better correlated at longer time delays. / Master of Science / Any fan or propeller that ingests any unsteady flow will produce noise. This is especially important in propeller aircraft and marine vehicles where turbulence is generated from appendages on the vehicle’s body. This self-generated turbulence travels downstream and is eventually drawn into the propeller and produces noise. The broad study that the present work is a part of is concerned with understanding this ingestion noise problem so that the interaction can be better modeled and the sound produced can be predicted. To predict the sound produced by a fan or propeller ingesting turbulence, detailed information about the inflow condition is needed. In the present study the turbulence structure of the wake shed by a circular cylinder at 20 meters per second. The two-point velocity correlation in the wake serves as the complete inflow condition for the turbulence ingestion problem. The structure of the cylinder wake inflow condition was compared with the structure of an airfoil wake to evaluate how the differences in the two flows would influence the sound produced by a rotor ingesting the two conditions. The two flows were found to be quite similar in the mean flow. The cylinder wake was found to be significantly more turbulent than the airfoil wake and was correlated over greater distances. This suggests that the structures in the cylinder wake are larger and remain coherent longer than those in the airfoil wake. The average instantaneous velocity fields were estimated in both wake flows and showed that the structures in the cylinder wake were significantly different from the structures in the airfoil wake. These flow structure comparisons show why the differences seen in the turbulence profiles and two-point correlations exist.

This thesis uses Proper Orthogonal Decomposition to model parametric and coupled systems. First, Proper Orthogonal Decomposition and its properties are introduced as well as how to numerically compute the decomposition. Next, a test case was used to show how well POD can be used to simulate and control a system. Finally, techniques for modeling a parametric system over a given range and a coupled system split into subdomains were explored, as well as numerical results. / Master of Science

Proper Orthogonal Decomposition (POD), combined with the Method of Snapshots and Galerkin projection, is a popular method for the model order reduction of nonlinear PDEs. The POD requires the left singular vectors from the singular value decomposition (SVD) of an n-by-m "snapshot matrix" S, each column of which represents the computed state of the system at a given time. However, the direct computation of this decomposition can be computationally expensive, particularly for snapshot matrices that are too large to fit in memory. Hierarchical Approximate POD (HAPOD) (Himpe 2016) is a recent method for the approximate truncated SVD that requires only a single pass over S, is easily parallelizable, and can be computationally cheaper than direct SVD, all while guaranteeing the requested accuracy for the resulting basis. This method processes the columns of S in blocks based on a predefined rooted tree of processors, concatenating the outputs from each stage to form the inputs for the next. However, depending on the selected parameter values and the properties of S, the performance of HAPOD may be no better than that of direct SVD. In this work, we numerically explore the parameter values and snapshot matrix properties for which HAPOD is computationally advantageous over the full SVD and compare its performance to that of a parallelized incremental SVD method (Brand 2002, Brand 2003, and Arrighi2015). In particular, in addition to the two major processor tree structures detailed in the initial publication of HAPOD (Himpe2016), we explore the viability of a new structure designed with an MPI implementation in mind. / Master of Science / Singular Value Decomposition (SVD) provides a way to represent numeric data that breaks the data up into its most important components, as well as measuring how significant each part is. This decomposition is widely used to assist in finding patterns in data and making decisions accordingly, or to obtain simple, yet accurate, representations of complex physical processes. Examples of useful data to decompose include the velocity of water flowing past an obstacle in a river, a large collection of images, or user ratings for a large number of movies. However, computing the SVD directly can be computationally expensive, and usually requires repeated access to the entire dataset. As these data sets can be very large, up to hundreds of gigabytes or even several terabytes, storing all of the data in memory at once may be infeasible. Thus, repeated access to the entire dataset requires that the files be read repeatedly from the hard disk, which can make the required computations exceptionally slow. Fortunately, for many applications, only the most important parts of the data are needed, and the rest can be discarded. As a result, several methods have surfaced that can pick out the most important parts of the data while accessing the original data only once, piece by piece, and can be much faster than computing the SVD directly. In addition, the recent bottleneck in individual computer processor speeds has motivated a need for methods that can efficiently run on a large number of processors in parallel. Hierarchical Approximate POD (HAPOD) [1] is a recently-developed method that can efficiently pick out the most important parts of the data while only accessing the original data once, and which is very easy to run in parallel. However, depending on a user-defined algorithm parameter (weight), HAPOD may return more information than is needed to satisfy the requested accuracy, which determines how much data can be discarded. It turns out that the input weights that result in less extra data also result in slower computations and the eventual need for more data to be stored in memory at once. This thesis explores how to choose this input weight to best balance the amount of extra information used with the speed of the method, and also explores how the properties of the data, such as the size of the data or the distribution of levels of significance of each part, impact the effectiveness of HAPOD.

With the maturity of conventional industrial robots, there has been increasing interest in designing robots that emulate realistic animal motions. This discipline requires careful and systematic investigation of a wide range of animal motions from biped, to quadruped, and even to serpentine motion of centipedes, millipedes, and snakes. Collecting optical motion capture data of such complex animal motions can be complicated for several reasons. Often there is the need to use many high-quality cameras for detailed subject tracking, and self-occlusion, loss of focus, and contrast variations challenge any imaging experiment. The problem of self-occlusion is especially pronounced for animals. In this thesis, we walk through the process of collecting motion capture data of a running lizard. In our collected raw video footage, it is difficult to make temporal correspondences using interpolation methods because of prolonged blurriness, occlusion, or the limited field of vision of our cameras. To work around this, we first make a model data set by making our best guess of the points' locations through these corruptions. Then, we randomly eclipse the data, use Gappy POD to repair the data and then see how closely it resembles the initial set, culminating in a test case where we simulate the actual corruptions we see in the raw video footage. / Master of Science / There has been increasing interest over the past few years in designing robots that emulate realistic animal motions. To make these designs as accurate as possible requires thorough analysis of animal motion. This is done by recording video and then converting it into numerical data, which can be analyzed in a rigorous way. But this conversion cannot be made when the raw video footage is ambiguous, for instance, when the footage is blurry, the shot is too dark or too light, the subject (or parts of the subject) are out of view of the camera, etc. In this thesis, we walk through the process of collecting video footage of a lizard running and then converting it into data. Ambiguities in the video footage result in an incomplete translation into numerical data and we use a mathematical technique called the Gappy Proper Orthogonal Decomposition to fill in this incompleteness in an intelligible way. And in the process, we lay your hands on the fundamental drivers of the animal’s motion.

The thesis deals with the flow in pipe bends and radial turbines geometries that are commonly found in an Internal Combustion Engine (ICE). The development phase of internal combustion engines relies more and more on simulations as an important complement to experiments. This is partly because of the reduction in development cost and the shortening of the development time. This is one of the reasons for the need of more accurate and predictive simulations. By using more complex computational methods the accuracy and predictive capabilities are increased. The disadvantage of using more sophisticated tools is that the computational time is increasing, making such tools less attractive for standard design purposes. Hence, one of the goals of the work has been to contribute to assess and improve the predictive capability of the simpler methods used by the industry. By comparing results from experiments, Reynolds Averaged Navier-Stokes (RANS) computations, and Large Eddy Simulations (LES) the accuracy of the different computational methods can be established. The advantages of using LES over RANS for the flows under consideration stems from the unsteadiness of the flow in the engine manifold. When such unsteadiness overlaps the natural turbulence the model lacks a rational foundation. The thesis considers the effect of the cyclic flow on the chosen numerical models. The LES calculations have proven to be able to predict the mean field and the fluctuations very well when compared to the experimental data. Also the effects of pulsatile exhaust flow on the performance of the turbine of a turbocharging system is assessed. Both steady and pulsating inlet conditions are considered for the turbine case, where the latter is a more realistic representation of the real flow situation inside the exhaust manifold and turbine. The results have been analysed using different methods: single point Fast Fourier Transforms (FFT), probe line means and statistics, area and volume based Proper Orthogonal Decomposition (POD) and Dynamic Mode Decomposition (DMD). / Denna avhandling behandlar flödet i rörkrökar och radiella turbiner som vanligtvis återfinns i en förbränningsmotor. Utvecklingsfasen av förbränningsmotorer bygger mer och mer på att simuleringar är ett viktigt komplement till experiment. Detta beror delvis på minskade utvecklingskostnader men även på kortare utevklningstider. Detta är en av anledningarna till att man behöver mer exakta och prediktiva simuleringsmetoder. Genom att använda mer komplexa beräkningsmetoder så kan både nogrannheten och prediktiviteten öka. Nackdelen med att använda mer sofistikerade metoder är att beräkningstiden ökar, vilket medför att sådana verktyg är mindre attraktiva för standardiserade design ändamål. Härav, ett av målen med projektet har varit att bidra med att bedöma och förbättra de enklare metodernas prediktionsförmåga som används utav industrin. Genom att jämföra resultat från experiment, Reynolds Averaged Navier-Stokes (RANS) och Large Eddy Simulations (LES) så kan nogrannheten hos de olika simuleringsmetoderna fastställas. Fördelarna med att använda LES istället för RANS när det gäller de undersökta flödena kommer ifrån det instationära flödet i grenröret. När denna instationäritet överlappar den naturligt förekommande turbulensen så saknar modellen en rationell grund. Denna avhandling behandlar effekten av de cykliska flöderna på de valda numeriska modellerna. LES beräkningarna har bevisats kunna förutsäga medelfältet och fluktuationerna väldigt väl när man jämför med experimentell data. Effekterna som den pulserande avgasströmning har på turboladdarens turbin prestanda har också kunnat fastställas. Både konstant och pulserande inlopps randvillkor har används för turbinfallet, där det senare är ett mer realistiskt representation av den riktiga strömningsbilden innuti avgasgrenröret och turbinen. Resultaten har analyserats på flera olika sätt: snabba Fourier transformer (FFT) i enskilda punkter, medelvärden och statistik på problinjer, area och volumsbaserade metoder så som Proper Orthogonal Decomposition (POD) samt Dynamic Mode Decomposition (DMD). / <p>QC 20140919</p>

In many areas of aerodynamics the technique of Large Eddy Simulation (LES) has proved a practical way of modelling the unsteady phenomena in numerical simulations. Few applications are as dependent on such an approach as the prediction of flow within a gas turbine combustor. Like any form of Computational Fluid Dynamics (CFD), LES requires specification of the velocity field at the inflow boundary, with much evidence suggesting the specification of inlet turbulence can be critical to the resultant accuracy of the prediction. While a database of time-resolved velocity data may be obtained from a precursor LES calculation, this technique is prohibitively expensive for complex geometries. An alternative is to use synthetic inlet conditions obtained from experimental data High-speed Particle Image Velocimetry (PIV) is used here to provide planar velocity data at up to 1kHz temporal resolution in two test cases representative of gas turbine combustor flows (a vortex generator in a duct and an idealised combustor). As the data sampling rate is approaching a typical LES time-step it introduces the possibility of applying instantaneous experimental data directly as an inlet condition. However, as typical solution domain inlet regions for gas turbine combustor geometries cannot be adequately captured in a single field of PIV data, it is necessary to consider a method by which a synchronous velocity field may be obtained from multiple PIV fields that were not captured concurrently. A method is proposed that attempts to achieve this by a combined process of Linear Stochastic Estimation and high-pass filtering. The method developed can be generally applied without a priori assumptions of the flow and is demonstrated to produce a velocity field that matches very closely that of the original PIV, with no discontinuities in the velocity correlations. The fidelity and computational cost of the method compares favourably to several existing inlet condition generation methods. Finally, the proposed and existing methods for synthetic inlet condition generation are applied to LES predictions of the two test cases. There is shown to be significant differences in the resulting flow, with the proposed method showing a marked ii reduction in the adjustment period that is required to establish turbulent equilibrium downstream of the inlet. However, it is noted the presence of downstream turbulence generating features can mask any differences in the inlet condition, to the extent that the flow in the core of the combustor test case is found to be insensitive to the inlet condition applied at the entry to the feed annulus for the test conditions applied here.

Assessing the impact of inlet flow distortion in turbomachinery is desired early in the design cycle. This thesis introduces and validates the use of methods based on the Proper Orthogonal Decomposition (POD) to analyze clean and 1/rev static pressure distortion simulation results at design and near stall operating condition. The value of POD comes in its ability to efficiently extract both quantitative and qualitative information about dominant spatial flow structures as well as information about temporal fluctuations in flow properties. Observation of the modes allowed qualitative identification of shock waves as well as quantification of their location and range of motion. Modal coefficients revealed the location of the passage shock at a given angular location. Distortion amplification and attenuation between rotors was also identified. A relationship was identified between how distortion manifests itself based on downstream conditions. POD provides an efficient means for extracting the most meaningful information from large CFD simulation data. Static pressure and axial velocity were analyzed to explore the flow physics of 3 rotors of a compressor with a distorted inlet. Based on the results of the analysis of static pressure using the POD modes, it was concluded that there was a decreased range of motion in passage shock oscillation. Analysis of axial velocity POD modes revealed the presence of a separated region on the low pressure surface of the blade which was most dynamic in rotor 1. The thickness of this structure decreased in the near stall operating condition. The general conclusion is made that as the fan approaches stall the apparent effects of distortion are lessened which leads to less variation in the operating condition. This is due to the change in operating condition placing the fan at a different position on the speedline such that distortion effects are less pronounced. POD modes of entropy flux were used to identify three distinct levels of entropy flux in the blade row passage. The separated region was the region with the highest entropy due to the irreversibilities associated with separation.