圆角化对方柱气动性能影响的流场机理

方形截面柱体的圆角化处理是常用的流动控制方法,但其流场作用机理尚未被澄清.采用大涡模拟方法,在雷诺数为2.2$\times$10$^{4}$时,考虑风攻角的影响,对均匀流作用下的标准方柱和圆角方柱的气动性能和流场特性进行了研究,定量分析了圆角化气动措施和风攻角变化对分离泡特性的影响规律,从流场角度澄清了圆角化气动措施对方柱气动性能的影响机理.研究表明:与标准方柱相比,圆角方柱的表面风压、气动力和涡脱强度呈整体下降的趋势,但圆角方柱的斯特劳哈尔数更高;圆角方柱的"分离泡流态'发生在更小的风攻角范围内,分离泡的出现会进一步造成方柱的尾流变窄,涡脱强度减弱;随着风攻角的增大,分离泡的长度会逐渐减小直至消失,分离泡的中心会逐渐向方柱前角(迎风向)和方柱壁面移动;与标准方柱相比,圆角方柱的气流发生初次分离的位置向下游移动,分离后的剪切层更贴近方柱,因而更易发生再附现象;方柱尾流宽度的减小和涡脱强度的减弱是导致圆角方柱气动力减小和斯特劳哈尔数增大的主要原因.      

熔喷双槽形喷嘴气体射流流场初探

在熔喷非织造布加工中,气体射流作为工作介质使聚合物熔体实现拉伸,气体射流流场的研究对熔喷气流拉抻数学模型研究非常重要。熔喷双槽形喷嘴形成的流场可以看作两股平面射流的合成。从单个点涡的性质出发,研究了涡偶的性质和涡偶代替射流的可行性。研究表明,在喷丝孔轴线附近,涡偶和射流的速度分布趋势相同,且有比较相近的速度分布,从而说明以涡偶代替射流是可行的。在此基础上,用两个涡偶分别代替两股射流,然后进行合成,推导出两股射流合成后速度分布的理论公式,该公式的计算结果与实验结果吻合较好。将该公式引入熔喷气流拉伸数学模型,预测出的纤维直径与采用经验公式时的预测结果几乎完全相同。结果表明,应用涡偶代替射流推导出的气流速度分布公式能够较好地描述熔喷双槽形喷嘴气体射流流场,可以用于完善熔喷气流拉伸数学模型。     

Wingtip vortex control via the use of a reverse half-delta wing

The effect of a 65° sweep reverse half-delta wing (RHDW), mounted at the squared tip of a rectangular NACA 0012 wing, on the tip vortex was investigated experimentally at Re?=?2.45?×?10⁵. The RHDW was found to produce a weaker tip vortex with a lower vorticity level and, more importantly, a reduced lift-induced drag compared to the baseline wing. In addition to the lift increment, the RHDW also produced a large separated wake flow and subsequently an increased profile drag. The reduction in lift-induced drag, however, outperformed the increase in profile drag and resulted in a virtually unchanged total drag in comparison with the baseline wing. Physical mechanisms responsible for the RHDW-induced appealing aerodynamics and vortex flow modifications were discussed.     

高亚声速运输机后体减阻措施研究

为了降低高亚声速运输机的阻力,改善飞机的气动性能,本文在某一典型运输机机身外形的基础上,研究了后体参数和涡流发生器对机身后体流场及阻力的影响。通过数值模拟计算,得出了不同外形的后体流态分布及阻力特性数据。结果表明,小的上翘角和收缩比、大的扁平度和长细比能够显著地减小压差阻力,而涡流发生器能够有效地减小摩擦阻力。该结果可为后续的型号研制提供参考。     

Evolution of an initially columnar vortex terminating normal to a no-slip wall

