Reynolds数对大攻角细长旋成体绕流滚转角效应影响的数值研究

使用低耗散的Roe格式,数值模拟了Reynolds数(Re)对大攻角细长旋成体绕流滚转角效应的影响.模型头部加了几何小扰动块以引发流场的不对称.在较大的Re数(Re=10 5)下,本文的计算结果与实验是相符的,此时细长体的滚转会导致双稳态、双周期现象,即侧向力随滚转角呈现类似方波形式的双周期变化,方波中侧向力基本保持不变的状态对应于流场的正则态,且两个正则态的侧向力方向相反,方波中侧向力基本保持不变的状态对应于流场的正则态,且两个正则态的侧向力方向相反;而在较小的Re数(Re=4 000)下,如果扰动足够大,细长体的滚转将导致不同的双稳态现象,此时两个正则态的侧向力方向相同,而在较小扰动下双稳态现象不再出现;Re数更小时(Re=1 000),即使在较大的扰动下,双稳态现象也不再出现,侧向力随滚动角仍是连续变化的.本文的计算结果表明,Re数越小,流场对头部扰动的感受性越弱.     

Flow past a square cylinder at low Reynolds numbers

Results are presented for the flow past a stationary square cylinder at zero incidence for Reynolds number, Re ? 150. A stabilized finite‐element formulation is employed to discretize the equations of incompressible fluid flow in two‐dimensions. For the first time, values of the laminar separation Reynolds number, Re_s, and separation angle, θ_s, at Re_s are predicted. Also, the variation of θ_s with Re is presented. It is found that the steady separation initiates at Re = 1.15. Contrary to the popular belief that separation originates at the rear sharp corners, it is found to originate from the base point, i.e. θ_s=180^° at Re = Re_s. For Re > 5, θ_s approaches the limit of 135 ^°. The length of the separation bubble increases approximately linearly with increasing Re. The drag coefficient varies as Re^?0.66. Flow characteristics at Re ? 40 are also presented for elliptical cylinders of aspect ratios 0.2, 0.5, 0.8 and 1 (circle) having the same characteristic dimension as the square and major axis oriented normal to the free‐stream. Compared with a circular cylinder, the flow separates at a much lower Re from a square cylinder leading to the formation of a bigger wake (larger bubble length and width). Consequently, at a given Re, the drag on a square cylinder is more than the drag of a circular cylinder. This suggests that a cylinder with square section is more bluff than the one with circular section. Among all the cylinder shapes studied, the square cylinder with sharp corners generates the largest amount of drag. Copyright © 2010 John Wiley & Sons, Ltd.     

Non-Local Interactions and Feedback Instability in a High Reynolds Number Flow

The question of non-locality is considered for a model supersonic flow at high Reynolds number in a channel formed between two parallel plates of different length, using the channel length as a control parameter. Examples are given of time-periodic stable and unstable flows forced by a disturbance positioned in the middle of the channel. It is shown that in certain parameter ranges the flow in a channel of ever increasing length is not approximated by the solutions obtained for infinitely long channels. This is interpreted in terms of a feedback interaction between the flow near the channel ends and the disturbance source. Feedback is shown to result from a slow upstream decay of disturbances coupled with a relatively fast downstream growth of instability waves. For a free (non-forced) flow, the feedback is found to lead to a form of global or resonant instability. Examples of growth rate calculations for the feedback modes are given.     

Flow of a viscous fluid past a flexible membrane at low Reynolds numbers

The numerical calculation of a steady two-dimensional viscous flow past a flexible membrane is treated. Both edges of the membrane are fixed in the flow and its chord is set normal to the flow. The Navier-Stokes equation in terms of the stream function and the vorticity is transformed to the body fitted coordinate system. The numerical calculations, based on a finite difference method and relaxation method, are carried out for several values of the membranes tension for cases when the Reynolds numbers are 5, 10 and 20. It is found that two different shapes of the membranes are possible at a given value of tension and Reynolds number: one with a small deformation, and the other with a large deformation. Two vortices appear in the concave region of the membrane if its deformation increases beyond a certain extent.     

Near-Wake Turbulence Properties around a Circular Cylinder at High Reynolds Number

