Simulation of Separating Flows over a Backward-Facing Step

Simulation results are reported for plane two-dimensional viscous incompressible flow in a channel with an abrupt expansion. The mathematical model is provided by the quasi-hydrodynamic equations in the incompressible fluid approximation. The computations are carried out in a range of Reynolds numbers including both laminar and turbulent flow. As the Reynolds number increases, the solution bifurcates and the steady laminar flow changes to time-dependent flow. The computation results are consistent with known experimental data. Turbulence models were not used for large Reynolds number computations.     

Interaction between inner and outer layer in drag-reduced turbulent flows

Effects of wall-based skin-friction drag reduction strategies on the statistical properties of large-scale motions in moderate Reynolds number turbulent flows have been investigated by exploiting Direct Numerical Simulation of turbulent channels. To educe large scales, a new efficient parallel distributed memory algorithm has been implemented which delivers data-driven modes of increasing characteristic lengthscales: the Fast and Adaptive Bidimensional Empirical Mode Decomposition (FABEMD). The influence of wall-based skin friction reduction on large scales is studied by comparing single point statistics, such as r.m.s. fluctuations, and two-point statistics, as cross-correlation functions in controlled and uncontrolled channel flow fields at constant friction Reynolds number. The traditional way of observing large-scale footprinting at the wall, as cross-correlation of the streamwise velocity components at different wall distances, has been found to be unreliable when comparing drag-reduced flows, due to the arbitrary choice of a reference plane in the logarithmic layer. A more sound way of observing the footprinting via the correlation of the streamwise velocity with the friction velocity is addressed and shows an increase of the footprinting in drag-reduced flows. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

SOME REMARKS AID PROBLEMS CONCERNING THE GEOGRAPHY OF FANO 4-FOLDS OF INDEX AND PICARDNUMBER ONE

The two-dimensional steady flow of a viscous incompressible fluid in a diverging symmetrical channel is examined. The paper exploits a new series summation and improvement technique (i.e. Drazin and Tourigny, 1996). The solutions are expanded into Taylor series with respect to the corresponding Reynolds number and the bifurcation study is perfomed. Parameter ranges for the Reynolds number, where no, one or two solutions of the given type exist, are computed.     

Influence of Mass Loading on Particle-Laden Turbulent Pipe Flow

A CFD code in the framework of OpenFOAM was validated for simulations of particle-laden pipe and channel flows at low to intermediate mass loadings. The code is based on an Eulerian two-fluid approach with Reynolds-averaged conservation equations, including turbulence modeling and four-way coupling. Pipe flow simulations of particles in air against gravity were conducted at Reynolds numbers up to 50000. The particle mass loading was varied and its effect on the mean velocities and turbulent fluctuations of the two phases was studied. Special attention was paid to the influence of mass loading on the centerline velocity and the wall shear velocity of the fluid phase for various flow parameters and particle properties. Empirical correlations were established between these two quantities and the flow Reynolds number, particle Reynolds number, Stokes number and particle to fluid density ratio for a range of particle mass loadings. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical investigation of vortical evolution in a backward-facing step expansion flow

A numerical investigation of laminar flow over a backward-facing step is presented for the Reynolds number in the range of 50Re2500. The objective of this numerical investigation is to add to the existing knowledge of the backward-facing step flow to deepen our understanding of the expansion flow structure. We proceed with the analysis by verifying the computer code through the Pearson vortex problem. We then perform a parametric study by varying the Reynolds number, with the aim of determining whether or not there exists a critical Reynolds number, above which reattachment length on the channel floor decreases. We also concentrate on subjects that have been little explored in the flow, examples of which are the onset of a single vortex in the primary eddy and how the recirculating bubble containing flow reversals is torn into smaller eddies. Eddy distortion, leading to mobile saddle points, and the merging of eddies are also discussed in this study.     

Prandtl's Mixing Length Model - Revisited

The paper is concerned with a modification of Prandtl's mixing length model of Reynolds stresses in fully developed turbulent channel flows. Here it is a well established fact that Prandtl's model falls short to describe the Reynolds stresses correctly very close to the wall. Furthermore, the resulting solution of the time averaged velocity fails to describe the channel flow correctly from the wall to its center. To overcome these shortcomings, the only characteristic mixing length in Prandtl's model is replaced by separate mixing lengths for velocity fluctuations parallel to the wall and normal to the wall, respectively. The modified model describes the mean velocity, all Reynolds stresses, and the functional dependence between Reynolds number based on the mean velocity and the one based on the friction velocity. For all Reynolds numbers these results - and corresponding results for the production terms of Reynolds stresses and the energy balance of the mean flow - agree quantitatively with experimental data and with data obtained by numerical simulations. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Effect of Deborah number and phase difference on peristaltic transport of a third-order fluid in an asymmetric channel

