A LES study of the flow interference between tandem square cylinder pairs

The results of an investigation on the interference effects of the tandem square cylinders exposed to a uniform flow are presented in this paper. Time-dependent and three-dimensional flow simulations are carried out using large eddy simulation with a one-equation subgrid model. An incompressible three-dimensional finite volume code with a collocated grid arrangement is used for solving filtered Navier–Stokes equations. These equations are solved with an implicit fractional two-step method. Simulations are conducted with different Reynolds numbers between 10³ and 10⁵. The longitudinal spacing between the cylinders is selected 4D for the chosen Reynolds numbers, where D is the side of the cylinders. Also the effect of the spacing between cylinders, ranging from 1D to 12D, is studied for the selected Reynolds numbers. The instantaneous flow field is studied by analyzing the vortices, pressure, streamlines and Q-criterion to assist understanding of the various flow patterns, vortical structures and Kelvin–Helmholtz vortices in the separating shear layers. The hysteresis is observed in a certain range of the gap spacing, which this range depends on the selected Reynolds number. The global results are also computed and compared with available experimental results. The results indicate that there is a satisfactory agreement between the predictions and available experimental data considering the fine grid adopted.     

Computation of high Reynolds number flow around a circular cylinder with surface roughness

Incompressible high-Reynolds-number flows around a circular cylinder are analyzed by direct integration of the Navier-Stokes equations using finite-difference method. A generalized coordinate system is used so that a sufficient number of grid points are distributed in the boundary layer and the wake. A numerical scheme which suppresses non-linear instability for calculations of high-Reynolds-number flows is developed. The computation of an impulsively started flow at Re = 1200 is compared with corresponding experimental observations, and excellent agreements are obtained.A series of computations are carried out on the flow around a circular cylinder with surface roughness. The height of the roughness in these computations is 0.5% of the diameter. The range of Reynolds numbers is from 10³ to 10⁵; no turbulence model is employed. Sharp reduction of drag coefficient is observed near Re = 2 × 10⁴, which indicates that the critical Reynolds number is captured in the present computation.     

Time-averaged topological flow patterns and their influence on vortex shedding of a square cylinder in crossflow at incidence

Flow characteristics around the square cylinder and their influence on the wake properties are studied. Time-averaged flow patterns on the surfaces of square cylinder in a cross-stream at incidence are experimentally probed by surface-oil flow technique and analyzed by flow topology for Reynolds numbers between 3.9×10⁴ and 9.4×10⁴ as the incidence angle changes from 0° to 45°. Vortex shedding characteristics are measured by a single-wire hot-wire anemometer for Reynolds numbers between 5×10³ and 1.2×10⁵. The effects of topological flow patterns on the wake properties then are revealed and discussed. Flows around the square cylinder are identified as three categories: the subcritical, supercritical, and wedge flows according to the prominently different features of the topological flow patterns. The Strouhal number of vortex shedding, turbulence in the wake, and wake width present drastically different behaviors in different characteristic flow regimes. A critical incidence angle of 15° separates the subcritical and supercritical regimes. At the critical incidence angle the wake width and shear-layer turbulence present minimum values. The minimum wake width appearing at the critical incidence angle, which leads to the maximum Strouhal number, is due to the reattachment of one of the separated boundary layer to the lateral face of the square cylinder. If the Strouhal numbers are calculated based on the wake width instead of the cross-stream projection width of cylinder, the data in the subcritical and supercritical regimes are well correlated into two groups, which would approach constants at high Reynolds numbers.     

Two‐dimensional prediction of time dependent,turbulent flow around a square cylinder confined in a channel

