方案详情文
智能文字提取功能测试中
Exp Fluids (2017)58:83DOI 10.1007/s00348-017-2331-0CrossMark 83 Page 2 of 14Exp Fluids (2017) 58:83 Aerodynamic drag ofa transiting sphere by large-scaletomographic-PIV W. Terral· A. Sciacchitanol·F. Scarano AbstracttAA method is introduced to measure the aero-dynamic drag of moving objects such as ground vehiclesor athletes in speed sports. Experiments are conductedas proof-of-concept that yield the aerodynamic drag ofa sphere towed through a square duct in stagnant air. Thedrag force is evaluated using large-scale tomographic PIVand invoking the time-average momentum equation withina control volume in a frame of reference moving with theobject. The sphere with 0.1 m diameter moves at a veloc-ity of 1.45 m/s, corresponding to a Reynolds number of10,000. The measurements in the wake of the sphere areconducted at a rate of 500 Hz within a thin volume ofapproximately 3×40×40 cubic centimeters. Neutrallybuoyant helium-filled soap bubbles are used as flow trac-ers. The terms composing the drag are related to the flowmomentum, the pressure and the velocity fluctuations andthey are separately evaluated. The momentum and pressureterms dominate the momentum budget in the near wakeup to 1.3 diameters downstream of the model. The pres-sure term decays rapidly and vanishes within 5 diameters.The term due to velocity fluctuations contributes up to10% to the drag. The measurements yield a relatively con-stant value of the drag coefficient starting from 2 diametersdownstream of the sphere. At 7 diameters the measurementinterval terminates due to the finite length of the duct. Errorsources that need to be accounted for are the sphere support Electronic supplementary material The online version of thisarticle (doi:10.1007/s00348-017-2331-0) contains supplementarymaterial, which is available to authorized users. ( 区 W. Terraw.terra@tudelft.nl ) ( Aerospace Engineering Department, TU De l ft, Delft, The Netherlands ) wake and blockage effects. The above findings can providepractical criteria for the drag evaluation of generic bluffobjects with this measurement technique. Keywords Tomographic PIV ·Momentum equation .Aerodynamic drag·HFSB·Sphere· Wake 1 Introduction Aerodynamic forces result from the relative motion of anobject in air. These forces are relevant for a multitude ofapplications, such as aircraft systems, wind turbines andground transportation. Aerodynamic forces are typicallymeasured by means of wind tunnel experiments (Baconand Reid 1924; Zdravkovich 1990 among others), wherethe model is immersed in a uniform air stream within accu-rately controlled conditions. In other cases, aerodynamicstudies are conducted with the object moving in a quies-cent fluid. This approach is needed for instance to study theflow behind accelerating objects (Coutanceau and Bouard1977a) or the development of the wake of an aircraft overa large distance (Scarano et al. 2002; Von Carmer et al.2008). Furthermore, instructive flow visualizations havebeen obtained of high-speed flight of a bullet in stagnant air(van Dyke 1982). Aerodynamic investigations for groundvehicles and speed sports such as cycling, skating and ski-ing require specific arrangements of the wind tunnel (mov-ing floor, model supports) to achieve realistic conditions.Yet, experiments that simulate specific conditions suchas curved trajectories, acceleration, or with the athlete inmotion remain challenging. In contrast, experiments con-ducted with the object moving in the laboratory frame ofreference are relatively easy to realize and can mimic moreclosely real-life conditions. An aspect that has not been sufficiently covered fromtowed experiments, is the quantitative analysis of the flowfield, such to be able to evaluate the aerodynamic forces,the drag in particular. The present work focuses on the lat-ter question to explore the viability of drag estimation fromaerodynamic data obtained during towed experiments. Theresults are intended as a proof-of-concept for applicationsin automotive industry and speed sports in particular. In professional cycling the aerodynamic drag is rou-tinely estimated using the mechanical power generated bythe rider. These estimates, however, include other forcesdue to the friction between wheels and ground and othermechanical resistance. Therefore, assumptions need to bemade to extract the aerodynamic drag from the total resist-ance (Grappe et al. 1997). Similar assumptions are madefor ground vehicle testing that relies upon constant speedtorque measurement (Fontaras et al. 2014) and coast downtests (Howell et al. 2002) of cars and trucks. Instrumenta-tion of a vehicle with pressure taps to measure the overallpressure distribution requires an unrealistic multitude ofsensing points. For an athlete, such approach is practicallyimpossible Aerodynamic drag measurements based on wake survey are potentially of more general applicability as theydo not dependent upon assumptions as discussed aboveand do not require extensive instrumentation of the model.The approach deals with the determination of the momen-tum deficit of the flow past the object. The airflow veloc-ity is measured in the wake and the aerodynamic force isobtained from the velocity deficit compared to the incom-ing stream. This approach has been long practiced in windtunnels by means of wake rakes (multi-hole Pitot probes;Goett 1939; Guglielmo and Selig 1996). In contrast tousing force balances and pressure taps, a wake rake offersthe advantage of yielding additional whole-field velocityinformation in the model’s wake. The identification of thestreamwise vortical structures behind cars, for instance, hasbeen among the most significant insights in vehicle aerody-namics in the past decades (Hucho and Sovran 1993). Fur-thermore, relating changes in the leg orientation of a cyclistto the wake flow topology recently provided new insightsfor cycling aerodynamics (e.g.Crouch et al. 2014). Particle image velocimetry (PIV) has been demonstratedas a valid instrument to replace Pitot probes for a stationaryobject (Kurtulus et al. 2007; Van Oudheusden 2013) as wellas for moving objects like rotor blades (Ragni et al. 2011,2012). More recently, Neeteson et al. (2016) have extendedthe approach to estimate the drag of a sphere freely fallingin water trying to reconstruct the pressure distribution allover its surface. The present study aim at measuring the aerodynamicdrag of transiting objects using tomographic PIV meas-urements in the wake of the model and invoking the conservation of momentum in a control volume. Becausethe use of micro-sized droplets as tracers is limited to arelatively small measurement domain (Scarano 2013),using helium-filled soap bubbles (HFSB) is consideredessential for the current experiment. The high light scat-tering efficiency and tracing fidelity of the HFSB (Bos-bach et al. 2009; Scarano et al. 2015, among others) allowmeasurements on the meter scale (Kiihn et al. 2011),preluding upscaling applications towards field tests in theautomotive and speed sports. The detailed goal of the present work is to examinethe accuracy of a method that measures the aerodynamicdrag of a transiting object and to assess its potentialapplicability for full-scale conditions. For this demon-stration, a sphere is towed within a rectangular channel ata velocity of 1.45 m/s (Re=10,000). Despite the simplegeometry, the flow exhibits an unsteady, turbulent wakewith complex vortex interactions (e.g. Achenbach 1972;Brucker 2001) mimicking conditions also encounteredbehind ground vehicles (Ahmed et al. 1984) or cyclists(Crouch et al. 2014). The experiments determine thetime-average drag and a detailed comparison is madewith data available from literature. The terms contribut-ing to the overall drag are studied separately to identifya criterion for simplified measurement configurations.Finally, the experimental uncertainty related to the tomo-graphic PIV measurements and disturbing factors suchas non-homogeneous initial flow conditions, supportingstrut and blockage effects are discussed. 2 Working principle 2.1 Drag from a control volume approach The drag force acting on a body moving in a fluid canbe derived by application of the conservation of momen-tum in a control volume containing this body (Anderson Fig.1 Schematic description of the control volume approach withina wind-tunnel setting (stationary object) 1991), as visualized in Fig. 1. In the incompressible flowregime, the time dependent drag force acting on the bodycan be written as: where v is the velocity vector with components [u,y,w]along the coordinate directions [x,y,z] respectively, p is thefluid density, p the static pressure and t the viscous stresstensor. V is the control volume, with S its boundary and nis the outward pointing normal vector. It can be shown that in most cases of interest the vis-cous stress is negligible with respect to the other contributions when the control surface is sufficiently far away fromthe body surface and the boundary is not aligned with theflow shear (Kurtulus et al. 2007). Furthermore, the contourS is defined as the contour abcd, with segments ad and bcapproximating streamlines. When the segments ab, ad andbc are taken sufficiently far away from the model, expression(1) can be rewritten such that the only surface integral to beevaluated is that in the wake of the model Swake (segment cd): where Uand p are, the freestream velocity and pressure,respectively. 2.2 Time-average force on a stationary model The evaluation of the volume integral on the right hand sideof Eq. (2) poses the typical problems due to limited opticalaccess all around the object. Evaluating this integral can beavoided by considering the time-average drag instead of itsinstantaneous value. When decomposing the equation intothe Reynolds average components and averaging both sidesof the equation, the time-average drag force is obtainedwith the sole contribution of surface integrals: (3) where u is the time-average streamwise velocity and u’ thefluctuating streamwise velocity 2.3 Time-average forces on a moving model at constantvelocity According to the principle of Galilean invariance, Eq. (3)holds in any reference frame moving at constant velocity. Following Ragni et al. (2011), the time-average drag actingon a model at constant speed UM is expressed as: where u is the streamwise velocity measured in the lab-oratory frame of reference. In the moving frame of refer-ence the freestream velocity stems from the velocity of themodel relative to that of its environment Ueny: For perfectly quiescent air Ueny=0.During real experi-ments it may occur that Ueny is different from zero and can-not be neglected. The value of Ueny needs therefore to bemeasured prior to the passage of the model. Expression (4) allows deriving the time-average dragfrom velocity and pressure statistics in a cross section ofthe wake. The time-average pressure is evaluated from thevelocity measurements solving the Poisson equation forpressure, according to Van Oudheusden (2013). An accu-rate evaluation of the pressure in a highly three-dimensionalflow requires the estimation of the full velocity gradienttensor components (Ghaemi et al. 2012; Van Oudheusden2013), which justifies the adoption of the tomographic PIVtechnique instead of planar stereo-PIV. The measurementvolume must be large enough to encompass the full wakeand reach the region of steady potential flow at its edgeswhere Dirichlet boundary conditions can be applied basedon Bernoulli law. Neumann conditions are prescribed at theinflow and outflow boundaries of the domain. The result-ing approach yields the time averaged drag using only thevelocity measurements and no other information. 3 Experimental apparatus and procedure 3.1 Measurement system and conditions A schematic view of the system devised for the experi-ments is shown in Fig. 2, with a photograph of the setup inFig. 3. The apparatus consists of a 170 cm long duct witha squared cross section of 50×50 cm, where the spheremodel is towed. Part of the duct has transparent walls foroptical access (Fig. 3). The model is a smooth sphere of diameter D=10 cmtowed at a constant speed UM = 1.45 m/s. The model issupported by an aerodynamically shaped strut with 20 mmchord and 3 mm thickness. The strut is 20 cm long and isinstalled onto a carriage moving on a rail beneath the bot-tom wall of the duct. The carriage is pulled by a linen wire experimental setup front-view Fig. 3 Overview of the experimental system and measurement con-figuration connected to the shaft of a digitally controlled electricmotor (Maxon Motor RE35). Markers on the surface of thesphere track its position during the transit. The wake behind the sphere starting from rest needstime to reach a fully developed regime. For an impulsivelystarted cylinder at low Reynolds number (Re<100) it isreported that about 15 cylinder diameters (Coutanceau andBouard 1977b) are needed before the wake is developed.In the present work a conservative value is taken with themodel beginning 25 sphere diameters before the measure-ment region. 3.2 Tomographic system The time-resolved tomo-PIV measurements are conductedusing neutrally buoyant helium-filled soap bubbles (HFSB,300 um diameter) as tracer particles produced with an arrayof ten generators that yield a total of about 300,000 parti-cles per second. The air, helium and soap fluid flow rates are controlled by a fluid supply unit provided by LaVisionGmbH. The average time response of such tracer particlesis expected to be below 20 s (Scarano et al. 2015). Con-sidering the relevant flow time scale (D/UM=70 ms) thesmall value of the Stokes number of the tracers indicatestheir adequacy for the current experiments. The illumination is provided by a Quantronix DarwinDuo Nd:YLF laser (2×25 mJ/pulse at 1 kHz). The laserbeam is first shaped into an elliptical cross section and isthen cut into a rectangular one with light stops (Fig.3). Thesize of the measurement volume is 3 cm × 40 cm × 40 cmin x, y and zdirection, respectively, and is located 80 cmfrom the exit of the duct (Fig. 2). The tomographic imag-ing system consists of three Photron Fast CAM SAl cam-eras (CMOS, 1024×1024 pixels, pixel pitch of 20 um, 12bits). Each camera is equipped with a 60 mm Nikkor objec-tive set to f/8. The optical magnification is approximately0.07. In the present conditions, the seeding concentrationis approximately 3 particles/cm’and the imaging densityis 0.04 particles/pixel. PIV acquisition is performed withinLaVision Davis 8.3 delivering single-exposure frames at arate of 500 Hz. 3.3 Measurement procedure The tunnel entrance and exit are closed to confine theHFSB seeding before the transit of the sphere. The exitis closed by a porous curtain, which maintains the seed-ing tracers inside, but prevents the buildup of an over-pressure. The HFSB generators are operated for approxi-mately two minutes until the steady-state concentrationiSIreached. Approximately quiescent conditionsareachieved 30 s after the generators are switched off.The tunnel entrance wall is then opened and the modelis put in motion through the duct. Image recording Fig. 4 Particle streaks for three positions of the sphere passingthrough the measurement domain. The particle streaks are obtainedby averaging over ten consecutive frames.The sphere positions areseparated by time increments of 70 ms begins 1 s before the sphere passes through the meas-urement region and stops when it touches the exit cur-tain. Figure 4 illustrates the motion of the air by particlestreaks at three positions of the model; before entering(top), inside (middle) and after leaving the measure-ment domain (bottom). The measurements are separatedby 70 ms (35 frames) and the streaklines are obtainedby averaging ten consecutive frames. A relatively highvelocity is observed in the center of the measurementregion behind the model (bottom image) with quiescentair conditions at the edges. An animation of raw imagesof the passing sphere is available at the multimedia storeof the publisher (Supplementary material 1). The experi-ment comprises 35 repeated measurements to form a sta-tistical estimate of the flow properties and the associatedaerodynamic drag. 3.4 Data reduction The tomographic-PIV data analysis is performed with theLaVision Davis 8.3 software. Image pre-processing com-prises background subtraction and Gaussian smoothing. Fig.5 Illumination distribution along the measurement depth (evalu-ated over an area of 25×25cm) The volume reconstruction and velocity evaluation fol-lows the sequential MTE-MART algorithm (SMTE, Lynchand Scarano 2015), yielding a discretized object with1074×1050×72 voxels. The interrogation is based on spa-tial cross-correlation volumes of 32’ voxels with an over-lap of 75%. The resulting velocity vector field has a den-sity of 3 vectors/cm. Figure 5 illustrates the reconstructedintensity distribution along the depth. The SMTE algorithmreturns a high reconstruction signal-to-noise ratio, indicat-ing that the cross correlation result is not affected by ghostparticles effects (Elsinga et al. 2006). A Galilean transformation of the instantaneous velocityis performed to represent the measurement in a frame ofreference consistent with the object moving at a constantvelocity U M. In this frame of reference, the drag is evalu-ated using Eq. (4). The phase-average velocity field in thelaboratory frame of reference, v* is obtained from the 35repeated measurements yielding the mean velocity andits fluctuations. A three-dimensional spatial representa-tion of the velocity field in the sphere frame of reference,y is obtained by considering the linear relation between thestreamwise coordinate and the time elapsed after the pas-sage of the sphere: where its origin matches the center of the sphere and bothcoordinate systems coinciding at t=0. The procedureencompasses the streamwise direction downstream of thesphere from x/D=0.5 till x/D=9.5. This approach is illus-trated in Fig. 6. reconstructed velocity field (v) Fig. 6 Illustration of the time-average streamlines after streamwise reconstruction of the velocity field (v) (in the sphere frame of reference)from phase-average time-history velocity fields (v") (in the laboratory reference frame) The Poisson equation for pressure is solved prescrib-ing Neumann conditions on all boundaries. The obtainedpressure field is scaled afterwards by a constant value,prescribing the freestream pressure at its top, bottom andside boundaries in the spatial range in which the veloc-ity at these boundaries best matches stagnant conditions(x/D>4), and excluding a small segment (5 cm wide) inthe center of the bottom boundary, which is affected by thewake of the strut. 4 Results 4.1 Instantaneous flow field At a Reynolds number of 10,000 the wake of a sphere is inthe unsteady regime, exhibiting vortex shedding and com-plex vortex interactions (Bakic et al. 2006). A snapshot ofthe flow structure will typically yield an asymmetric pat-tern, while the time-average structure is known to be axi-symmetric. Figure 7 shows the instantaneous velocity fieldin the laboratory frame of reference in the center YZ-planeat four consecutive time instants. A supplementary anima-tion of the time-resolved velocity field is available as mul-timedia file (Supplementary material 2). Non-dimensionaltime is defined as t"=tUM/D. Each increment △t*=1 cor-responds to a translation in space of one sphere diameterin negative x-direction. At t *=0.5, a region of acceler-ated flow is visible at the periphery of the wake. Further-more, a negative peak of streamwise velocity is present inthe near wake of the sphere. The maximum velocity deficitdecays with time, consistently with the observations frompast investigations (Jang and Lee 2008; Constantinescu andSquires 2003). 4.2 Time-average flow structure The ensemble-statistics yield the time-average velocityfield, the fluctuating velocity and time-average pressuredistribution. These terms are inspected to understand howthe individual terms from Eq. 4 contribute to the aerody-namic drag. Figure 8 illustrates the streamwise velocity dis-tribution in the separated wake (x/D=0.85, top) and afterthe flow reattachment (x/D=3, bottom). The velocity fieldin the wake is close to the axi-symmetric condition, withsome slight deviations due to the supporting strut (5-10%velocity deficit). The latter will be accounted for in the sec-tion on drag derivation. Furthermore, the spatial velocitydistribution shows a radial velocity directed towards theflow symmetry axis, decreasing in magnitude at increasingdistance from the sphere, which is consistent with literature(e.g. Jang and Lee 2008). The expected flow reversal in thecenter of the wake is also captured in the present measure-ment (Fig. 8-top-right). The streamwise velocity contour at x/D=3 (Fig. 8 bot-tom-right) shows a slight asymmetry in the spatial veloc-ity distribution outside the wake. At the top of the domainthe non-dimensional streamwise velocity is about 0.98,while at the bottom it is 1.01. This asymmetry stems fromthe flow conditions prior to the transit of the sphere andis ascribed to the motion induced during injection of theHFSB tracers. In the derivation of the aerodynamic drag,the momentum term expresses a deficit in the wake, relativeto the fluid momentum prior to the passage of the sphere(Eq. 5). Therefore any residual motion before the passageof the sphere is accounted for the drag evaluation. The streamwise velocity distributioninintthecen-tral XY-plane is depicted in Fig. 9. in the spatial range0.5
还剩12页未读,是否继续阅读?
继续免费阅读全文产品配置单
北京欧兰科技发展有限公司为您提供《平移运动球体中空气动力学阻力,3D3C速度场检测方案(粒子图像测速)》,该方案主要用于其他中空气动力学阻力,3D3C速度场检测,参考标准《暂无》,《平移运动球体中空气动力学阻力,3D3C速度场检测方案(粒子图像测速)》用到的仪器有体视层析粒子成像测速系统(Tomo-PIV)、LaVision DaVis 智能成像软件平台。
我要纠错相关方案