The effect of a third-order fluid on the peristaltic transport in an asymmetric channel is studied. The wavelength of the peristaltic waves is assumed to be large compared to the varying channel width, whereas the wave amplitudes need not be small compared to the varying channel width. The channel asymmetry is produced by choosing the peristaltic wave train on the walls to have different amplitudes and phase. The flow is investigated in a wave frame of reference moving with velocity of the wave. The effects of Deborah number, phase difference, varying channel width and wave amplitudes on the pumping characteristics, streamline pattern and trapping phenomena are investigated. It is observed that the trapping regions increase as the channel becomes more and more symmetric and the trapped bolus volume decreases for increasing Deborah number, phase difference and varying channel width whereas it increases for increasing flow rate and wave amplitudes. Furthermore, the obtained results could also have applications to a range of peristaltic flows for a variety of non-Newtonian fluids such as aqueous solutions of high-molecular weight polyethylene oxide and polyacrylamide.     

On the Influence of Thermal Convection in Classical Taylor‐Couette Flows

The presented work deals with the instabilities that occur in the flow of a viscous fluid between axisymmetric cylinders with a rotating inner and stationary outer cylinder. The results of a numerical study of convective flows are presented. The inner cylinder is rotating and heated from within, while the outer cylinder is stationary and cooled outside. Stationary horizontal endplates are used to seal the annulus, forming an enclosure. The working fluid is silikon oil M3. The flow of oil was rendered visible by injecting aluminium powder. By increasing the Reynolds number with angular velocity of the driving inner cylinder, the flow bifurcates into different types of instabilities. Investigation was aimed to find the values of critical Reynolds and Rayleigh numbers corresponding to the critical speeds and temperature differences at which these instabilities set in. The three‐dimensional problem was modelled numerically using software package FLUENT in which discretization is performed by means of finite volume techniques. Computational grid was created in preprocessor Gambit. Numerical experiments are conducted to determine the interdependence between the heat transfer mechanism and the structure of secondary flows     

Self-similar turbulent boundary layer in pressure gradient

Self-similar flows in a turbulent boundary layer when the free-stream velocity is specified as a power function of longitudinal coordinate are investigated. The self-similar formulation not only simplifies solving of the problem by reducing the equations of motion to ordinary differential equations but also provides a mean for formulating closure conditions. It is shown that for the class of flows under consideration that depend on three governing parameters the dimensionless mixing length is a function of the normalised distance from the wall and the exponent in the law specifying the free-stream velocity distribution in the outer region and a universal function of local Reynolds number in the wall region, the latter corollary being true even when the skin friction vanishes. In calculations this function is set to be independent of pressure gradient, which gives the results very close to experimental data. There exist four different self-similar flow regimes. Each regime is related to its similarity parameter, one of which is the well-known Clauser equilibrium parameter and the other three are established for the first time. In case of adverse pressure gradient when the exponent lies within certain limits, which depend on Reynolds number, the problem has two solutions with different values of the boundary layer thickness and skin friction, which points out the possibility of hysteresis in near-separating flow. Separation occurs not at the minimum value of the exponent that corresponds to the strongest adverse pressure gradient but at a higher one whose dependence on Reynolds number is calculated in the paper. The results of the theory are in good agreement with experimental data. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Three-dimensional numerical simulation of nonisotropic thermal density flow in a strongly curved open channel

The results from a 3D nonisotropic algebraic stress/flux turbulence model are presented to investigate the structure of thermal density flow and the temperature distribution in a strongly curved open channel (180° bend). The numerically simulated results show that (i) several secondary flows take place at the bend cross-section 90° of the curved open channel, the feature which is not found for the isothermal flows and thermal density flow in a straight channel, and (ii) the thermocline in a curved channel is thicker than that in a straight channel due to the secondary flows-induced strong mixing process taking place in the former. Such features may be ascribed to the complex interaction of the buoyant force, the centrifugal force and the Reynolds stresses taking place only in curved channels. The simulated results are in good agreement with available experimental data, which indicates that the developed model can be applied for predicting the motion of the nonisotropic thermal density flow in the curved open channel.     

DNS of a Turbulent Channel Flow with Streamwise Rotation - Study of the Reverse Effect of the Cross Flow

