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

    Turbulent drag reduction by spanwise slot blowing pulsed plasma actuation

    2022-11-17 02:59:28BoruiZHENG鄭博睿YuanzhongJIN金元中MinghaoYU喻明浩YueqiangLI李躍強(qiáng)BinWU武斌andQuanlongCHEN陳全龍
    Plasma Science and Technology 2022年11期
    關(guān)鍵詞:金元

    Borui ZHENG(鄭博睿),Yuanzhong JIN(金元中),Minghao YU(喻明浩),Yueqiang LI(李躍強(qiáng)),Bin WU(武斌) and Quanlong CHEN(陳全龍)

    1 School of Automation and Information Engineering,Xi’an University of Technology,Xi’an 710048,People’s Republic of China

    2 Department of Aeronautics and Astronautics,Northwestern Polytechnical University,Xi’an 710072,People’s Republic of China

    3 The Green Aerotechnics Research Institute of Chongqing Jiaotong University,Chongqing 401120,People’s Republic of China

    Abstract This work studies the turbulent drag reduction(TDR)effect of a flat plate model using a spanwise slot blowing pulsed plasma actuator(SBP-PA).Wind tunnel experiments are carried out under a Reynolds number of 1.445×104.Using a hot-wire anemometer and an electrical data acquisition system,the influences of millisecond pulsed plasma actuation with different burst frequencies and duty cycles on the microscale coherent structures near the wall of the turbulent boundary layer(TBL)are studied.The experimental results show that the SBP-PA can effectively reduce the frictional drag of the TBL.When the duty cycle exceeds 30%,the TDR rate is greater than 11%,and the optimal drag reduction rate of 13.69% is obtained at a duty cycle of 50%.Furthermore,optimizing the electrical parameters reveals that increasing the burst frequency significantly reduces the velocity distribution in the logarithmic region of the TBL.When the normalized burst frequency reaches f+=2πfpd/U∞=7.196,the optimal TDR effectiveness is 16.97%,indicating a resonance phenomenon between the pulsed plasma actuation and the microscale coherent structures near the wall.Therefore,reasonably selecting the electrical parameters of the plasma actuator is expected to significantly improve the TDR effect.

    Keywords:turbulent boundary layer control,plasma flow control,hot-wire,turbulent frictional drag,turbulent drag reduction

    1.Introduction

    Plasma flow control,a relatively new active flow control concept based on plasma aerodynamic actuation,benefits from many technical advantages;for example,it is inexpensive to implement,imposes little burden in terms of weight,and has a simple structure and fast frequency response.To date,this technology has achieved important progress in many fields,such as anti-ici ng on aircraft,lift increase and drag reduction,and the performance improvement of unmanned aerial vehicles(UAVs)[1-4].Accordingly,plasma flow control is expected to greatly improve the aerodynamic and economic performance of aircraft.

    Dielectric barrier discharge(DBD)plasma actuation,a classic active flow control method,has been applied to control the turbulent boundary layer(TBL)for more than 20 years.The method of using a plasma actuator to modify the TBL was first proposed in 1994 by Roth[5].In this approach,a plasma actuator is excited by a high-frequency and highvoltage alternating current,which can alter the drag and momentum of the TBL.

    In recent years,the DBD plasma actuation technique has received substantial attention internationally,with considerable research being carried out in relevant academic fields.In 1994,Choi[6]carried out an active flow control experiment of plasma blowing in a channel and obtained a turbulent drag reduction(TDR)rate of 25%under a Reynolds number(Re)of 1800.Moreover,Dalvand[7]used the method of unsteady actuation to study the actuation characteristics of DBD,Pouryouussefi[8]used a neural network to optimize the discharge state of DBD,and after optimizing the geometric configuration and actuation parameters of the plasma actuator,low-speed wind tunnel experiments revealed that the wall oscillations during plasma actuation can suppress turbulent near-wall coherent structures and increase the local TDR rate to 45%[9,10].Jukes applied a spanwise DBD alternate actuation method to control the TBL,and hot-wire measurement results showed that the wall shear stress(WSS)was reduced by 45%[10].Similarly,Li[11]performed particle image velocimetry(PIV)experiments and found that pulsed plasma actuation can effectively reduce the WSS of a turbulent flow while stabilizing the streamwise vortices.In 2020,Corke[12]found that pulsed DC-DBD can generate a stable spanwise velocity,reduce the velocity fluctuation on the flat plate surface in the spanwise direction,hinder the ejection of low-velocity streaks near the wall,and inhibit the generation and evolution of coherent structures in the TBL,thereby reducing the WSS.In summary,the outcomes of previous wind tunnel experiments and numerical simulations confirm that pulsed DBD plasma actuation can effectively reduce turbulent friction drag;nevertheless,how to precisely control the TBL based on the evolution of microscale coherent structures therein is still worthy of in-depth study.

    To better understand the evolution of microscale coherent structures near the wall within the TBL and the corresponding control methods,scholars have conducted rigorous research on the development processes of near-wall streamwise vortices in the TBL.Hamilton[13]discovered that streamwise vortices near the wall of the TBL are self-sustaining and postulated both that the instability or transient growth of lowvelocity streaks leads to the generation of streamwise vortices and that low-velocity streaks represent the trajectories of selfsustaining streamwise vortices.Therefore,interrupting any growth stage of self-sustaining streamwise vortices is expected to suppress their further development and reduce the WSS of the turbulent flow.Baron[14]verified the above inference through numerical simulation and found that wall oscillations can suppress the ejection events of streamwise vortices and achieve a local TDR of 30%.Park found that local uniform blowing can uplift the whole trajectories of streamwise vortices and reduce the WSS generated by contact with the wall during the sweep events of streamwise vortices;however,the uplift of upstream streamwise vortices can enhance the turbulence intensity of streamwise vortices downstream and increase the downstream WSS[15].Therefore,reasonably choosing the control method and impact area are highly important.

    To interfere with self-sustaining streamwise vortices more effectively,scholars have compared the TDR effects of streamwise blowing layouts and spanwise blowing layouts through experiments and simulations.One analysis revealed that the transverse traveling waves generated by spanwise blowing can induce additional streamwise vortices,weaken the burst frequencies of ejection and sweep events of the original streamwise vortices,and suppress self-sustaining streamwise vortices near the wall,thereby stabilizing streaky structures;therefore,blowing in the spanwise direction can achieve greater TDR than blowing in the streamwise direction[16-21].Chen[18]conducted a TDR wind tunnel experiment with a pulsed slot jet and found that when the scale of blowing-induced vortices is similar to the thickness of the viscous sublayer,the flow control effectiveness of induced jet can act as a virtual barrier to streamwise vortices near the wall,which can isolate the interactions between the ejection and sweep events of streamwise vortices with the wall surface,reduce the turbulent WSS.However,it is worth noting that a powerful blowing intensity may aggravate the movements of microscale coherent structures near the wall,thereby increasing the WSS.Indeed,beyond a specific energy injection range,the TDR method may aggravate the cyclic regeneration of streamwise vortices near the wall.In addition,as the Reynolds number increases,the TDR effect gradually weakens,which may be due to the flow control with an improper intensity and impact area.As a result,the induced jets and vortices cannot penetrate the viscous sublayer and thus cannot affect the microscale coherent structures near the wall.

    As evidenced above,reasonably selecting the TDR control method,influence area,and control intensity according to the external environment is crucial[21-24].Therefore,there is an urgent need to develop a control method that can precisely apply pulsed plasma actuation to the internal structure of the TBL and control the energy therein appropriately without increasing turbulent drag or wasting energy.

    The authors assert that pulsed plasma actuation can effectively suppress the generation and evolution of streamwise vortices in the TBL by disturbing the self-sustainability of streamwise vortices near the wall,thereby reducing the WSS of the turbulent flow.In addition,spanwise blowing can act as a virtual barrier to near-wall streamwise vortices,as spanwise blowing isolates the interaction between the ejection and sweep events of streamwise vortices with the wall surface,and further enhances the TDR effect.Studies have shown that the current control method of slit blowing requires a gas source to generate airflow[14-18,21].Although DBD can operate without a gas source,DBD cannot generate jets perpendicular to the wall[7-9].Duong and Thomas found that DC-DBD can play a crucial role in turbulence drag reduction at a higher Reynolds number[12],but Xie found that DC-DBD also brings in substantial strong electromagnetic interference(EMI)during DBD operation,which will seriously interfere with the operation of data-acquisition system[25,26].Compared with DC-DBD,AC-DBD has lower power consumption with less EMI,and has the similar effect of reducing drag[1,7-9].

    Based on the above analysis,this work combines pulsed AC plasma actuation with spanwise slot blowing to analyze the TDR effect of the spanwise slot blowing pulsed plasma actuator(SBP-PA),and to obtain the optimal controlling parameters of the SBP-PA,and reveal the TDR control mechanism of plasma actuation.

    2.Experimental syste m

    A wind tunnel experiment was conducted in the low-speed wind tunnel of Northwestern Polytechnical University in China,as shown in figure 1.The dimensions of the experimental section of the wind tunnel are 675×720×200 mm3,and the length from the exit of the contraction section to the center of the test section is 2000 mm.The lower wall of the experimental section is an integrally processed organic glass panel,and the upper wall panel is a rugged steel plate.A rectangular measuring window with a length of 675 mm in the streamwise direction and a width of 220 mm in the spanwise direction is set in the upper section of the wind tunnel,and a three-dimensional measurement system comprising a hot-wire probe is installed in the upper part of the wind tunnel to acquire high-precision measurements(displacement precision:10 μm).A motor equipped with a frequency converter controls the flow speed of the wind tunnel in the range of 0-20 m s-1,and the turbulence intensity is less than 0.3%[27].

    This experiment investigated the TDR effect achieved by an SBP-PA on a plexiglass plate on the lower wall of the wind tunnel,and the incoming flow velocity was maintained at 9.6 m s-1(Re=1.445×104).The spanwise blowing device was placed 1820 mm downstream of the outlet of the contraction section of the wind tunnel,as shown in figure 2.To ensure that the TBL can fully develop in the test section of wind tunnel,a strip of 120-grit sandpaper with a width of 50 mm was pasted on the wall 150 mm from the exit of the contraction section as a transition roughness trip.In the experiment,the velocity distribution of the boundary layer was measured in the logarithmic region downstream of the plasma actuator to verify that the TBL fully developed in this region,and to verify the repeatability of the SBP-PA experiments.

    2.1.Slot blowing pulsed plasma actuation system

    The plasma actuation system included a sinusoidal AC highvoltage source(model CTP-2000K,CORONA Lab.),our SBP-PA,a high-voltage probe(type P6015A,Tektronix),a function generator(model AFG3021,Tektronix)and an oscilloscope(type TDSI1001B,Tektronix).The CTP-2000K power supply was used for the plasma discharge.The carrier frequency was fixed atF=7.6 kHz,which was the optimum impedance frequency.A digital pulse modulator was used to generate a control signal with a pulse frequency of 1-1000 Hz and a duty cycle of 1%-100%.The control signal time was adjusted with high precision to form the actuation waveform shown in figure 3,which was used to generate the pulsed plasma discharge.The pulse frequencyfpand the duty cycle τpare defined by equations(1)and(2)[28],respectively:

    The SBP-PA is composed of a DBD plasma actuator and a designed blowing device.The structure of the DBD plasma actuator is shown in figure 4.The length of the actuator is 70 mm.The upper and lower electrodes are each composed of a copper sheet with a width of 2 mm and 10 mm,respectively,and the actuator is underlain by an insulating layer composed of 0.24 mm thick Kapton insulating tape.The upper electrode is arranged 10 mm in front of the L-shaped module,and the distance between the adjacent electrodes is 25 mm.The DBD plasma actuator is pasted onto the L-shaped module on the lower wall of the wind tunnel,and the L-shaped module is fixed onto the flat plate by bolts to form the blowing module,as shown in figure 5.The blowing module is located 1820 mm downstream of the transition point of this Wind Tunnel and the hot-wire probe used to detect the flow control point is set at a distance of 50 mm downstream from the blowing module(the red five-pointed star in figure 2).In this work,the boundary layer velocity was recorded before and after plasma actuation.

    The SBP-PA is connected to a CTP-2000K power supply,and the pulse duty cycle and pulse frequency are adjusted through an AFG3021 function generator.When the actuator applies a sinusoidal AC voltage that exceeds a voltage threshold,the air on the upper surface of the SBP-PA is broken down,triggering an electronic avalanche,and generating nitrogen and oxygen ions.These ions are accelerated under the driving force of the electric field,thereby transferring momentum to neutral gas molecules through impacts and inducing jets from the positive electrode to the negative electrode,as shown in figure 4.As the actuation voltage continues to increase,the jet velocity gradually increases,causing the jet to pass through the slot in the spanwise direction and blow onto the turbulent near-wall surface,where the induced jet attempts to disrupt the self-sustaining streamwise vortices in the TBL by cutting off the sweep events of streamwise vortex,which can result in high friction near the wall,as shown in figure 6.In this way,the goal of TDR is achieved.

    2.2.Data acquisition system and method

    In this experiment,a 55P15 boundary layer probe and a Dantec StreamLine constant-temperature hot-wire anemometer were used to measure the TBL velocity profile.The hot-wire material of the 55P15 probe is tungsten,the length is 1.25 mm,the diameter is 5 μm,the resistance at 20 °C is 3.58 Ω,and the temperature coefficient of resistance(TCR)is 0.36% K-1.In this experiment,the measurement position of the hot-wire probe is along the centerline of the actuator and 50 mm downstream of the actuator,as shown in figure 7.The hot-wire anemometer mainly includes a host computer,a calibrator,and a hot-wire probe.Before the experiment,the automatic calibration unit matched with the hot-wire anemometer was used to perform a velocity-voltage calibration on the 55P15 probe within the calibration velocity range of 0.5-20 m s-1.According to the King formula and experimental conditions(environmental humidity,temperature,and air viscosity),the overheat ratio of the hot-wire anemometer was set to 1.8 to ensure the accuracy of the hot-wire measurement[29].During the experiment,because the noise generated by SBP-PA is much smaller than the turbulent background noise,the low-pass filter frequency of the hot wire is selected as 10 kHz to ensure the measurement accuracy.To measure the temperature change produced by the plasma actuator during operation,a 34460A Keysight multimeter was adopted to measure the resistance of the hot-wire material of the 55P15 probe;based on the relationship between the resistance and the ambient temperature,the temperature change was quantified,and the hot-wire temperature data were corrected.

    3.Experimental results

    3.1.Verification of TBL

    The frictional drag of the TBL can be calculated by the frictional velocity of the wall.Therefore,to quantify the WSS of the TBL,it is necessary to ensure that the velocity profile of the TBL is accurately measured.

    To indirectly obtain the wall WSS,the Clauser equation was utilized in this study to fit the velocity in the logarithmic region of the boundary layer by the least-squares method.The equation is as follows:

    whereu+=u/uτrepresents the dimensionless flow velocity,y+=yuτ/ν represents the dimensionless wall distance,νrepresents the kinematic viscosity in air,anduτis the wall friction velocity.Deriving formula(3)further leads to the following:

    The relationship between the friction velocityuτiand the measured time-averaged sequence(yi,ui)can be obtained,where different values oficorrespond to different normal distances:

    Equation(6)demonstrates that each time-averaged sequence(yi,ui)in the logarithmic layer corresponds to a certain friction velocityuτi,and only the values ofuτiare unknown.By applying least-squares regression to all the time-averaged sequences in the logarithmic region,the goodness of fit(R2)approaches 1,and the optimal friction velocityuτcan be obtained[29].

    Through equations(6)and(7),the skin friction coefficientCfand WSS τw=ρuτ2can be obtained:

    In equation(9),Rrepresents the TDR rate.WhenR>0,plasma actuation achieves a positive TDR;whenR<0,plasma actuation generates drag.τw-offrepresents the WSS at the measuring point without plasma actuation,whereas τw-onrepresents the WSS at the measuring point with plasma actuation.

    To ensure that the area containing the actuator is located within the TBL,when the plasma actuator was not working,the velocity of the boundary layer was measured 10 mm downstream of the actuator at a wind speed of 9.6 m s-1.Figure 8 shows the velocity distribution of the boundary layer without plasma actuation.The velocity distribution within the logarithmic region(50<y+<200)of the boundary layer downstream of the actuator conforms to a classic logarithmic distribution,indicating that the boundary layer in the measurement area is a fully developed TBL.Therefore,the placement of a plasma actuator on the lower wall of the wind tunnel does not affect the natural development of turbulence in the measurement area.

    Figure 1.Schematic diagram of the low-speed wind tunnel.

    Figure 2.Schematic diagram of the experimental layout.

    Figure 3.Pulse waveform diagram.

    Figure 4.Schematic diagram of the DBD plasma actuator.

    Figure 5.Schematic diagram of the installation of the blowing module.

    Figure 6.Schematic diagram of the blowing TDR strategy.

    Figure 7.The spatial location of the hot-wire anemometer.

    Figure 8.Influence of the arrangement of the plasma actuator on the TBL.

    Figure 9.The deviation of friction velocity from the average value of five measurements and the corresponding fitting residuals.

    3.2.Repeatability test

    To verify the repeatability of the SBP-PA experiment,we kept the incoming velocity 9.6 m s-1(Re=1.445×104),Vp-p=8 kV,τ=50%,fp=200 Hz,five repeated measurements were performed,and the corresponding friction speed was fitted.As shown in figure 9,the friction velocity deviationis measured five times,whereu τiis theith measurement value andis the average value of the five measurements.The confidence interval is from-0.0174 to+0.0174(m s-1),the root mean square(RMS)is 0.005 44(m s-1),and the RMS value is within the confidence interval,so the experiment is highly repeatable.

    3.3.TDR experimental results with slot blowing pulsed plasma actuation

    3.3.1.TDR effects under different duty cycles.Previous studies have found that a higher DBD actuation voltage corresponds to a stronger induced jet velocity and a longer jet distance;therefore,more energy can be injected into the TBL.This phenomenon is capable of interfering with self-sustaining near-wall streamwise vortices and coherent structures,and expected to further interfere with the structure of large-scale vortices in the outer region of the TBL,thereby regulating the influence of the structure of largescale vortices on the streamwise vortices within the inner boundary layer.However,the results of previous experiments indicate that more input energy to the plasma actuator does not correspond to a better outcome[18,30].First,an excessively high actuation voltage can cause the insulating layer of the plasma actuator to age prematurely;after the insulating layer ages,the equivalent resistance decreases,the internal current of the plasma actuator increases,and heat accumulates quickly on the surface of the insulating layer,which causes the plasma actuator to break down,catch fire,or begin arcing.Second,when the jet velocity induced by the plasma actuator is too high,after the jet flow merges with the incoming flow,a synthetic jet is generated at a certain angle to the main flow direction and easily causes the flow to separate from the wall.Third,the application of excessively strong plasma actuation on the near-wall surface may accelerate the growth of self-sustaining streamwise vortices,thereby increasing the WSS[18].Therefore,choosing the appropriate plasma actuation parameters according to external flow phenomena and the scale of the vortices is necessary for indepth study.

    To study the influence of the SBP-PA with different energies on the WSS of turbulent flow and evaluate the change in the WSS of the TBL under different pulse duty cycles,the experiment in this paper maintained an incoming flow velocityU∞=9.6 ms-1,a peak-to-peak plasma actuation voltageVp-p=8kV,a pulse frequencyfp=100 Hz,and a carrier frequencyF=7.6 kHz.Figure 10 shows the TDR achieved by spanwise blowing under different duty cycles(τ=15%,30%,50%,70%,and 90%),and the error distribution of friction velocity fitting at 95%confidence,it can be found that the friction speed fittingR2under different duty ratios is close to 1,which indicates that the fitting accuracy is good and the results are accurate.It can also be concluded from the changing trend that the TDR rateRis not positively correlated with the duty cycle τ.When the duty cycle increases from 15% to 30%,the TDR rate greatly increases by 216.8%.When the duty cycle increases from 30%to 90%,the TDR rate first increases slightly and then gradually decreases.The highest TDR rate(Rmax=13.67%)appears at a duty cycle of 50%,while the lowest TDR rate(Rmin=11.78%)appears at a duty cycle of 90%.In the initial stage of flow control,the plasma actuation energy needs to be strong enough to meet the basic requirements of the flow control energy.However,when the plasma actuation energy is higher than a certain threshold(50%),the injection of more energy fails to improve the TDR effect;in fact,the outcome is a reduction in the TDR effect.Hence,it is essential to match the physical parameters between the energy being injected into the TBL and the scale of the flow being controlled.

    Figure 10.The variation in the TDR rate under different duty cycles.

    The plasma flow control mechanism closely relates to the vortex structure,and the duty cycle technique can produce a series of periodic vortex structures that enhance the flow field’s vorticity and momentum.Studies have shown that lowfrequency pulse actuation can generate a series of large-scale shedding vortices,leading to momentum mixing between the mainstream and the shear layer.In comparison,highfrequency pulse actuation causes small-scale vortices to promote the transition to mainstream[31].A new vortex is formed in the pulsed actuation mode at each pulse cycle.The vortices induced by pulsed actuation propagate faster than that of continuous actuation,and the vortex core velocity of pulsed actuation increases with the value of the duty cycle.

    In this experiment,the duty cycle value represents the length of the discharge time in each pulse period,the period of duty-off is the development time of the induced vortices,and the period of duty-on is the evolution time of the induced vortices.When the duty cycle increases gradually,the duty-on time of the plasma actuator increases,the plasma density becomes more significant,and the discharge energy becomes more intensive.More energy will be injected into the flow field by the momentum transfer,which further promotes the formation of the induced vortex.

    However,when the duty cycle increases beyond a certain threshold,the vorticity of the induced vortices decreases,the distance of adjacent vortices decreases,and the propagation mode gradually evolves from intermittent pulse to continuous mode[32,33].So,when the duty ratio is smaller,the interaction between the primary vortex and the residual vortexof the last cycle is negligible.The induced vortex can resist the residual vortex’s impact,maintain a complete vortex configuration,and control the boundary layer in a pulsed manner.On the contrary,when the duty ratio is larger,the interaction between the primary vortex and the residual vortex of the last cycle increases,and the vorticity of the primary vortex is not enough to resist the impact of the residual vortex.Thus,it is difficult to form a complete vortex,and the scale of the induced vortex cannot match the thickness of the viscous bottom layer,resulting in an efficiency decrease in drag reduction.

    When the pulse duty cycle is greater than 30%,an increase in energy has little effect on the TDR effect of plasma actuation.Therefore,with the flow conditions remaining unchanged,a smaller duty cycle can be selected to improve the energy efficiency,which is more conducive to conserving energy and is therefore valuable for environmental protection applications.

    Table 1 compares the changes in the boundary layer parameters under different duty ratios.When the duty ratio is less than 50%,uτi,WSS,andCfgradually decrease,and the TDR effect gradually increases;however,when the duty cycle

    Table 1.Comparison of the related calculation parameters from different duty cycles at the measuring point.

    τ≥50%,the TDR effect of the SBP-PA begins to gradually decrease.This analysis shows that when the duty cycle tends toward approximately 50%,the scale of vortices induced by the SBP-PA jet is close to the size of the turbulent near-wall streamwise vortices,and the blowing height is similar to the thickness of the viscous sublayer;thus,this process has a better inhibitory effect on self-sustaining near-wall streamwise vortices.However,when the duty cycle deviates from the optimal control parameter(τ≈50%),the energy produced by the SBP-PA and the scale of the induced vortices deviate from the energy and spatial scales of the current mainstream streamwise vortices,and the control effect is slightly reduced.In a future study,a high-frequency and high-spatial resolution flow visualization experiment will be carried out to further verify this inference.

    3.4.TDR effects at different pulse frequencies

    Previous plasma experiments discovered that a resonance-like phenomenon occurs when there is a coupling relationship between the pulse frequency of plasma actuation and the natural frequency of the incoming flow,and this phenomenon greatly enhances the plasma flow control effect[15].To obtain the optimal of plasma parameters,improve the effectiveness of the plasma actuator,and expand the application scope of plasma flow control engineering,based on the experimental results of the duty cycle described in our previous section,the influences of different pulse frequencies on the frictional drag of turbulent flow were studied in detail by maintaining an incoming flow velocityU∞=10 m s-1,with a peak-to-peak voltageVp-p=8 kV,and a duty cycle τ=50%.This analysis essentially demonstrates that the TDR effect of the plasma actuator on the TBL varies at different pulse frequencies(fp=50 Hz,100 Hz,and 200 Hz).

    The movement of the hot-wire probe in the boundary layer is controlled by the three-dimensional measurement system,and the change in the velocity of the boundary layer before and after plasma actuation is measured 10 mm downstream from the center of the actuator(Δx+≈1253).Figures 11(a)-(c)compare the boundary layer velocity distributions measured by applying the SBP-PA with different pulse frequencies.According to the internal dimensions of the boundary layer(friction velocityuτand viscous length

    lv=v/uτ,wherevis the kinematic viscosity),when the SBPPA is not applied,the velocity distribution is dimensionless.However,the velocity distribution of the boundary layer changes after applying different pulse frequencies,and the velocities of the buffer layer and the logarithmic region decrease significantly,resulting in a velocity reduction.After applying the SBP-PA,the blowing jet interacts with the flow near the wall,and the reverse vortex induced by the SBP-PA forms a kind of‘virtual barrier’effect near the wall;consequently,the streamwise vortices and coherent structures near the wall of the turbulent flow are lifted away from the wall,and the frictional drag induced by the interaction between the sweep events of the streamwise vortices is attenuated.At the same time,the regeneration of the near-wall streamwise vortices is suppressed,decreasing the velocity in the buffer layer and the logarithmic region,and ultimately reaching the goal of TDR.

    Figure 11(d)compares the changes in the boundary layer velocity distribution under different pulse frequencies.When

    Figure 11.Comparison of the boundary layer velocity distributions at different pulse frequencies:(a)fp=50 Hz;(b)fp=100 Hz;(c)fp=200 Hz;(d)velocity distributions at different pulse frequencies.

    y+<125,applying the SBP-PA with different pulse frequencies triggers varying downward deviations in the velocity profile.As the pulse frequency continues to increase,the deviations in the velocities of the buffer layer and logarithmic region gradually increase;in other words,the viscous effect in the near-wall area gradually plays the dominant role in the inner area of the TBL.When 125<y+<400,the effects of different pulse frequencies are only marginally different,and the velocity distributions under the three pulse frequencies tend to be consistent.Wheny+>400,the velocity profiles of the boundary layer with and without plasma actuation essentially coincide,which indicates that the control effect of the SBP-PA is basically lost.The following conclusions can be drawn.(1)At different pulse frequencies,the TDR effect of the SBP-PA varies considerably,and there is an optimal pulse frequency.(2)The SBP-PA has a limited impact area,and the current actuator configuration affects mainly the structure and viscous effect of streamwise vortices near the wall in the viscous sublayer,buffer layer,and logarithmic layer.With increasing distance from the wall,the kinetic energy effect of the SBP-PA gradually weakens,making it difficult to affect the structure of large-scale vortices in the outer region of the turbulent flow.

    Table 2 shows a comprehensive comparison of the distributions of various parameters of the boundary layer under different pulse frequencies.As the pulse frequency gradually increases,uτigradually decreases,leading to a gradual increase in the TDR rate.When the pulse frequencyfp=200 Hz(f+=2πfpd/U∞=7.196),the optimal local TDR rateRmax=16.97%.This also means that,whenf+=7.196,the jet frequency induced by the SBP-PA is closer to the natural frequency of the external mainstream;that is,the reverse vortex generated by a single pulse of the jet blowing through the slot is close to the spatial scale of the streamwise vortices near the wall,thereby producing the optimal TDR effect.It is worth noting that in this experiment,to prevent electromagnetic interference and protect the measurement and control equipment,the measuring position of the hot-wire probe in the boundary layer isΔx+≈1253 mm downstream of the plasma actuator.It can be inferred that the actual TDR rate measured in the near-wall area above the SBP-PA or at a distance less thanΔx+downstream of the SBP-PA will be greater than the TDR rate measured in the current position.Future experiments will employ noncontact technology to carry out near-wall turbulence measurements.

    Table 2.Parameters measured under different actuation frequencies at the measuring point.

    Figures 12(a)-(c)compare the velocity fluctuation distributions in the TBL measured by applying the SBP-PA with different pulse frequencies.Using the external dimensions of the boundary layer without plasma actuation,the velocity

    Figure 12.Distributions of the flow direction velocity fluctuations at different pulse frequencies:(a)fp=50 Hz;(b)fp=100 Hz;(c)fp=200 Hz;(d)velocity fluctuations with different fp.

    fluctuation is dimensionless.In the range of 0<y/δ<0.0014,the three pulse frequencies all reduce the velocity fluctuation to varying degrees,which may be due to the SBP-PA uplifting the turbulent near-wall streamwise vortices and coherent structures,thereby weakening the interaction between the sweep events of streamwise vortices and the wall and reducing the velocity fluctuation in the near-wall region(buffer layer and logarithmic region).However,when 0.0014<y/δ<0.56,the velocity fluctuation increases to varying degrees,and the pulse frequencyfp=200 Hz greatly increases the velocity fluctuation.This may be because the SBP-PA control effect increases with rising pulse frequency,and the number of vortices generated by the coupling between the SBP-PA and the main flow increases,which enhance both the velocity fluctuation and the intensity of turbulence.Wheny/δ>0.56,the velocity fluctuation distribution basically tends to be the same;this trend is similar to the experimental results in the flow direction arrangement of the slot blowing TDR strategy[18].This comparison also shows that the pulse frequencyfp=200 Hz(f+=2πfpd/U∞=7.196)corresponds to the optimal TDR effect for the current flow.With increasing distance from the wall,the control effect of the SBP-PA gradually weakens,which further confirms that the effect of the SBP-PA is limited to only the buffer layer and logarithmic layer near the wall,whereas the plasma actuator is still unable to exert an influence on the structure of large-scale vortices in the outer region of turbulence flow.Figure 12(d)comprehensively compares the distributions of the velocity fluctuation in the boundary layer under different pulse frequencies,indicating that within the range of 0.0014<y/δ<0.56,a greater pulse frequency corresponds to stronger pulse enhancement.The slot blowing strategy increases the amplitude of the velocity fluctuation in the outer region but reduces the velocity fluctuations near the wall.This result is similar to the outcomes reported in the literature[15].

    The increase of the pulse frequency will lead to more collisions between ions and electrons,which will increase the interaction between the induced vortices and the main flow,which may improve the flow control effectiveness.In this experiment,the drag reduction effect is also enhanced with the increase of the pulse frequency.Still,the ratio of TDR is not simply positively correlated with the pulse frequency.Reference[32]pointed out that DBD will reduce the vortex scale when the pulse frequency exceeds 200 Hz,and the vortex scale has a particular relationship with the TBL control.The control effect will be improved when the vortex scale is equivalent to the viscous bottom layer.Therefore,when the pulse frequency exceeds 200 Hz,the size of the vortex-induced by SBP-PA may be reduced,and then the drag reduction effect may be weakened.

    From the view of time-averaged analysis,the low-and high-velocity streak structure in the near-wall region of the TBL is distributed randomly,which increases the velocity gradient in the near-wall region and generates high friction drag.In this paper,a method of TDR using SBP-PA was proposed by suppressing the streak structure in the near-wall region.When the SBP-PA is actuated,a series of induced vortices perpendicular to the wall is generated.They evolve into large-scale streamwise vortices(LSSV)in the near-wall region of the boundary layer.

    The induced vortices perpendicular to the wall can inhibit the spanwise motion of the streak structure near the wall,thereby stabilizing the streak structure to act like riblets[34].Such riblets-like drag reduction effect can be called virtual riblets,which can suppress the interaction with the wall by inhibiting the development of the streamwise vortex and streak structure near the wall.Also,the virtual riblets can reduce the burst events of the vortex structure in the near-wall region of the TBL.

    It is worth noting that,to make the virtual riblets play a more significant role in inhibiting streamwise vortices and low-velocity streaks,the scale of induced vortex structure should be on the same level as the thickness of the viscous bottom layer.If the scale of the induced vortex structure is larger than the thickness of the viscous bottom layer,the flow effect of SBP-PA is similar to rendering a series of large-scale streamwise vortex instead of virtual riblets.At this time,the LSSV generated by SBP-PA will be dominant,attenuating the drag reduction effect and increasing the turbulence intensity.Conversely,when the virtual riblets effect dominates,the drag reduction effect will increase,and the turbulence intensity can be reduced.In the future,our group will analyze the near-wall turbulent structures to verify the drag reduction mechanism of SBP-PA by using high-speed PIV and microimaging technology.

    4.Conclusion

    In this experiment,a plasma flow control strategy with a slot blowing configuration is adopted,and the flow control effect is measured by means of a hot-wire anemometer and a data acquisition system.Experimental research is carried out on the TDR in a low-speed wind tunnel to study the influence of an SBP-PA on the streamwise vortices and coherent structures near the wall of the TBL.The analysis reveals the following conclusions:

    (1)Implementing the SBP-PA under different duty cycles indicates that the TDR effect is not positively related to the duty cycle.When the plasma actuation energy is higher than a certain threshold(50%),the injection of additional energy fails to improve the TDR effect and instead reduces the effect.

    (2)At different pulse frequencies,the TDR effect varies considerably,and there is an optimal pulse frequencyfp=200 Hz(f+=7.196).The SBP-PA causes the jet flow to interact with the flow near the wall,induces a reverse vortex,uplifts the streamwise vortices and coherent structures near the wall of the turbulent flow away from the wall,and alleviates the frictional effect induced by the interaction between the sweep events of the streamwise vortices and the wall.At the same time,the regeneration of the near-wall streamwise vortices is suppressed,and the goal of TDR is finally achieved.

    (3)The area over which plasma actuation exerts an influence is limited.The proposed SBP-PA affects mainly the structure of streamwise vortices near the wall in the buffer layer,viscous sublayer,and logarithmic region.With increasing distance from the wall,the kinetic effect of the plasma actuator gradually weakens,and it is difficult to affect the structure of large-scale vortices in the outer region of the turbulence flow.

    This work systematically studies the TDR effects of an SBP-PA under different electrical parameters,the outcomes of which are pertinent to the ongoing development of adaptive intelligent control.In future work,our research group will combine microscopy and stereoscopic PIV(SPIV)with an SBP-PA to study the influences of streamwise vortices in the inner area of the TBL near the wall and large-scale coherent structures in the outer region.Combined with the results of wind tunnel experiments and computational fluid dynamics(CFD)simulations,the velocity fluctuation quadrant decomposition method will be used to study the influences of the position,arrangement,pulse frequency,and working phase of the plasma actuator on the coherent structures of turbulent flow under different incoming flows,especially the control laws of sudden phenomena such as the ejection and sweep events of streamwise vortices.Furthermore,plasma actuator configurations that can more effectively control large-scale coherent structures in the outer region of turbulence will be explored,and the regulation mechanism between the flow structures in the inner and outer regions of turbulence flow will be obtained.Ultimately,a method to optimize the TDR effect,the mechanism by which the TDR effect will be generated,and the parameter influences will be scientifically explained.

    Acknowledgments

    The present study is supported by National Natural Science Foundation of China(Nos.61971345 and 12175177)and the Foundation for Key Laboratories of National Defense Science and Technology of China(No.614220120030810).

    猜你喜歡
    金元
    Flow control performance evaluation of a tri-electrode sliding discharge plasma actuator
    Electrical and aerodynamic characteristics of sliding discharge based on a microsecond pulsed plasma supply
    梅樊璟 金元奕 陳傲宇
    金代衛(wèi)紹王朝進(jìn)士輯錄——兼談其在金元文化傳承中的地位
    金元明時期北京地區(qū)古橋研究概述
    活力(2019年21期)2019-04-01 12:17:18
    金元浦:文創(chuàng)讓城市煥發(fā)新活力
    商周刊(2019年1期)2019-01-31 02:36:32
    轉(zhuǎn)彎
    遼宋西夏金元日常生活史研究概述
    金元餃
    漫畫腦力快車
    免费观看精品视频网站| 午夜精品在线福利| 免费在线观看影片大全网站| 欧美日韩国产mv在线观看视频| 男女之事视频高清在线观看| 老熟女久久久| 亚洲一区二区三区不卡视频| 伦理电影免费视频| 人成视频在线观看免费观看| 999久久久国产精品视频| 50天的宝宝边吃奶边哭怎么回事| 黄片大片在线免费观看| 啦啦啦 在线观看视频| 亚洲专区中文字幕在线| 女人爽到高潮嗷嗷叫在线视频| 老司机午夜福利在线观看视频| 桃红色精品国产亚洲av| 中文字幕另类日韩欧美亚洲嫩草| 日韩欧美三级三区| 国产野战对白在线观看| 国产成人免费观看mmmm| 久久久水蜜桃国产精品网| 嫩草影视91久久| a级毛片黄视频| 精品第一国产精品| 波多野结衣一区麻豆| 国产成人精品在线电影| 精品一区二区三卡| 国产视频一区二区在线看| 精品少妇久久久久久888优播| 热99re8久久精品国产| 亚洲av第一区精品v没综合| 精品乱码久久久久久99久播| 18禁观看日本| 欧美乱色亚洲激情| av一本久久久久| 99热国产这里只有精品6| 国产亚洲精品久久久久5区| 丰满的人妻完整版| 91字幕亚洲| 国产精品国产高清国产av | 高清在线国产一区| 成人黄色视频免费在线看| 亚洲精品国产精品久久久不卡| 两个人免费观看高清视频| 91成人精品电影| 成人av一区二区三区在线看| 自线自在国产av| 久久国产亚洲av麻豆专区| 国产成人免费观看mmmm| 欧美日韩亚洲综合一区二区三区_| 国产欧美亚洲国产| 午夜福利乱码中文字幕| 国产精品1区2区在线观看. | av国产精品久久久久影院| 中文欧美无线码| 日韩大码丰满熟妇| 交换朋友夫妻互换小说| 19禁男女啪啪无遮挡网站| 一a级毛片在线观看| 久久香蕉精品热| 不卡av一区二区三区| 人妻丰满熟妇av一区二区三区 | 老司机深夜福利视频在线观看| 黑人欧美特级aaaaaa片| 大型av网站在线播放| 人妻一区二区av| 亚洲欧美色中文字幕在线| 国产又色又爽无遮挡免费看| 国产高清国产精品国产三级| 女性生殖器流出的白浆| 亚洲第一av免费看| 国产成+人综合+亚洲专区| 国产亚洲欧美98| 午夜福利,免费看| 搡老岳熟女国产| 大型黄色视频在线免费观看| 国产乱人伦免费视频| 免费在线观看视频国产中文字幕亚洲| 久久久国产成人精品二区 | 久久精品人人爽人人爽视色| 国产精品国产av在线观看| 18禁观看日本| 男人舔女人的私密视频| 动漫黄色视频在线观看| 国产精品亚洲一级av第二区| 精品国产超薄肉色丝袜足j| 亚洲五月色婷婷综合| 亚洲欧美激情在线| 亚洲欧洲精品一区二区精品久久久| 午夜精品久久久久久毛片777| 久久天躁狠狠躁夜夜2o2o| 51午夜福利影视在线观看| 亚洲久久久国产精品| 亚洲美女黄片视频| 黄网站色视频无遮挡免费观看| 9191精品国产免费久久| 欧美久久黑人一区二区| 多毛熟女@视频| 夜夜夜夜夜久久久久| 久久精品国产a三级三级三级| 亚洲欧美日韩高清在线视频| 欧美日韩亚洲国产一区二区在线观看 | 久久香蕉激情| 亚洲精品久久成人aⅴ小说| 中文字幕制服av| 久久香蕉精品热| 国产视频一区二区在线看| 美国免费a级毛片| av一本久久久久| 日本wwww免费看| 欧美成狂野欧美在线观看| 亚洲av欧美aⅴ国产| 天天添夜夜摸| 99在线人妻在线中文字幕 | 国产成人免费观看mmmm| 91国产中文字幕| av超薄肉色丝袜交足视频| 精品人妻在线不人妻| 国产亚洲精品第一综合不卡| 身体一侧抽搐| 国产精品1区2区在线观看. | 999久久久精品免费观看国产| 下体分泌物呈黄色| 香蕉久久夜色| 热re99久久国产66热| av国产精品久久久久影院| 最新美女视频免费是黄的| 亚洲avbb在线观看| 久久久久视频综合| 午夜老司机福利片| 久9热在线精品视频| 18在线观看网站| 中文字幕另类日韩欧美亚洲嫩草| 欧美激情极品国产一区二区三区| 热99久久久久精品小说推荐| 久久国产乱子伦精品免费另类| 天天躁夜夜躁狠狠躁躁| 热re99久久精品国产66热6| 亚洲欧美精品综合一区二区三区| 欧美国产精品va在线观看不卡| 看黄色毛片网站| 看免费av毛片| 久久人人爽av亚洲精品天堂| 国产日韩一区二区三区精品不卡| 欧美成人免费av一区二区三区 | 久久国产亚洲av麻豆专区| 免费不卡黄色视频| x7x7x7水蜜桃| 亚洲成国产人片在线观看| 欧美中文综合在线视频| 久久国产精品人妻蜜桃| 国产成人av激情在线播放| 亚洲 欧美一区二区三区| 久久热在线av| 久久久久国产一级毛片高清牌| 欧美成人免费av一区二区三区 | 91精品国产国语对白视频| av国产精品久久久久影院| 欧美另类亚洲清纯唯美| 老司机午夜十八禁免费视频| 国产区一区二久久| 国产色视频综合| 视频区图区小说| 无人区码免费观看不卡| 国产精华一区二区三区| 国产黄色免费在线视频| 欧美在线一区亚洲| 国产成人精品无人区| 90打野战视频偷拍视频| 亚洲av日韩精品久久久久久密| 日韩三级视频一区二区三区| 国产精品二区激情视频| 精品久久久久久,| 少妇粗大呻吟视频| 丝瓜视频免费看黄片| 精品国产一区二区三区四区第35| 极品少妇高潮喷水抽搐| 男人舔女人的私密视频| 中文字幕制服av| 好男人电影高清在线观看| 欧美 日韩 精品 国产| 成人av一区二区三区在线看| 欧美 日韩 精品 国产| 男人舔女人的私密视频| 人人妻,人人澡人人爽秒播| 十八禁高潮呻吟视频| 国产又爽黄色视频| 午夜福利欧美成人| 亚洲,欧美精品.| 成年人午夜在线观看视频| 免费一级毛片在线播放高清视频 | 99久久人妻综合| 韩国精品一区二区三区| 老汉色∧v一级毛片| 久久婷婷成人综合色麻豆| 免费在线观看影片大全网站| 纯流量卡能插随身wifi吗| ponron亚洲| 韩国精品一区二区三区| av中文乱码字幕在线| 精品亚洲成a人片在线观看| 无遮挡黄片免费观看| 久久亚洲精品不卡| 一进一出好大好爽视频| 亚洲av成人av| 国产一区二区激情短视频| 最近最新中文字幕大全免费视频| 超色免费av| 很黄的视频免费| 欧美黄色淫秽网站| 99热网站在线观看| 欧美丝袜亚洲另类 | 亚洲自偷自拍图片 自拍| 一本大道久久a久久精品| 国产1区2区3区精品| 午夜福利影视在线免费观看| 婷婷成人精品国产| 在线观看一区二区三区激情| 80岁老熟妇乱子伦牲交| 日本黄色日本黄色录像| 久久人妻熟女aⅴ| 电影成人av| 成人手机av| 欧美日韩福利视频一区二区| 一级片'在线观看视频| av国产精品久久久久影院| 天天躁日日躁夜夜躁夜夜| 9191精品国产免费久久| 不卡av一区二区三区| 老司机深夜福利视频在线观看| 久久久精品区二区三区| 天天躁夜夜躁狠狠躁躁| 欧美 日韩 精品 国产| 久久人人97超碰香蕉20202| 精品亚洲成国产av| 免费少妇av软件| 欧美激情 高清一区二区三区| 国产欧美亚洲国产| 免费女性裸体啪啪无遮挡网站| 亚洲五月天丁香| 美女 人体艺术 gogo| 多毛熟女@视频| 久久国产亚洲av麻豆专区| 亚洲国产精品合色在线| 9热在线视频观看99| 日韩中文字幕欧美一区二区| 三上悠亚av全集在线观看| 亚洲av成人不卡在线观看播放网| 丁香欧美五月| 欧美午夜高清在线| 日韩熟女老妇一区二区性免费视频| 欧美日韩乱码在线| 久久久久国内视频| 18禁黄网站禁片午夜丰满| 大片电影免费在线观看免费| 久久久久久久久免费视频了| 色综合婷婷激情| 亚洲一区中文字幕在线| 精品一区二区三区视频在线观看免费 | 女性生殖器流出的白浆| 叶爱在线成人免费视频播放| 1024视频免费在线观看| 精品一区二区三区四区五区乱码| e午夜精品久久久久久久| 一级片'在线观看视频| 黄色a级毛片大全视频| 免费在线观看视频国产中文字幕亚洲| 精品电影一区二区在线| 久热这里只有精品99| 中文字幕人妻丝袜制服| 老司机午夜十八禁免费视频| 免费在线观看亚洲国产| 一进一出抽搐动态| 中文字幕色久视频| 国产激情欧美一区二区| 国产精品亚洲一级av第二区| 日韩成人在线观看一区二区三区| 91麻豆精品激情在线观看国产 | 精品亚洲成a人片在线观看| 午夜福利一区二区在线看| 狠狠狠狠99中文字幕| 极品人妻少妇av视频| 看片在线看免费视频| 免费日韩欧美在线观看| 精品国产超薄肉色丝袜足j| 在线观看午夜福利视频| 最新的欧美精品一区二区| 中文欧美无线码| www.精华液| 内射极品少妇av片p| 国产精品综合久久久久久久免费| 五月伊人婷婷丁香| 黄色视频,在线免费观看| 欧美色视频一区免费| 日本三级黄在线观看| 亚洲一区二区三区不卡视频| 成人特级黄色片久久久久久久| 成人三级黄色视频| 国产精品一区二区三区四区久久| 99久久无色码亚洲精品果冻| 法律面前人人平等表现在哪些方面| 黑人欧美特级aaaaaa片| 色老头精品视频在线观看| ponron亚洲| 国产主播在线观看一区二区| a级毛片a级免费在线| 桃红色精品国产亚洲av| 国产v大片淫在线免费观看| 美女被艹到高潮喷水动态| 国产精品98久久久久久宅男小说| 午夜福利在线在线| 午夜福利欧美成人| 国产精品女同一区二区软件 | 在线国产一区二区在线| 动漫黄色视频在线观看| 国产一区在线观看成人免费| 亚洲国产欧美人成| 欧美xxxx黑人xx丫x性爽| 国产精品亚洲av一区麻豆| 看免费av毛片| 精品不卡国产一区二区三区| 亚洲美女黄片视频| 国产成人福利小说| 国产成人影院久久av| 欧美成人免费av一区二区三区| 精品久久久久久久人妻蜜臀av| 老汉色∧v一级毛片| avwww免费| 在线免费观看的www视频| 亚洲五月天丁香| 亚洲第一电影网av| 亚洲av二区三区四区| 美女免费视频网站| 国产伦人伦偷精品视频| 国产成人欧美在线观看| 脱女人内裤的视频| 搡老熟女国产l中国老女人| 身体一侧抽搐| 黄色女人牲交| 日韩欧美精品免费久久 | 桃色一区二区三区在线观看| 国产探花在线观看一区二区| 黄色日韩在线| 此物有八面人人有两片| 亚洲熟妇中文字幕五十中出| 午夜免费激情av| 日本三级黄在线观看| 老熟妇乱子伦视频在线观看| 久久婷婷人人爽人人干人人爱| 国产高清有码在线观看视频| 欧美中文日本在线观看视频| 91九色精品人成在线观看| 亚洲成人中文字幕在线播放| 免费看美女性在线毛片视频| 1024手机看黄色片| 亚洲av成人精品一区久久| 1000部很黄的大片| 综合色av麻豆| 欧美黄色淫秽网站| netflix在线观看网站| 亚洲国产日韩欧美精品在线观看 | 淫秽高清视频在线观看| 欧美色视频一区免费| 中文字幕精品亚洲无线码一区| 高潮久久久久久久久久久不卡| 9191精品国产免费久久| 欧美区成人在线视频| 午夜a级毛片| 欧美成人性av电影在线观看| 久久久国产成人精品二区| netflix在线观看网站| 伊人久久大香线蕉亚洲五| 国产精品三级大全| 精品熟女少妇八av免费久了| 免费观看精品视频网站| 身体一侧抽搐| 久久久久久国产a免费观看| 中文字幕精品亚洲无线码一区| 中文字幕高清在线视频| 波多野结衣高清无吗| 久久久国产精品麻豆| 久9热在线精品视频| 夜夜躁狠狠躁天天躁| 国产极品精品免费视频能看的| 18禁国产床啪视频网站| 免费看十八禁软件| 国产精品 欧美亚洲| 亚洲av二区三区四区| 国产伦精品一区二区三区视频9 | 日韩亚洲欧美综合| 久久久久久久亚洲中文字幕 | 一个人看视频在线观看www免费 | 1024手机看黄色片| 怎么达到女性高潮| 99久久久亚洲精品蜜臀av| www.999成人在线观看| 久久九九热精品免费| 国产色爽女视频免费观看| 精品国产亚洲在线| 国产精品一区二区三区四区久久| 国产亚洲精品久久久com| 香蕉久久夜色| 午夜精品久久久久久毛片777| 亚洲欧美日韩高清在线视频| 身体一侧抽搐| 国产一区二区激情短视频| 狂野欧美白嫩少妇大欣赏| 在线国产一区二区在线| 免费看a级黄色片| 欧美一区二区亚洲| 国产v大片淫在线免费观看| 成人国产综合亚洲| 国产高清激情床上av| 在线观看午夜福利视频| 九色国产91popny在线| 我的老师免费观看完整版| 国产成人福利小说| 亚洲人成伊人成综合网2020| 好看av亚洲va欧美ⅴa在| 91av网一区二区| 久久性视频一级片| 在线国产一区二区在线| 亚洲美女黄片视频| 搡老妇女老女人老熟妇| 欧美极品一区二区三区四区| 亚洲人成网站在线播放欧美日韩| 一二三四社区在线视频社区8| aaaaa片日本免费| h日本视频在线播放| 国产精品一区二区三区四区久久| 日本撒尿小便嘘嘘汇集6| 午夜两性在线视频| 亚洲熟妇熟女久久| 99热6这里只有精品| 在线十欧美十亚洲十日本专区| 国产乱人视频| 国产三级中文精品| 欧美成人性av电影在线观看| xxx96com| 婷婷精品国产亚洲av在线| 夜夜爽天天搞| 啦啦啦韩国在线观看视频| 在线视频色国产色| 国产精品久久视频播放| 看片在线看免费视频| 国产一区二区在线观看日韩 | 波野结衣二区三区在线 | 亚洲黑人精品在线| 制服人妻中文乱码| 亚洲 国产 在线| 免费看美女性在线毛片视频| 成人18禁在线播放| 丰满人妻熟妇乱又伦精品不卡| 给我免费播放毛片高清在线观看| 丁香六月欧美| 国产男靠女视频免费网站| 少妇裸体淫交视频免费看高清| tocl精华| 成人一区二区视频在线观看| 国产午夜精品久久久久久一区二区三区 | 欧美最黄视频在线播放免费| 欧美午夜高清在线| 国产中年淑女户外野战色| 色老头精品视频在线观看| 久久精品影院6| 亚洲avbb在线观看| 日本免费一区二区三区高清不卡| 免费人成视频x8x8入口观看| 国产成人av激情在线播放| 日本五十路高清| 老鸭窝网址在线观看| 成人一区二区视频在线观看| 亚洲成人精品中文字幕电影| 久久久久国内视频| 中文字幕高清在线视频| 少妇丰满av| 别揉我奶头~嗯~啊~动态视频| 国产午夜精品论理片| 高清日韩中文字幕在线| x7x7x7水蜜桃| 国产精品98久久久久久宅男小说| 变态另类丝袜制服| 亚洲欧美精品综合久久99| 亚洲精品在线美女| 日韩欧美免费精品| 九色国产91popny在线| 国产野战对白在线观看| 午夜老司机福利剧场| netflix在线观看网站| 久久精品国产亚洲av香蕉五月| 日韩有码中文字幕| 国产精品一区二区三区四区免费观看 | 欧美日韩瑟瑟在线播放| 十八禁网站免费在线| 亚洲第一电影网av| 久久人人精品亚洲av| 成年人黄色毛片网站| 欧美不卡视频在线免费观看| 国产一区二区三区在线臀色熟女| 99热这里只有是精品50| 国产不卡一卡二| 热99re8久久精品国产| 在线天堂最新版资源| 久久亚洲真实| 在线十欧美十亚洲十日本专区| 欧美日韩中文字幕国产精品一区二区三区| 香蕉久久夜色| 亚洲美女黄片视频| or卡值多少钱| 可以在线观看毛片的网站| 一区福利在线观看| 桃色一区二区三区在线观看| 老熟妇仑乱视频hdxx| 99riav亚洲国产免费| 岛国在线观看网站| 亚洲无线在线观看| 久久久久久久亚洲中文字幕 | 日本黄色视频三级网站网址| 午夜亚洲福利在线播放| 欧美最黄视频在线播放免费| 国产精品久久久人人做人人爽| 国产精品亚洲av一区麻豆| 午夜福利高清视频| 亚洲av免费高清在线观看| 亚洲av二区三区四区| 狠狠狠狠99中文字幕| 搡老熟女国产l中国老女人| 88av欧美| 法律面前人人平等表现在哪些方面| 欧美又色又爽又黄视频| 禁无遮挡网站| 国产黄片美女视频| 香蕉丝袜av| 一区二区三区国产精品乱码| 中文字幕人妻丝袜一区二区| 亚洲av电影不卡..在线观看| 男女午夜视频在线观看| 国产成人aa在线观看| 亚洲精品在线美女| 精品久久久久久久末码| 亚洲一区二区三区不卡视频| 亚洲国产欧美人成| 国产精品一区二区免费欧美| 日韩欧美国产一区二区入口| 欧美乱妇无乱码| 亚洲av二区三区四区| 极品教师在线免费播放| 成人一区二区视频在线观看| 国产精品,欧美在线| 熟女少妇亚洲综合色aaa.| 日本与韩国留学比较| 一进一出抽搐动态| 亚洲精品在线美女| 日韩高清综合在线| 又紧又爽又黄一区二区| 国产精品久久视频播放| av福利片在线观看| 最近最新中文字幕大全免费视频| 中文资源天堂在线| 久久久精品欧美日韩精品| 丁香六月欧美| 国产国拍精品亚洲av在线观看 | 欧美一级a爱片免费观看看| 婷婷六月久久综合丁香| 99久久综合精品五月天人人| 亚洲真实伦在线观看| 国产麻豆成人av免费视频| 不卡一级毛片| 一本久久中文字幕| av天堂中文字幕网| 亚洲av熟女| 国产探花极品一区二区| 国产精品自产拍在线观看55亚洲| 啦啦啦韩国在线观看视频| 国产成年人精品一区二区| 丰满人妻熟妇乱又伦精品不卡| 国产真人三级小视频在线观看| 日本黄色片子视频| 3wmmmm亚洲av在线观看| 日本撒尿小便嘘嘘汇集6| 精品电影一区二区在线| 久久性视频一级片| 国产视频一区二区在线看| 欧美成人一区二区免费高清观看| 日韩精品中文字幕看吧| 亚洲欧美一区二区三区黑人| 国产午夜精品论理片| 俄罗斯特黄特色一大片| 人人妻人人看人人澡| 琪琪午夜伦伦电影理论片6080| 精品久久久久久久久久免费视频| 色综合亚洲欧美另类图片| 色尼玛亚洲综合影院| 国产精品香港三级国产av潘金莲| 成年女人毛片免费观看观看9| 午夜福利18| 性欧美人与动物交配| 村上凉子中文字幕在线| 亚洲欧美激情综合另类| 欧美日本亚洲视频在线播放| 成人av在线播放网站| 91久久精品电影网| 亚洲最大成人手机在线| 日本黄大片高清| 好男人在线观看高清免费视频| 老汉色∧v一级毛片| 国产欧美日韩一区二区精品| 人妻夜夜爽99麻豆av| 精品欧美国产一区二区三| 亚洲成人久久性| 欧美日韩福利视频一区二区| 久久国产精品影院| av在线天堂中文字幕| 可以在线观看毛片的网站|