细长体背风面横向非对称分离的研究

用数值模拟方法研究了超声速情况下，无限长细长体背风面的涡结构。数值模拟的出发方程和计算格式分别为全Ｎ－Ｓ方程和二阶空间精度的ＴＶＤ格式。数值结果给出了圆锥、半球柱体和椭圆锥在不同攻角下的流场结果。结果表明圆锥在攻角α＝１５°，２０°和２５°时背风面呈现明显的稳定非对称横向分离，而半球柱体和椭圆锥在３２．５°和２５°时背风面均未出现非对称的横向分离结构。     

Supersonic viscous perfect-gas flow past a slender sharp elliptic cone

The supersonic M_∞ = 5 flow past slender elliptic cones with the semi-vertex angle in the plane of the major semi-axis ? _c = 4° and an isothermal surface is investigated under the assumption of the flow symmetry. Calculations on the basis of the time-dependent three-dimensional Navier-Stokes and Reynolds equations are carried out on the Reynolds number and angle of attack ranges 10⁴ ≤ Re ≤ 108 and 0 ≤ α ≤ 15° for cones with ellipticity coefficients 1/32 ≤ δ= b/a ≤ 1. The effect of the relevant parameters of the problem on the flowfield structure and the aerodynamic characteristics of the cones is demonstrated.     

Heat transfer and flow behavior around four staggered elliptic cylinders in cross flow

Experimental and numerical studies were carried out to investigate forced convection heat transfer and flow features around the downstream elliptic cylinder in four staggered cylinders in cross flow. The elliptic cylinders examined had an axis ratio (b/c) of 1:2, and they were arranged with zero angle of attack to the upstream flow. The present heat transfer measurements were obtained by heating only the downstream elliptic cylinder (test cylinder) under the condition of constant heat flux. The testing fluid was air and the Reynolds number based on the major axis length (c) was ranged from 4,000 to 45,570. The tested longitudinal spacing ratio (S_x/c) and the transversal spacing ratio (S_y/b) were in the ranges of 1.5 ≤ S_x/c ≤ 4.0 and 1.5 ≤ S_y/b ≤ 4.0, respectively. The air flow pattern and temperature fields around the four staggered elliptic cylinders were predicted by using CFD software package. Also, a flow visualization study was made to show the flow features around the elliptic cylinders. It was observed that Nu_m of the downstream elliptic cylinder in four staggered cylinders was higher than that of three in-line cylinders for all tested spacing ratios and Reynolds numbers except for Re = 4,000. It was clear that, at lower Reynolds number values (Re < 14,100), the average Nusselt number of the downstream elliptic cylinder in three staggered arrangement was higher than that of the downstream cylinder in four staggered arrangement for all tested spacing ratios. On the other hand, at Re > 14,100, the tested elliptic cylinder in four staggered arrangement had the higher values of the average Nusselt number. Moreover, in four staggered arrangement, the maximum average Nusselt number enhancement ratio (average Nusselt number of the tested downstream cylinder/average Nusselt number of a single elliptic cylinder) was found to be about 2.0, and was obtained for spacing ratios of S_x/c = 2.5, S_y/b = 2.5 and at Re = 32,000. Finally, the average Nusselt number of the tested cylinder in four staggered arrangement was correlated in terms of Reynolds number and cylinder spacing ratios.     

Numerical simulation of high‐Reynolds number flow around circular cylinders by a three‐step FEM–BEM model

An innovative computational model, developed to simulate high‐Reynolds number flow past circular cylinders in two‐dimensional incompressible viscous flows in external flow fields is described in this paper. The model, based on transient Navier–Stokes equations, can solve the infinite boundary value problems by extracting the boundary effects on a specified finite computational domain, using the projection method. The pressure is assumed to be zero at infinite boundary and the external flow field is simulated using a direct boundary element method (BEM) by solving a pressure Poisson equation. A three‐step finite element method (FEM) is used to solve the momentum equations of the flow. The present model is applied to simulate high‐Reynolds number flow past a single circular cylinder and flow past two cylinders in which one acts as a control cylinder. The simulation results are compared with experimental data and other numerical models and are found to be feasible and satisfactory. Copyright © 2001 John Wiley & Sons, Ltd.     

Transient numerical simulation of buoyancy driven flow adjacent to an elliptic tube