This paper presents two‐dimensional and unsteady RANS computations of time dependent, periodic, turbulent flow around a square block. Two turbulence models are used: the Launder–Sharma low‐Reynolds number k–ε model and a non‐linear extension sensitive to the anisotropy of turbulence. The Reynolds number based on the free stream velocity and obstacle side is Re=2.2×10⁴. The present numerical results have been obtained using a finite volume code that solves the governing equations in a vertical plane, located at the lateral mid‐point of the channel. The pressure field is obtained with the SIMPLE algorithm. A bounded version of the third‐order QUICK scheme is used for the convective terms. Comparisons of the numerical results with the experimental data indicate that a preliminary steady solution of the governing equations using the linear k–ε does not lead to correct flow field predictions in the wake region downstream of the square cylinder. Consequently, the time derivatives of dependent variables are included in the transport equations and are discretized using the second‐order Crank–Nicolson scheme. The unsteady computations using the linear and non‐linear k–ε models significantly improve the velocity field predictions. However, the linear k–ε shows a number of predictive deficiencies, even in unsteady flow computations, especially in the prediction of the turbulence field. The introduction of a non‐linear k–ε model brings the two‐dimensional unsteady predictions of the time‐averaged velocity and turbulence fields and also the predicted values of the global parameters such as the Strouhal number and the drag coefficient to close agreement with the data. Copyright © 2009 John Wiley & Sons, Ltd.     

Experimental investigation of the separation of flow around a sphere

The results of an experimental investigation of the flow around a sphere over a broad range of Mach numbers M=0.3–3 and Reynolds numbers Re=3·10⁴–3·10⁷ are presented. The experiments were carried out on a ballistic test stand and in a wind tunnel. Flow patterns and pressure distributions were obtained. In particular, the effect of the Mach and Reynolds numbers on the position of the separation point and the edge shock was studied; the pressure distribution on the sphere was measured; and a nonmonotonic displacement of the flow separation point upon passage through the speed of sound was established.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 152–156, January–February, 1991.     

Anisotropic Organised Eddy Simulation for the prediction of non-equilibrium turbulent flows around bodies

The unsteady turbulent flow around bodies at high Reynolds number is predicted by an anisotropic eddy-viscosity model in the context of the Organised Eddy Simulation (OES). A tensorial eddy-viscosity concept is developed to reinforce turbulent stress anisotropy, that is a crucial characteristic of non-equilibrium turbulence in the near-region. The theoretical aspects of the modelling are investigated by means of a phase-averaged PIV in the flow around a circular cylinder at Reynolds number 1.4×10⁵. A pronounced stress–strain misalignment is quantified in the near-wake region of the detached flow, that is well captured by a tensorial eddy-viscosity concept. This is achieved by modelling the turbulence stress anisotropy tensor by its projection onto the principal matrices of the strain-rate tensor. Additional transport equations for the projection coefficients are derived from a second-order moment closure scheme. The modification of the turbulence length scale yielded by OES is used in the Detached Eddy Simulation hybrid approach. The detached turbulent flows around a NACA0012 airfoil (2-D) and a circular cylinder (3-D) are studied at Reynolds numbers 10⁵ and 1.4×10⁵, respectively. The results compared to experimental ones emphasise the predictive capabilities of the OES approach concerning the flow physics capture for turbulent unsteady flows around bodies at high Reynolds numbers.     

Finite element analysis of vortex shedding using equal order interpolations

An operator splitting and element‐by‐element conjugated gradient solver, and equal order interpolations are applied for solving time dependent Navier–Stokes (NS) equations to simulate flow induced vortex shedding in the present study. In addition, the convection term is corrected by balanced tensor diffusivity, which can stabilize the numerical simulation and overcome the numerical oscillations. The evolution of the interested flowing properties with time is analyzed by using spectral analysis. The developed code has been validated by the application of two examples: a driven cavity flow and a flow induced vortex vibration. Results from the first example for Reynolds number Re=10³ and Re=10⁴ are compared with other numerical simulations. Results from the second example, uniform flow past a square rod over a wide range of high Reynolds numbers from Re=10³～10⁵, are compared with experimental data and other numerical studies. Copyright © 2002 John Wiley & Sons, Ltd.     

A least-squares finite element method for time-dependent incompressible flows with thermal convection

The time-dependent Navier–Stokes equations and the energy balance equation for an incompressible, constant property fluid in the Boussinesq approximation are solved by a least-squares finite element method based on a velocity–pressure–vorticity–temperature–heat-flux ( u –P–ω–T– q ) formulation discretized by backward finite differencing in time. The discretization scheme leads to the minimization of the residual in the l²-norm for each time step. Isoparametric bilinear quadrilateral elements and reduced integration are employed. Three examples, thermally driven cavity flow at Rayleigh numbers up to 10⁶, lid-driven cavity flow at Reynolds numbers up to 10⁴ and flow over a square obstacle at Reynolds number 200, are presented to validate the method.     

