F. Durst

Benchmark Computations of Laminar Flow Around a Cylinder

An overview of benchmark computations for 2D and 3D laminar flows around a cylinder is given, which have been defined for a comparison of different solution approaches for the incompressible Navier-Stokes equations developed within the Priority Research Programme. The exact definitions of the benchmarks are recapitulated and the numerical schemes and computers employed by the various participating groups are summarized. A detailed evaluation of the results provided is given, also including a comparison with a reference experiment. The principal purpose of the benchmarks is discussed and some general conclusions which can be drawn from the results are formulated.

Low-Reynolds-number flow around an oscillating circular cylinder at low Keulegan-Carpenter numbers

Time-averaged LDA measurements and time-resolved numerical flow predictions were performed to investigate the laminar flow induced by the harmonic in-line oscillation of a circular cylinder in water at rest. The key parameters, Reynolds number Re and Keulegan–Carpenter number KC , were varied to study three parameter combinations in detail. Good agreement was observed for Re =100 and KC =5 between measurements and predictions comparing phase-averaged velocity vectors. For Re =200 and KC =10 weakly stable and non-periodic flow patterns occurred, which made repeatable time-averaged measurements impossible. Nevertheless, the experimentally visualized vortex dynamics was reproduced by the two-dimensional computations. For the third combination, Re =210 and KC =6, which refers to a totally different flow regime, the computations again resulted in the correct fluid behaviour. Applying the widely used model of Morison et al . (1950) to the computed in-line force history, the drag and the added-mass coefficients were calculated and compared for different grid levels and time steps. Using these to reproduce the force functions revealed deviations from those originally computed as already noted in previous studies. They were found to be much higher than the deviations for the coarsest computational grid or the largest time step. The comparison of several in-line force coefficients with results obtained experimentally by Kuhtz (1996) for β=35 confirmed that force predictions could also be reliably obtained by the computations.

LDA measurements in the near-wall region of a turbulent pipe flow

This paper presents laser-Doppler measurements of the mean velocity and statistical moments of turbulent velocity fluctuations in the near-wall region of a fully developed pipe flow at low Reynolds numbers. A refractive-index-matched fluid was used in a Duran-glass test section to permit access to the near-wall region without distortion of the laser beams. All measurements were corrected for the influence of the finite size of measuring control volume. Measurements of long-time statistical averages of all three fluctuating velocity components in the near-wall region are presented. It is shown that the turbulence intensities in the wall region do not scale with inner variables. However, the limiting behaviour of the intensity components very close to the wall show only small variations with the Reynolds number. Measurements of higher-order statistical moments, the skewness and flatness factors, of axial and tangential velocity components confirm the limiting behaviour of these quantities obtained from direct numerical simulations of turbulent channel flow. The comparison of measured data with those obtained from direct numerical simulations reveals that noticeable discrepancies exist between them only with regard to the flatness factor of the radial velocity component near the wall. The measured v ’ flatness factor does not show the steep rise close to the wall indicated by numerical simulations. Analysis of the measured data in the near-wall region reveals significant discrepancies between the present LDA measurements and experimental results obtained using the hot-wire anemometry.

A PARALLEL BLOCK-STRUCTURED MULTIGRID METHOD FOR THE PREDICTION OF INCOMPRESSIBLE FLOWS

In this paper a parallel multigrid finite volume solver for the prediction of steady and unsteady flows in complex geometries is presented. For the handling of the complexity of the geometry and for the parallelization a unified approach connected with the concept of block-structured grids is employed. The parallel implementation is based on grid partitioning with automatic load balancing and follows the message-passing concept, ensuring a high degree of portability. A high numerical efficiency is obtained by a non-linear multigrid method with a pressure correction scheme as smoother. By a number of numerical experiments on various parallel computers the method is investigated with respect to its numerical and parallel efficiency. The results illustrate that the high performance of the underlying sequential multigrid algorithm can largely be retained in the parallel implementation and that the proposed method is well suited for solving complex flow problems on parallel computers with high efficiency.

Momentum and heat transfer from cylinders in laminar crossflow at 10-4 ≤ Re ≤ 200

Abstract A summary of results of numerical investigations of the two-dimensional flow around a heated circular cylinder located in a laminar crossflow is presented. Numerical investigations were carried out for the Reynolds number range 10 −4 ⩽ Re ⩽ 200 and for temperature loadings of 1.003–1.5. The computations yield information on Nu and C D variation with Reynolds number. The temperature dependence of the fluid properties (air) was taken into account and this resulted in a temperature dependence of the Nu – Re and C D – Re results. Information is also provided on the Strouhal number dependence on the Re number and on the critical Re number where vortex shedding starts.

The plane symmetric sudden-expansion flow at low Reynolds numbers

