No description available.

The following study investigates hydrodynamic stability for two-dimensional, incompressible flow past a cylinder and compares it alongside four different variations of a wave-like ground introduced within the wake region of the cylinder wake. These different variations include changing the distance of the cylinder both horizontally from the wave-like structure and vertically from the ground. The geometry and meshes were initially constructed using GMSH and imported into Nektar++. The baseflows were then obtained in Nektar++ using the Velocity Correction Scheme, continuous Galerkin method, and Unsteady Navier Stokes solver. Then, the Implicitly Restarted Arnoldi Method driver was used to retrieve the various eigenvalues/eigenmodes and growth rates. Finally, the results were visualized in Paraview which allowed clear comparisons between the stability of the flow between each case. The findings obtained show a clear effect on stability when considering different cases, for a plain cylinder and for each case there are observations to be made in how the various eigenmodes varied in terms of magnitude and shape, other observations were made in the differing critical Reynolds number and frequencies among the cases. This study is relevant to various natural environments where a blunt object may come in range of a bumpy or wavy ground. In these scenarios it can be important to monitor how instabilities propagate and cause effects such as turbulence or drag. Additionally, investigation like these can detail how to effectively avoid undesirable characteristics of instability.

Wavlet analysis and its advantages in determining time-varying characteristics are discussed. The Morlet wavelet is defined and procedures for choosing its parameters are described. The recovery of modulation characteristics using the Morlet wavelet is demonstrated. Hydrodynamic linear stability is reviewed and its application to steady and unsteady mixing layers is discussed. Modulation effects are demonstrated by using the magnitude and phase of the wavelet coefficients. The time-varying characteristics of the most unstable modes are determined using the real part of the wavelet coefficients. It is found that mean flow unsteadiness increases the amplitude and phase modulation of the mixing layers. Synchronized variations of the two most unstable modes, the fundamental and the subharmonic, are also observed in the region of subharmonic growth. In a second application of wavelet analysis, the phase lag of the wavelet coefficients is used to determine the phase relation between the fundamental and the subharmonic in acoustically forced mixing layers. The results show that selective forcing affects the time-variations of the phase relation. In a third application, the magnitude and phase of the wavelet coefficients are used to decompose propagating waves measured at a single location. / Master of Science

The steady problem of free convective heat transfer from an isothermal inclined elliptic cylinder and its stability is investigated. The cylinder is inclined at an arbitrary angle with the horizontal and immersed in an unbounded, viscous, incompressible fluid. It is assumed that the flow is laminar and two-dimensional and that the Boussinesq approximation is valid. The full steady Navier-Stokes and thermal energy equations are transformed to elliptical co-ordinates and an asymptotic analysis is used to find appropriate far-field conditions. A numerical scheme based on finite differences is then used to obtain numerical solutions. Results are found for small to moderate Grashof and Prandtl numbers, and varying ellipse inclinations and aspect ratios. <br /><br /> A linear stability analysis is performed to determine the critical Grashof number at which the flow loses stability. Comparisons are made with long-time unsteady solutions.

The steady problem of free convective heat transfer from an isothermal inclined elliptic cylinder and its stability is investigated. The cylinder is inclined at an arbitrary angle with the horizontal and immersed in an unbounded, viscous, incompressible fluid. It is assumed that the flow is laminar and two-dimensional and that the Boussinesq approximation is valid. The full steady Navier-Stokes and thermal energy equations are transformed to elliptical co-ordinates and an asymptotic analysis is used to find appropriate far-field conditions. A numerical scheme based on finite differences is then used to obtain numerical solutions. Results are found for small to moderate Grashof and Prandtl numbers, and varying ellipse inclinations and aspect ratios. <br /><br /> A linear stability analysis is performed to determine the critical Grashof number at which the flow loses stability. Comparisons are made with long-time unsteady solutions.