The present study investigates the turbulent properties of the flow around a circular cylinder in the near-wake and in the near-wall upstream region at the Reynolds number 140,000. A detailed cartography of the mean and turbulent velocity fields using a moderate blockage and aspect ratio is provided in order to use the present results for direct comparisons with realisable 3D Navier-Stokes computations. The flow structure is analysed by means of two experiments using respectively the LDV and the PIV techniques, both providing a refined grid of measurement points. The dynamics of the separation region, the growth and decay of turbulence in the near wake, as well as the spatial growth of the organised mode are analysed.     

Flow of a viscous fluid past a heterogeneous porous sphere at low Reynolds numbers

A flow past a heterogeneous porous sphere is investigated by using the perturbation theory. The flow through the sphere is divided into two zones, which are fully saturated with the viscous fluid, and the flow in these zones is governed by the Brinkman equation. The space outside the sphere, where a clear fluid flows, is also divided into two zones: the Navier–Stokes zone and the Oseen flow zone. The solutions on the interface inside the sphere are matched with the condition proposed by Merrikh and Mohammad. The stream function in the Navier–Stokes zone is matched with that on the sphere surface by the condition proposed by Ochoa-Tapia and Whitaker. It is found that the drag on the spherical shell decreases as the permeability toward the sphere boundary increases.     

Flow past a transversely oscillating square cylinder in free stream at low Reynolds numbers

This paper reports simulation results for free‐stream flow past an oscillating square cylinder at Re=100 and 150, for oscillating‐to‐natural‐shedding frequency ratios of 0.5?f_r?3.0 at a fixed oscillation amplitude of 0.2 of the cylinder width. The transformed governing equations are solved in a non‐inertial frame of reference using the finite volume technique. The ‘lock‐in’ phenomena, where the vortex shedding becomes one with the oscillation frequency, is observed near the natural shedding frequency (f_r≈1). Beyond the synchronization band, downstream recovery of the wake to its stationary (natural) state (frequency) is observed in cross‐stream velocity spectra. At higher forcing frequencies, a phase lag between the immediate and the far wake results in a shear layer having multi‐polar vortices. A ‘Vortex‐switch’ accompanied by a change in the direction of energy transfer is identified at the ‘lock‐in’ boundaries. The variation of aerodynamic forces is noticed to be different in the lock‐in regime. The velocity phase portrait in the far wake revealed a chaotic state of flow at higher excitation though a single (natural) frequency appears in the spectra. Copyright © 2008 John Wiley & Sons, Ltd.     

Unsteady Separation Processes at High Reynolds Number and Their Control

Physical situations where a viscous boundary layer breaks down and interacts strongly with an effectively inviscid external flow are common place. For large Reynolds numbers, viscous effects are normally confined to thin boundary layers on all solid surfaces for the majority of any observation time. In most practical situations, exposure of such layers to an adverse pressure gradient is inevitable and in this circumstance, a sequence of events commences near the wall that culminates in an eruption and a strong viscous-inviscid interaction with the external flow. The events leading up to eruption are known as the Van Dommelen–Shen process and the eruption itself is referred to as boundary-layer separation; here the term ‘separation’ denotes the first process of interaction between a hitherto thin boundary layer and the external flow. The event is sufficiently complicated that extraordinary measures are needed to compute its evolution. In most situations, the onset of separation is subtle and hard to detect and thus development of rational control procedures is a challenging task. Here recent calculations of unsteady separation events are discussed for two- and three-dimensional flows. The phenomena involved are generic but leading-edge separation on airfoils and rotorcraft blades is emphasized. Recent studies on various control mechanisms are described, which are found to have the effect of slowing down and/or weakening the separation process. For some control processes, it has proved possible to eliminate separation entirely.     

Flow Rate Prediction for a Semi-permeable Membrane at Low Reynolds Number in a Circular Pipe

Transport in Porous Media - A concise and accurate prediction method is required for membrane permeability in chemical engineering and biological fields. As a preliminary study on this topic, we...     

Approximate Wall Boundary Conditions in the Large-Eddy Simulation of High Reynolds Number Flow

