• <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)彎
    遼宋西夏金元日常生活史研究概述
    金元餃
    漫畫腦力快車
    天堂√8在线中文| 精品不卡国产一区二区三区| 香蕉丝袜av| 身体一侧抽搐| 久久久久久久久中文| www国产在线视频色| 俺也久久电影网| 18禁美女被吸乳视频| 老汉色av国产亚洲站长工具| 久久国产精品人妻蜜桃| 色播亚洲综合网| 好男人电影高清在线观看| av黄色大香蕉| 国产亚洲av嫩草精品影院| 亚洲国产精品sss在线观看| 国产激情偷乱视频一区二区| 欧美成人一区二区免费高清观看 | 色噜噜av男人的天堂激情| 国产极品精品免费视频能看的| 成人高潮视频无遮挡免费网站| 精品国内亚洲2022精品成人| 啦啦啦观看免费观看视频高清| 啦啦啦韩国在线观看视频| 长腿黑丝高跟| 国产一级毛片七仙女欲春2| 亚洲在线自拍视频| 后天国语完整版免费观看| 琪琪午夜伦伦电影理论片6080| 男女午夜视频在线观看| 熟妇人妻久久中文字幕3abv| 久久精品亚洲精品国产色婷小说| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利视频1000在线观看| 午夜福利免费观看在线| 精品国产三级普通话版| 夜夜躁狠狠躁天天躁| 亚洲五月天丁香| 国产精品香港三级国产av潘金莲| 亚洲专区中文字幕在线| 成年版毛片免费区| 色老头精品视频在线观看| 欧美黄色淫秽网站| 久久这里只有精品中国| 国产精品乱码一区二三区的特点| 国产成人啪精品午夜网站| 国产真实乱freesex| 国产精品一区二区三区四区久久| 午夜a级毛片| 神马国产精品三级电影在线观看| 精品福利观看| 麻豆一二三区av精品| 亚洲天堂国产精品一区在线| 亚洲avbb在线观看| 欧美成人免费av一区二区三区| 久久99热这里只有精品18| 老熟妇乱子伦视频在线观看| 欧美最黄视频在线播放免费| 国产成人aa在线观看| 日韩欧美三级三区| 九九在线视频观看精品| 黄色日韩在线| 欧美最黄视频在线播放免费| 天天躁日日操中文字幕| 亚洲国产欧洲综合997久久,| 亚洲精品一区av在线观看| 全区人妻精品视频| 一区二区三区国产精品乱码| 黑人欧美特级aaaaaa片| 亚洲欧洲精品一区二区精品久久久| 中文字幕精品亚洲无线码一区| 黄片大片在线免费观看| 麻豆国产av国片精品| 五月玫瑰六月丁香| 亚洲真实伦在线观看| 无限看片的www在线观看| 99久国产av精品| 偷拍熟女少妇极品色| 又粗又爽又猛毛片免费看| 网址你懂的国产日韩在线| 日日摸夜夜添夜夜添小说| 色av中文字幕| 国产成人av教育| 日韩高清综合在线| 欧美日韩乱码在线| 日韩 欧美 亚洲 中文字幕| 亚洲 国产 在线| 国语自产精品视频在线第100页| 午夜激情福利司机影院| 九九热线精品视视频播放| 亚洲自拍偷在线| 精品人妻1区二区| 午夜影院日韩av| 12—13女人毛片做爰片一| 久久久久久人人人人人| 一本精品99久久精品77| 一个人免费在线观看电影 | 久久亚洲真实| 欧美成人一区二区免费高清观看 | 亚洲第一欧美日韩一区二区三区| 中文字幕精品亚洲无线码一区| 别揉我奶头~嗯~啊~动态视频| 精品一区二区三区视频在线 | 国产黄a三级三级三级人| 久久精品影院6| 国产午夜精品久久久久久| 麻豆成人av在线观看| 搡老熟女国产l中国老女人| 欧美极品一区二区三区四区| 国产黄色小视频在线观看| 午夜a级毛片| 亚洲国产看品久久| 曰老女人黄片| 人妻丰满熟妇av一区二区三区| 国产成人欧美在线观看| 精品国产乱码久久久久久男人| 一级a爱片免费观看的视频| 啦啦啦免费观看视频1| 国产主播在线观看一区二区| 欧美午夜高清在线| 99热只有精品国产| 亚洲精品美女久久av网站| av在线天堂中文字幕| 久久久久久大精品| 九九久久精品国产亚洲av麻豆 | 欧美精品啪啪一区二区三区| 国产成人欧美在线观看| 色精品久久人妻99蜜桃| 亚洲av第一区精品v没综合| 久久久久久大精品| 女生性感内裤真人,穿戴方法视频| 久久中文看片网| 国产精品久久久久久精品电影| www日本黄色视频网| 久久精品夜夜夜夜夜久久蜜豆| 18禁黄网站禁片免费观看直播| 精品人妻1区二区| 黄色 视频免费看| 女同久久另类99精品国产91| 亚洲av成人一区二区三| 欧美一级毛片孕妇| 国产成人av教育| 欧美3d第一页| 蜜桃久久精品国产亚洲av| 亚洲国产精品久久男人天堂| 亚洲天堂国产精品一区在线| 国产精品久久电影中文字幕| 久久九九热精品免费| 51午夜福利影视在线观看| 90打野战视频偷拍视频| 午夜福利欧美成人| 无遮挡黄片免费观看| 精华霜和精华液先用哪个| 男女午夜视频在线观看| 我的老师免费观看完整版| 99国产极品粉嫩在线观看| 少妇人妻一区二区三区视频| 热99在线观看视频| 国产欧美日韩精品亚洲av| or卡值多少钱| 99国产精品一区二区蜜桃av| 日韩高清综合在线| 欧美日韩一级在线毛片| 哪里可以看免费的av片| 中文字幕人妻丝袜一区二区| 18禁黄网站禁片免费观看直播| ponron亚洲| av在线蜜桃| 亚洲自偷自拍图片 自拍| 小蜜桃在线观看免费完整版高清| 中文字幕久久专区| 99久久精品国产亚洲精品| 午夜福利视频1000在线观看| 熟女人妻精品中文字幕| 国产高清有码在线观看视频| 男插女下体视频免费在线播放| 午夜精品久久久久久毛片777| 夜夜夜夜夜久久久久| 97超级碰碰碰精品色视频在线观看| 熟女少妇亚洲综合色aaa.| 国产综合懂色| 国产探花在线观看一区二区| 亚洲在线观看片| 九色国产91popny在线| 中出人妻视频一区二区| 可以在线观看的亚洲视频| 我要搜黄色片| 国产一区二区在线av高清观看| 长腿黑丝高跟| 综合色av麻豆| 久久中文字幕人妻熟女| 日日干狠狠操夜夜爽| 欧美在线一区亚洲| 国产成人影院久久av| 久久天躁狠狠躁夜夜2o2o| 色播亚洲综合网| 高清毛片免费观看视频网站| 午夜福利在线观看吧| 日韩欧美一区二区三区在线观看| 欧美日韩黄片免| 欧美最黄视频在线播放免费| 免费在线观看日本一区| 欧美乱色亚洲激情| 成人特级黄色片久久久久久久| 此物有八面人人有两片| 欧美日本亚洲视频在线播放| 三级国产精品欧美在线观看 | 亚洲国产精品999在线| 欧美又色又爽又黄视频| 国产伦一二天堂av在线观看| 久久久久国内视频| 国产精品久久久久久人妻精品电影| 男插女下体视频免费在线播放| 黄色片一级片一级黄色片| 欧美极品一区二区三区四区| 国产人伦9x9x在线观看| 中文字幕高清在线视频| 国产精品99久久99久久久不卡| 成人国产一区最新在线观看| 欧美av亚洲av综合av国产av| 网址你懂的国产日韩在线| 国产精品1区2区在线观看.| 黑人操中国人逼视频| 九色国产91popny在线| 制服人妻中文乱码| 狂野欧美激情性xxxx| 国产精品美女特级片免费视频播放器 | 91av网一区二区| 在线视频色国产色| 一级a爱片免费观看的视频| 久久人人精品亚洲av| 午夜免费观看网址| 国产午夜精品久久久久久| av视频在线观看入口| 亚洲色图 男人天堂 中文字幕| 久久久久性生活片| 国产在线精品亚洲第一网站| 亚洲无线观看免费| 久久午夜综合久久蜜桃| 日韩欧美在线二视频| 色av中文字幕| 在线看三级毛片| 一本久久中文字幕| 欧美日韩一级在线毛片| 日本黄色片子视频| 亚洲av片天天在线观看| 黄色成人免费大全| 国产成人系列免费观看| 色哟哟哟哟哟哟| 人人妻,人人澡人人爽秒播| 男女午夜视频在线观看| 国产主播在线观看一区二区| 亚洲熟妇中文字幕五十中出| 成在线人永久免费视频| 亚洲国产欧美一区二区综合| 亚洲色图 男人天堂 中文字幕| 国产日本99.免费观看| www.精华液| 久久草成人影院| 可以在线观看的亚洲视频| 欧美色视频一区免费| 亚洲精品乱码久久久v下载方式 | 成人无遮挡网站| 桃红色精品国产亚洲av| 又爽又黄无遮挡网站| 久久久精品大字幕| 一个人看的www免费观看视频| 国产极品精品免费视频能看的| 国产一区二区三区视频了| 亚洲精品一区av在线观看| 日韩精品青青久久久久久| 一进一出好大好爽视频| 欧美性猛交╳xxx乱大交人| 亚洲五月天丁香| av女优亚洲男人天堂 | 在线观看66精品国产| 婷婷精品国产亚洲av在线| 国产精品免费一区二区三区在线| 露出奶头的视频| 国产亚洲av嫩草精品影院| 国产精品永久免费网站| 又黄又粗又硬又大视频| 动漫黄色视频在线观看| 欧美日韩综合久久久久久 | 高清毛片免费观看视频网站| 99久久精品国产亚洲精品| 非洲黑人性xxxx精品又粗又长| 欧美成人免费av一区二区三区| 婷婷丁香在线五月| 老司机福利观看| 日日摸夜夜添夜夜添小说| 欧美高清成人免费视频www| 一个人看的www免费观看视频| 国产乱人伦免费视频| av天堂在线播放| АⅤ资源中文在线天堂| 国产精品久久视频播放| av片东京热男人的天堂| 最近最新中文字幕大全电影3| 亚洲午夜精品一区,二区,三区| 精品久久久久久,| 久久中文看片网| 舔av片在线| tocl精华| 国产成人一区二区三区免费视频网站| 高清在线国产一区| 亚洲 欧美一区二区三区| 少妇的逼水好多| 少妇熟女aⅴ在线视频| 两个人的视频大全免费| 一级a爱片免费观看的视频| 中文资源天堂在线| 国产午夜福利久久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av成人av| 日本 欧美在线| aaaaa片日本免费| 亚洲在线自拍视频| 国产乱人伦免费视频| 黄色 视频免费看| 在线观看一区二区三区| 国产伦人伦偷精品视频| 狂野欧美白嫩少妇大欣赏| 亚洲在线观看片| 欧美日韩瑟瑟在线播放| 色哟哟哟哟哟哟| 丰满的人妻完整版| 99国产综合亚洲精品| 男女视频在线观看网站免费| 特大巨黑吊av在线直播| 国产伦一二天堂av在线观看| 久久久久亚洲av毛片大全| 日韩三级视频一区二区三区| 日韩 欧美 亚洲 中文字幕| 女生性感内裤真人,穿戴方法视频| 俄罗斯特黄特色一大片| 黄片大片在线免费观看| 成人亚洲精品av一区二区| 亚洲自偷自拍图片 自拍| 国产极品精品免费视频能看的| 日本免费一区二区三区高清不卡| 美女黄网站色视频| 18禁裸乳无遮挡免费网站照片| 国产精品永久免费网站| 黄色日韩在线| 一区福利在线观看| 男女之事视频高清在线观看| 中文字幕精品亚洲无线码一区| 99久久精品一区二区三区| 日本黄大片高清| 国产一区二区在线av高清观看| 久久精品国产清高在天天线| 久久中文字幕人妻熟女| 成人无遮挡网站| 国产精品永久免费网站| 国产精品亚洲美女久久久| 国产野战对白在线观看| 精品人妻1区二区| 国产精品久久久久久精品电影| 国产精品一区二区三区四区久久| 中亚洲国语对白在线视频| 三级毛片av免费| 国产人伦9x9x在线观看| 99热精品在线国产| 国产毛片a区久久久久| 久久天堂一区二区三区四区| 美女高潮喷水抽搐中文字幕| 久久久精品欧美日韩精品| 免费在线观看成人毛片| 麻豆一二三区av精品| 国产精品久久久久久亚洲av鲁大| 亚洲性夜色夜夜综合| 国产97色在线日韩免费| 国产熟女xx| 此物有八面人人有两片| 亚洲18禁久久av| 亚洲成av人片免费观看| 亚洲在线观看片| 天天一区二区日本电影三级| 欧美3d第一页| 在线免费观看不下载黄p国产 | 亚洲av电影在线进入| 别揉我奶头~嗯~啊~动态视频| 亚洲人成网站高清观看| 精品不卡国产一区二区三区| 美女大奶头视频| 午夜亚洲福利在线播放| 免费在线观看成人毛片| 波多野结衣巨乳人妻| 麻豆久久精品国产亚洲av| 黄片小视频在线播放| 国产成人影院久久av| svipshipincom国产片| 手机成人av网站| 美女大奶头视频| 两个人视频免费观看高清| 国产精品久久视频播放| 久久久久国产一级毛片高清牌| 麻豆一二三区av精品| 免费无遮挡裸体视频| 五月伊人婷婷丁香| 90打野战视频偷拍视频| 久久九九热精品免费| 亚洲av中文字字幕乱码综合| 人人妻人人澡欧美一区二区| 1024香蕉在线观看| 欧美日韩乱码在线| 久久久国产成人精品二区| 成人永久免费在线观看视频| 亚洲自拍偷在线| av在线天堂中文字幕| 日韩有码中文字幕| 亚洲国产欧美网| 全区人妻精品视频| 99精品欧美一区二区三区四区| 久久性视频一级片| 午夜精品一区二区三区免费看| 嫩草影视91久久| 亚洲性夜色夜夜综合| 久久久久国产一级毛片高清牌| 国产精品国产高清国产av| 草草在线视频免费看| 国产欧美日韩一区二区精品| 亚洲中文日韩欧美视频| 2021天堂中文幕一二区在线观| 亚洲aⅴ乱码一区二区在线播放| 亚洲成a人片在线一区二区| 99久久精品一区二区三区| 精品国产亚洲在线| 国产欧美日韩一区二区三| 中国美女看黄片| 国产高清激情床上av| 国产精品爽爽va在线观看网站| 一级作爱视频免费观看| 99国产精品一区二区三区| 欧美日本亚洲视频在线播放| 岛国在线免费视频观看| 久久久久国内视频| 国产精品久久久久久亚洲av鲁大| 色综合站精品国产| 成人鲁丝片一二三区免费| 久久久久国产精品人妻aⅴ院| 中文字幕人成人乱码亚洲影| 啪啪无遮挡十八禁网站| 国产人伦9x9x在线观看| 亚洲中文字幕一区二区三区有码在线看 | 久久天堂一区二区三区四区| 男人和女人高潮做爰伦理| 国产伦在线观看视频一区| 一a级毛片在线观看| 欧美乱妇无乱码| 亚洲 欧美一区二区三区| 亚洲美女黄片视频| 啪啪无遮挡十八禁网站| 国产探花在线观看一区二区| 日韩欧美精品v在线| 免费观看的影片在线观看| 成人精品一区二区免费| 无遮挡黄片免费观看| 精品熟女少妇八av免费久了| 日韩精品中文字幕看吧| 亚洲av日韩精品久久久久久密| 国产aⅴ精品一区二区三区波| 巨乳人妻的诱惑在线观看| 男人和女人高潮做爰伦理| 亚洲国产日韩欧美精品在线观看 | 不卡一级毛片| 国产欧美日韩精品亚洲av| 成人鲁丝片一二三区免费| 99视频精品全部免费 在线 | 级片在线观看| 久久久久亚洲av毛片大全| 桃红色精品国产亚洲av| 国产熟女xx| 少妇丰满av| 在线观看免费午夜福利视频| 亚洲欧美日韩东京热| 丰满的人妻完整版| 国产成人一区二区三区免费视频网站| 老熟妇乱子伦视频在线观看| 久久香蕉国产精品| 亚洲av电影不卡..在线观看| 精品国产超薄肉色丝袜足j| 亚洲狠狠婷婷综合久久图片| 久久伊人香网站| 久久久久亚洲av毛片大全| 色精品久久人妻99蜜桃| 真人做人爱边吃奶动态| 久久精品aⅴ一区二区三区四区| 好男人在线观看高清免费视频| 免费看光身美女| 国产精品一区二区三区四区久久| 免费在线观看亚洲国产| 国产亚洲精品久久久久久毛片| 视频区欧美日本亚洲| 色尼玛亚洲综合影院| 午夜影院日韩av| 搡老妇女老女人老熟妇| 亚洲精品456在线播放app | 久久草成人影院| 动漫黄色视频在线观看| 久久久久久国产a免费观看| 少妇的丰满在线观看| 欧美日韩精品网址| 深夜精品福利| 亚洲欧美激情综合另类| 久久久久国产精品人妻aⅴ院| 美女黄网站色视频| 无人区码免费观看不卡| 小说图片视频综合网站| 久久精品综合一区二区三区| www国产在线视频色| 12—13女人毛片做爰片一| 精品久久久久久久末码| 亚洲中文日韩欧美视频| 黄片小视频在线播放| 婷婷丁香在线五月| 人妻夜夜爽99麻豆av| 老鸭窝网址在线观看| 欧美乱码精品一区二区三区| av国产免费在线观看| 国产欧美日韩一区二区精品| 在线十欧美十亚洲十日本专区| 一二三四社区在线视频社区8| 亚洲人成电影免费在线| 久久久久国产精品人妻aⅴ院| 琪琪午夜伦伦电影理论片6080| 午夜成年电影在线免费观看| 久久欧美精品欧美久久欧美| 69av精品久久久久久| 一区二区三区国产精品乱码| 亚洲片人在线观看| av女优亚洲男人天堂 | 国产高清视频在线观看网站| 婷婷丁香在线五月| 成人特级黄色片久久久久久久| 欧美性猛交╳xxx乱大交人| xxx96com| 亚洲成人精品中文字幕电影| 好男人在线观看高清免费视频| 亚洲专区国产一区二区| 露出奶头的视频| 两个人看的免费小视频| 神马国产精品三级电影在线观看| 成人无遮挡网站| 无限看片的www在线观看| 麻豆国产97在线/欧美| 欧美国产日韩亚洲一区| 免费看a级黄色片| 真人一进一出gif抽搐免费| 国产亚洲精品一区二区www| 男女床上黄色一级片免费看| 亚洲国产精品久久男人天堂| 亚洲国产欧美人成| 国产av一区在线观看免费| 亚洲美女黄片视频| 国产精品美女特级片免费视频播放器 | 久久久久国产精品人妻aⅴ院| 亚洲精品粉嫩美女一区| 精品国产乱子伦一区二区三区| 午夜a级毛片| 精品无人区乱码1区二区| 国产极品精品免费视频能看的| 国产高清有码在线观看视频| 黑人操中国人逼视频| 中文字幕人妻丝袜一区二区| 免费看日本二区| 97碰自拍视频| 免费看十八禁软件| 亚洲一区二区三区不卡视频| 琪琪午夜伦伦电影理论片6080| 黑人欧美特级aaaaaa片| 色综合亚洲欧美另类图片| 99久国产av精品| 久久欧美精品欧美久久欧美| 国产精华一区二区三区| 国产黄片美女视频| 亚洲七黄色美女视频| 欧美日韩国产亚洲二区| 成年女人永久免费观看视频| 欧美一级a爱片免费观看看| 一进一出好大好爽视频| 精品久久久久久成人av| 此物有八面人人有两片| 亚洲国产色片| 九九在线视频观看精品| 男女视频在线观看网站免费| 欧美绝顶高潮抽搐喷水| 99在线人妻在线中文字幕| 国产高清videossex| 国产成人系列免费观看| 国产精品美女特级片免费视频播放器 | 超碰成人久久| 成在线人永久免费视频| 天堂网av新在线| 成人国产一区最新在线观看| 男女视频在线观看网站免费| 黑人操中国人逼视频| 床上黄色一级片| 九色国产91popny在线| 一a级毛片在线观看| 国产97色在线日韩免费| 不卡一级毛片| 毛片女人毛片| 三级国产精品欧美在线观看 | 欧美乱色亚洲激情| av在线天堂中文字幕| 老司机深夜福利视频在线观看| 最近最新中文字幕大全免费视频| 超碰成人久久| 91久久精品国产一区二区成人 | 国产aⅴ精品一区二区三区波| 亚洲国产欧洲综合997久久,| 丁香六月欧美| 精品久久久久久久毛片微露脸| 搡老熟女国产l中国老女人|