圆角化对受迫振动方柱绕流特性影响机理分析

为分析圆角化对低雷诺数下受迫振动方柱绕流特性的影响机理, 对Ansys Fluent软件进行二次开发, 即通过用户自定义函数中的DEFINE_ CG_MOTION宏对柱体周期性受迫振动的函数进行编程, 并对流场计算域进行区域划分以便利用动网格技术中动态层法实现柱体受迫振动, 从而实现对受迫振动柱体绕流流场的流固耦合模拟.在雷诺数Re = 200时, 考虑方柱截面不同圆角的影响, 对均匀流作用下5种圆角化r/D = 1/2, 1/4, 1/5, 1/8和0受迫振动方柱的绕流进行数值模拟, 分析了这5种参数下受迫振动方柱的升阻力系数、尾流涡量和锁定区间的变化规律, 澄清了圆角化对受迫振动方柱稳定性的影响机理.研究表明: 与尖角方柱相比, 圆角化方柱升阻力系数有了明显的减小, 且升力、阻力系数随圆角增大而减小; 低振幅比下圆角方柱的涡旋脱落模式均为2S模态, 涡旋尾迹变窄; 锁定区间范围基本关于F = 1对称, 锁定区间的变化趋势与圆柱类似.      

Flow past a cylinder: shear layer instability and drag crisis

Flow past a circular cylinder for Re=100 to 10⁷ is studied numerically by solving the unsteady incompressible two‐dimensional Navier–Stokes equations via a stabilized finite element formulation. It is well known that beyond Re ～ 200 the flow develops significant three‐dimensional features. Therefore, two‐dimensional computations are expected to fall well short of predicting the flow accurately at high Re. It is fairly well accepted that the shear layer instability is primarily a two‐dimensional phenomenon. The frequency of the shear layer vortices, from the present computations, agree quite well with the Re^0.67 variation observed by other researchers from experimental measurements. The main objective of this paper is to investigate a possible relationship between the drag crisis (sudden loss of drag at Re ～ 2 × 10⁵) and the instability of the separated shear layer. As Re is increased the transition point of shear layer, beyond which it is unstable, moves upstream. At the critical Reynolds number the transition point is located very close to the point of flow separation. As a result, the shear layer eddies cause mixing of the flow in the boundary layer. This energizes the boundary layer and leads to its reattachment. The delay in flow separation is associated with narrowing of wake, increase in Reynolds shear stress near the shoulder of the cylinder and a significant reduction in the drag and base suction coefficients. The spatial and temporal power spectra for the kinetic energy of the Re=10⁶ flow are computed. As in two‐dimensional isotropic turbulence, E(k) varies as k^?5/3 for wavenumbers higher than energy injection scale and as k^?3 for lower wavenumbers. The present computations suggest that the shear layer vortices play a major role in the transition of boundary layer from laminar to turbulent state. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical simulation and hydrodynamic visualization of transient viscous flow around an oscillating aerofoil

Unsteady viscous flow around a large-amplitude and high-frequency oscillating aerofoil is examined in this paper by numerical simulation and experimental visualization. The numerical method is based on the combination of a fourth-order Hermitian finite difference scheme for the stream function equation and a classical second-order scheme to solve the vorticity transport equation. Experiments are carried out by a traditional visualization method using solid tracers suspended in water. The comparison between numerical and experimental results is found to be satisfactory. Time evolutions of the flow structure are presented for Reynolds numbers of 3 × 10³ and 10⁴. The influence of the amplitude and frequency of the oscillating motion on the dynamic stall is analysed.     

Numerical Simulation of Aerodynamic Characteristics of Two Rectangular Cylinders in Side-by-side Arrangement

We present numerical results of flows around two rectangular cylinders in side-by-side arrangement by three-dimensional computations. The two rectangular cylinders are arranged with various distances between the cylinders. The three-dimensional flow structures around two rectangular cylinders denote significant features depending on the distance between the cylinders. In our computations, the Reynolds number Re is set as 10,000, and both the aspect ratio of section of two rectangular cylinders and the distance between the cylinders are considered as parameters to calculate the flow around two rectangular cylinders. The obtained numerical results are also compared with experimental data.     