In this paper the problem of laminar, transient, two-dimensional free convective heat transfer from the surface of a horizontal elliptic tube is considered. The tube, whose surface is suddenly subjected to uniform heat flux, is placed in a quiescent Boussinesq Newtonian fluid with its major axis horizontal. The details of both flow and thermal fields are obtained by solving the full governing Navier–Stokes and energy equations. These equations, expressed in terms of stream function, vorticity and temperature, are numerically solved using an implicit spectral finite difference procedure. The parameters involved are the modified Rayleigh number, Prandtl number and axis-ratio. The investigation covers a Rayleigh number range up to 10⁷. The minor–major axis ratio of elliptic cylinder ranges between 0.05 and 0.998 and Prandtl number ranges between 0.1 and 10. The effects of these parameters on the surface temperature distribution and heat transfer coefficients are determined and the different aspects of the results are discussed for some selected cases.     

Three-dimensional simulations of flow past two circular cylinders in side-by-side arrangements at right and oblique attacks

The flow past two identical circular cylinders in side-by-side arrangements at right and oblique attack angles is numerically investigated by solving the three-dimensional Navier–Stokes equations using the Petrov–Galerkin finite element method. The study is focused on the effect of flow attack angle and gap ratio between the two cylinders on the vortex shedding flow and the hydrodynamic forces of the cylinders. For an oblique flow attack angle, the Reynolds number based on the velocity component perpendicular to the cylinder span is defined as the normal Reynolds number Re_N and that based on the total velocity is defined as the total Reynolds number Re_T. Simulations are conducted for two Reynolds numbers of Re_N=500 and Re_T=500, two flow attack angles of α=0° and 45° and four gap ratios of G/D=0.5, 1, 3 and 5. The biased gap flow for G/D=0.5 and 1 and the flip-flopping bistable gap flow for G/D=1 are observed for both α=0° and 45°. For a constant normal Reynolds number of Re_N=500, the mean drag and lift coefficients at α=0° are very close to those at α=45°. The difference between the root mean square (RMS) lift coefficient at α=0° and that at α=45° is about 20% for large gap ratios of 3 and 5. From small gap ratios of 0.5 and 1, the RMS lift coefficients at α=0° and 45° are similar to each other. The present simulations show that the agreement in the force coefficients between the 0° and 45° flow attack angles for a constant normal Reynolds number is better than that for a constant total Reynolds number. This indicates that the normal Reynolds number should be used in the implementation of the independence principle (i.e., the independence of the force coefficients on the flow attack angle). The effect of Reynolds number on the bistable gap flow is investigated by simulating the flow for Re_N=100–600, α=0° and 45° and G/D=1. Flow for G/D=1 is found to be two-dimensional at Re_N=100 and weak three-dimensional at Re_N=200. While well defined biased flow can be identified for Re_N=300–600, the gap flow for Re_N=100 and 200 changes its biased direction too frequently to allow stable biased flow to develop.     

Fluid-dynamic forces and wake frequencies on a tilted rectangular cylinder near a solid wall

We investigate the impact of different boundary conditions on the flow field developing around a tilted rectangular cylinder with two different values of the aspect ratio (l/s=3 and 4). We are mainly interested in analyzing the changes in force coefficients and in the vortex shedding Strouhal number when the cylinder is placed at various distances from a bottom wall and different values of attack angle. The angle of attack ranges between −30° and +30° and the cylinder elevation above the bottom wall is varied between almost zero and 5 times the thickness of the cylinder. A large body of experimental results is related to the small elevation conditions at different attack angles, where the presence of the wall has a non-negligible effect on the behavior of the force coefficients and Strouhal number of the vortex shedding.     

棱柱绕流流动的数值模拟

利用数值方法对长宽比为1／3, 1和3的棱柱绕流在雷诺数为100的非稳态流动特性进行了分析和研究。采用有限体积法对棱柱绕流的二维流动N-S方程进行离散求解,分析和研究了非稳态的棱柱绕流流场,升力系数,阻力系数和涡动特性,数值模拟的结果与相关文献的数据比较吻合。通过上述研究能够为了解棱柱绕流的非稳态流动特性提供有力的帮助。而对棱柱三维流动的模拟分析和对雷诺数的变化对棱柱流动特性的影响进行研究,将为掌握棱柱绕流的工程特性打下基础。     