The spatial stability of hypersonic conical boundary layers is investigated utilizing different numerical techniques. First, the development and verification of a Linearized Compressible Navier-Stokes solver (LinCS) is presented, followed by an investigation of different effects that affect the stability of the flow in free-flight/ground tests, such as: high-enthalpy effects, wall-temperature ratio, and three-dimensionality (i.e. angle-of-attack). A temporally/spatially high-order of accuracy parallelized Linearized Compressible Navier-Stokes solver in disturbance formulation was developed, verified and employed in stability investigations. Herein, the solver was applied and verified against LST, PSE and DNS, for different hypersonic boundary-layer flows over several geometries (e.g. flat plate - M=5.35 & 10; straight cone - M=5.32, 6 & 7.95; flared cone - M=6; straight cone at AoA = 6 deg - M=6). The stability of a high-enthalpy flow was investigated utilizing LST, LinCS and DNS of the experiments performed for a 5 deg sharp cone in the T5 tunnel at Caltech. The results from axisymmetric and 3D wave-packet investigations in the linear, weakly, and strongly nonlinear regimes using DNS are presented. High-order spectral analysis was employed in order to elucidate the presence of nonlinear couplings, and the fundamental breakdown of second mode waves was investigated using parametric studies. The three-dimensionality of the flow over the Purdue 7 deg sharp cone at M=6 and AoA =6 deg was also investigated. The development of the crossflow instability was investigated utilizing suction/blowing at the wall in the LinCS/DNS framework. Results show good agreement with previous computational investigations, and that the proper basic flow computation/formation of the vortices is very sensitive to grid resolution.

This thesis concerns the stability of plane Couette flow and pipe Poiseuille flow in three space dimensions. The mathematical model for both flows is the incompressible Navier--Stokes equations. Both analytical and numerical techniques are used. We present new results for the resolvent corresponding to both flows. For plane Couette flow, analytical bounds on the resolvent have previously been derived in parts of the unstable half-plane. In the remaining part, only bounds based on numerical computations in an infinite parameter domain are available. Due to the need for truncation of this infinite parameter domain, these results are mathematically insufficient. We obtain a new analytical bound on the resolvent at s=0 in all but a compact subset of the parameter domain. This is done by deriving approximate solutions of the Orr--Sommerfeld equation and bounding the errors made by the approximations. In the remaining compact set, we use standard numerical techniques to obtain a bound. Hence, this part of the proof is not rigorous in the mathematical sense. In the thesis, we present a way of making also the numerical part of the proof rigorous. By using analytical techniques, we reduce the remaining compact subset of the parameter domain to a finite set of parameter values. In this set, we need to compute bounds on the solution of a boundary value problem. By using a validated numerical method, such bounds can be obtained. In the thesis, we investigate a validated numerical method for enclosing the solutions of boundary value problems. For pipe Poiseuille flow, only numerical bounds on the resolvent have previously been derived. We present analytical bounds in parts of the unstable half-plane. Also, we derive a bound on the resolvent for certain perturbations. Especially, the bound is valid for the perturbation which numerical computations indicate to be the perturbation which exhibits largest transient growth. The bound is valid in the entire unstable half-plane. We also investigate the stability of pipe Poiseuille flow by direct numerical simulations. Especially, we consider a disturbance which experiments have shown is efficient in triggering turbulence. The disturbance is in the form of blowing and suction in two small holes. Our results show the formation of hairpin vortices shortly after the disturbance. Initially, the hairpins form a localized packet of hairpins as they are advected downstream. After approximately $10$ pipe diameters from the disturbance origin, the flow becomes severely disordered. Our results show good agreement with the experimental results. In order to perform direct numerical simulations of disturbances which are highly localized in space, parallel computers must be used. Also, direct numerical simulations require the use of numerical methods of high order of accuracy. Many such methods have a global data dependency, making parallelization difficult. In this thesis, we also present the process of parallelizing a code for direct numerical simulations of pipe Poiseuille flow for a distributed memory computer. / QC 20100825