Drag and Separation Flow Past Solid Sphere with Porous Shell at Moderate Reynolds Number

Axisymmetric viscous, two-dimensional steady and incompressible fluid flow past a solid sphere with porous shell at moderate Reynolds numbers is investigated numerically. There are two dimensionless parameters that govern the flow in this study: the Reynolds number based on the free stream fluid velocity and the diameter of the solid core, and the ratio of the porous shell thickness to the square root of its permeability. The flow in the free fluid region outside the shell is governed by the Navier–Stokes equation. The flow within the porous annulus region of the shell is governed by a Darcy model. Using a commercially available computational fluid dynamics (CFD) package, drag coefficient and separation angle have been computed for flow past a solid sphere with a porous shell for Reynolds numbers of 50, 100, and 200, and for porous parameter in the range of 0.025–2.5. In all simulation cases, the ratio of b/a was fixed at 1.5; i.e., the ratio of outer shell radius to the inner core radius. A parametric equation relating the drag coefficient and separation point with the Reynolds number and porosity parameter were obtained by multiple linear regression. In the limit of very high permeability, the computed drag coefficient as well as the separation angle approaches that for a solid sphere of radius a, as expected. In the limit of very low permeability, the computed total drag coefficient approaches that for a solid sphere of radius b, as expected. The simulation results are presented in terms of viscous drag coefficient, separation angles and total drag coefficient. It was found that the total drag coefficient around the solid sphere as well as the separation angle are strongly governed by the porous shell permeability as well as the Reynolds number. The separation point shifts toward the rear stagnation point as the shell permeability is increased. Separation angle and drag coefficient for the special case of a solid sphere of radius r = a was found to be in good agreement with previous experimental results and with the standard drag curve.     

Experimental investigations of local heat fluxes on a sphere and on the spherical blunt end of an axisymmetric body

This article presents the results of experimental investigations of local heat transfer on a sphere and on the spherical blunt end of axisymmetric bodies in a flow of low-density gas. The data were obtained in the region of Mach numbers from 3.3 to 8, at Reynolds numbers from 5 to 2.6·10³.     

Fluid-dynamic interaction between two spheres

Experiment of fluid-dynamic interaction between two spheres was conducted to obtain basic information concerning the two-phase flow, especially in dense phase. Two or three spheres were set up in a water tunnel in the longitudinal or transverse direction with Reynolds numbers less than 10³. The flow behind the sphere was visualized by the use of condense milk and change in vortex structure due to the interaction was observed in detail. Additionally, drag force on the sphere was measured by a pendulum method which was developed to detect small drag, and the range of distance in which the drag is affected by the interaction was shown.     

Flow in the base region of flat bodies

A considerable number of theoretical and experimental papers has been devoted in recent years to studying the flow details in the base region which were not included in the Chapman-Korst scheme. The authors obtained new data supplementing the known facts, or altering the conception about the substance of certain phenomena (the origination of the edge shock, for example). Results are presented below of experimental investigations of the flow configuration in the base region of flat bodies, and the fundamental parameters governing the difference between the flow behind pointed and blunt bodies are examined. The experiments were conducted in the wind tunnel of the Moscow State University Institute of Mechanics in the M=0.3–3.8 Mach number range. The Reynolds numbers, referred to 0.1 m and computed by means of the free-stream parameters, varied between 2·10⁶ and 3·10⁶. The flow behind a wedge (with a 45° vertex angle and 0.025 m altitude) mounted on a plate 0.2 m wide and 0.4 m long was investigated. The flow was studied in the near wake on wedges withapex half-angles =15, 20, and 30°, width 0.2, altitude of the rear section 0.05 m, and a cylinder of the same width and 0.05 m diameter. The models were mounted in the center of the working section by using slender lateral pylons fastened to the perforated walls of the tunnel. A strip of emery paper was glued to turbulize the boundary layer on the models. The pressure in the base domain was measured by using total and static pressure detectors (1.2 mm outer diameter) mounted on a traversing gear.Institute of Mechanics, Moscow State University. Moscow. Translated from Izvestiya Akademii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 6, pp. 61–70, November–December, 1972.     

Explicit Friction Factor Accuracy and Computational Efficiency for Turbulent Flow in Pipes