The early evolution of an initially columnar vortex normal to a solid wall was examined. The vortex was generated by a pair of flaps in a water tank. Detrimental effects from the wall during the vortex generation were avoided by producing the vortex normal to a free surface and subsequently bringing a horizontal plate into contact with the surface. Digital particle image velocimetry (DPIV) measurements of the velocity and vorticity, together with laser induced fluorescence (LIF) visualizations, in a meridional plane revealed a toroidal structure with the appearance of an axisymmetric vortex breakdown bubble. Agreement was found between the measurements and numerical simulations of the axisymmetric Navier–Stokes equations. The results show that the flow in the effusive corner region is dominated by a Bödewadt-type spatially oscillatory boundary layer within the core region and a potential-like vortex boundary layer at large radii. The toroidal structure results from the interaction between these two boundary layers, leading to the roll up of a dominant shear layer within the Bödewadt structure, and does not develop from the columnar vortex itself.     

Wall shear stress characterization of a 3D bluff-body separated flow

Efficient flow control strategies aimed at reducing the aerodynamic drag of road vehicles require a detailed knowledge of the reference flow. In this work, the flow around the rear slanted window of a generic car model was experimentally studied through wall shear stress measurements using an electrochemical method. The mean and fluctuating wall shear stress within the wall impact regions of the recirculation bubble and the main longitudinal vortex structures which develop above the rear window are presented. Correlations allow a more detailed characterization of the recirculation phenomenon within the separation bubble. In the model symmetry plane the recirculation structure compares well with simpler 2D configurations; specific lengths, flapping motion and shedding of large-scale vortices are observed, these similarities diminish when leaving the middle plane due to the strong three-dimensionality of the flow. A specific attention is paid to the convection processes occurring within the recirculation: a downstream convection velocity is observed, in accordance with 2D recirculations from the literature, and an upstream convection is highlighted along the entire bubble length which has not been underlined in some previous canonical configurations.     

Upwind sail aerodynamics: A RANS numerical investigation validated with wind tunnel pressure measurements

The aerodynamics of a sailing yacht with different sail trims are presented, derived from simulations performed using computational fluid dynamics. A Reynolds-averaged Navier–Stokes approach was used to model sixteen sail trims first tested in a wind tunnel, where the pressure distributions on the sails were measured. An original approach was employed by using two successive simulations: the first one on a large domain to model the blockage due to the wind tunnel walls and the sails model, and a second one on a smaller domain to model the flow around the sails model. A verification and validation of the computed aerodynamic forces and pressure distributions was performed. The computed pressure distribution is shown to agree well with the measured pressures. The sail surface pressure was correlated with the increase of turbulent viscosity in the laminar separation bubble, the flow reattachment and the trailing edge separation. The drive force distribution on both sails showed that the fore part of the genoa (fore sail) provides the majority of the drive force and that the effect of the aft sail is mostly to produce an upwash effect on the genoa. An aerodynamic model based on potential flow theory and a viscous correction is proposed. This model, with one free parameter to be determined, is shown to fit the results better than the usual form drag and induced drag only, even if no friction drag is explicitly considered.     

Drag reduction and thrust generation by tangential surface motion in flow past a cylinder

Sensitivity of drag to tangential surface motion is calculated in flow past a circular cylinder in both two- and three-dimensional conditions at Reynolds number \(\textit{Re} \le 1000\). The magnitude of the sensitivity maximises in the region slightly upstream of the separation points where the contour lines of spanwise vorticity are normal to the cylinder surface. A control to reduce drag can be obtained by (negatively) scaling the sensitivity. The high correlation of sensitivities of controlled and uncontrolled flow indicates that the scaled sensitivity is a good approximation of the nonlinear optimal control. It is validated through direct numerical simulations that the linear range of the steady control is much higher than the unsteady control, which synchronises the vortex shedding and induces lock-in effects. The steady control injects angular momentum into the separating boundary layer, stabilises the flow and increases the base pressure significantly. At \(\textit{Re}=100\), when the maximum tangential motion reaches 50% of the free-stream velocity, the vortex shedding, boundary-layer separation and recirculation bubbles are eliminated and 32% of the drag is reduced. When the maximum tangential motion reaches 2.5 times of the free-stream velocity, thrust is generated and the power savings ratio, defined as the ratio of the reduced drag power to the control input power, reaches 19.6. The mechanism of drag reduction is attributed to the change of the radial gradient of spanwise vorticity (\(\partial _{r} \hat{\zeta }\)) and the subsequent accelerated pressure recovery from the uncontrolled separation points to the rear stagnation point.     

