No description available.

This thesis deals with numerical simulation of turbulentflows in geometrically simple cases. Both plane and cylindricalgeometries are used. The simplicity of the geometry allows theuse of spectral methods which yield a very high accuracy usingrelatively few grid points. A spectral method for planegeometries is implemented on a parallel computer. Thetransitional Reynolds number for plane Couette flow is verifiedto be about 360, in accordance with earlier findings. TurbulentCouette flow at twice the transitional Reynolds number isstudied and the findings of large scale structures in earlierstudies of Couette flow are substantiated. These largestructures are shown to be of limited extent and give anintegral length scale of six half channel heights, or abouteight times larger than in pressure-driven channel flow.Despite this, they contain only about 10 \% of the turbulentenergy. This is demonstrated by applying a very smallstabilising rotation, which almost eliminates the largestructures. A comparison of the Reynolds stress budget is madewith a boundary layer flow, and it is shown that the near-wallvalues in Couette flow are comparable with high-Reynolds numberboundary layer flow. A new spectrally accurate algorithm isdeveloped and implemented for cylindrical geometries andverified by studying the evolution of eigenmodes for both pipeflow and annular pipe flow. This algorithm is a generalisationof the algorithm used in the plane channel geometry. It usesFourier transforms in two homogeneous directions and Chebyshevpolynomials in the third, wall-normal, direction. TheNavier--Stokes equations are solved with a velocity-vorticityformulation, thereby avoiding the difficulty of solving for thepressure. The time advancement scheme used is a mixedimplicit/explicit second order scheme. The coupling between twovelocity components, arising from the cylindrical coordinates,is treated by introducing two new components and solving forthem, instead of the original velocity components. TheChebyshev integration method and the Chebyshev tau method isboth implemented and compared for the pipe flow case.

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

<p>This thesis deals with numerical simulation of turbulentflows in geometrically simple cases. Both plane and cylindricalgeometries are used. The simplicity of the geometry allows theuse of spectral methods which yield a very high accuracy usingrelatively few grid points. A spectral method for planegeometries is implemented on a parallel computer. Thetransitional Reynolds number for plane Couette flow is verifiedto be about 360, in accordance with earlier findings. TurbulentCouette flow at twice the transitional Reynolds number isstudied and the findings of large scale structures in earlierstudies of Couette flow are substantiated. These largestructures are shown to be of limited extent and give anintegral length scale of six half channel heights, or abouteight times larger than in pressure-driven channel flow.Despite this, they contain only about 10 \% of the turbulentenergy. This is demonstrated by applying a very smallstabilising rotation, which almost eliminates the largestructures. A comparison of the Reynolds stress budget is madewith a boundary layer flow, and it is shown that the near-wallvalues in Couette flow are comparable with high-Reynolds numberboundary layer flow. A new spectrally accurate algorithm isdeveloped and implemented for cylindrical geometries andverified by studying the evolution of eigenmodes for both pipeflow and annular pipe flow. This algorithm is a generalisationof the algorithm used in the plane channel geometry. It usesFourier transforms in two homogeneous directions and Chebyshevpolynomials in the third, wall-normal, direction. TheNavier--Stokes equations are solved with a velocity-vorticityformulation, thereby avoiding the difficulty of solving for thepressure. The time advancement scheme used is a mixedimplicit/explicit second order scheme. The coupling between twovelocity components, arising from the cylindrical coordinates,is treated by introducing two new components and solving forthem, instead of the original velocity components. TheChebyshev integration method and the Chebyshev tau method isboth implemented and compared for the pipe flow case.</p>

The study of turbulence has been dominated historically by a bottom-up approach, with a much stronger emphasis on the physical structure of flows than on that of the dynam- ical state space. Turbulence has traditionally been described in terms of various visually recognizable physical features, such as waves and vortices. Thanks to recent theoretical as well as experimental advancements, it is now possible to take a more top-down approach to turbulence. Recent work has uncovered non-trivial equilibria as well as relative periodic orbits in several turbulent systems. Furthermore, it is now possible to verify theoretical results at a high degree of precision, thanks to an experimental technique known as Particle Image Velocimetry. These results squarely frame moderate Reynolds number Re turbulence in boundary shear flows as a tractable dynamical systems problem. In this thesis, I intend to elucidate the finer structure of the state space of moderate Re wall-bounded turbulent flows in hope of providing a more accurate and precise description of this complex phenomenon. Computation of new undiscovered equilibria, relative equilibria, and their heteroclinic connections provide a skeleton upon which a numerically accurate description of turbulence can be framed. The behavior of the equilibria under variation of Reynolds number and cell aspect ratios is also examined. It is hoped that this description of the state space will provide new avenues for research into nonlinear control systems for shear flows as well as quantitative predictions of transport properties of moderate Re fluid flows.

No description available.