The implicit Colebrook–White equation is the accepted method for accurately estimating the friction factor for turbulent flow in pipes. This study reviews 28 explicit equations for approximating the friction factor to integrate both the accuracy to the implicit Colebrook–White equation and the relative computational efficiency of the explicit equations. A range of 901 Reynolds numbers were selected for the review between Re?≥?4 ×10³ and? ≤?4 × 10⁸ and 20 relative pipe roughness values were selected between $\varepsilon \mathord{\left/ {\vphantom {\varepsilon D}} \right.}D\ge 10^{-6}\le 10^{-1}$ , thus producing a matrix of 18,020 points for each explicit equation, covering all the values to be encountered in pipeline hydraulic analysis for turbulent flow. The accuracy of the estimation of the friction factor using the explicit equations to the value obtained using the implicit Colebrook–White equation were calculated and reported as absolute, relative percentage and mean square errors. To determine the relative computational efficiency, 300 million friction factor calculations were performed using randomly generated values for the Reynolds number and the relative pipe roughness values between the limits specified for each of the explicit equations and compared to the time taken by the Colebrook–White equation. Finally, 2D and 3D contour models were generated showing both the range and magnitude of the relative percentage accuracy across the complete range of 18,020 points for each explicit equation.     

The determination of pipe contraction pressure loss coefficients for incompressible turbulent flow

The accurate prediction of pipe contraction presure loss is important in the design of pipe systems, such as heat exchangers, particularly when close control of the flow distribution in a network of pipes is required. The prediction of the contraction pressure loss depends heavily on experimental data. Large discrepancies in these predictions are evident in the literature. New experimental results giving pressure loss coefficients for range of Reynolds numbers of 4×10⁴ to 2×10⁵ and area ratios 0.13 to 0.7 are presented and compared with those of other workers and with predictions from a method that allows for velocity profile variation through the contraction. The results show a Reynolds number dependence. The effects of small-bore pipe inlet geometry on the loss coefficients are also examined.     

Mixed convection cooling of a heated circular cylinder by laminar upward-directed slot jet impingement

An experimental and numerical study has been carried out to investigate the heat transfer characteristics of a horizontal circular cylinder exposed to a slot jet impingement of air. A square-edged nozzle is mounted parallel with the cylinder axis and jet flow impinges on the bottom of the cylinder. The study is focused on low Reynolds numbers ranging from 120 to 1,210, Grashof numbers up to Gr = 10Re ² and slot-to-cylinder spacing from 2 to 8 of the slot width. The flow field is greatly influenced by the slot exit velocity and the buoyancy force due to density change. A Mach–Zehnder Interferometer is used for measurement of local Nusselt number around the cylinder at 10° interval. It is observed that the average Nusselt number decreases with increasing the jet spacing and increases with rising the Reynolds number. A finite volume method utilizing a curvilinear coordinate transformation is used for numerical modeling. The numerical results show good agreement with the experimental results. The flow and thermal field are seen to be stable and symmetric around the cylinder over the range of parameters studied.     

Time‐dependent turbulent cavitating flow computations with interfacial transport and filter‐based models

Turbulent cavitating flow computations need to address both cavitation and turbulence modelling issues. A recently developed interfacial dynamics‐based cavitation model (IDCM) incorporates the interfacial transport into the computational modelling of cavitation dynamics. For time‐dependent flows, it is known that the engineering turbulence closure such as the original k–ε model often over‐predicts the eddy viscosity values reducing the unsteadiness. A recently proposed filter‐based modification has shown that it can effectively modulate the eddy viscosity, rendering better simulation capabilities for time‐dependent flow computations in term of the unsteady characteristics. In the present study, the IDCM along with the filter‐based k–ε turbulence model is adopted to simulate 2‐D cavitating flows over the Clark‐Y airfoil. The chord Reynolds number is Re=7.0 × 10⁵. Two angles‐of‐attack of 5 and 8° associated with several cavitation numbers covering different flow regimes are conducted. The simulation results are assessed with the experimental data including lift, drag and velocity profiles. The interplay between cavitation and turbulence models reveals substantial differences in time‐dependent flow results even though the time‐averaged characteristics are similar. Copyright © 2005 John Wiley & Sons, Ltd.