Topology of vortex creation in the cylinder wake

We analyze the topology of the two-dimensional flow around a circular cylinder at moderate Reynolds numbers in the regime where the vortex wake is created. A normal form for the stream function close to the cylinder is presented and used to predict the streamline pattern both in the steady and the periodic regime, where two different vortex shedding scenarios are identified. The theoretical predictions are verified numerically. For the vorticity, a very different topology occurs with infinite nested sequences of iso-curves moving downstream. General equations of motion for critical points are derived.     

Heat convection from a cylinder performing steady rotation or rotary oscillation – Part I: Steady rotation

In this paper, the problem of laminar, two dimensional heat convection from a circular cylinder performing steady rotation is investigated. The cylinder is␣placed with its axis horizontal in a quiescent fluid of infinite extent. Because of viscous dissipation, the flow process is confined to the region adjacent to the cylinder and is mainly driven by shear and buoyancy forces. The study is based on the solution of the full conservation equations of mass, momentum and energy for Rayleigh numbers up to 10⁴ and Reynolds numbers (based on surface velocity) up to 400 while Prandtl number ranges between 0.7 and 7.0. For the range of parameters considered, the study revealed that the rate of heat transfer increases with the increase of Rayleigh number and decreases with the increase of speed of rotation. The increase of Prandtl number resulted in an appreciable increase in the average Nusselt number only at low Reynolds numbers. The effect of Prandtl number at high Reynolds number is negligibly small. The resulting flow field in all cases is steady with no vortex shedding. The streamlines and isotherms are plotted for a number of cases to show the details of the velocity and thermal fields. Received on 15 December 1997     

FLOW CHARACTERISTICS AROUND AN INCLINED ELLIPTIC CYLINDER IN A TURBULENT BOUNDARY LAYER

The flow characteristics around an inclined elliptic cylinder located near a flat plate were investigated experimentally. The axis ratio of the elliptic cylinder was AR=2. The pressure distributions along the surface of the cylinder and the flat plate were measured by varying the angle of attack of the elliptic cylinder. The velocity profiles behind the cylinder were measured using hot-wire anemometry. When the angle of attack varies, the peak pressure location on the windward cylinder surface moves towards the rear edge of the cylinder, while that on the leeward surface moves towards the front edge of the cylinder. The vortex-shedding frequency also gradually decreases, defining a critical angle of attack for each gap ratio. The location of the minimum pressure on the flat plate surface moves downstream for positive angles of attack, while it moves upstream for negative angles of attack. Negative angles of attack cause a greater disturbance in the boundary layer near the wall compared to positive angles of attack. This shows that the separated wall shear layer from the boundary layer and the lower shear layer of the cylinder wake are strongly merged compared to other cases.     

NUMERICAL SOLUTION OF THE SLOW TRANSLATION OF A SPHERE MOVING ALONG THE AXIS OF A ROTATING VISCOUS FLUID

The motion of a sphere along the axis of rotation of an incompressible viscous fluid that is rotating as a solid mass is investigated by means of numerical methods for small values of Reynolds numbers and moderate values of Taylor numbers. The Navier-Stokes equations governing the steady, axisymmetric, viscous flow can be written as three coupled, nonlinear, elliptic partial differential equations for the stream function, vorticity and rotational velocity component. Finite difference method is used for solving the governing equations. Second order derivatives are approximated by central differences and nonlinear terms are approximated by upwind differences. Results are presented mostly in the form of graphs of the streamlines and vorticity lines. When 1/ Ro > 2.2, separation occurs and reverse flow is obtained.     

The distortion of a non-aligned magnetic field by the flow of an inviscid conducting fluid past a circular cylinder

The distortion of a magnetic field by the flow of a conducting fluid past a cylinder of the same permeability is found for small and infinite values of the magnetic Reynolds number. For small values good agreement is obtained with the results of Seebass and Tamada when the flow is aligned with the field at large distances from the body.For infinite magnetic Reynolds number, all the lines of force are dragged into the cylinder and upstream and downstream wake regions are present on the axis of the flow.     

Effect of Large Scale 3-D Structures on the Flow Around a Heated Cylinder at Low Reynolds Number