Modeling of rotating turbulent flows is a major issue in engineering applications. In this work a turbulent channel flow rotating about the streamwise direction is presented. The theory is based on the investigations of [3] employing Lie group analysis. It was found that a cross flow in spanwise direction is induced. A series of direct numerical simulations (DNS) has been conducted for both different rotation rates and different Reynolds numbers to validate the cross flow. In addition some new interesting effects were observed. The averaged profile ū₃ of the cross flow is formed like a ‘S’ that means it exhibits a triple zero-crossing which denotes regions of reverse flow. Alaso a reverse effect is seen which means that for small rotation rates up to Ro=10 the spanwise mean velocity profiles increase and at rotation number Ro=14 this effect appears to reverse. Both effects were observed at two different Reynolds numbers Re = 180 and Re = 270. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Lower Branch Modes of the Compressible Boundary Layer Flow Due to a Rotating-Disk

In this article, a theoretical study is pursued to investigate the structure of the lower branch neutral stability modes of three-dimensional small disturbances imposed on the compressible boundary layer flow due to a rotating-disk. Special attention is focused on to the short-wavelength stationary/nonstationary compressible crossflow vortex modes at sufficiently high Reynolds numbers with reasonably small scaled frequencies. Following closely the asymptotic framework introduced in [ 1 ] for the incompressible stationary modes, it is demonstrated here that the compressible modes having sufficiently long time scale can also be described by an asymptotic expansion procedure based on the triple-deck approach. Making use of this rational asymptotic technique, which rigorously takes into account the nonparallel effects, the asymptotic structure of the nonstationary modes is shown to be adjusted by a balance between viscous and Coriolis forces, and resulted from the fact of vanishing shear stress at the disk surface, as in the incompressible Von Karman's flow. As a consequence of matching successive regions in the asymptotic procedure, it is found that the wavenumber and the orientation of the compressible lower branch modes are governed by an eigenrelation, which is akin to the one obtained previously in [ 1 ] for the incompressible stationary mode and in [ 2 ] for the compressible stationary modes. The nonparallel influences are toward destabilizing all the modes, though the wall insulation and heating are relatively stabilizing for the modes in the vicinity of the stationary mode, unlike the wall cooling. The asymptotic compressible data obtained at high Reynolds number limit compares fairly well with the numerical results generated directly solving the linearized compressible system with usual parallel flow approximation.     

A numerical study on statistical temporal scales in inertia particle dispersion

The statistical temporal scales involved in inertia particle dispersion are analyzed numerically. The numerical method of large eddy simulation, solving a filtered Navier-Stokes equation, is utilized to calculate fully developed turbulent channel flows with Reynolds numbers of 180 and 640, and the particle Lagrangian trajectory method is employed to track inertia particles released into the flow fields. The Lagrangian and Eulerian temporal scales are obtained statistically for fluid tracer particles and three different inertia particles with Stokes numbers of 1, 10 and 100. The Eulerian temporal scales, decreasing with the velocity of advection from the wall to the channel central plane, are smaller than the Lagrangian ones. The Lagrangian temporal scales of inertia particles increase with the particle Stokes number. The Lagrangian temporal scales of the fluid phase ‘seen’ by inertia particles are separate from those of the fluid phase, where inertia particles travel in turbulent vortices, due to the particle inertia and particle trajectory crossing effects. The effects of the Reynolds number on the integral temporal scales are also discussed. The results are worthy of use in examining and developing engineering prediction models of particle dispersion.     

On the Generation of Mean Flows by the Interaction of Görtler Vortices and Tollmien-Schlichting Waves in Curved Channel Flows

There are many fluid flows where the onset of transition can be caused by different instability mechanisms which compete in the nonlinear regime. Here the interaction of a centrifugal instability mechanism with the viscous mechanism which causes Tollmien-Schlichting waves is discussed. The interaction between these modes can be strong enough to drive the mean state; here the interaction is investigated in the context of curved channel flows so as to avoid difficulties associated with boundary layer growth. Essentially it is found that the mean state adjusts itself so that any modes present are neutrally stable even at finite amplitude. In the first instance the mean state driven by a vortex of short wavelength in the absence of a Tollmien-Schlichting wave is considered. It is shown that for a given channel curvature and vortex wavelength there is an upper limit to the mass flow rate which the channel can support as the pressure gradient is increased. When Tollmien-Schlichting waves are present then the nonlinear differential equation to determine the mean state is modified. At sufficiently high Tollmien-Schlichting amplitudes it is found that the vortex flows are destroyed, but there is a range of amplitudes where a fully nonlinear mixed vortex-wave state exists and indeed drives a mean state having little similarity with the flow which occurs without the instability modes. The vortex and Tollmien-Schlichting wave structure in the nonlinear regime has viscous wall layers and internal shear layers; the thickness of the internal layers is found to be a function of the Tollmien-Schlichting wave amplitude.     

Effect of constriction height on flow separation in a two-dimensional channel