The near-wall regions of high Reynolds numbers turbulent flows must be modelled to treat many practical engineering and aeronautical applications. In this review we examine results from simulations of both attached and separated flows on coarse grids in which the near-wall regions are not resolved and are instead represented by approximate wall boundary conditions. The simulations use the dynamic Smagorinsky subgrid-scale model and a second-order finite-difference method. Typical results are found to be mixed, with acceptable results found in many cases in the core of the flow far from the walls, provided there is adequate numerical resolution, but with poorer results generally found near the wall. Deficiencies in this approach are caused in part by both inaccuracies in subgrid-scale modelling and numerical errors in the low-order finite-difference method on coarse near-wall grids, which should be taken into account when constructing models and performing large-eddy simulation on coarse grids. A promising new method for developing wall models from optimal control theory is also discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Hydrodynamic Interactions at Low Reynolds Number

We consider the hydrodynamic interactions of low Reynolds number microswimmers, presenting a review of recent work based upon models of linked sphere swimmers. Particular attention is paid to those aspects that are generic, applicable to all microswimmers and not only to the simple models considered. The importance of the relative phase in swimmer–swimmer interactions is emphasised, as is the role of simple symmetry arguments in both understanding and constraining the hydrodynamic properties of microswimmers.     

雷诺数对三角翼绕流的影响

应用染色液流动显示技术研究了雷诺数对６０°尖前缘三角翼前缘涡破裂位置、背风面流动结构等的影响，并详细分析了背风面流动随攻角的变化．     

Statistical Structure at the Wall of the High Reynolds Number Turbulent Boundary Layer