This work presents results of flow around a heated circular cylinder in mixed convection regime and demonstrates that Prandtl number and angle of attack of the incoming flow have a large influence on the characterisation of the flow transition from 2-D to 3-D. Previous studies show that heat transfer can enhance the formation of large 3-D structures in the wake of the cylinder for Reynolds numbers between 75 and 127 and a Richardson number larger than 0.35. This transitional mode is generally identified as “mode E”. In this work, we compare the results for water-based flow (large Prandtl number) with the ones for air-based flows (low Prandtl number). The comparison is carried out at two Reynolds numbers (100 and 150) and at a fixed Richardson number of 1. It shows that at the low Reynolds number of 100 the low Prandtl number flow does not enter into transition. This is caused by the impairment of the baroclinic vorticity production provoked by the spanwise temperature gradient. At low Prandtl number temperature gradients are less steep. For an air-based flow at Reynolds number 150, several Richardson numbers have been simulated. In this situation, the flow enters into transition and exhibits the characteristics of “mode E”, with the development of Λ-shaped structures in the near wake and mushroom-like structures in the far wake. It is also observed that the transition is delayed at Richardson number of 0.5. Simulations are also carried to investigate the effect of the angle of attack on the incoming flow on the development of large coherent structures. When the angle of attack is positive, the development of the wake tends to return to a more bi-dimensional configuration, where large scale coherent structures are impaired. In contrast, when the angle of attack is negative, large scale tri-dimensional structures dominate the flow in the wake, but with a very chaotic behaviour and the regular pattern of zero angle of attack is destroyed. The different behaviour of the flow with the variation of the angle of attack is also related to the baroclinic vorticity production, where new terms appear in the equations, leading to a positive effect of the vorticity production in case of a negative angle of attack and the opposite for a positive angle of attack.     

Effects of an axial flow component on the Honji instability

This paper is concerned with a numerical study of the three-dimensional Honji instability that can arise in an oscillatory flow impinging on a circular cylinder. It is well known that when the fluid motion far from the cylinder is perpendicular to its axis then the flow is liable to a three-dimensional breakdown via this instability which initially appears as an axially periodic mushroom-like structure attached to the surface of the cylinder. Here the focus is on examining the Honji instability under an oblique inflow. The obliqueness of the free stream is represented by an angle of attack through introducing an axial flow component. It is found that the Honji mode is suppressed by increasing the axial flow component, and when this component is sufficiently large the instability mechanism is no longer operative so that all that remains is a featureless two-dimensional columnar flow. At smaller values of angle of attack, though the Honji structure remains, it is deformed by the axial flow component. The developed two-layer near-cylinder vortical structures can be related to the energy and momentum transfer between the two layers.     

On the effect of flow direction on mixed convection from a horizontal cylinder

The influence of free stream direction on mixed (natural and forced) convective heat transfer from a circular cylinder is investigated. The cylinder, which has an isothermal surface, is placed with its axis horizontal and normal to the oncoming flow. The free stream direction varies between the vertically upward (parallel flow) and the vertically downward (contraflow) directions. The investigation is based on the time integration of the unsteady, two-dimensional equations of motion and energy until reaching steady conditions. The study is limited to Reynolds numbers up to Re = 40 and Grashoff numbers of Gr = Re². The results are compared with the available experimental data and the agreement is satisfactory.     

Boundary layer flow of Maxwell fluid due to torsional motion of cylinder: modeling and simulation

This paper investigates the boundary layer flow of the Maxwell fluid around a stretchable horizontal rotating cylinder under the influence of a transverse magnetic field. The constitutive flow equations for the current physical problem are modeled and analyzed for the first time in the literature. The torsional motion of the cylinder is considered with the constant azimuthal velocity E. The partial differential equations (PDEs) governing the torsional motion of the Maxwell fluid together with energy transport are simplified with the boundary layer concept. The current analysis is valid only for a certain range of the positive Reynolds numbers. However, for very large Reynolds numbers, the flow becomes turbulent. Thus, the governing similarity equations are simplified through suitable transformations for the analysis of the large Reynolds numbers. The numerical simulations for the flow, heat, and mass transport phenomena are carried out in view of the bvp4c scheme in MATLAB. The outcomes reveal that the velocity decays exponentially faster and reduces for higher values of the Reynolds numbers and the flow penetrates shallower into the free stream fluid. It is also noted that the phenomenon of stress relaxation, described by the Deborah number, causes to decline the flow fields and enhance the thermal and solutal energy transport during the fluid motion. The penetration depth decreases for the transport of heat and mass in the fluid with the higher Reynolds numbers. An excellent validation of the numerical results is assured through tabular data with the existing literature.     