A study of the aerodynamics of a generic container freight wagon using Large-Eddy Simulation

In this work simulations using the Large Eddy Simulation technique have been made of the flow around a generic container freight wagon model. The model consists of one 11.8 m standard length container placed on a wagon. Details of the undercarriage such as wheels are included, but the container is generic and smoothed in comparison to a real freight wagon. The Reynolds number of the flow is 10⁵ based on the container width of 2.354 m. Two cases have been considered in the study, one case where the wagon is standing alone and one case where it is submerged into a train set with wagons ahead and behind the wagon. The latter case is simulated using periodic boundary condition. Both the time-averaged and the instantaneous flow around the wagon for the two cases are described. For the single wagon case, it is found that the separation bubble formed on the roof of the container oscillates back and forth in the streamwise direction and that this oscillation is in phase with oscillations found in the upper shear layer of the ring vortex in the wake. The mechanism that is causing the synchronization of the oscillations of the separation bubble at the front and the upper shear layers in the wake is found to be waves of vorticity being shed from the separation bubble. The time-averaged ring vortex in the near wake of the single wagon is found to be inclined due to the disturbance of the undercarriage details on flow in the lower shear layer. The lower center of the ring vortex is located closer to the base face than the upper center. The drag coefficient of the wagon in the periodic case was found to be only 10% of that of the single wagon case. This is due to two symmetrical counter-rotating vortices found in the gaps which make the train set appear as a single body to the oncoming flow and shielding the wagon from any direct impingement of the flow. The counter-rotating vortices in the gap are found to inhibit periodic oscillations in the lateral direction. These oscillations cause vortical structures to form by the air that is pushed out from the gap and these flow structures cause a dominating oscillation of non-dimensional frequency St=0.12 in the side force signal.     

Early stages of an impulsively started unsteady laminar flow past tapered trapezoidal cylinders

Characteristics of the developing recirculation region behind a tapered trapezoidal cylinder and its interaction with the separating shear layer from the leading edges were studied numerically for an impulsively started laminar flow. An unsteady stream function–vorticity formulation was used. The Reynolds numbers considered range from 25 to 1000. Pressure contours, surface pressure coefficient, wake length and drag coefficient were studied through the streamline flow field. Main flow and subflow regimes were identified by an analysis of the evolution of the flow characteristics. It was found that typically, for a given trapezoidal cylinder, flow starts with no separation. As time advances, the symmetrical standing zone of recirculation develops aft of the trapezoidal cylinder. The rate of growth in width, length and structure of the aft end eddies depends on the Reynolds number. In time, separated flow from the leading edges of the trapezoidal cylinder also develops and forms growing separation bubbles on the upper and lower inclined surfaces of the trapezoidal cylinder. As time advances, the separation bubbles on the upper and lower inclined surfaces of the cylinder grow towards the downstream regions and eventually merge with the swelling symmetrical eddies aft of the cylinder. This merging of the flows creates a complex flow regime with a disturbed tertiary flow zone near the merging junction. Eventually, depending on the Reynolds number and the tapered angle of the trapezoidal cylinder, the flow develops into a specific category of symmetrical standing recirculatory flow with its own distinct characteristics. Comparisons with the available results of other investigators showed very good agreement. © 1998 John Wiley & Sons, Ltd.     

Distributed forcing of the flow past a blunt-based axisymmetric bluff body

In this paper, we address the influence of a blowing-/suction-type distributed forcing on the flow past a blunt-based axisymmetric bluff body by means of direct numerical simulations. The forcing is applied via consecutive blowing and suction slots azimuthally distributed along the trailing edge of the bluff body. We examine the impact of the forcing wavelength, amplitude and waveform on the drag experienced by the bluff body and on the occurrence of the reflectional symmetry preserving and reflectional symmetry breaking wake modes, for Reynolds numbers 800 and 1,000. We show that forcing the flow at wavelengths inherent to the unforced flow drastically damps drag oscillations associated with the vortex shedding and vorticity bursts, up to their complete suppression. The overall parameter analysis suggests that this damping results from the surplus of streamwise vorticity provided by the forcing that tends to stabilize the ternary vorticity lobes observed at the aft part of the bluff body. In addition, conversely to a blowing-type or suction-type forcing, the blowing-/suction-type forcing involves strong nonlinear interactions between locally decelerated and accelerated regions, severely affecting both the mean drag and the frequencies representative of the vortex shedding and vorticity bursts.     

