• <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)彎
    遼宋西夏金元日常生活史研究概述
    金元餃
    漫畫腦力快車
    亚洲精品成人久久久久久| 天天一区二区日本电影三级| 亚洲欧美激情综合另类| 亚洲精华国产精华精| 国内毛片毛片毛片毛片毛片| 亚洲精品国产精品久久久不卡| 久久亚洲精品不卡| 久久精品人妻少妇| 亚洲专区国产一区二区| 亚洲精品乱码久久久v下载方式 | 日本免费a在线| 免费在线观看影片大全网站| 51国产日韩欧美| 美女黄网站色视频| 搡老熟女国产l中国老女人| 免费搜索国产男女视频| 国产极品精品免费视频能看的| 91麻豆av在线| 日日干狠狠操夜夜爽| 一区二区三区高清视频在线| 午夜久久久久精精品| 成年免费大片在线观看| 国产亚洲欧美98| 男人和女人高潮做爰伦理| 国产精品98久久久久久宅男小说| 色老头精品视频在线观看| 婷婷精品国产亚洲av| 亚洲精品在线美女| 国产爱豆传媒在线观看| 亚洲精品成人久久久久久| 精品乱码久久久久久99久播| 久久久久久国产a免费观看| 乱人视频在线观看| 精品乱码久久久久久99久播| 又紧又爽又黄一区二区| 色老头精品视频在线观看| 在线观看66精品国产| 91久久精品电影网| 嫩草影视91久久| 老熟妇仑乱视频hdxx| 免费看日本二区| 在线观看午夜福利视频| 真人做人爱边吃奶动态| 黄色女人牲交| 国产精品 国内视频| www.999成人在线观看| 免费高清视频大片| 天堂动漫精品| 免费电影在线观看免费观看| 国产精品久久久久久精品电影| 国产不卡一卡二| 美女免费视频网站| 精华霜和精华液先用哪个| 搡老岳熟女国产| avwww免费| 精品国产三级普通话版| 国产高清激情床上av| 国内揄拍国产精品人妻在线| 国产蜜桃级精品一区二区三区| 黄色女人牲交| 国产91精品成人一区二区三区| 久久久成人免费电影| 亚洲精品成人久久久久久| 少妇丰满av| 男女做爰动态图高潮gif福利片| eeuss影院久久| 国产伦精品一区二区三区视频9 | 看片在线看免费视频| 黑人欧美特级aaaaaa片| 丁香六月欧美| 午夜视频国产福利| 人人妻人人看人人澡| 国产av麻豆久久久久久久| 99久久精品一区二区三区| 少妇的逼水好多| 欧美日本视频| 成人高潮视频无遮挡免费网站| 少妇熟女aⅴ在线视频| 国内精品一区二区在线观看| 老司机深夜福利视频在线观看| 国产精品99久久久久久久久| 岛国在线观看网站| 免费无遮挡裸体视频| 一个人看的www免费观看视频| 国产乱人伦免费视频| av天堂中文字幕网| 欧美日韩黄片免| 悠悠久久av| av在线蜜桃| 国产午夜福利久久久久久| 两个人的视频大全免费| 日韩精品青青久久久久久| 最好的美女福利视频网| av国产免费在线观看| 久久亚洲精品不卡| 国产欧美日韩一区二区三| 国产精品一及| 欧美日韩中文字幕国产精品一区二区三区| 国产一区二区在线av高清观看| 九九热线精品视视频播放| 一夜夜www| www.www免费av| 在线观看免费午夜福利视频| 成人亚洲精品av一区二区| av黄色大香蕉| 99久久成人亚洲精品观看| 香蕉久久夜色| 波多野结衣巨乳人妻| 网址你懂的国产日韩在线| 亚洲成人精品中文字幕电影| 欧美日韩中文字幕国产精品一区二区三区| 黄片大片在线免费观看| 国内久久婷婷六月综合欲色啪| 伊人久久精品亚洲午夜| 欧美av亚洲av综合av国产av| 欧美一区二区国产精品久久精品| 男女视频在线观看网站免费| 国产97色在线日韩免费| 国产色婷婷99| 国产免费男女视频| 中文亚洲av片在线观看爽| 亚洲 国产 在线| eeuss影院久久| 午夜免费观看网址| 麻豆国产97在线/欧美| 亚洲七黄色美女视频| 成人性生交大片免费视频hd| 法律面前人人平等表现在哪些方面| 亚洲欧美日韩高清在线视频| 亚洲国产日韩欧美精品在线观看 | 叶爱在线成人免费视频播放| 男人和女人高潮做爰伦理| aaaaa片日本免费| 国产91精品成人一区二区三区| 免费高清视频大片| av欧美777| 他把我摸到了高潮在线观看| 他把我摸到了高潮在线观看| 国产真实伦视频高清在线观看 | 亚洲欧美一区二区三区黑人| 国产欧美日韩一区二区精品| 麻豆成人午夜福利视频| 又紧又爽又黄一区二区| 国产高清有码在线观看视频| 搡老妇女老女人老熟妇| 午夜激情欧美在线| 国产精品久久电影中文字幕| 国产亚洲精品综合一区在线观看| 亚洲av五月六月丁香网| 亚洲狠狠婷婷综合久久图片| 精品不卡国产一区二区三区| 无限看片的www在线观看| 国产精品美女特级片免费视频播放器| 51午夜福利影视在线观看| 日韩欧美一区二区三区在线观看| 久久这里只有精品中国| 综合色av麻豆| 国产日本99.免费观看| 亚洲第一电影网av| 国产成人a区在线观看| 国产探花极品一区二区| 最近在线观看免费完整版| 天堂√8在线中文| 丰满乱子伦码专区| 久久天躁狠狠躁夜夜2o2o| 中文字幕精品亚洲无线码一区| 老司机福利观看| 两个人看的免费小视频| 99久久九九国产精品国产免费| 国产乱人伦免费视频| 国产色婷婷99| 又粗又爽又猛毛片免费看| 欧美黑人欧美精品刺激| 91麻豆精品激情在线观看国产| av天堂中文字幕网| 别揉我奶头~嗯~啊~动态视频| 亚洲黑人精品在线| 一进一出抽搐动态| 国产精品日韩av在线免费观看| 老司机深夜福利视频在线观看| 少妇熟女aⅴ在线视频| 手机成人av网站| 成熟少妇高潮喷水视频| 精品久久久久久久末码| 精品乱码久久久久久99久播| 亚洲专区中文字幕在线| 中文字幕久久专区| 黄片大片在线免费观看| 国产精品亚洲美女久久久| 亚洲性夜色夜夜综合| 欧美乱码精品一区二区三区| 亚洲精品一区av在线观看| 国产午夜福利久久久久久| 国产精品亚洲av一区麻豆| 久久国产精品影院| 国产精品野战在线观看| 日韩亚洲欧美综合| 久久精品影院6| 午夜福利高清视频| 国产伦精品一区二区三区四那| 高清日韩中文字幕在线| 麻豆久久精品国产亚洲av| 精品一区二区三区视频在线 | 欧美成狂野欧美在线观看| 国产爱豆传媒在线观看| 国产成人av激情在线播放| 床上黄色一级片| 久久精品国产自在天天线| 又紧又爽又黄一区二区| 久久久久精品国产欧美久久久| 2021天堂中文幕一二区在线观| 午夜免费男女啪啪视频观看 | 久久国产精品人妻蜜桃| 国产精品三级大全| 一级毛片女人18水好多| 亚洲成a人片在线一区二区| a级毛片a级免费在线| 国产蜜桃级精品一区二区三区| 长腿黑丝高跟| 大型黄色视频在线免费观看| 亚洲av二区三区四区| 最新在线观看一区二区三区| 久久久久久久午夜电影| 亚洲无线在线观看| 国产黄色小视频在线观看| 中亚洲国语对白在线视频| 国产精华一区二区三区| 最近在线观看免费完整版| 一本一本综合久久| 久久久久亚洲av毛片大全| 一个人看视频在线观看www免费 | 18禁在线播放成人免费| 女人被狂操c到高潮| 亚洲国产高清在线一区二区三| 精品乱码久久久久久99久播| 少妇人妻精品综合一区二区 | 亚洲人与动物交配视频| 一本综合久久免费| avwww免费| 久久6这里有精品| 亚洲欧美日韩高清在线视频| 亚洲无线观看免费| 午夜福利欧美成人| 色播亚洲综合网| 久久久久精品国产欧美久久久| 99视频精品全部免费 在线| 少妇丰满av| 别揉我奶头~嗯~啊~动态视频| 欧美成人免费av一区二区三区| 国产在视频线在精品| 成人性生交大片免费视频hd| 神马国产精品三级电影在线观看| 怎么达到女性高潮| 男插女下体视频免费在线播放| 色吧在线观看| 99精品欧美一区二区三区四区| 亚洲美女视频黄频| 人人妻,人人澡人人爽秒播| 手机成人av网站| 欧美乱码精品一区二区三区| 免费在线观看日本一区| 精品人妻1区二区| 久久精品国产亚洲av涩爱 | 老汉色av国产亚洲站长工具| 男人舔奶头视频| 老司机午夜福利在线观看视频| 精品国产美女av久久久久小说| 国产高清有码在线观看视频| 国产97色在线日韩免费| 在线十欧美十亚洲十日本专区| 老熟妇仑乱视频hdxx| 亚洲美女黄片视频| 在线免费观看不下载黄p国产 | 国产毛片a区久久久久| 五月玫瑰六月丁香| 久久伊人香网站| 国产男靠女视频免费网站| 亚洲五月婷婷丁香| 成人国产综合亚洲| 国产精品女同一区二区软件 | 久久草成人影院| 91久久精品电影网| 搡老妇女老女人老熟妇| 欧美日本视频| 国产精品三级大全| 精品一区二区三区av网在线观看| 婷婷六月久久综合丁香| 亚洲人成伊人成综合网2020| 在线观看av片永久免费下载| 日韩高清综合在线| 亚洲色图av天堂| 深爱激情五月婷婷| 高潮久久久久久久久久久不卡| 少妇丰满av| 亚洲欧美激情综合另类| 国产高清视频在线播放一区| 久久午夜亚洲精品久久| 麻豆成人av在线观看| 免费看美女性在线毛片视频| 日日摸夜夜添夜夜添小说| 一区福利在线观看| 中出人妻视频一区二区| 欧美又色又爽又黄视频| 91字幕亚洲| 亚洲精华国产精华精| 亚洲一区高清亚洲精品| 老司机在亚洲福利影院| 床上黄色一级片| 久久人妻av系列| 日韩 欧美 亚洲 中文字幕| 亚洲人成伊人成综合网2020| 成人国产一区最新在线观看| 麻豆成人午夜福利视频| 在线天堂最新版资源| 国模一区二区三区四区视频| 黑人欧美特级aaaaaa片| 久久久成人免费电影| 婷婷丁香在线五月| 午夜福利高清视频| 欧美日韩一级在线毛片| 久久国产精品人妻蜜桃| 国产私拍福利视频在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲精品一区av在线观看| а√天堂www在线а√下载| 18+在线观看网站| 久久欧美精品欧美久久欧美| 欧美日韩一级在线毛片| 无限看片的www在线观看| 宅男免费午夜| 亚洲欧美精品综合久久99| 三级男女做爰猛烈吃奶摸视频| 日韩免费av在线播放| av天堂中文字幕网| 色尼玛亚洲综合影院| 欧洲精品卡2卡3卡4卡5卡区| 人人妻人人看人人澡| 韩国av一区二区三区四区| 尤物成人国产欧美一区二区三区| 又爽又黄无遮挡网站| 一本综合久久免费| 成人国产一区最新在线观看| 宅男免费午夜| 国产精品久久视频播放| 欧美中文日本在线观看视频| 看片在线看免费视频| 日本黄大片高清| 亚洲精品色激情综合| 亚洲欧美日韩高清专用| 观看免费一级毛片| 国产精品99久久99久久久不卡| 国产成年人精品一区二区| 91在线精品国自产拍蜜月 | 国产成人啪精品午夜网站| 欧洲精品卡2卡3卡4卡5卡区| 国内精品美女久久久久久| 中文字幕人成人乱码亚洲影| 亚洲电影在线观看av| 久久精品综合一区二区三区| 在线观看日韩欧美| 国产精品亚洲av一区麻豆| 91在线观看av| 国产伦在线观看视频一区| 天堂av国产一区二区熟女人妻| 精品一区二区三区视频在线 | 亚洲第一电影网av| 国产熟女xx| 老熟妇仑乱视频hdxx| 国内少妇人妻偷人精品xxx网站| 亚洲精品在线观看二区| 日韩免费av在线播放| 亚洲,欧美精品.| 性色av乱码一区二区三区2| 在线观看美女被高潮喷水网站 | 身体一侧抽搐| 国产精品综合久久久久久久免费| 日韩精品青青久久久久久| 神马国产精品三级电影在线观看| 变态另类丝袜制服| 久久久色成人| 国内精品一区二区在线观看| 欧美最黄视频在线播放免费| 在线免费观看不下载黄p国产 | 人人妻,人人澡人人爽秒播| 亚洲最大成人中文| 国产免费男女视频| 麻豆久久精品国产亚洲av| 亚洲国产欧洲综合997久久,| 成人无遮挡网站| 成人无遮挡网站| 婷婷丁香在线五月| 91九色精品人成在线观看| 欧美又色又爽又黄视频| 变态另类丝袜制服| 在线观看免费午夜福利视频| 国产成人啪精品午夜网站| 韩国av一区二区三区四区| 成人av一区二区三区在线看| 婷婷精品国产亚洲av在线| 午夜福利在线观看吧| 在线国产一区二区在线| 欧美黄色片欧美黄色片| 久久性视频一级片| 亚洲国产欧美网| 久久精品国产清高在天天线| 国产精华一区二区三区| 身体一侧抽搐| 国产精品久久久久久精品电影| 欧美xxxx黑人xx丫x性爽| 精品福利观看| 99精品欧美一区二区三区四区| 天堂网av新在线| 亚洲av成人精品一二三区| 极品教师在线视频| 国产极品天堂在线| 国产精品伦人一区二区| av在线观看视频网站免费| 国产视频内射| 国产爱豆传媒在线观看| 51国产日韩欧美| 插逼视频在线观看| 亚洲四区av| 亚洲色图av天堂| 丰满少妇做爰视频| 国产精品.久久久| 国产精品久久久久久av不卡| 久久精品国产亚洲av涩爱| 亚洲欧美日韩卡通动漫| 国产午夜精品论理片| 观看美女的网站| 91精品国产九色| 精品国产三级普通话版| 久久人人爽人人片av| 国产真实伦视频高清在线观看| 97精品久久久久久久久久精品| 男女边摸边吃奶| 91精品伊人久久大香线蕉| 黄色日韩在线| 国产黄片视频在线免费观看| 成人亚洲精品一区在线观看 | 国产在视频线精品| 欧美+日韩+精品| 久久人人爽人人爽人人片va| 国产精品一二三区在线看| 晚上一个人看的免费电影| 美女脱内裤让男人舔精品视频| 久久久久久久久久久免费av| 欧美丝袜亚洲另类| 观看免费一级毛片| 欧美bdsm另类| 国国产精品蜜臀av免费| 国产单亲对白刺激| 黄色一级大片看看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品第二区| 特大巨黑吊av在线直播| 天堂影院成人在线观看| 国产精品综合久久久久久久免费| 综合色av麻豆| videos熟女内射| 成年女人看的毛片在线观看| 床上黄色一级片| 99热这里只有精品一区| av网站免费在线观看视频 | 亚洲国产精品sss在线观看| 久久久久久久午夜电影| 别揉我奶头 嗯啊视频| 日韩精品青青久久久久久| 亚洲国产色片| 内射极品少妇av片p| 搡老乐熟女国产| 久久久a久久爽久久v久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日本色播在线视频| 秋霞在线观看毛片| 天堂中文最新版在线下载 | 有码 亚洲区| 日韩视频在线欧美| 久久久国产一区二区| 亚洲av不卡在线观看| 国产成人精品久久久久久| 淫秽高清视频在线观看| 精品久久久久久成人av| av网站免费在线观看视频 | 国产精品久久久久久久电影| 十八禁国产超污无遮挡网站| 久久鲁丝午夜福利片| 天天躁夜夜躁狠狠久久av| 成年女人在线观看亚洲视频 | 亚洲av成人精品一二三区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲四区av| 色尼玛亚洲综合影院| 久久久亚洲精品成人影院| 舔av片在线| 麻豆av噜噜一区二区三区| 大陆偷拍与自拍| 精品久久久久久久久av| 五月天丁香电影| 人妻制服诱惑在线中文字幕| 爱豆传媒免费全集在线观看| 国产乱来视频区| 国产日韩欧美在线精品| av.在线天堂| 免费播放大片免费观看视频在线观看| 看黄色毛片网站| 伦理电影大哥的女人| 久久精品夜夜夜夜夜久久蜜豆| 国产精品女同一区二区软件| 亚洲美女搞黄在线观看| 亚洲精华国产精华液的使用体验| 精品国产露脸久久av麻豆 | 老女人水多毛片| 亚洲最大成人中文| 久久久久久九九精品二区国产| 黄色一级大片看看| 天天躁夜夜躁狠狠久久av| 91在线精品国自产拍蜜月| 欧美高清性xxxxhd video| 亚洲国产日韩欧美精品在线观看| 久久精品人妻少妇| 99久久精品国产国产毛片| 国产伦在线观看视频一区| 国产不卡一卡二| av在线播放精品| 小蜜桃在线观看免费完整版高清| 亚州av有码| 欧美日韩亚洲高清精品| 少妇被粗大猛烈的视频| 丝瓜视频免费看黄片| 99久久精品热视频| 亚洲一区高清亚洲精品| 国产有黄有色有爽视频| 欧美xxxx黑人xx丫x性爽| 国产精品一区二区三区四区免费观看| 综合色av麻豆| 日韩在线高清观看一区二区三区| 爱豆传媒免费全集在线观看| 国产成人精品一,二区| 国产精品国产三级专区第一集| 非洲黑人性xxxx精品又粗又长| 22中文网久久字幕| 午夜激情欧美在线| av在线观看视频网站免费| 免费观看的影片在线观看| 亚洲欧美一区二区三区黑人 | 99久久精品国产国产毛片| 国产成人精品福利久久| 日韩av在线大香蕉| 日本欧美国产在线视频| 免费看不卡的av| 亚洲国产欧美在线一区| 十八禁国产超污无遮挡网站| 成人毛片60女人毛片免费| 毛片女人毛片| 国产白丝娇喘喷水9色精品| 3wmmmm亚洲av在线观看| 亚洲精品国产av成人精品| 欧美 日韩 精品 国产| 一级av片app| 十八禁网站网址无遮挡 | 久久久久九九精品影院| 午夜福利成人在线免费观看| 国产一区二区三区综合在线观看 | 色综合亚洲欧美另类图片| 成人毛片60女人毛片免费| 久久人人爽人人爽人人片va| 观看美女的网站| 尤物成人国产欧美一区二区三区| 三级男女做爰猛烈吃奶摸视频| 不卡视频在线观看欧美| 国产免费又黄又爽又色| 啦啦啦韩国在线观看视频| 纵有疾风起免费观看全集完整版 | 乱系列少妇在线播放| 禁无遮挡网站| 国产成人a∨麻豆精品| 免费人成在线观看视频色| 国产男人的电影天堂91| 在线免费观看不下载黄p国产| 日本熟妇午夜| 亚洲18禁久久av| 国产69精品久久久久777片| 国产一级毛片七仙女欲春2| 日日干狠狠操夜夜爽| 99久久精品一区二区三区| 亚洲自拍偷在线| 日本欧美国产在线视频| 久久久久国产网址| 有码 亚洲区| 别揉我奶头 嗯啊视频| 欧美xxⅹ黑人| 国产黄色小视频在线观看| 黄色配什么色好看| 国产精品爽爽va在线观看网站| 久久久久久久国产电影| 国产高潮美女av| 免费在线观看成人毛片| 精品不卡国产一区二区三区| 深爱激情五月婷婷| 国产在视频线在精品| 91aial.com中文字幕在线观看| 国产91av在线免费观看| av播播在线观看一区| 亚洲成人久久爱视频| 九九在线视频观看精品| 日韩成人av中文字幕在线观看| 看十八女毛片水多多多| 丰满乱子伦码专区| 免费黄频网站在线观看国产| 国产真实伦视频高清在线观看| 国产精品久久视频播放| 大话2 男鬼变身卡| 女人十人毛片免费观看3o分钟| 成人av在线播放网站|