In this thesis the orientation of macro-sized fibres in turbulent flows is studied, as well as the effect of nano-sized fibrils on hydrodynamic stability. The focus lies on enabling processes for new materials where cellulose is the main constituent. When fibres (or any elongated particles) are added to a fluid, the complexity of the flow-problem increases. The fluid flow will influence the rotation of the fibres, and therefore also effect the overall fibre orientation. Exactly how the fibres rotate depends to a large extent on the mean velocity gradient in the flow. In addition, when fibres are added to a suspending fluid, the total stress in the suspension will increase, resulting in an increased apparent viscosity. The increase in stress is related to the direction of deformation in relation to the orientation of the particle, i.e. whether the deformation happens along the long or short axis of the fibre. The increase in stress, which in most cases is not constant neither in time nor space, will in turn influence the flow. This thesis starts off with the orientation and spatial distribution of fibres in the turbulent flow down an inclined plate. By varying fibre and flow parameters it is discovered that the main parameter controlling the orientation distribution is the aspect ratio of the fibres, with only minor influences from the other parameters. Moreover, the fibres are found to agglomerate into streamwise streaks. A new method to quantify this agglomeration is developed, taking care of the problems that arise due to the low concentration in the experiments. It is found that streakiness, i.e. the tendency to agglomerate in streaks, varies with Reynolds number. Going from fibre orientation to flow dynamics of fibre suspensions, the influence of cellulose nanofibrils (CNF) on laminar/turbulent transition is investigated in three different setups, namely plane channel flow, curved-rotating channel flow, and the flow in a flow focusing device. This last flow case is selected since it is can be used for assembly of CNF based materials. In the plane channel flow, the addition of CNF delays the transition more than predicted from measured viscosities while in the curved-rotating channel the opposite effect is discovered. This is qualitatively confirmed by linear stability analyses. Moreover, a transient growth analysis in the plane channel reveals an increase in streamwise wavenumber with increasing concentration of CNF. In the flow focusing device, i.e. at the intersection of three inlets and one outlet, the transition is found to mainly depend on the Reynolds number of the side flow. Recirculation zones forming downstream of two sharp corners are hypothesised to be the cause of the transition. With that in mind, the two corners are given a larger radius in an attempt to stabilise the flow. However, if anything, the flow seems to become unstable at a smaller Reynolds number, indicating that the separation bubble is not the sole cause of the transition. The choice of fluid in the core flow is found to have no effect on the stability, neither when using fluids with different viscosities nor when a non-Newtonian CNF dispersion was used. Thus, Newtonian model fluids can be used when studying the flow dynamics in this type of device. As a proof of concept, a flow focusing device is used to produce a continuous film from CNF. The fibrils are believed to be aligned due to the extensional flow created in the setup, resulting in a transparent film, with an estimated thickness of 1 um. / <p>QC 20141003</p>

Shallow flows are those whose width is significantly larger than their depth. In these types of flows, two dimensional coherent structures can be generated and can influence the flow greatly by the lateral transfer of mass and momentum. The development of coherent structures as a result of flow instabilities has been a topic of interest for environmental fluid mechanics for decades. Studies on the use of linear modal stability analysis is commonly found in literature. However, the relatively recent development in the field of hydrodynamic stability suggests that the traditional linear modal stability analysis does not describe the behaviour of the perturbations in finite time. The discrepancy between asymptotic behaviour and finite time behaviour is particularly large in shear driven flows and it is most likely to be the case for shallow flows. This study aims to provide a better understanding of finite time growth of perturbation energy in shallow flows. The three cases of shallow flows evaluated are the mixing layer, jet and wake. The critical cases are obtained through the linear modal analysis and nonmodal analysis was conducted to show the transient behaviour in finite time for what is so-called marginally stable. Finally, the thesis concludes by generalizing the finite time energy growth in the S-k space.

Detonation is a supersonic mode of combustion that is modeled by a system of conservation laws of compressible fluid mechanics coupled with the equations describing thermodynamic and chemical properties of the fluid. Mathematically, these governing equations admit steady-state travelling-wave solutions consisting of a leading shock wave followed by a reaction zone. However, such solutions are often unstable to perturbations and rarely observed in laboratory experiments. The goal of this work is to study the stability of travelling-wave solutions of detonation models by the following novel approach. We linearize the governing equations about a base travelling-wave solution and solve the resultant linearized problem using high-order numerical methods. The results of these computations are postprocessed using dynamic mode decomposition to extract growth rates and frequencies of the perturbations and predict stability of travelling-wave solutions to infinitesimal perturbations. We apply this approach to two models based on the reactive Euler equations for perfect gases. For the first model with a one-step reaction mechanism, we find agreement of our results with the results of normal-mode analysis. For the second model with a two-step mechanism, we find that both types of admissible travelling-wave solutions exhibit the same stability spectra. Then we investigate the Fickett’s detonation analogue coupled with a particular reaction-rate expression. In addition to the linear stability analysis of this model, we demonstrate that it exhibits rich nonlinear dynamics with multiple bifurcations and chaotic behavior.