The one and two-point statistical structure of very high Reynolds number turbulence in the surface layer near a rigid `wall' is analysed. The essential mechanisms for turbulent eddies impinging on the wall are studied using linearised rapid distortion theory, which show how the mean shear and blocking actions of the surface act first independently and then, over the life time of the eddy, interactively. Previous analytical results are reinterpreted and some new results are derived to show how the integral length scales, cross correlations and spectra of the different components of the turbulence are distorted depending on the form of the spectra of eddies above the surface layer and how they are related to motions of characteristic eddy structures near the surface. These results are applied to derive some quantitative and qualitative predictions in the surface layers (SL), where the eddies are affected by local shear dynamics, and in the `eddy surface layer' (ESL) where quasi independents loping elongated eddies interact directly with the wall, and where there is a large range of wave number within which the spectra of the horizontal velocity components are proportional to k ⁻¹. The longest eddies in the boundary layer occur near the wall. Field experiments agree with the theoretical model predictions for the quite different forms for the spectra, cospectra and cross correlations for the vertical and horizontal components of the velocity field. By showing that in SL the energy exchange between the large and small scale eddies is local(`staircase') energy cascade, whereas in ESL there is a direct nonlocal (`elevator-like')energy transfer to the small scales, it is shown why the thickness of the ESL increases over rougher surfaces and as the Reynolds number decreases. This revised version was published online in August 2006 with corrections to the Cover Date.     

Large-Eddy Simulation of the Flow Over a Circular Cylinder at Reynolds Number 3900 Using the OpenFOAM Toolbox

The flow over a circular cylinder at Reynolds number 3900 and Mach number 0.2 was predicted numerically using the technique of large-eddy simulation. The computations were carried out with an O-type curvilinear grid of size of 300 × 300 × 64. The numerical simulations were performed using a second-order finite-volume method with central-difference schemes for the approximation of convective terms. A conventional Smagorinsky and a dynamic k-equation eddy viscosity sub-grid scale models were applied. The integration time interval for data sampling was extended up to 150 vortex shedding periods for the purpose of obtaining a fully converged mean flow field. The present numerical results were found to be in good agreement with existing experimental data and previously obtained large-eddy simulation results. This gives an indication on the adequacy and accuracy of the selected large-eddy simulation technique implemented in the OpenFOAM toolbox.     

Reynolds Number Dependence of Velocity Structure Functions in a Turbulent Pipe Flow

Low-order moments of the increments δu andδv where u and v are the axial and radial velocity fluctuations respectively, have been obtained using single and X-hot wires mainly on the axis of a fully developed pipe flow for different values of the Taylor microscale Reynolds numberR _λ. The mean energy dissipation rate〉ε〈 was inferred from the uspectrum after the latter was corrected for the spatial resolution of the hot-wire probes. The corrected Kolmogorov-normalized second-order structure functions show a continuous evolution withR _λ. In particular, the scaling exponentζ _v , corresponding to the v structure function, continues to increase with R _λ in contrast to the nearly unchanged value of ζ _u . The Kolmogorov constant for δu shows a smaller rate of increase with R _λ than that forδv. The level of agreement with local isotropy is examined in the context of the competing influences ofR _λ and the mean shear. There is close but not perfect agreement between the present results on the pipe axis and those on the centreline of a fully developed channel flow. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

Direct Numerical Simulation of Turbulent Pipe Flow at Moderately High Reynolds Numbers

Fully resolved direct numerical simulations (DNSs) have been performed with a high-order spectral element method to study the flow of an incompressible viscous fluid in a smooth circular pipe of radius R and axial length 25R in the turbulent flow regime at four different friction Reynolds numbers Re _τ ?=?180, 360, 550 and $1\text{,}000$ . The new set of data is put into perspective with other simulation data sets, obtained in pipe, channel and boundary layer geometry. In particular, differences between different pipe DNS are highlighted. It turns out that the pressure is the variable which differs the most between pipes, channels and boundary layers, leading to significantly different mean and pressure fluctuations, potentially linked to a stronger wake region. In the buffer layer, the variation with Reynolds number of the inner peak of axial velocity fluctuation intensity is similar between channel and boundary layer flows, but lower for the pipe, while the inner peak of the pressure fluctuations show negligible differences between pipe and channel flows but is clearly lower than that for the boundary layer, which is the same behaviour as for the fluctuating wall shear stress. Finally, turbulent kinetic energy budgets are almost indistinguishable between the canonical flows close to the wall (up to y ^?+??≈?100), while substantial differences are observed in production and dissipation in the outer layer. A clear Reynolds number dependency is documented for the three flow configurations.     

Formula for the Flow Resistance Factor in a Pipe with a Sudden Expansion at Small Reynolds Numbers

A formula for the flow resistance factors in a pipe with a sudden expansion of the cross section at Reynolds numbers of 0.2 to 10 is obtained by numerical solution of the complete Navier–Stokes equations for incompressible fluids. The flow resistance factors obtained using the derived formula are compared to those found by numerical solution of the Navier–Stokes equations.     

On the Lower Critical Reynolds Number for Flow in a Circular Pipe

Flow in a circular pipe is investigated experimentally at Reynolds numbers higher than that at which the resistance coefficients calculated from the Blasius formula for laminar flow and from the Prandtl formula for turbulent flow are equal. The corresponding Reynolds number based on the mean-flow velocity and the pipe diameter is about 1000. The experiments were performed at a high level of inlet pulsations produced by feeding gas into the pipe through a hole with a diameter several times smaller than the pipe diameter. In our experiments the critical Reynolds number was determined as the value, independent of the distance from the inlet, at which the ratio of the axial to the mean-flow velocity as a function of the Reynolds number deviated from 2. At the maximum ratio of the pipe cross-sectional area to the area of the hole through which the gas entered the pipe, equal to 26, the critical Reynolds number was about 2300. After a fivefold increase in the hole area the critical Reynolds number increased by approximately 4%.At Reynolds numbers below 2000, after at a high level of the inlet pulsations an almost laminar flow had developed in the pipe, a perturbation was introduced by inserting a diametrically oriented cylindrical rod with a diameter 10–20 times smaller than the pipe diameter. In these experiments, at Reynolds numbers higher than 1000, at a distance from the rod equal to 50 pipe diameters the axial to mean-flow velocity ratio was less than 2, approaching this value again at large distances from the rod. The insertion of the rod led to a decrease in the critical Reynolds number by approximately 12%.     

A Unified Introduction to Fluid Mechanics of Flying and Swimming at High Reynolds Number

Modeling the flow around a deformable and moving surface is required to calculate the forces exerted by a swimming or flying animal on the surrounding fluid. Assuming that viscosity plays a minor role, linear potential models can be used. These models derived from unsteady airfoil theory are usually divided in two categories depending on the aspect ratio of the moving surface: for small aspect ratios, slender-body theory applies while for large aspect ratios two-dimensional or lifting-line theory is used. This paper aims at presenting these models with a unified approach. These potential models being analytical, they allow fast computations and can therefore be used for optimization or control.