We studied numerically the effect of the constriction height on viscous flow separation past a two-dimensional channel with locally symmetric constrictions. A numerically stable scheme in primitive variables (velocity and pressure) for the solution of two-dimensional incompressible time-dependent Navier–Stokes equations is employed using finite-difference approximation in staggered grid. The wall shear stresses at different heights of the constriction are computed and presented graphically. It is noticed that the maximum stress and the length of the recirculating region associated with two shear layers of the constriction increase with the increase of the area reduction of the constriction. The critical Reynolds number for symmetry breaking bifurcation for the 50%, 60% and 70% area reduction are obtained numerically. The flow field separates after the symmetry breaking bifurcation and the symmetry of the flow depends on the Reynolds number and the height of the constriction.     

微通道周期流动电位势及电粘性效应

求解了双电层的Poisson-Boltzmann方程和流体运动的Navier-Stokes方程,得到在周期压差作用下,二维微通道的周期流动电位势,流动诱导电场和液体流动速度的解析解．量纲分析表明,流体电粘性力与以下3个参数有关:1) 电粘性数,它表示定常流动时,通道最大电粘性力与压力梯度的比;2) 形状函数,它表示电粘性力在通道横截面的分布形态; 3) 耦合系数,它表示电粘性力的振幅衰减特征和相位差．分析结果表明,微通道周期流动诱导电场、流动速度与频率Reynolds数有关．在频率Reynolds数小于1时,流动诱导电场随频率Reynolds数变化很慢．在频率Reynolds数大于1时,流动诱导电场随频率Reynolds数的增加快速衰减．在通道宽度与双电层厚度比值较小情况下,电粘性效应对周期流动速度和流动诱导电场有重要影响．     

On Flow Past Bodies in Aligned Magnetic Fields

In inviscid, incompressible flow past a body in a strong alignedmagnetic field at very low magnetic Reynolds number, the flowfar from the body may be rectilinear but non-uniform. The paperdiscusses the evolution of small perturbations towards or awayfrom such a rectilinear flow and concludes that such evolutioncan only occur gradually (i.e. over many diameters) ahead ofthe body. Downstream the decay must occur close to the body.Two dimensional and axisymmetric flows are considered.     

Numerical computation of pulsatile flow through a locally constricted channel

This paper deals with the numerical solution of a pulsatile laminar flow through a locally constricted channel. A finite difference technique has been employed to solve the governing equations. The effects of the flow parameters such as Reynolds number, flow pulsation in terms of Strouhal number, constriction height and length on the flow behaviour have been studied. It is found that the peak value of the wall shear stress has significantly changed with the variation of Reynolds numbers and constriction heights. It is also noted that the Strouhal number and constriction length have little effect on the peak value of the wall shear stress. The flow computation reveals that the peak value of the wall shear stress at maximum flow rate time in pulsatile flow situation is much larger than that due to steady flow. The constriction and the flow pulsation produce flow disturbances at the vicinity of the constriction of the channel in the downstream direction.     

Rotation-Symmetric Referencebodys For Energy Efficient Flow Around Bodies

Topology optimization is used to optimize problems of flow around bodies and problems of guided flow. Within the context of research, optimization criteria are developed to increase the energy efficiency of these problems [1] [2] [3] [4] [5]. In order to evaluate the new criteria in respect to the increasing of energy efficiency, reference bodies for different Reynolds numbers in combination with given design space limitations are needed. Therefore, an optimal body at Reynolds number against 0 was analytically determined by Bourot [6]. At higher Reynolds numbers, in the range of laminar and turbulent flows, no analytical solution is known. Accordingly, reference bodies are calculated by CFD calculations at three technical relevant Reynolds numbers (1.000, 32.000, 100.000) in combination with parameter optimization. The cross section of the bodies is described by a parameterized model. To get the optimal body, a parameter optimization based on a “brute force”; algorithm is used to optimize with regard to the friction loss and pressure loss in order to minimize the total loss (c_d-value). The result is an optimal parameter constellation, depending on the Reynolds number. Within the results, it is possible to develop the optimal geometries. The identified characteristics of the flow field around these bodies are used as base for new optimization criteria. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Construction of circle bifurcations of a two-dimensional spatially periodic flow

The study by Yudovich [V.I. Yudovich, Example of the generation of a secondary stationary or periodic flow when there is loss of stability of the laminar flow of a viscous incompressible fluid, J. Math. Mech. 29 (1965) 587-603] on spatially periodic flows forced by a single Fourier mode proved the existence of two-dimensional spectral spaces and each space gives rise to a bifurcating steady-state solution. The investigation discussed herein provides a structure of secondary steady-state flows. It is constructed explicitly by an expansion that when the Reynolds number increases across each of its critical values, a unique steady-state solution bifurcates from the basic flow along each normal vector of the two-dimensional spectral space. Thus, at a single Reynolds number supercritical value, the bifurcating steady-state solutions arising from the basic solution form a circle.