PIV investigation of flow behind surface mounted permeable ribs

The flow behind surface mounted permeable rib geometries, i.e. solid, slit, split-slit and inclined split-slit ribs have been studied using flow visualization and PIV (2-C and 3-C) technique in streamwise and cross-stream measurement planes. The objective behind this study is to understand the flow structures responsible for heat transfer/mixing enhancement with simultaneous pressure penalty reduction by permeable rib geometries. The Reynolds number based on the rib height has been set equal to 5,538 and the open area ratio of permeable ribs is equal to 20%. The permeable rib geometries have shorter reattachment length in comparison to the solid rib. The maximum 41% reduction in reattachment length is observed for the inclined split-slit rib. The splitter mounted inside the slit leads to two corner vortices behind it. The corner vortices drag the flow from the primary recirculation bubble region towards the rib resulting in drop of the reattachment length. Two horseshoe vortices are present in the flow through the slit at both sides of the splitter due to the upstream flow separation. The slit inclination moves these horseshoe vortices closer to the bottom wall. A film like flow through the slit is present near the downstream corner of the inclined split-slit rib. The spanwise velocity gradient due to the splitter leads to vorticity and turbulence enhancement by vortex stretching. The inclination of the slit and the use of a splitter inside the slit are two important design parameters responsible in generation of near-wall longitudinal vortices. The flow field behind permeable ribs is dominated by vortical structures with definable critical flow patterns, i.e. node, saddle and foci. These predominant swirling flow motions contribute to the mixing enhancement behind permeable rib geometries. On leave from Mechanical Engineering Department, IIT Kanpur, U.P. 208016, India     

Drag and lift reduction of a 3D bluff-body using active vortex generators

In this study, a passive flow control experiment on a 3D bluff-body using vortex generators (VGs) is presented. The bluff-body is a modified Ahmed body (Ahmed in J Fluids Eng 105:429–434 1983) with a curved rear part, instead of a slanted one, so that the location of the flow separation is no longer forced by the geometry. The influence of a line of non-conventional trapezoïdal VGs on the aerodynamic forces (drag and lift) induced on the bluff-body is investigated. The high sensitivity to many geometric (angle between the trapezoïdal element and the wall, spanwise spacing between the VGs, longitudinal location on the curved surface) and physical (freestream velocity) parameters is clearly demonstrated. The maximum drag reduction is ?12%, while the maximum global lift reduction can reach more than ?60%, with a strong dependency on the freestream velocity. For some configurations, the lift on the rear axle of the model can be inverted (?104%). It is also shown that the VGs are still efficient even downstream of the natural separation line. Finally, a dynamic parameter is chosen and a new set-up with motorized vortex generators is proposed. Thanks to this active device. The optimal configurations depending on two parameters are found more easily, and a significant drag and lift reduction (up to ?14% drag reduction) can be reached for different freestream velocities. These results are then analyzed through wall pressure and velocity measurements in the near-wake of the bluff-body with and without control. It appears that the largest drag and lift reduction is clearly associated to a strong increase of the size of the recirculation bubble over the rear slant. Investigation of the velocity field in a cross-section downstream the model reveals that, in the same time, the intensity of the longitudinal trailing vortices is strongly reduced, suggesting that the drag reduction is due to the breakdown of the balance between the separation bubble and the longitudinal vortices. It demonstrates that for low aspect ratio 3D bluff-bodies, like road vehicles, the flow control strategy is much different from the one used on airfoils: an early separation of the boundary layer can lead to a significant drag reduction if the circulation of the trailing vortices is reduced.     

Flow development upstream of a fence