Reynolds-number-effects in flow around a rectangular cylinder with aspect ratio 1:5

The paper reports on experiments carried out over a wide range of Reynolds numbers in a high pressure wind tunnel. The model was a sharp-edged rectangular cylinder with aspect ratio height/width 1:5 (width/span ratio 1:10.8), which was investigated in both basic orientations, lengthwise (4×10³<Re<4×10⁵) and perpendicular to the flow (2.7×10⁴<Re<6.4×10⁵). The Reynolds number is based on the height of the model normal to the flow. Steady and unsteady forces were measured with a piezoelectric balance. Thus along with steady (i.e. time averaged values) including the base pressure coefficient, also power spectra and probability density functions were measured yielding for example Strouhal numbers, higher statistical moments, etc. A response diagram for the vortex resonance phenomenon was taken for the natural bending motion of the slender model. If lift coefficient for constant angle of attack is plotted against Reynolds number, a significant Reynolds number effect is seen. For α=4°, the curve shows an inflection point and the lift varies between 0.3 and 0.6. For α=6° and 2° there are similar variations shifted to lower and higher values of Re, respectively. Probably the shapes of separation bubbles that depend on the Reynolds number are responsible for these effects. No Reynolds number effects were observed when the long side was normal to the flow, an orientation where reattachment at the side walls is not possible. Comparing both basic cases (α=0° and 90°), the interpretation of the probability distributions of lift force leads to the conclusion that the possibility of reattachment (α=0°) seems to enhance the degree of order in the vortex shedding process.     

A high-order compact scheme for solving the 2D steady incompressible Navier-Stokes equations in general curvilinear coordinates

In this paper, a high-order compact finite difference algorithm is established for the stream function-velocity formulation of the two-dimensional steady incompressible Navier-Stokes equations in general curvilinear coordinates. Different from the previous work, not only the stream function and its first-order partial derivatives but also the second-order mixed partial derivative is treated as unknown variable in this work. Numerical examples, including a test problem with an analytical solution, three types of lid-driven cavity flow problems with unusual shapes and steady flow past a circular cylinder as well as an elliptic cylinder with angle of attack, are solved numerically by the newly proposed scheme. For two types of the lid-driven trapezoidal cavity flow, we provide the detailed data using the fine grid sizes, which can be considered the benchmark solutions. The results obtained prove that the present numerical method has the ability to solve the incompressible flow for complex geometry in engineering applications, especially by using a nonorthogonal coordinate transformation, with high accuracy.     

Numerical simulation with a TVD–FVM method for circular cylinder wake control by a fairing

A water drop-shaped fairing is applied to control the wake behind a circular cylinder and to suppress the formation of Karman vortex street in this paper. The results are evaluated using high resolution CFD technique. A finite-volume total variation diminishing (TVD) approach based upon the recently proposed elemental velocity vector transformation (EVVT) method, which aims at solving the incompressible turbulent flow for irregular boundary conditions with renormalization group (RNG) turbulence model, is used to simulate the flow field around circular cylinder systems. The calculations are carried out with cylinder systems with and without fairings, while the fairings have different top shape angles within the range of $30°~90°$. The Reynolds number ranges from 1000 to 50 000. It is shown that the simulation results of present numerical method reaches good agreement with the available experimental and numerical simulation data of typical circular cylinder flow and a fixed fairing cylinder system flow. Compared with bare cylinder, the faired bluff structures can obviously reduce the lift and drag forces and alter the vortex shedding frequency. Overall, the mean drag coefficient can be reduced up to about (10–31)% and the RMS lift coefficient can be reduced up to (30–99)% for all faired systems at given Reynolds numbers. The influence of Reynolds number and attack angles on the flow field characters of bare cylinder and faired cylinders is also discussed. The faired structures with shape angles within $30°~45°$under zero-attack-angle-inflow case are considered as the optimal structures, with which the mean drag coefficient and the RMS lift coefficient can be reduced up to (26–31)% and (98–99)%, respectively. Considering the influence of attack angles on lift and drag coefficients reduction, 75° shaped faired structure may be taken as a proper option.