亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        Numerical Simulation on Flow Control for Drag Reduction of Revolution Body Using Dimpled Surface

        2011-03-09 11:57:32WANGJing王晶ZHANGChengchun張成春RENLuquan任露泉HANZhiwu韓志武
        Defence Technology 2011年1期

        WANG Jing(王晶),ZHANG Cheng-chun(張成春),REN Lu-quan(任露泉),HAN Zhi-wu(韓志武)

        (1.College of Agriculture,Jilin University,Changchun 130062 Jilin,China;2.MOE Key Laboratory of Bionics Engineering,Jilin University,Changchun 130022 Jilin,China)

        Introduction

        The drag reduction effect and methods,especially the flow control using non-smooth surfaces have been paid much more attention in recent years.Walsh[1]found that riblet surfaces can develop a lower shear stress than that of smooth surfaces on the basis of fluiddynamic reasoning.Similar studies performed by Reif in Germany were motivated by observation on shark skin.Subsequently,many researchers studied the various riblet surfaces,including a shark skin replica,and confirmed that the riblet surfaces can reduce the skin friction drag[2-7].In recent years,many researchers have paid more attention to engineering of non-smooth surface from the perspective of bionics.A typical application is the swimming suit based on the hydrodynamics of a shark skin.The investigators realize that the dimpled surface can also have the fluid dynamic effects.Bearman and Harvey[8]found that,for Reynolds number of 4×104to 3×105,the dimpled circular cylinder has a lower drag coefficient than a smooth cylinder.It was measured that,at transonic speed,the revolution bodies with dimples can reduce the total drag by 3%[9].

        In our past researches,we examined the drag reduction effect of bionic dimpled surface by wind tunnel tests.However,it is difficult to obtain the mechanism of drag reduction of the dimpled surface by experiments.Therefore,now,we focus our attention on the numerical analyses for the external flow of smooth and dimpled surfaces to obtain the drag reduction mechanism of the dimpled surface.

        1 Modeling Approach

        1.1 Governing Equation

        The governing equation can be writtenas where u is the velocity vector,and the components of u in x,y and z are u,v and w;φ is common variable;Γφis general diffusion coefficient;Sφis general source item.

        Table 1 shows the concrete terms of Eq.(1)for a compressible steady state flow.

        Table 1 Concrete terms of governing equation of flow

        1.2 Turbulence Modeling

        Some turbulence models are now widely used in engineering application,as they offer a good compromise between numerical effort and computational accuracy.SST(shear-stress transport)k-ω model developed by Menter[9]can effectively blend the robust and accurate formulation of the k-ω model in the near-wall region with the free-stream independence of the k-ε model in the far field was used in these fully turbulent simulations.The transport equations governing k and ω take the following form without regarding to the buoyancy

        where τijis the turbulent shear stress.The production term of τijin Eq.(2)and(3)is

        The blending functions F1and F2are defined as

        where

        The turbulent eddy viscosity can be calculated by using following formula

        The model constants can be calculated by using F1blending function,

        where φ1represents a generic constant in the k-ω equations and φ2the same constant in the k-ε equations.

        The k-ω model constants are given by

        The value of k-ε model constants are

        2 Computational Details

        2.1 Computational Case

        The revolution body with length of 517.5 mm and maximum diameter of 62.5 mm is shown in Fig.1.Fig.2 shows the position and dimensions of the dimples.The radius of the dimple is 1 mm;the depth of the dimple is 1 mm;the space between two adjacent dimples in axial direction is 4 mm.

        Fig.1 Main dimensions of model of revolution body

        2.2 Grids and Boundary Conditions

        The computational region is meshed by unstructured hybrid grids with tetrahedral elements and triangular prism elements.The distance of the first layer near the wall is determined according to

        Fig.2 Dimensional parameters of dimpled surface

        where uτis the wall friction velocity given by

        By solving Eq.(12),we can obtain

        where y+is the dimensionless distance of the first layer to the wall.The appropriate range of y+suiting SST kω turbulence model is y+<5,in these simulations,y+=4.U∞is the flow velocity calculated by using

        where Ma is the Mach number of flow,k the specific heat ratio of air,R the gas constant,and T the absolute temperature.k and R are 1.4 and 287 m2/s2K respectively.

        The average friction coefficientcan be approximately calculated by using

        where ReLis the Reynolds number based on the length of revolution body,which can be defined as

        where L is the length of the revolution body,and v the kinetic viscosity coefficient defined as

        where T=299 K,ρ=1.185 kg/m3,C=110.4 K.

        Thus,the distance between the first layer grids and the wall is Δy≈0.01 mm,and the boundary layer thickness δ is

        Figure 3 shows the unstructured Cartesian hybrid grids with prism and tetrahedral elements.There are 10 layers of prisms over the surface of the revolution body.The dimpled region is provided with high-density meshes to observe the variation of the vortexes forming in the dimples.The computational region has 1.706 ×106cells composed of tetrahedral and prisms approximately.

        Fig.3 Grids in boundary layer dimpled region

        As an axially symmetric body,the revolution body can be simulated in its quarter to reduce the amount of calculation.Thus,the computational region can be defined as a quarter of a cylinder.It is 1/4 of 20Rm×10L,where Rmis the maximum radius of the revolution body,L the length.

        The pressure far-field boundary condition is used to model the compressible free-stream at Mach number of 0.4,the static pressure of 100 750.4 Pa,the temperature of 299 K.The angle of attack is α =0°.

        3 Calculation Results and Analyses

        3.1 Calculation Results

        The drag reduction ratio of the dimpled surface can be defined as

        where Csmoothis the drag coefficient of the smooth revolution body,and Cbionicthe counterpart of the dimpled one.According to Eq.(20),the drag reduction ratios of the viscous drag,the pressure drag and the total drag are all obtained.The viscous drag of the revolution body decreases by 8.05%,the pressure drag by 1.9%and the total drag by 6.24%.

        3.2 Mechanism of Friction Drag Reduction

        The viscous fiction stress of revolution body includes viscous shear stress and turbulent Reynolds stress and can be defined as

        where μ is the dynamic viscosity coefficient,μtthe turbulent viscosity coefficient.The relationship of the instantaneous velocity vx,the time averaged velocityvxand the fluctuation velocity v'is vx=vx+v'.The friction drag reduction mechanism can be explained on the basis of Eq.(21).

        Figure 4 shows the wall shear stresses of the dimpled and smooth revolution bodies.Evidently,the wall shear stress of the downstream of the dimples is reduced observably.Fig.5 compares the velocity to that of smooth body in the same location near the wall of the dimpled body,and it shows the wall shear stress reduction of the dimpled revolution body.It can be seen from Fig.5 that the flow speed in the dimples and downstream is visibly lower than that of the smooth body.In Fig.5(b),δ+is defined as

        where Δy is the distance from monitor point to the wall of the model,δ is the boundary layer thickness.Form Fig.5(b),we can deduce that the dimples reduce the flow field velocity gradient near the dimpled surface,and increase the viscous sub-layer of the boundary layer.Therefore,the viscous shear stress decreases.

        According to the physical factor of Reynolds stress,we know that the speed fluctuation can reflect the scale of Reynolds stress.Moreover,the ratio of its square root to the time averaged velocity is the turbulent intensity.Fig.6 compares the turbulent intensity and the turbulent kinetic energy of the dimpled body to the smooth one.From Fig.6,we can draw that the Reynolds stress is reduced,because the turbulent intensity near the dimpled surface is lower than that of the smooth one.The turbulent kinetic energy of the dimpled model,which is consumed by the turbulent fluctuation,is lower than that of the smooth one.

        Fig.4 Wall shear stress comparison between dimpled and smooth revolution bodies

        Fig.5 Velocity comparison between dimpled and smooth revolution bodies

        In summary,the bionic dimpled surface reduces the viscous shear stress and the Reynolds stress simultaneously.Therefore,its friction drag is reduced evidently.

        Fig.6 Turbulent comparison between dimpled and smooth revolution bodies

        3.3 Mechanism of Pressure Drag Reduction

        The pressure drag of the revolution body includes the base drag and the shock wave drag.The base drag can beexpressed as

        where Sbis the base area,and SMis the maximum cross section.Fig.7 compares the base pressures of the dimpled revolution bode to that of the smooth one.It indicates that the base pressure near the bottom of the dimpled body is higher than that of the smooth one.Namely,the dimples reduce the pressure difference.Accordingly,the pressure drag is reduced.

        3.4 Analysis of Boundary Layer Control Using Dimpled Surface

        Figure 8 shows the flow pattern in the dimple and the velocity vector on the dimpled surface.Evidently,there are a low-speed rotating vortexes forming in the dimples.Therefore,it is deduced that the dimples produce two effects.First,the low speed rotating vortexes in the dimples can be entitled an air cushion effect.Second,it produces friction drag against to the other area and results in a driving effect.

        Fig.7 Base static pressure comparison between dimpled and smooth revolution bodies

        Fig.8 Flow pattern of interior dimples and velocity vector of dimpled surface

        Figure 9 compares the streamlines on the dimpled surface to that on the smooth one.Obviously,the streamlines of the dimpled revolution body are rarer than those of the smooth one.It indicates that the flow does not entry into the dimples,but passes over the air cushion and further explains the existence of the air cushion effect.

        Fig.9 Streamline comparison between dimpled and smooth revolution bodies

        4 Conclusions

        Following conclusions can be drawn from the numerical simulation with SST k-ω model.

        1)At the Mach number of 0.4,the dimpled surface arranged on the rearward of the revolution body reduces the viscous drag by 8.05%,the pressure by 1.9%and the total drag by 6.24%.

        2)The dimples simultaneously reduce the wall shear stress and the Reynolds stress,then,the skin friction drag decreases.The dimples can also weaken the pumping action on the dead water region behind the revolution body caused by the external flow,which is the reason why the pressure drag decreases.

        3)The basic reason of the skin friction drag reduction is that the dimples control and correct the boundary layer.The low speed rotating airflows in the dimples hinder the exterior airflow and the wall of the dimples to contact directly.It can be regarded as an air cushion effect.Secondly,the low speed rotating vortexes forming in the dimples can produce a friction drag against to the other area and results in a driving effect.

        [1]Walsh M J.Riblets as a viscous drag reduction technique[J].AIAA Journal,1983,21(4):485 -486.

        [2]Reif W E,Dinkelacher A.Hydrodynamics of the squamation in fast swimming sharks[J].Neues Jahrbuch für Geologieund Palaeontologie,(Abhandlungen),1982,164:184-187.

        [3]Bechert D W,Bartenwerfer M,Hoppe G,et al.Drag reduction mechanisms derived from shark skin[C].Proceedings of the 15th ICAS Congress,London,1986:1044-1068.

        [4]Choi K S.Near-wall structures of a turbulent boundary layer with riblets[J].J Fluid Mech,1989,208:417 -458.

        [5]Walsh M J,Sellers W L,McGinley C B.Riblet drag at flight conditions[J].Journal of Aircraft,1989,26(6):570-575.

        [6]Bechert D W,Bruse M,Hage W.Experiments with three dimensional riblets as an idealized model of shark skin[J].Exp Fluids,2000,28:403 -412.

        [7]Bechert D W,Bruse M,Hage W,et al.Fluid mechanics of biological surfaces and their technological application[J].Naturwissenschaften,2000,87:157-171.

        [8]Bearman P W,Harvey J K.Control of circular cylinder flow by the use of dimples[J].AIAA Journal,1993,31:1753-1756.

        [9]REN Lu-quan,ZHANG Cheng-chun,TIAN Li-mei.Experiment study on drag reduction for bodies of revolution using bionic non-smoothness[J].Journal of Jilin University(Eng and Techn Ed),2005,35(4):431-436.

        [10]Menter F R.Two-equation eddy-viscosity models for engineering applications[J].AIAA Journal,1994,32:1598-1605.

        亚洲国产精品中文字幕久久| 中文字幕一区二区三区四区久久| 亚洲va在线∨a天堂va欧美va| 日韩人妻精品无码一区二区三区| 无遮挡十八禁在线视频国产制服网站 | 亚洲中文字幕第一页免费| 美国黄色av一区二区| 日韩一区二区三区人妻免费观看| 久久伊人这里都是精品| 色老板美国在线观看| 国产特级毛片aaaaaa视频| 亚洲中文字幕久在线| 国产无遮挡又黄又爽又色| 就去吻亚洲精品欧美日韩在线| 亚洲AV一二三四区四色婷婷| 在线无码精品秘 在线观看| 中文字幕一区二区三区.| 人妻露脸国语对白字幕| 亚洲国产综合人成综合网站| 亚洲av高清一区二区三| 精品国产粉嫩内射白浆内射双马尾| 亚洲av精品一区二区三区| 国产国拍亚洲精品mv在线观看| 人人妻人人澡人人爽人人精品| 五月中文字幕| 免费国产调教视频在线观看| 亚洲女同av一区二区在线观看| 人妻少妇中文字幕在线观看| 精品成在人线av无码免费看| 麻麻张开腿让我爽了一夜| 国模精品无码一区二区二区| 区一区一日本高清视频在线观看| 免费av一区男人的天堂| 99精品国产成人一区二区 | 男女视频一区二区三区在线观看| 亚洲高清三区二区一区 | 免费av片在线观看网址| 国产无遮挡又爽又刺激的视频老师| 公粗挺进了我的密道在线播放贝壳| 妺妺窝人体色www聚色窝韩国| 日本一区二区三区资源视频|