The effect of Reynolds number on the flow development upstream of a rigid, non-porous, static fence is investigated experimentally. The flow field is measured using time-resolved, two-component particle image velocimetry at Reynolds numbers based on fence height of 18000, 36000, and 54000. The results show that a laminar separation bubble forms upstream of the junction vortex at the base of the fence. The mean extent of the bubble decreases with increasing Reynolds number, with mean separation moving downstream and mean reattachment moving upstream. In the aft portion of the bubble, shear layer vortices form and are shed at scaled frequencies and wavelengths that are comparable to laminar separation bubble shedding in low Reynolds number airfoils and flat plates with an imposed adverse pressure gradient. The strong periodicity of the associated coherent structures and the proximity of shear layer roll-up relative to the fence should be taken into consideration in the relevant designs due to potential implications to structural loading. A simple flow separation prediction model combining inviscid fence flow solution with Thwaites’ method is introduced and shows good agreement with the experimental results for the Reynolds number range considered.     

NUMERICAL PREDICTION OF FORCE ON RECTANGULAR CYLINDERS IN OSCILLATING VISCOUS FLOW

Numerical results are presented for an oscillating viscous flow past a square cylinder with square and rounded corners and a diamond cylinder with square corners at Keulegan–Carpenter numbers up to 5. This unsteady flow problem is formulated by the two-dimensional Navier–Stokes equations in vorticity and stream-function form on body-fitted coordinates and solved by a finite-difference method. Second-order Adams-Bashforth and central-difference schemes are used to discretize the vorticity transport equation while a third-order upwinding scheme is incorporated to represent the nonlinear convective terms. Since the vorticity distribution has a mathematical singularity at a sharp corner and since the force coefficients are found in experiments to be sensitive to the corner radius of rectangular cylinders, a grid-generation technique is applied to provide an efficient mesh system for this complex flow. Local grid concentration near the sharp corners, instead of any artificial treatment of the sharp corners being introduced, is used in order to obtain high numerical resolution. The elliptic partial differential equation for stream function and vorticity in the transformed plane is solved by a multigrid iteration method. For an oscillating flow past a rectangular cylinder, vortex detachment occurs at irregular high frequency modes at KC numbers larger than 3 for a square cylinder, larger than 1 for a diamond cylinder and larger than 3 for a square cylinder with rounded corners. The calculated drag and inertia coefficients are in very good agreement with the experimental data. The calculated vortex patterns are used to explain some of the force coefficient behavior.     

鼓泡流化床中流动特性的多尺度数值模拟

鼓泡流化床因其较高的传热特性以及较好的相间接触已经被广泛应用于工业生产中,而对鼓泡流态化气固流动特性的充分认知是鼓泡流化床设计的关键.在鼓泡流化床中,气泡相和乳化相的同时存在使得床中呈现非均匀流动结构,而这种非均匀结构给鼓泡流化床的数值模拟造成了很大的误差.基于此,以气泡作为介尺度结构,建立了多尺度曳力消耗能量最小的稳定性条件,构建了适用于鼓泡流化床的多尺度气固相间曳力模型.结合双流体模型,对A类和B类颗粒的鼓泡流化床中气固流动特性进行了模拟研究,分析了气泡速度、气泡直径等参数的变化规律.研究表明,与传统的曳力模型相比,考虑气泡影响的多尺度气固相间曳力模型给出的曳力系数与颗粒浓度的关系是一条分布带,建立了控制体内曳力系数与局部结构参数之间的关系.通过模拟得到的颗粒浓度和速度与实验的比较可以发现,考虑气泡影响的多尺度曳力模型可以更好地再现实验结果.通过A类和B类颗粒的鼓泡床模拟研究发现,A类颗粒的鼓泡床模拟受多尺度曳力模型的影响更为显著.      

Flow behavior around a square cylinder subject to modulation of a planar jet issued from upstream surface