Detailed velocity measurements and numerical predictions are presented for the flow through a plane nominally two-dimensional duct with a Symmetric sudden expansion of area ratio 1:2. Both the experiments and the predictions confirm a symmetry-breaking bifurcation of the flow leading to one long and one short Separation zone for channel Reynolds numbers above 125, based on the upstream channel height and the maximum flow velocity upstream. With increasing Reynolds numbers above this value, the short separated region remains approximately constant in length whereas the long region increases in length. The experimental data were obtained using a one-component laser-Doppler anemometer at many Reynolds number values, with more extensive measurements being performed for the three Reynolds numbers 70, 300 and 610. Predictions were made using a finite volume method and an explicit quadratic Leith type of temporal discretization. In general, good agreement was found between measured and predicted velocity profiles for all Reynolds numbers investigated.

The Development Lengths of Laminar Pipe and Channel Flows

The authors’ research work into fully developed pulsating and oscillating laminar pipe and channel flows raised questions regarding the development length of the corresponding steady flow. For this development length, i.e., the distance from the entrance of the pipe to the axial position where the flow reaches the parabolic velocity profile of the Hagen-Poiseuille flow, a wide range of contradictory data exists. This is shown through a short review of the existing literature. Superimposed diffusion and convection, together with order of magnitude considerations, suggest that the normalized development length can be expressed as L∕D=C0+C1Re and for Re→0 one obtains C0=0.619, whereas for Re→∞ one obtains C1=0.0567. This relationship is given only once in the literature and it is presumed to be valid for all Reynolds numbers. Numerical studies show that it is only valid for Re→0 and Re→∞. The development length of laminar, plane channel flow was also investigated. The authors obtained similar results to those for the pipe flow: L∕D=C0′+C1′; Re, where C0′=0.631 and C1′=0.044. Finally, correlations are given to express L∕D analytically for the entire Re range for both laminar pipe and channel flows.

Backward-Facing Step Flows for Various Expansion Ratios at Low and Moderate Reynolds Numbers

This paper is concerned with the behavior of flows over a backward-facing step geometry for various expansion ratios H/h=1.9423, 2.5 and 3.0. A literature survey was carried out and it was found that the flow shows a strong two-dimensional behavior, on the plane of symmetry, for Reynolds numbers ReD=ρUbD/μ below approximately 400 (Ub= bulk velocity and D= hydraulic diameter). In this Reynolds number range, two-dimensional predictions were carried out to provide information on the general integral properties of backward-facing step flows, on mean velocity distributions and streamlines. Information on characteristic flow patterns is provided for a wide Reynolds number range, 10−4≤ReD≤800. In the limiting case of ReD→0, a sequence of Moffatt eddies of decreasing size and intensity is verified to exist in the concave corner also at ReD=1. The irreversible pressure losses are determined for various Reynolds numbers as a function of the expansion ratio. The two-dimensional simulations are known to underpredict the primary reattachment length for Reynolds numbers beyond which the actual flow is observed to be three-dimensional. The spatial evolution of jet-like flows in both the streamwise and the spanwise direction and transition to three-dimensionality were studied at a Reynolds number ReD=648. This three-dimensional analysis with the same geometry and flow conditions as reported by Armaly et al. (1983) reveals the formation of wall jets at the side wall within the separating shear layer. The wall jets formed by the spanwise component of the velocity move towards the symmetry plane of the channel. A self-similar wall-jet profile emerges at different spanwise locations starting with the vicinity of the side wall. These results complement information on backward-facing step flows that is available in the literature.

Experiments on the rise of air bubbles in clean viscous liquids

A number of new experiments have been performed on the rise of air bubbles in clean mixtures of distilled water and pure, reagent grade, glycerine covering a range of the relevant parameter, the Morton number, Mo = gv 4 ρ 3 /σ 3 , of 10 13 . Here g is the acceleration due to gravity, v the kinematic viscosity, ρ the density and σ the surface tension of the mixture. In these careful measurements several scaling regimes have been found that have not been discussed before in the extensive literature on the subject. The transitions between these regimes have been delineated and attempts made to discuss the dynamical processes that might be important in each of them.

Further contributions on the two-dimensional flow in a sudden expansion

Two-dimensional laminar flow of an incompressible viscous fluid through a channel with a sudden expansion is investigated. A continuation method is applied to study the bifurcation structure of the discretized governing equations. The stability of the different solution branches is determined by an Arnoldi-based iterative method for calculating the most unstable eigenmodes of the linearized equations for the perturbation quantities. The bifurcation picture is extended by computing additional solution branches and bifurcation points. The behaviour of the critical eigenvalues in the neighbourhood of these bifurcation points is studied. Limiting cases for the geometrical and flow parameters are considered and numerical results are compared with analytical solutions for these cases.