No description available.

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.

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.

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

Mon travail de thèse a porté sur la compréhension «systèmes dynamiques de la dynamique à grande échelle dans l’écoulement pleinement développé de cisaillement turbulent. Dans le plan écoulement de Couette, simulation des grandes échelles (LES) est utilisée pour modéliser petits mouvements d’échelle et de ne résoudre mouvements à grande échelle afin de calculer non linéaire ondes progressives (SNT) et orbites périodiques relatives (RPO). Artificiel sur-amortissement a été utilisé pour étancher une gamme croissante de petite échelle motions et prouvent que les motions grande échelle sont auto-entretenue. Les solutions d’onde inférieure branche itinérantes qui se trouvent sur le bassin laminaire turbulent limite sont obtenues pour ces simulation sur-amortie et continue encore dans l’espace de paramètre à des solutions de branche supérieure. Cette approche ne aurait pas été possible si, comme supposé dans certains enquêtes précédentes, les mouvements à grande échelle dans le mur bornées flux de cisaillement sont forcée par un mécanisme fondé sur l’existence de structures actives à plus petite échelle. En flux Poseuille, orbites périodiques relatives à décalage réflexion symétrie sur la limite du bassin laminaire turbulent sont calculés en utilisant DNS. Nous montrons que le RPO trouvé sont connectés à la paire de voyager vague (TW) solution via bifurcation mondiale (noeud-col-période infinie bifurcation). La branche inférieure de cette solution TW évoluer dans un état de l’envergure localisée lorsque le domaine de l’envergure est augmentée. La solution de branche supérieure développe plusieurs stries avec un espacement de l’envergure compatible avec des mouvements à grande échelle en régime turbulent. / My thesis work focused on ‘dynamical systems’ understanding of the large-scale dynamics in fully developed turbulent shear flow. In plane Couette flow, large-eddy simulation (L.E.S) is used to model small scale motions and to only resolve large-scale motions in order to compute nonlinear traveling waves (NTW) and relative periodic orbits (RPO). Artificial over-damping has been used to quench an increasing range of small-scale motions and prove that the motions in large-scale are self-sustained. The lower-branch traveling wave solutions that lie on laminar-turbulent basin boundary are obtained for these over-damped simulation and further continued in parameter space to upper branch solutions. This approach would not have been possible if, as conjectured in some previous investigations, large-scale motions in wall bounded shear flows are forced by mechanism based on the existence of active structures at smaller scales. In Poseuille flow, relative periodic orbits with shift-reflection symmetry on the laminar-turbulent basin boundary are computed using DNS. We show that the found RPO are connected to the pair of traveling wave (TW) solution via global bifurcation (saddle-node-infinite period bifurcation). The lower branch of this TW solution evolve into a spanwise localized state when the spanwise domain is increased. The upper branch solution develops multiple streaks with spanwise spacing consistent with large-scale motions in turbulent regime.

We study behavior of incompressible non-Newtonian fluids with a relation be- tween the shear stress and the shear rate given by a non-monotone S-shaped curve. These fluids are described with a special class of implicit constitutive relations that may be derived in a thermodynamically consistent manner us- ing the entropy production maximization principle or gradient dynamics. In the latter approach, the constitutive relation is given as the derivative of a non-convex dissipation potential. The concept of dissipation potential allows us to discuss stability of the constitutive relation and explain the experimen- tally observed response discontinuities. We are also concerned with hydrody- namic stability of flows of implicitly constituted fluids. Finally, we propose a numerical scheme for simulation of transient flows of fluids with a specific non-monotone constitutive relation. We employ the numerical scheme in a simulation of two-dimensional Taylor-Couette flow and the numerical results confirm our theoretical observations concerning the admissible flow states.