The flow behavior in the up- and downstream regions of a square cylinder subject to the modulation of a planar jet issued from the cylinder׳s front surface was studied using the laser-assisted smoke flow visualization method and hot-wire anemometer measurement. Reynolds numbers were from 1628 to 13 000. The drag force experienced by the square cylinder was obtained by measuring the surface pressures on the up- and downstream faces. The temporally evolving smoke flow patterns in the up- and downstream regions were synchronously revealed through the smoke flow visualization. The frequency characteristics of the instability waves in the up- and downstream regions were synchronously detected by the two hot-wire anemometers. Four characteristic flow modes were observed within the different ranges of the injection ratios. At the low injection ratios (IR<1), the ‘swinging jet’ appeared. The jet swung periodically leftward and rightward and formed a fluid bubble on the front surface. The fluid bubble contained a pair of counter-rotating vortices and presented a periodic variation in its height. At moderately low injection ratios (1<IR<4.3), the ‘deflected oscillating jet’ appeared. The jet was deflected in either the left or the right direction and wrapped around one of the edges of the square cylinder. Both the swinging and oscillating motions of the jet in the swinging jet and deflected oscillating jet modes were induced by the periodic feedback pressure signals generated by the vortex shedding in the wake. At the moderately high (4.3<IR<8.3) and high (IR>8.3) injection ratios, the ‘deflection jet’ and ‘penetrating jet’ appeared. The jet detached from the cylinder׳s front surface and penetrated a long distance into the upstream region due to large jet momentum. Neither periodic jet oscillation in the upstream region nor vortex shedding in the wake was observed. The drag coefficient was found to be decreasing quickly with increasing the injection ratio.     

Effect of vortex shedding in unsteady aerodynamic forces for a low Reynolds number stationary wing at low angle of attack

This study elucidates the relation between wake vortex shedding and aerodynamic force fluctuations for a low Reynolds number wing from time resolved particle image velocimetry (TR-PIV) experimental measurements. The results reveal a periodic lift and drag variation within the shedding cycle and resolve the frequencies of those fluctuations from a proper orthogonal decomposition (POD) and power spectral density (PSD) analysis. To show the effect of vortex shedding on the body force fluctuations, the evolution of instantaneous aerodynamic forces is compared to the pressure field of the fluid flow and to the vortical structures in the wake of the airfoil. A six step model describing the vortex-force relation is proposed. It shows that changes in lift such as maximum lift and minimum lift are associated with the detachment of a vortex. It also shows that the minimum or local minimum drag value is obtained at the onset formation of a vortex on the airfoil wake. Similarly, the maximum or local maximum drag is obtained at the onset formation of the saddle on the airfoil wake. The model further explains the asymmetry observed in the unsteady drag force evolution. The model can be used to optimize flow control and fluid-structure interaction applications.     

Interaction of a line vortex with a round parachute canopy

The interaction of a rectilinear vortex with an inflated round parachute canopy model was studied experimentally in a water tunnel where the vortex core was aligned with the axis of the canopy. Three different canopy diameters were used, and the canopy model was attached to a streamlined forebody. Dye flow visualization indicated that vortex breakdown was present when the core trajectory was within the canopy opening. Vortex breakdown occurred about one to two canopy diameters upstream of the canopy opening. The vortex core completely disintegrated when it interacted with the forebody near the canopy centerline. The vortex breakdown and disintegration caused unsteady, asymmetric deformations on the canopy surface. A reduction in the time-averaged drag and an increase in the fluctuating drag was observed when the vortex core was within the canopy opening. The disintegration of the vortex core near the canopy centerline lessened the drag reduction brought on by the presence of the core.     

EMMS mixture model with size distribution for two-fluid simulation of riser flows

Further development of an energy-minimization multiscale modeling approach to simulating two-phase flow under turbulent conditions that considers the size distribution of mesoscale structures, i.e. bubbles and clusters, is presented. User-defined values of minimum and maximum cluster or bubble diameters were specified. A uniform size distribution was first considered as a test case, in which the drag force comprised contributions from each size group. The mathematical form of the objective function describing the energy for suspension and transport was not altered. The heterogeneity index of this new drag modification was then used to simulate pilot-scale circulating fluidized-bed risers involving Geldart group A particles. The results were validated against available experimental data. The model is capable of capturing both axial and radial profiles of flow-field variables.