• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Effect of an end plate on surface pressure distributions of two swept wings

    2017-11-17 08:31:43MohmmdRezSOLTANIMehrnMASDARIMohmmdRsoulTIRANDAZ
    CHINESE JOURNAL OF AERONAUTICS 2017年5期

    Mohmmd Rez SOLTANI,Mehrn MASDARI,Mohmmd Rsoul TIRANDAZ

    aDepartment of Aerospace Engineering,Sharif University of Technology,Tehran 19166,Iran

    bFaculty of New Sciences and Technologies,University of Tehran,Tehran,Iran

    Effect of an end plate on surface pressure distributions of two swept wings

    Mohammad Reza SOLTANIa,Mehran MASDARIb,*,Mohammad Rasoul TIRANDAZa

    aDepartment of Aerospace Engineering,Sharif University of Technology,Tehran 19166,Iran

    bFaculty of New Sciences and Technologies,University of Tehran,Tehran,Iran

    A series of wind tunnel tests was conducted to examine how an end plate affects the pressure distributions of two wings with leading edge(LE)sweep angles of 23°and 40°.All the experiments were carried out at a midchord Reynolds number of 8×105,covering an angle of attack(AOA)rangefrom-2°to 14°.Static pressure distribution measurements were acquired over the upper surfaces of the wings along three chordwise rows and one spanwise direction at the wing quarter-chord line.The results of the tests confirm that at a particular AOA,increasing the sweep angle causes a noticeable decrease in the upper-surface suction pressure.Furthermore,as the sweep angle increases,the development of a laminar separation bubble near the LEs of the wings takes place at higher AOAs.On the other hand,spanwise pressure measurements show that increasing the wing sweep angle results in forming a stronger vortex on the quarter-chord line which has lower sensitivity to AOA variation and remains substantially attached to the wing surfacefor higher AOAs than that can be achieved in the case of a lower sweep angle.In addition,data obtained indicate that installing an end platefurthe r reinforces the spanwiseflow over the wing surface,thus affecting the pressure distribution.

    1.Introduction

    Calculation of pressure distribution over airfoils operating in inviscid flow,is one of the classic problems of the fluid motion the ory.Consequently,it has regained significance in designing swept wing aircraft.

    The reasons for the application of swept wings to highsubsonic-speed aircraft are now well-established.Namely,a swept wing has a greater critical Mach number than that of a corresponding unswept wing.In addition,the variation of pressure coefficient with Mach number for a given swept wing is smaller compared to that of an unswept wing with the same airfoil shape.Moreover,on a swept wing,the critical value of pressure coefficient is reached at a higher freestream Mach number.Furthermore,at low Mach numbers,the perturbation velocities are small compared with those of an unswept wing.

    The tendency of a sweptback wing to tip stall is due to the induced spanwiseflow of the boundary layer from root to tip.Since the tip of a sweptback wing is located aft part of the wing(behind the center of lift),it has the effect of moving the center of lift forward,causing a sharp pitch-up movement and thus making the stall worse than it was.The tendency to tip stall is the greatest when wing sweep and taper are combined.

    There are many studies containing experimental pressure distribution measurements.For example,in 1952,a series of low-speed wind tunnel tests was conducted to determine boundary layer and surface pressure distribution characteristics of a family of sweptback wings.As a result,it was found that the experimental data would have had application in giving a better understanding of the viscous flow phenomena on sweptback wings,particularly in relation to the stall phenomena.1Likewise,in 1964,an experimental study was conducted on sweptback wings at the Ames Research Center to investigate the stall characteristics of swept wings.2In addition,in 1990,a comprehensive experimental study at the NASA Langley Research Center was conducted on F-14A variable swept wings modified for laminar flow,to assess the transition location under various conditions.Accordingly,left-and rightwing pressure distribution measurements at various Mach numbers and two different altitudes were used to give an explanation for the increased response amplitude of the right wing.3Furthermore,in 1997,a detailed investigation was conducted on F-111 aircraft modified with a variable-sweep supercritical mission adaptive wing(MAW)at the Dryden Flight Research Center,to improve the aerodynamic characteristics.4Moreover,in 2011,a series of experimental studies was conducted to examine the flow field over a swept wing under various conditions at Shar if University of Technology.Subsequently,measurements of pressure distribution were used to predict the transition point at each chordwise section.5Furthe r,surface pressure distribution over a swept wing was exerted by Yen and Hsu6to investigate the effects of angle of attack α and chord Reynolds number Recon vortex shedding and aerodynamic coefficients.

    On all subsonic wings the re is a tendency for a secondary flow to develop from the high-pressure region below the wing round the wingtip to the relatively low-pressure region on the upper surface.7At a positive lift condition,the flow close to the tip of a three-dimensional wing tends to curl around the tip,because the flow is forced from the lower surface region(high-pressure side)to the upper surface region(low-pressure side)in a circular fashion(Fig.1).This tendency for the flow to leak around the tip causes a circular three-dimensional flow pattern named vortex,which features a low-pressure core.Viewed from behind,the port wing tip vortex rotates clockwise and the starboard vortex rotates counterclockwise.Aircraft that encounter the wing tip vortices of large preceding aircraft can experience severe upward or downward loads as well as an overpowering rolling moment.

    Fig.1 Wing tip vortex.

    The research on the wing tip vortices began in the fifties and it is still an ongoing study because ofits wide range of applications and because of the big number of fluid dynamic features that it implies:vorticity sheet rolling up,multiple vortex system,high turbulence level and relaminarization,diffusion,decay,collapse,instability and merging of vortices are just a few numbers of the aspects that coexist in this phenomena.8

    Significant efforts have been utilized to investigate the formation and decay process of a wing tip vortex.Some of those efforts have been summarized by Spalart.9Winglets,end plates,tip sails,wing-grids,and blowing jets are typical examples of mechanisms for wing tip vortex control.They are effective in decreasing the strength of the trailing wing tip vortex and reducing the induced drag.10

    The effect of eliminating wing tip vortices on aerodynamic characteristics has not been fully investigated until now.Only few experimental investigations on the effect of an end plate were conducted by Fink and Lastinger11and Carter.12Fink and Lastinger performed a series of wind-tunnel investigations to determine the ground effect on the aerodynamic characteristics of thick and highly cambered wings with various aspect ratios.Furthe r,in 2008,a numerical analysis was performed to investigate the aerodynamic characteristics and static height stability of an end plate on an aspect-ratio-one wing in ground effect.13

    End plates are commonly installed on aerodynamic models to mitigate three-dimensional end effects.Furthermore,experimental studies on flow development over non-circular geometries have similarly found that improved mean spanwise uniformity can be achieved by installing end plates on the se models.14Moreover,end plates have been used in a number of experimental studies on airfoils operating at low Reynolds numbers.15–18

    It should furthe r be noted that a preceding work done by Soltani et al.19has successfully investigated the influence of Reynolds number on the surface pressure distribution over the same wing models without an end plate.In addition,Pelletier and Mueller18investigated the effects of end plates on force balance measurements for a two-dimensional Eppler 61 airfoil model at chord Reynolds numbers between 40×103and 100×103.As a result,the y demonstrated that the presence of end plates has a strong effect on lift and drag coefficients at low Reynolds numbers.

    The purpose of this investigation is to assess the effects of end plates on the surface pressure distributions of two sweptback wings.Static pressure distributions along three chordwise rows and similar pressure measurements along one spanwise direction are obtained at a midchord Reynolds number of 0.8×105and at angles of attack(AOAs)ranging from α =-2°to α =14°.

    2.Experimental apparatus

    All experiments were conducted in a closed-loop subsonic wind tunnel with a test section size of 80 cm×80 cm×200 cm in height,width,and depth,respectively.The tunnel can be operated over a speed range of 10–100 m/s.Thefree-stream velocity V∞was measured using a pitot-static tube.Flow conditioning was accomplished by a section offour large,anti-turbulence screens and aluminum honeycomb installed upstream at a 7:1 contraction section,producing a freestream turbulence intensity of less than 0.1%.In addition,the non-uniformity of the average velocity across the cross-section was less than 0.5 percent.

    Fig.2 Wing model and test setup schematic,and locations of rows and surface pressure tabs at the midchord.

    In this investigation,two semi-span 1:2.5 scaled models,which are made of aluminum alloy,were used as test cases for which the effects of experimental setup on surface pressure distribution were evaluated.The wing models,which have airfoil sections similar to that of NACA 6-series,have a leading edge sweep Λ of 23°and 40°.Furthermore,the root chord length Cris 394 mm,the tip chord length Ctis 172 mm,and the wing semi span(b/2)is 516 mm,which yield a semi-span aspect ratio of 0.43(Fig.2).A polished aluminum alloy plate with a rounded leading edge(LE)and trailing edge(TE)was installed over the test section floor to decrease the boundary layer effect of the test section on the model.In addition,to eliminate the formation of wing tip vortices,a Plexiglas flat plate with a rounded LE and TE and a size of 100 cm×100 cm×1 cm in width,length,and height respectively,was installed on the tip of a wing.The wing models were mounted on a support and the n inserted through a boundarylayer thickness controlling plate.In Fig.2(a),the general arrangement of each model used for this investigation when installed in the wind tunnel is illustrated.

    In this study,a pressure instrumentation was located on the right-wing model of the real airplane.A total of ninety static pressure taps of 0.4 mm in diameter were distributed along the model surface in three chordwise rows,i.e.,the inboard station of the wing(Section 1,2y/b=0.2),the middle station(Section 2,2y/b=0.44),and the outboard station(Section 3,2y/b=0.78),as well as in one spanwise row along the quarterchord line(Fig.2(b)).

    The pressure taps distributed on the upper surface of the wing were connected by stainless steel tubing to accurate pressure transducers.In most cases,the lengths of the pressure lines from the orifices to the transducers were limited to less than 1.5 m,thus minimizing pressure lag effects.Furthermore,each transducer’s data were collected via a terminal board and transformed to a computer through a 64-channel,12-bit analog-to-digital(A/D)board capable of an acquisition rate of up to 500 kHz.The raw data were the n digitally filtered using a low-pass filtering routine.During the filtering process,the cutofffrequency was calculated from eithe r the power spectrum estimation orfrequency-domain analysis.The advantage of this method is that the noise,which may dominate the signal in the time domain,appears only as a single peak or spike in the frequency-domain.20

    3.Presentation of results

    Static pressure distributions over the upper surfaces of the wings,in both chordwise and spanwise directions,were measured at different angles of attack.Chordwise pressure data will be presented first,followed by those of spanwise measurements,and finally a comparison between the measurements obtained in this experiment(wing with an end plate,WWE)and the results of the case without an end plate(WOE)will be given.Each figure of the pressure plots is associated with a sketch to illustrate the sweep angles of the wings and the relative locations of the rows of pressure taps.

    3.1.Chordwise pressure distribution results

    3.1.1.Effect of angle of attack on chordwise Cp

    Middle-section static pressure coefficient Cpdistributions for angles of attack from-2°to 6°are presented in Fig.3 for both wing models.As illustrated in this figure,at low AOAs,-2°≤ α ≤ 2°,the wing section airfoil has a thickness distribution that leads to a favorable pressure gradient over a large part of the airfoil.Furthermore,over a significant part of the wing chordwise Section,the flow is attached to the surface.It should furthe r be noted that in the region over the wing sur-face where dCp/dx<0,called favorable pressure gradient region,the value of Cpincreases negatively,creating a stronger suction force over the wing surface;conversely in the region over the wing surface where dCp/dx>0,called adverse pressure gradient,the magnitude of Cpdecreases,but the value of Cpincreases positively.

    Fig.3 Effect of AOA on Cpover the middle section,for-2° ≤ α ≤ 6°.

    As the AOA increases,a suction peak in the chordwise pressure distribution can be seen that moves upstream(Fig.3).For example,on Section 2 of the 40°-sweep(40DS)wing,at an angle of attack of 0°,the pressure peak is detected at x/C?0.47,C is airfoil chord,but as the angle of attack increases to 2°and 4°,the point of Cp,minmoves to x/C?0.44 and x/C?0.03,respectively(Fig.3(b)).The increment of the suction peak,which is an indication ofincreasing the speed of the air moving past the wing and generating stronger suction pressure on its surface,is clearly visible when increasing the angle of attack(Fig.3).

    The results show that when the 23°-sweep(23DS)wing is at angles of attack greater than 4°,the adverse pressure gradient aft of the point of the suction peak on the middle station becomes stronger(Figs.3(a)and 4).A high suction peak near the wing leading edge indicates that the flow becomes turbulent close to the leading edge where x/C>0.Static surface pressure measurements demonstrate that for angles of attack of 6°and 8°,the flow will definitely become turbulent in a region of 0≤x/C≤0.17.On the contrary,the flow over the wing surface of this station(Section 2)seems to be completely attached for an angle of attack range of-2°≤ α ≤ 8°.It should be pointed out that due to practical reasons,no pressure tap was installed at x/C>0.8;the refore,results could not be obtained close to the wing trailing edge(Fig.2(b)).

    Fig.4 Effect of AOA on Cpover the middle section,for Λ =23°and 8° ≤ α ≤ 14°.

    With an increasing angle of attack to 10°,the pressure distribution over the surface of the 23DS wing represents a significant decrease in the magnitude of Cp,minand a slight peak at the LE,x/C=0,followed by a region of relatively constant pressurefrom x/C=0.04 aft to x/C=0.1,and a lack of complete pressure recovery near the TE of the wing(Fig.4).The plateau of the pressure distribution is a characteristic of a laminar separation bubble.21

    For Reynolds numbers from about 0.8×106to 1.2×106,laminar separation on an airfoil occurs when an adverse pressure gradient of sufficient strength is imposed on an approaching laminar boundary layer and the flow near the surface is forced to zero velocity at the point of separation.When sufficient pressure recovery is achieved with some distance downstream,the flow reattaches,forming a separation bubble.Furthermore,flow separation leads to the formation of a shear layer,which divides the recirculating flow within the bubble and the external freestream flow.Before the laminar shear layer transits to turbulent flow,the static pressure in the bubble remains fairly constant,and the speed of the reverseflow is very slow.The beginning and end of this pressure plateau,which is commonly termed as the ‘dead-air” region,mark the locations of laminar boundary layer separation and separated shear layer transition,respectively.Entrainment of high-energy external flow by the turbulent shear layer allows for adequate pressure recovery and results in flow reattachment.22

    Considering the data obtained in this investigation,local separation of the flow first appears near the wing LE,followed by reattachment of the flow farthe r aft and the n recovery to essentially freestream pressure.From Fig.4,it can be seen that with an increasing angle of attack from 8°to 14°,the suction peak value continues to increasefrom-3.23 to-1.37, |Cp|increases near the TE from 0.22 to 0.67,the pressure plateau lengthe ns from 0.04≤x/C≤0.1 for α=10°to 0.04 ≤ x/C ≤ 0.27 for α =14°,and the point of reattachment moves downstream from x/C=0.1 to x/C=0.27.

    3.1.2.Effect of spanwise location on chordwise Cp

    Fig.5 shows a comparison between the measured pressure distributions at three spanwise locations of the 23DS wing for angles of attack of 0°,8°,10°,12°,and 14°.From the data obtained,a slight amount of longitudinal shift in the pressure peak can be detected,as the distancefrom the wing root to the wing tip is increased(Fig.5(a)).

    The pressure distribution measurements show a sudden drop in the magnitude of Cpnear the TE of Section 3 for the 23DS wing(Fig.5(a)).Moreover,the data obtained indicate that with an increasing angle of attack, |Cp|drops more sharply in this area of the wing surface(Fig.5(a)–(e).As will be discussed in Section 3.3 of this article,with an increasing AOA,a vortex-shaped flow forms on the surface of the end plate,to which an increasing angle of attack causes its size and strength to increase,thus affecting a larger area of the flow field near the tip of the wing surface and resulting in a stronger spanwise flow.As a result,this vortex affects the airstream velocity along the outboard section near the TE,and causes a sudden drop in the pressure distribution near the trailing edge of the outboard Section.

    Fig.5 Effect of spanwise location on chordwise pressure distribution for Λ =23°.

    From Fig.5(b),it can be noted that when the 23DS wing is at an 8°angle of attack,at a location downstream from the suction peak,the rate of pressure recovery on Section 3 decreases significantly,creating a pressure plateau in the pressure distribution.As it was previously mentioned,this plateau of the pressure distribution which coversa region of 0.04

    As the angle of attack increases to 1°,the length of the area of constant pressure on Section 2 increases in size and covers a region of 0.07≤x/C≤0.2.At the same time,increasing the minimum value of the outboard-Section pressure distribution,lack of pressure recovery near the trailing edge,and a constant Cpfrom x/C=0.05 to x/C=0.77 all indicate complete separation of the flow over the entire chord of Section 3(Fig.5(d)).From Fig.5(e),it can be observed that at a 14°angle of attack,the re exists a pressure plateau near the leading edge of Section 1 covering a region of 0.03≤x/C≤0.09.As noted before,this area of constant pressure is a characteristic of a laminar separation bubble which indicates local separation.As the distancefrom the wing root to the wing tip is increased,the length of the pressure plateau increases in size and covers a region of 0.04≤x/C≤0.27 on the mean aerodynamic chord C(Section 2);however,the re is some pressure recovery after this area of constant pressure,indicating a tendency toward reattachment of the flow.On Section 3,Cpis nearly constant in a region of 0.05≤ x/C≤ 0.77,indicating that the flow is separated in this region of the wing.

    It is thought that this laminar bubble type of separation leads to formation of a continuous vortex streak,which streams inboard and aft when increasing the angle of attack from 8°(initial appearance near the LE of Section 3)to 14°(initial appearance near the LE of Section 1).When the 23°-sweep wing is at a 14°angle of attack,a spreading-out of the separation streak over a larger chordwise distance can be noted as the distancefrom the origin of the streak is increased(Fig.5(e).It may be that the complete separation noted above is the result of this spreading-out process.On the other hand,with respect to,Ref.23because of the usual sectional lift coefficient distribution on a swept wing,the bubble appears first at the tip(for an angle of attack of 8°,Fig.5(b))and spreads slowly toward the root,as the wing angle of attack is increased(comparing Fig.5(b)–(e)).In this case,before the leading-edge bubble spreads to the root,the section near the wing tip(Section 3)has complete separation,and the bubble has curved back to leave the wing at the farthe st inboard point where separation has reached the trailing edge(Fig.5(e).As the angle of attack is increased,both the origin of the bubble and the point at which it leaves the wing move inboard(Fig.5(b)–(e)).

    3.1.3.Effect ofincreasing sweep angle on chordwise Cp

    Fig.6 illustrates the chordwise pressure distributions for the two wing models at various angles of attack.In Fig.6(a),the pressure distributions for the middle station(Section 2)are compared for the two sweep angles at an angle of attack of 2°.It can be seen that as the sweep angle increases,the value of pressure coefficient increases negatively,creating less suction over the wing surface.

    Results show that as the sweep angle increases,the formation of a laminar separation bubble takes place at higher angles of attack.Furthermore,with an increasing sweep angle,at a specified angle of attack,the length of the pressure plateau decreases.The results in Fig.6(b)show that at a 12°angle of attack,the re is no indication of local separation or sectional stall in the inboard section pressure distribution,for both sweep angles.As the angle of attack increases to 14°,a pressure plateau close to the LE of the inboard section of the 23DS wing can be seen, which covers a region of 0.03≤x/C≤0.09,while on the inboard station(Section 1)of the 40DS wing,the re is no sign of a pressure plateau(Fig.6(c)).As discussed before,this area of constant pressure is a characteristic of a laminar separation bubble and indicates local separation of the flow.

    Theformation of a laminar separation bubble near the middle-station LE of the 23DS wing,which covers a region of 0.04≤x/C≤0.1,takes place at an angle of attack of 10°.This constant pressure area is followed by flow reattachment farthe r aft and recovery to essentially free-stream pressure(Fig.6(d)).From Fig.6(d),it can be observed that the laminar separation bubble on the middle station of the 40DS wing is formed at an angle of attack of 10°,covering a region of 0≤x/C≤0.03.It should be noted that,however,the separation bubbles on the middle sections of the two models are detected at the same AOA,and it is seen that as the sweep angle increases,the separation bubble’s covering region decreases from 0.041≤x/C≤0.1 in the case of the 23DS wing to 0≤x/C≤0.03 in the case of the 40DS wing.Moreover,Fig.6(d)shows that with an increasing sweep angle,the slope of the pressure recovery aft of the separation bubble increases.Pressure measurements show that as the angle of attack increases,the lengths of the laminar separation bubbles on Section 2 of the models increase in size,but in the range of angle of attack under experiments,the re is no sign of completeflow separation over the entire chord of Section 2(Fig.6(e)and(f)).

    Fig.6 Effect of sweep angle on chordwise pressure distributions of wing models at various AOAs.

    The separation bubble in the outboard section pressure distribution of the 23DS wing is observed at an angle of attack of 8°(Fig.6(g)),which covers a region of 0.05

    It can be noted that at a 10°angle of attack,the re is an indication of separation over the entire chord length on the outboard station of the 23DS wing(Fig.6(h)),while in the outboard section pressure distribution of the 40DS wing,characteristics of a stalled airfoil can be seen at an AOA of 14°(Fig.6(i)).

    As illustrated in Fig.6,at a particular angle of attack,the lengths of the laminar separation bubbles over the two wings’chordwise sections decrease with an increasing sweep angle.Furthermore,increasing the sweep angle causes the stall phenomena on the chordwise sections to delay.It can be seen that at an angle of attack of 10°,the flow over the outboard section of the 23DS wing is completely separated,and characteristics of a stalled airfoil are clearly visible;neverthe less,over the outboard section of the 40DS wing,the pressure measurements show alocalized separation which coversa region of 0.05

    Fig.7 Spanwise pressure distribution along quarter-chord line.

    3.2.Spanwise pressure distribution results

    3.2.1.Effect of angle of attack on spanwise Cp

    Fig.7(a)and(b)shows the measured spanwise pressure distributions in terms of Cpversus 2y/b along the quarter-chord lines of the models at angles of attack ranging from-2°to 14°.As illustrated in Fig.7(a)and(b),for low angles of attack,the values of Cpslightly decrease along the span from the wing root to the wing tip,which indicate a stronger suction force dominating the outboard region of the wing.It is found that the magnitudes of pressure coefficients are the lowest near the wing root and increase outboard,which of course is a characteristic of swept wings.Furthermore,from the pressure results,it can be seen that with an increasing AOA from-2°to 4°,the slope of the pressure distribution increases,which suggests that the suction pressure near the wing tip is getting stronger more rapidly,compared to that of the inboard region of the wing.The rate of change of the Cpcurve slopefor AOAs ranging from-2°to 4°is presented in Fig.8 for both wing models.

    The results obtained suggest that at a specific AOA,a distortion in the pressure distribution near the wing tip area can be seen.Furthermore,with an increasing angle of attack,a large portion of the spanwise section near the wing tip is dominated by a very-low-pressure area.Generally,the distortion observed in this area,which is due to the formation of a vortex,and the low pressure associated with the centrifugal force caused by this vortex lead to the formation of a lowerpressure area on the wing surface(Fig.7(b)and(d)).

    Fig.8 Variation of spanwise pressure distribution curve slope.

    Spanwise pressure results show that this vortex due to the sweep angle,forms near the tip of the 23DS wing at an 8°AOA,covering an area of 0.54<2y/b<0.97(Fig.7(b)).As it is shown in Fig.7(b),the point of the minimum pressure,which shows the vortex core location,pushes inward from the wing tip to the wing root,increasing the covered area over the quarter-chord line,when increasing the AOA.At a 12°angle of attack,a slight drop in the minimum value of the pressure coefficient suggests that the core of the vortex just starts to get farthe r from the surface.In addition,it can be easily seen that with increasing the AOA to 14°,the strength of the vortex decreases significantly and the pressure peak drops from-1.21 to-1.19,due to a greater distance between the vortex core and the wing surface.Moreover,it is found that as the AOA increases from 8°to 10°,a sudden drop in the spanwise Cpnear the tip of the wing indicates that this region of the wing is experiencing an incipient stall due to the vortex breakdown.In addition,the results show that the separation area moves from the wing tip to the wing root,when increasing the angle of attack from 8°to 14°.

    In the case of the 40DS wing,the distortion observed in the pressure distribution,which leads to forming a low-pressure area near the wing tip,is seen at an angle of attack of 6°.Like the case of the 23DS wing,as the AOA gets higher,the vortex increases in size and gets stronger,furthe r reducing the pressure over the quarter-chord line near the wing tip(Fig.7(c)and(d)).Results show that in the range of angles of attack under experiments,the re is no sign of a spanwise-pressurepeak drop,which means that the core of the vortex is successively close to the 40DS wing surface;however,as it is illustrated in Fig.8,the curve slope ofincreasing the value of spanwise Cp,mindecreases,when increasing the AOA from 8°to 14°.Moreover,like the 23DS wing,with an increasing AOA from 10°to 12°,a sudden pressure drop near the tip of the 40DS wing shows an incipient stall region due to the vortex breakdown,which moves toward the wing root when increasing the angle of attack(Fig.7(d)).

    3.2.2.Effect ofincreasing sweep angle on spanwise Cp

    Fig.9 shows the effect ofincreasing the sweep angle on the spanwise pressure distribution at a 2°angle of attack.As it is shown,with an increasing sweep angle,a longitudinal shift in the pressure distribution curve can be observed,which represents a weaker suction force over the quarter-chord line of the 40DS wing.Furthermore,it is clearly seen that as the sweep angle increases,the curve slope of the spanwise Cpdecreases.

    Fig.9 Effect of sweep angle on spanwise pressure distribution at α =2°.

    Fig.10 Variation of spanwise Cp,min.

    The results suggest that at a specific AOA,as the wing sweep angle increases,the region covered by the vortex generated on the quarter-chord line near the wing tip decreases in size.As can be seen in Fig.7(b)and(d),at an AOA of 8°,this area of low pressure decreases from 0.54<2y/b<0.97 for the 23DS wing to 0.54<2y/b<0.74 for the 40DS wing.

    Fig.10 illustrates the effects of the AOA and increasing the wing sweep angle on the minimum value of the spanwise pressure distribution.As it can be seen in the figure,at an 8°AOA,|Cp,min|decreases when increasing the sweep anglefrom 23°to 40°.In other words,at an 8°angle of attack,with an increasing sweep angle,the strength of the vortex formed on the quarterchord line diminishes.However,it is found that with increasing the AOA to higher values,the magnitude of spanwise Cp,minover the 23DS wing decreases when compared to that of the 40DS wing(Fig.10).As it is shown,for angles of attack of 12°and 14°,the values of spanwise Cp,minfor the wings go in opposite directions;in other words,the vortex formed on the quarter-chord line of the 23DS wing weakens in strength,and the magnitude ofits Cp,mindrops.On the other hand,the vortex formed on the 40°-sweep wing quarter-chord line gets stronger,increasing |Cp,min.|.Although in the range of AOAs under experiments,the re is no sign of a sudden drop in spanwise Cp,minfor the 40DS wing,it can be predicted that the same process will take place at a higher AOA.

    Results show that at a specific angle of attack,increasing the sweep angle results in generating a stronger vortex on the quarter-chord line,which has lower sensitivity to angle of attack variations.In other words,as the sweep angle increases,the vortex remains attached to the wing surfacefor higher angles of attack,causing vortex breakdown and Sectional stall to occur at higher angles of attack.

    3.3.Effect of end plate on chordwise Cp

    Figs.11 and 12 show the results of this investigation(pressure distribution over two swept WWE)and those of19(pressure distribution over two swept WOE)at angles of attack of 0°,2°,4°,6°,and 8°.In addition,for more accurate deduction,corresponding streamline tracks over the surfaces of the wings are shown in the sefigures.

    Fig.11 Comparison of chordwise Cpfor Λ =23°.

    A comparison of the results shows that in the case of the 23DS wing,at an AOA of 0°,the end plate has almost no effect on the chordwise Cp,and the predicted data are in very good agreement with those of19;however,the small difference between the two pressure data in the outboard station is within the measurement accuracy range.Moreover,comparing the CFD results of WWE to those of WOE in Fig.12(a)shows no remarkable difference in streamline traces on the two wing surfaces.As the AOA increases to 4°,data obtained in this investigation present a larger|Cp|near the LE of the chordwise sections than that in the case of WOE.However,the qualitativefeatures of the pressure distributions remain the same as the end plate interference is included(Fig.11(b)).

    The end plate causes a blockage effect on the flow and additional viscous effects,especially near the end plate24.It should be pointed out that the use of an end plate which acts as a barrier to the spanwiseflow along the outboard portion of a wing and the refore eliminates generation of wing tip vortices,causes the streamlines in the neighborhood of a model to curve as the y pass over the wing surface.Moreover,near the wing tip area,the streamlines curve more sharply near the LE and close to the end plate.These spanwise streamlines encountering the end plate result in forming a vortex-shaped flow on the surface of the end plate.It seems that installing an end plate on the tip of a wing model causes the magnitude of pressure coefficient near the LE of the chordwise sections to increase,especially in the outboard-section pressure distribution,when compared to that of WOE,but as the distancefrom the se regions is increased,the comparison shows good agreement specially in the middle and inner sections of the wing(Fig.11(b)).

    Fig.11(c)and(d)shows that in the case of WWE,increasing the AOA causes the streamlines curve more sharply.Furthermore,the vortex formed on the surface of the end plate increases in size and moves downstream to the end-plate TE,so that at an 8°AOA,its trace amplitude moves toward the middle section of the wing.

    Fig.12 Comparison of chordwise Cpfor Λ =40°.

    A comparison between Fig.11(c)and 11(d)shows that as the AOA increases from 6°to 8°,in Sections 1 and 2 pressure distributions,only the areas near the LE are affected by the disturbed flow created on the surface of the end plate.On the contrary,it seems that curving the streamlines more sharply and increasing the size of the end-plate-surface vortex have a greater effect on the wing’s outboard-section pressure distribution.As a result,it causes a rapid drop in Cp,minnear the outboard-section LE,followed by a lack of pressure recovery toward the TE.As can be seen in Fig.11(d),in the case of WWE,the pressure data show an area of constant pressure near the LE of Section 3 which covers a region of 0.05≤x/C≤0.14,while in the case of WOE,no indication of a constant pressure area is seen.As previously discussed,the plateau of the pressure distribution curve is a characteristic of a laminar separation bubble,which takes place at lower angles of attack for the case of WWE compared to those of WOE.Moreover,the results show a sudden drop in the magnitude of Cpnear the TE of WWE compared to that of WOE.

    It is found that by increasing the strength and size of the end-plate-surface vortex,its effect on Cpnear the TE increases,which results in a decrease of the airstream velocity near the outboard-section TE.As it is shown in Fig.11(a)–(d),with an increasing AOA, |Cp|in this area of the wing drops more sharply.

    A comparison between Figs.11 and 12 shows that the end plate has a greater effect on Cpnear the LE and TE of the 40DS wing.Moreover,Comparing the CFD results of the two studies solidly confirms that the streamline traces on the 40DS wing surface are affected more greatly by the end plate than in the case of the 23DS wing.

    As it is presented in Fig.12(a),eliminating wing tip vortices causes more gradual acceleration of the flow and more gradual pressure recovery on all chordwise sections of the 40DS wing.This gradual acceleration of the flow from the LE to the minimum pressure point is clearly visible on Sections 1 and 2.Furthermore,the most significant effect of the end plate on pressure recovery occurs in the outboard-Section pressure distribution of the 40DS wing,which is in the vicinity of the vortex formed on the end plate(comparing pressure data of WWE to those of WOE in Fig.11(a)).

    CFD results show that in the case of the 40DS wing with an end plate,increasing the AOA causes the streamlines to curve more sharply toward the wing tip in comparison with those of WOE(Fig.12(a)–(d)).Furthermore,as the AOA increases,the differences between the pressure results in Sections 1 and 2 remain the same,but in Section 3,the suction force near the LE gets stronger compared to the results of WOE(Fig.12(b)–(d)).Moreover,like the case of the 23DS wing with an end plate,the results show a sudden drop in |Cp|near the TE of the 40DS wing with an end plate compared to that of WOE.As discussed before,with an increasing AOA,the end-plate-surface vortex gets stronger,which results in a stronger spanwiseflow,thus affecting the airstream velocity near the wing’s outboard-Section TE and causing a sudden drop in the pressure distribution.As for the case of the 23DS wing with an end plate,Fig.12(a)–(d)show that with an increasing AOA,|Cp|drops more sharply in this region of the 40DS wing.According to Ref.25,the presence of the lateral boundaries around the test-section(the walls of the test-section)produces a lateral constraint to the flow pattern around a body known as solid blocking.For closed throats,solid blocking is the same as an increase in dynamic pressure,increasing all forces and moments at a given angle of attack.Furthermore,the lateral boundaries produce an alteration to the local angle of attack along the span.

    It is found that the use of an end plat decreases the test section’s cross-sectional area,thus increasing the effect of lateral boundaries.Moreover,according to Ref.24,the end plate causes a blockage effect on the flow and additional viscous effect,especially near the end plate.These combined effects significantly increase the magnitudes of Cpnear the LE of the wing,especially in the outboard Section,and reinforce the spanwiseflow over the surface of the wing due to the formation of a vortex-shaped flow on the surface of the end plate.The results show that the use of an end plate to eliminate wing tip vortices have a greater effect on the 40DS wing.

    4.Conclusions

    (1)The results show that localized regions of separation followed by reattachment of the flow in the chordwise-Se ction pressure distributions occur at angles of attack somewhat below the stall.The region of separation starts near the leading edge of the outboard station(Section 3)and runs inboard and aft in a continuous bubble streak.

    (2)As the sweep angle increases,the formation of a laminar separation bubble takes place at higher angles of attack,and at a specified angle of attack,with an increasing sweep angle,the length of the area of constant pressure is decreased.Furthermore,when increasing the sweep angle,the origin of this vortex streak near the leading edge of a wing moves to the outboard Section of the wing.

    (3)Increasing the sweep angle causes the stall phenomena on the chordwise Sections of a swept wing to delay.

    (4)The results obtained show that with an increasing angle of attack,a vortex forms near the wing tip on the quartet-chord line,which pushes inward from the wing tip to the wing root and increases in size when increasing the angle of attack.

    (5)At a specific angle of attack,as the sweep angle increases,the vortex formed on the quarter-chord line gets stronger,furthe r increasing the negative pressure in this region of the wing.In addition,as the sweep angle increases,the vortex remains attached to the wing surfacefor higher angles of attack,causing vortex breakdown and Sectional stall to delay.

    (6)The results show that an end plate causes a blockage effect on the flow and extra viscous effect specifically in the vicinity of the end plate.These joined effects of the end plate significantly result in forming a vortexshaped flow on the surface of the end plate,and despite our expectations,reinforcement of the spanwiseflow over the surface of the wing,thus affecting the pressure distribution,when compared to a wing without an end plate.

    1.Nay HO.Low speed boundary layer and pressure distribution tests on a family of sweep back wings[dissertation].Pasadena(CA):California Institute of Technology;1952.

    2.Harper CW,Maki RL.A review of the stall characteristics of swept wings.Moffett Field(CA):NASA Ames Research Center;1964.Report No.:64N24822.

    3.Anderson BT,Mayer RR.Effects of wing sweep on boundarylayer transition for a smooth F-14A wing at Mach numbers from 0.700 to 0.825.Edwards(CA):NASA Hugh L.Dryden Flight Research Facility;1990.Report No.:91N24556.

    4.Wing F,Supercritical VC,Powers SG,Webb LD.Flight wing surface pressure and boundary-layer data report from the F-111 smooth variable-camber supercritical mission adaptive wing.Edwards(CA):NASA Dryden Flight Research Center;1997.Report No.:91N24556.

    5.Soltani MR,Ghorbanian K,Masdari M.Investigation of the pressure distribution and transition point over a swept wing.Sci Iran 2011;18(6):1277–86.

    6.Yen SC,Hsu CM.Investigation on vortex shedding of a sweptback wing.Exp Therm Fluid Sci 2007;31(8):849–55.

    7.Houghton EL,Carpenter PW.Aerodynamics for engineering students.5th ed.Oxford:Butterworth-Heinemann;2003.p.522–3.

    8.Giuni M,Benard E,Green RB.SPIV experiments on wing trailing vortices.Engineering graduate school conference;2010.

    9.Spalart PR.Airplane trailing vortices.Annu Rev Fluid Mech 1998;30(1):107–38.

    10.Sohn MH,Chang JW.Visualization and PIV study of wingtip vortices for three different tip configurations.Aerosp Sci Technol 2012;16(1):40–6.

    11.Fink MP,Lastinger JL.Aerodynamic characteristics of lowaspect-ratio wings in close proximity to the ground.Hampton(VA):NASA Langley Research Center;1961.Report No.:NASATN-D-926,L-1367.

    12.Carter AW.Effect of ground proximity on the aerodynamic characteristics of aspect-ratio-1 airfoils with and without end plate.Hampton(VA):NASA Langley Research Center;1961.Report No.:NASA-TN-D-970,L-1693.

    13.Lee J.Influence of end plate on aerodynamic characteristics of low-aspect-ratio wing in ground effect.J Mech Sci Technol 2008;22(12):2578–89.

    14.Boutilier MSH,Yarusevych S.Effect of end plates and blockage on a low Reynolds number airfoil experiment.41st AIAA fluid dynamics conference and exhibit.Reston:AIAA;2011.

    15.Alam MM,Zhou Y,Yang HX.The ultra-low Reynolds number airfoil wake.Exp Fluids 2010;48(1):81–103.

    16.Lee T,Gerontakos P.Investigation offlow over an oscillating airfoil.J Fluid Mech 2004;512:313–41.

    17.Pelletier A,Mueller TJ.Low Reynolds number aerodynamics of low-Aspect-Ratio,thin/flat/cambered-plate wings.J Aircraft 2000;37(5):825–32.

    18.Pelletier A,Mueller TJ.Effect of endplates on two-dimensional airfoil testing at low Reynolds numbers.J Aircraft 2001;38(6):1056–9.

    19.Soltani MR,Masdari M,Seidjafari M,Rasi F.Experimental study of the flow field over two swept wings.13th annual and 2nd international fluid dynamics conference;2010.

    20.Howard RM.Principles of random signal analysis and low noise design.1st ed.Australia:John Wiley&Sons Inc.;2002.p.76–8.

    21.Davis RL,Carter JE.Analysis of transitional separation bubble on infinite swept wing.AIAA J 1987;25(3):421–8.

    22.Lee S.Effect of supercooled large-droplet icing on airfoil aerodynamics[dissertation].Champaign(IL):University of Illinois;2001.

    23.White EB,Saric WS,Gladden RD,Gabet PM.Stage of swept wing transition.39th aerospace sciences meeting&exhibit.Reston:AIAA;2001.

    24.Jung JH,Kim MJ,Yoon HS,Hung PA,Chun HH,Park DW.Endplate effect on aerodynamic characteristics of three dimensional wings in closefree surface proximity.Int J Nav Arch Ocean 2012;4(4):477–87.

    25.Rathakrishnan E.Instrumentation,measurements and experiments in fluids.1st ed.New York:CRC Press;2007.p.141–2.

    19 July 2016;revised 19 December 2016;accepted 24 January 2017

    Available online 22 August 2017

    Chordwise;

    Flow field;

    Pressure distribution;

    Swept wing;

    Wind tunnel

    *Corresponding author.

    E-mail address:m.masdari@ut.ac.ir(M.MASDARI).

    Peer review under responsibility of Editorial Committee of CJA.

    Production and hosting by Elsevier

    http://dx.doi.org/10.1016/j.cja.2017.07.008

    1000-9361?2017 Production and hosting by Elsevier Ltd.on behalf of Chinese Society of Aeronautics and Astronautics.This is an open access article under the CC BY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    ?2017 Production and hosting by Elsevier Ltd.on behalf of Chinese Society of Aeronautics and Astronautics.This is an open access article under the CC BY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    精品久久蜜臀av无| 99热精品在线国产| 一级黄色大片毛片| 在线观看一区二区三区| 极品教师在线免费播放| 别揉我奶头~嗯~啊~动态视频| 欧美中文综合在线视频| 老司机深夜福利视频在线观看| 神马国产精品三级电影在线观看| 亚洲中文av在线| 国产亚洲精品久久久com| 亚洲精品乱码久久久v下载方式 | 日韩国内少妇激情av| 波多野结衣高清无吗| 变态另类丝袜制服| 精品久久久久久久末码| 国产av一区在线观看免费| 欧美日韩瑟瑟在线播放| 亚洲真实伦在线观看| 国产精品日韩av在线免费观看| 成年人黄色毛片网站| 伦理电影免费视频| 一区二区三区激情视频| 欧美黄色片欧美黄色片| 久9热在线精品视频| 国产 一区 欧美 日韩| 亚洲专区国产一区二区| 国产精品久久久av美女十八| 亚洲av熟女| 久久草成人影院| 成人永久免费在线观看视频| 亚洲欧美日韩无卡精品| 免费大片18禁| 99久久综合精品五月天人人| 欧美激情久久久久久爽电影| 精品电影一区二区在线| 国产精品久久久人人做人人爽| 国产蜜桃级精品一区二区三区| 久久婷婷人人爽人人干人人爱| 欧美+亚洲+日韩+国产| 国产一区在线观看成人免费| 日韩免费av在线播放| 午夜福利高清视频| 午夜视频精品福利| 人人妻人人澡欧美一区二区| 国产高清激情床上av| 国产午夜精品久久久久久| 俄罗斯特黄特色一大片| 99riav亚洲国产免费| 狂野欧美白嫩少妇大欣赏| 亚洲国产中文字幕在线视频| 国产亚洲精品综合一区在线观看| 亚洲成人免费电影在线观看| 午夜影院日韩av| 久久精品夜夜夜夜夜久久蜜豆| 两个人视频免费观看高清| 欧美3d第一页| 好看av亚洲va欧美ⅴa在| 人人妻人人澡欧美一区二区| 欧美在线一区亚洲| 丰满人妻熟妇乱又伦精品不卡| 国产乱人伦免费视频| 97人妻精品一区二区三区麻豆| 一个人免费在线观看电影 | 亚洲成人中文字幕在线播放| а√天堂www在线а√下载| 国产69精品久久久久777片 | 又大又爽又粗| 久久久精品大字幕| 此物有八面人人有两片| 日本a在线网址| 婷婷丁香在线五月| 欧美不卡视频在线免费观看| 热99在线观看视频| 午夜精品在线福利| 国产在线精品亚洲第一网站| av在线蜜桃| 亚洲av免费在线观看| 午夜福利视频1000在线观看| 波多野结衣巨乳人妻| 日本免费一区二区三区高清不卡| 黄色视频,在线免费观看| 桃色一区二区三区在线观看| 国产亚洲精品一区二区www| 村上凉子中文字幕在线| 男女视频在线观看网站免费| 天堂网av新在线| 老汉色∧v一级毛片| 国产黄a三级三级三级人| 在线观看日韩欧美| 一区二区三区国产精品乱码| 国产野战对白在线观看| а√天堂www在线а√下载| 美女扒开内裤让男人捅视频| 精品一区二区三区视频在线观看免费| 九九热线精品视视频播放| 亚洲成a人片在线一区二区| 好看av亚洲va欧美ⅴa在| 国产成人系列免费观看| 夜夜爽天天搞| 日韩大尺度精品在线看网址| 午夜影院日韩av| 两个人视频免费观看高清| 亚洲精品乱码久久久v下载方式 | 男女床上黄色一级片免费看| 国产真实乱freesex| 一区二区三区国产精品乱码| 日韩欧美国产一区二区入口| 欧美日韩乱码在线| 九色国产91popny在线| 在线a可以看的网站| 男女下面进入的视频免费午夜| 亚洲国产欧洲综合997久久,| 老熟妇乱子伦视频在线观看| 一进一出抽搐gif免费好疼| 97超视频在线观看视频| 日本免费a在线| 国产亚洲av嫩草精品影院| 国产综合懂色| 一个人看视频在线观看www免费 | 亚洲国产日韩欧美精品在线观看 | 一进一出抽搐gif免费好疼| 国产午夜福利久久久久久| 精品久久久久久,| 欧美av亚洲av综合av国产av| 国产亚洲欧美在线一区二区| 午夜福利视频1000在线观看| 国产精品永久免费网站| h日本视频在线播放| 亚洲欧美日韩无卡精品| 久久中文字幕一级| 亚洲人成伊人成综合网2020| 成人精品一区二区免费| 999久久久国产精品视频| 国产精品一区二区三区四区久久| 免费大片18禁| 免费在线观看成人毛片| 男女那种视频在线观看| 无限看片的www在线观看| 亚洲精品在线美女| 男人舔女人下体高潮全视频| 日韩欧美免费精品| 18禁黄网站禁片午夜丰满| 国产av麻豆久久久久久久| 久久久久久久精品吃奶| 国产精品久久久久久人妻精品电影| 麻豆成人午夜福利视频| 国内久久婷婷六月综合欲色啪| 99国产精品99久久久久| 变态另类成人亚洲欧美熟女| 欧美又色又爽又黄视频| АⅤ资源中文在线天堂| 久久天躁狠狠躁夜夜2o2o| 美女高潮的动态| 成人精品一区二区免费| 色哟哟哟哟哟哟| 99久久精品热视频| 香蕉丝袜av| 热99re8久久精品国产| 成人一区二区视频在线观看| 国产黄片美女视频| 深夜精品福利| 亚洲欧美一区二区三区黑人| 久久久久性生活片| 国产精品野战在线观看| 亚洲va日本ⅴa欧美va伊人久久| 成年免费大片在线观看| 欧洲精品卡2卡3卡4卡5卡区| 人妻丰满熟妇av一区二区三区| 一二三四在线观看免费中文在| 少妇丰满av| 一个人免费在线观看电影 | 看片在线看免费视频| 人妻夜夜爽99麻豆av| 美女高潮喷水抽搐中文字幕| 日本 欧美在线| 成人性生交大片免费视频hd| 日韩欧美国产在线观看| 成人高潮视频无遮挡免费网站| av天堂在线播放| 国产精品日韩av在线免费观看| 嫁个100分男人电影在线观看| 精品久久久久久成人av| 熟妇人妻久久中文字幕3abv| 91字幕亚洲| 久久精品综合一区二区三区| 美女 人体艺术 gogo| 婷婷亚洲欧美| 嫩草影院精品99| 国产久久久一区二区三区| 一区二区三区高清视频在线| 精品免费久久久久久久清纯| 国产又黄又爽又无遮挡在线| 免费看a级黄色片| 国产精品久久视频播放| 亚洲自偷自拍图片 自拍| 在线免费观看的www视频| 国产视频内射| 日本黄大片高清| 午夜亚洲福利在线播放| 国产精品野战在线观看| 男女下面进入的视频免费午夜| 久9热在线精品视频| 天天一区二区日本电影三级| 亚洲中文av在线| 精品久久久久久成人av| 午夜影院日韩av| 男人的好看免费观看在线视频| 黑人操中国人逼视频| 国产精品久久久久久精品电影| 黄色片一级片一级黄色片| 婷婷亚洲欧美| 国产精品自产拍在线观看55亚洲| 在线播放国产精品三级| 成人国产一区最新在线观看| 男插女下体视频免费在线播放| 国产午夜精品论理片| 黄片小视频在线播放| 两人在一起打扑克的视频| 亚洲成人免费电影在线观看| 欧美在线一区亚洲| 亚洲一区二区三区不卡视频| 五月玫瑰六月丁香| 91麻豆av在线| 久久精品人妻少妇| 男人的好看免费观看在线视频| 国产成+人综合+亚洲专区| 97超视频在线观看视频| 国产精品 欧美亚洲| 日韩精品青青久久久久久| 久久天躁狠狠躁夜夜2o2o| 9191精品国产免费久久| 亚洲精品国产精品久久久不卡| 亚洲精品乱码久久久v下载方式 | 99国产综合亚洲精品| 久久午夜综合久久蜜桃| 国产主播在线观看一区二区| 久久久国产成人精品二区| 亚洲精品在线美女| 国产精品久久视频播放| 久久久久国产精品人妻aⅴ院| 亚洲熟妇熟女久久| 黑人欧美特级aaaaaa片| 最新中文字幕久久久久 | 国产一级毛片七仙女欲春2| 亚洲欧美一区二区三区黑人| 欧美乱码精品一区二区三区| 国产三级黄色录像| 夜夜躁狠狠躁天天躁| 亚洲国产看品久久| 欧美丝袜亚洲另类 | 亚洲一区高清亚洲精品| 伊人久久大香线蕉亚洲五| 性色avwww在线观看| 亚洲专区字幕在线| 19禁男女啪啪无遮挡网站| 99热这里只有精品一区 | tocl精华| 国产亚洲精品久久久久久毛片| 国产黄色小视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 久久精品亚洲精品国产色婷小说| 后天国语完整版免费观看| 欧美色欧美亚洲另类二区| 久久精品国产综合久久久| 97人妻精品一区二区三区麻豆| 老鸭窝网址在线观看| 日本一本二区三区精品| av在线蜜桃| 一个人免费在线观看的高清视频| 精品熟女少妇八av免费久了| 中文字幕久久专区| av在线蜜桃| 成人特级av手机在线观看| 最新在线观看一区二区三区| 亚洲一区二区三区不卡视频| 五月伊人婷婷丁香| 18美女黄网站色大片免费观看| 国产激情久久老熟女| 999久久久精品免费观看国产| 色综合站精品国产| 又爽又黄无遮挡网站| 最近在线观看免费完整版| 久久精品国产99精品国产亚洲性色| 婷婷亚洲欧美| 亚洲熟妇中文字幕五十中出| 国产成人aa在线观看| 看黄色毛片网站| 国产亚洲精品久久久久久毛片| 国产伦在线观看视频一区| 日韩欧美国产一区二区入口| 五月玫瑰六月丁香| 可以在线观看的亚洲视频| 欧美日韩瑟瑟在线播放| 男女床上黄色一级片免费看| 亚洲欧美激情综合另类| 精品福利观看| 国产黄a三级三级三级人| 桃色一区二区三区在线观看| 又紧又爽又黄一区二区| 国产午夜精品久久久久久| 成人鲁丝片一二三区免费| 国产乱人伦免费视频| АⅤ资源中文在线天堂| 少妇的丰满在线观看| 亚洲人成网站在线播放欧美日韩| 久久香蕉国产精品| 久久草成人影院| 亚洲 欧美 日韩 在线 免费| 久久久久久久午夜电影| 老司机午夜福利在线观看视频| 99视频精品全部免费 在线 | 国产高清视频在线观看网站| 色在线成人网| 亚洲av日韩精品久久久久久密| 久久草成人影院| 日本黄色视频三级网站网址| 国产精品日韩av在线免费观看| 99久国产av精品| 亚洲成人免费电影在线观看| 最近视频中文字幕2019在线8| 国产激情久久老熟女| 999精品在线视频| 少妇裸体淫交视频免费看高清| 黄色丝袜av网址大全| 999久久久国产精品视频| av视频在线观看入口| 精品国内亚洲2022精品成人| 日韩欧美精品v在线| 精品欧美国产一区二区三| 亚洲专区字幕在线| 亚洲av美国av| www国产在线视频色| 国产一级毛片七仙女欲春2| 欧美成人一区二区免费高清观看 | 日韩欧美免费精品| 免费在线观看视频国产中文字幕亚洲| 国产伦人伦偷精品视频| 超碰成人久久| 51午夜福利影视在线观看| 精品久久久久久久人妻蜜臀av| 日韩欧美精品v在线| 亚洲色图 男人天堂 中文字幕| 日韩三级视频一区二区三区| 日本 欧美在线| 国产成人啪精品午夜网站| 麻豆久久精品国产亚洲av| 国产一区二区三区在线臀色熟女| 色av中文字幕| 久久久久国产精品人妻aⅴ院| 国产伦精品一区二区三区视频9 | 免费电影在线观看免费观看| xxxwww97欧美| 在线a可以看的网站| 99国产综合亚洲精品| 亚洲天堂国产精品一区在线| 97碰自拍视频| 这个男人来自地球电影免费观看| 在线国产一区二区在线| 草草在线视频免费看| 18禁黄网站禁片午夜丰满| 国产精品99久久99久久久不卡| 久久久久久大精品| 日韩欧美国产在线观看| av中文乱码字幕在线| 青草久久国产| 亚洲最大成人中文| 亚洲中文av在线| 久久久久精品国产欧美久久久| www.自偷自拍.com| 欧美成人一区二区免费高清观看 | 视频区欧美日本亚洲| 久久久国产成人精品二区| 美女扒开内裤让男人捅视频| 亚洲国产中文字幕在线视频| 久久精品人妻少妇| 欧美乱妇无乱码| 岛国在线观看网站| 最近视频中文字幕2019在线8| 99热精品在线国产| 国产蜜桃级精品一区二区三区| 久久久国产成人精品二区| 国产精品久久久久久亚洲av鲁大| 青草久久国产| 国产精品国产高清国产av| a级毛片a级免费在线| 真实男女啪啪啪动态图| 日韩 欧美 亚洲 中文字幕| 狂野欧美白嫩少妇大欣赏| 黄色 视频免费看| 高清毛片免费观看视频网站| av天堂中文字幕网| av在线天堂中文字幕| 国产aⅴ精品一区二区三区波| 美女午夜性视频免费| 在线看三级毛片| 日韩欧美 国产精品| 一本一本综合久久| 欧美日韩精品网址| 日本免费a在线| 免费观看精品视频网站| 日韩欧美 国产精品| 两个人视频免费观看高清| netflix在线观看网站| 亚洲18禁久久av| 欧美一级毛片孕妇| 又黄又粗又硬又大视频| 婷婷六月久久综合丁香| 18禁黄网站禁片午夜丰满| 午夜日韩欧美国产| 99re在线观看精品视频| 国产精品亚洲av一区麻豆| 亚洲欧美精品综合一区二区三区| 99精品欧美一区二区三区四区| 亚洲精品456在线播放app | 亚洲欧美日韩无卡精品| 久久九九热精品免费| 特级一级黄色大片| 变态另类丝袜制服| 在线播放国产精品三级| 欧美日本亚洲视频在线播放| 真人一进一出gif抽搐免费| 成年版毛片免费区| 精品国产乱码久久久久久男人| 精品一区二区三区av网在线观看| 啦啦啦免费观看视频1| 三级毛片av免费| 亚洲一区二区三区色噜噜| 老鸭窝网址在线观看| 日韩国内少妇激情av| 露出奶头的视频| 怎么达到女性高潮| 国产欧美日韩一区二区三| 欧美极品一区二区三区四区| 国产乱人视频| 毛片女人毛片| 亚洲成人中文字幕在线播放| 国产一区二区在线观看日韩 | 久久精品aⅴ一区二区三区四区| bbb黄色大片| 在线观看免费视频日本深夜| 日韩欧美免费精品| 亚洲av片天天在线观看| 亚洲av免费在线观看| 伦理电影免费视频| 一区福利在线观看| av片东京热男人的天堂| 国产高清视频在线播放一区| 成人无遮挡网站| 亚洲欧美日韩东京热| 久久香蕉国产精品| 精品国产乱子伦一区二区三区| 香蕉av资源在线| 一卡2卡三卡四卡精品乱码亚洲| av中文乱码字幕在线| 女警被强在线播放| 身体一侧抽搐| 午夜免费激情av| 黄色成人免费大全| 两个人视频免费观看高清| 91九色精品人成在线观看| 欧美性猛交黑人性爽| 免费观看精品视频网站| 国产成年人精品一区二区| 国产av一区在线观看免费| 男女下面进入的视频免费午夜| 国产成人欧美在线观看| 国产av在哪里看| 午夜激情欧美在线| av在线蜜桃| 两性夫妻黄色片| 成年女人看的毛片在线观看| 给我免费播放毛片高清在线观看| 女同久久另类99精品国产91| 啪啪无遮挡十八禁网站| 亚洲人成网站在线播放欧美日韩| 国产成人一区二区三区免费视频网站| 亚洲国产精品999在线| 国产高清videossex| 午夜免费观看网址| 欧洲精品卡2卡3卡4卡5卡区| 久久99热这里只有精品18| 中出人妻视频一区二区| 欧美在线黄色| www.精华液| 99久久精品热视频| 欧美另类亚洲清纯唯美| 巨乳人妻的诱惑在线观看| 国产97色在线日韩免费| av天堂在线播放| 日本撒尿小便嘘嘘汇集6| 日韩人妻高清精品专区| 精品熟女少妇八av免费久了| 国产高潮美女av| 真人一进一出gif抽搐免费| 亚洲一区二区三区色噜噜| 美女cb高潮喷水在线观看 | 夜夜爽天天搞| 欧美日韩精品网址| 老鸭窝网址在线观看| 成年版毛片免费区| 久久亚洲真实| 国产精品99久久99久久久不卡| 欧美黑人巨大hd| 国产精品99久久99久久久不卡| 最近最新免费中文字幕在线| 黄色日韩在线| 夜夜爽天天搞| 黄色成人免费大全| 免费看十八禁软件| 国产亚洲精品久久久久久毛片| 久久精品影院6| 中亚洲国语对白在线视频| 特级一级黄色大片| 欧美一区二区国产精品久久精品| 精品国内亚洲2022精品成人| 视频区欧美日本亚洲| 成人永久免费在线观看视频| 成人午夜高清在线视频| 欧美日韩精品网址| 亚洲黑人精品在线| 老司机深夜福利视频在线观看| 天天添夜夜摸| 91九色精品人成在线观看| 成年女人永久免费观看视频| 日韩精品中文字幕看吧| 欧美日韩一级在线毛片| 国产不卡一卡二| 免费在线观看日本一区| 亚洲欧美日韩卡通动漫| 亚洲成人久久性| 色综合欧美亚洲国产小说| 一二三四社区在线视频社区8| 日韩国内少妇激情av| 香蕉国产在线看| 亚洲国产精品sss在线观看| 日本一本二区三区精品| 国产黄色小视频在线观看| 国产一级毛片七仙女欲春2| 久久久国产成人精品二区| 亚洲av美国av| 老司机深夜福利视频在线观看| 午夜精品在线福利| 一级a爱片免费观看的视频| 人妻丰满熟妇av一区二区三区| 亚洲欧美一区二区三区黑人| 12—13女人毛片做爰片一| 欧美日韩亚洲国产一区二区在线观看| 国产成人福利小说| 久久久久国产精品人妻aⅴ院| 国产真实乱freesex| 最新中文字幕久久久久 | 97碰自拍视频| 国产人伦9x9x在线观看| 国产精品自产拍在线观看55亚洲| 亚洲熟女毛片儿| 欧美3d第一页| 国产欧美日韩精品亚洲av| 国产精品乱码一区二三区的特点| 午夜精品久久久久久毛片777| 最新中文字幕久久久久 | 亚洲熟妇中文字幕五十中出| 国内精品一区二区在线观看| 国产淫片久久久久久久久 | 亚洲av五月六月丁香网| 一进一出抽搐动态| 男女做爰动态图高潮gif福利片| 亚洲av片天天在线观看| 久久久久九九精品影院| 国产精品久久久av美女十八| 亚洲 欧美一区二区三区| 男女视频在线观看网站免费| 国产精品影院久久| 久久人妻av系列| 精品免费久久久久久久清纯| 欧美一区二区精品小视频在线| 日本在线视频免费播放| 俄罗斯特黄特色一大片| 身体一侧抽搐| www.精华液| 国产精品自产拍在线观看55亚洲| 精品一区二区三区视频在线观看免费| 特大巨黑吊av在线直播| 变态另类丝袜制服| av天堂中文字幕网| 首页视频小说图片口味搜索| 亚洲精品美女久久av网站| 国产黄片美女视频| 国产视频内射| 日韩欧美一区二区三区在线观看| 一级毛片女人18水好多| 免费大片18禁| 老汉色∧v一级毛片| 国产精品 欧美亚洲| 国产一区二区激情短视频| 国产免费男女视频| 成人无遮挡网站| www.自偷自拍.com| 欧美成人免费av一区二区三区| 91麻豆精品激情在线观看国产| 精品国内亚洲2022精品成人| 亚洲无线在线观看| 国产精品乱码一区二三区的特点| 性色avwww在线观看| 老熟妇仑乱视频hdxx| www国产在线视频色| 国产私拍福利视频在线观看| 国产激情欧美一区二区| 在线观看美女被高潮喷水网站 | 偷拍熟女少妇极品色| 一个人看的www免费观看视频| 午夜精品久久久久久毛片777| 久久久国产欧美日韩av| 亚洲第一电影网av| 国产精品日韩av在线免费观看|