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

    Influence of SPS casing treatment on axial flow compressor subjected to radial pressure distortion

    2017-11-20 12:07:21LiFnyuLiJunDongXuSunDkunSunXiofeng
    CHINESE JOURNAL OF AERONAUTICS 2017年2期

    Li Fnyu,Li Jun,Dong Xu,Sun Dkun,*,Sun Xiofeng

    aScience and Technology on Plasma Dynamics Laboratory,Air Force Engineering University,Xi’an 710038,ChinabCo-Innovation Center for Advanced Aero-Engine,Beihang University,Beijing 100083,China

    Influence of SPS casing treatment on axial flow compressor subjected to radial pressure distortion

    Li Fanyua,b,Li Juna,Dong Xub,Sun Dakunb,*,Sun Xiaofengb

    aScience and Technology on Plasma Dynamics Laboratory,Air Force Engineering University,Xi’an 710038,ChinabCo-Innovation Center for Advanced Aero-Engine,Beihang University,Beijing 100083,China

    Axial compressor;Casing treatment;Fan stability;Inlet distortion;Rotating stall

    The present work is about the stall margin enhancement ability of a kind of stall precursor-suppressed(SPS)casing treatment when fan/compressor suffers from a radial total pressure inlet distortion.Experimental researches are conducted on a low-speed compressor with and without SPS casing treatment under radial distorted inlet flow of different levels as well as uniform inlet flow.The distorted flow fields of different levels are generated by annular distortion flow generators of different heights.The characteristic curves under these conditions are measured and analyzed.The results show that the radial inlet distortion could cause a stall margin loss from 2%to 30%under different distorted levels.The SPS casing treatment could remedy this stall margin loss under small distortion level and only partly make up the stall margin loss caused by distortion in large level without leading to perceptible additional efficiency loss and obvious change of characteristic curves.The pre-stall behavior of the compressor is investigated to reveal the mechanism of this stall margin improvement ability of the SPS casing treatment.The results do show that this casing treatment delays the occurrence of rotating stall by weakening the pressure perturbations and suppressing the nonlinear amplification of the stall precursor waves in the compression system.

    1.Introduction

    Modern aircrafts,especially the military fighters,are often subjected to complex flight conditions.Their propulsion systems using air-breathing gas turbine engine often suffer from distorted inlet conditions.For example,the flow separation at the nacelle lip when planes take off or fly with large angle of attack can lead to a pressure inlet distortion.The crosswinds at the inlet and the growth of the boundary layer in the inlet duct can also cause a pressure loss inside the inlet.In addition,the S-duct inlet can cause a couple of counter rotating vortexes at the interface of the duct and the engine,which will finally result in a distorted inlet flow.In the multispool engine,the rotating stall occurring in the upstream compressor stages seems to be a rotating pressure inlet distortion to the downstream stages.For fighters,in the formation flying,the separation flow and exhaust gases from the former flight may cause non-uniform temperature and total pressure inletcondition to the later flight.When launching missiles,the fighter may ingest the exhaust gases from the weapon,which can also lead to non-uniform temperature pro files.These aero-engine distorted inlet conditions directly in fluence the performance and stability of its compression system,such as fan/compressor,and then reduce the operation function,safety and fuel economy of the engine.Thus,the investigation of various inlet distortions and their in fluences on the performance and stability of the compression system are of great signi ficance for researchers and aero-engine companies.

    Among these different types of inlet distortion,the total pressure distortions are more common.According to whether the distortion area is transient or not,the total pressure distortion can be classified as the steady type containing radial inlet distortion,circumferential inlet distortion and their combination,and the unsteady type including rotating inlet distortion and swirl inlet distortion.Lots of works have been carried out to study the impact of the circumferential and radial total pressure distortion,1–8the rotating inlet distortion9–11and the swirl inlet distortion12–14on the performance and stability of the fan/compressor.Summarily,the circumferential and radical total pressure distortion can often lead to an impairment in the performance and a decrease of stall margin.The rotating and swirl type inlet distortion can also cause adverse impacts under some conditions.While in other special circumstances,they can even bring benefits to the performance and stall margin.All these studies on different types of inlet distortion and their impacts on the performance and stall margin of fan/compressor in aero-engine are very important and meaningful.Furthermore,the research works on how to decrease or even overcome these negative influences of the inlet distortion have more significance in practice.

    As we all know,there are many techniques to increase the stall margin of fan/compressor,such as the active rotating stall/surge control and casing treatment technique.The idea of compressor instability active control was originally presented by Epstein et al.15and first realized in laboratory by Ffowcs Williams and Huang on a centrifugal compressor.16After that,many active rotating stall control researches were carried out.17–22Some of them were conducted to remedy the stall margin loss caused by inlet distortion.Van Schalkwyk et al.20used an active feedback control approach to stabilize the circumferential inlet distortion induced instability associated with rotating stall onset on a three-stage low-speed compressor.The results showed that the stabilization allowed stallfree operation to be extended below the natural(distorted flow)stall point by up to 3.7 percent for a 0.8 dynamic head distortion.For a 1.9 dynamic head distortion,40 percent of the mass flow range lost due to inlet distortion was regained through active control.Spakovszky et al.attempted to stabilize rotating stall in a single-stage transonic axial flow compressor with radial21or circumferential22inlet distortion by active feedback control using 12 jet-injectors annularly arrayed upstream of the rotor tip.The results showed that under radial inlet distortion,steady axisymmetric injection of 4%of the compressor mass flow could make a reduction in stalling mass flow of 9.7%and unsteady non-axisymmetric injection could achieve a further 7.5%reduction.21Under circumferential total pressure distortion of about one dynamic head and a 120-degree extent,steady axisymmetric,unsteady asymmetric and unsteady non-axisymmetric injection resulted in 6.2%,9%and 16.4%reduction of stalling mass flow,respectively.22

    On the other hand,casing treatment technique is well known as an effective approach to extend stall margin.This approach was first found by Koch23in experiments accidentally.Since then,research works about the stall margin enhancement of fan/compressor with this technique and the mechanism behind it were conducted.24–26Some works also have been carried out to study its ability to make up the stall margin loss caused by the distorted inlet flow.Langston27studied some techniques to minimize the adverse effects of inlet distortion and found that the use of honeycomb blade tip shrouds could reduce the sensitivity of compressor to distortion.Madden and West28experimentally and numerically investigated the effects of the rotor casing treatment on a single stage military fan with tip-low total pressure radial inlet distortion.The results showed that the stability line can be successfully extended by using the rotor casing treatment.However,the designs of these kinds of casing treatment are based on trial and error method by experiments.Its perforated ratio is often more than 50%.Using this technique to extend the stall margin loss will bring additional efficiency loss with respect to the efficiency decrease caused by inlet distortion.24,28

    Recently,a kind of stall precursor-suppressed(SPS)casing treatment29–35has been proposed to overcome the weakness of conventional casing devices.Its mechanism is based on the interaction of the pressure wave and the vortex wave,which can add some damping on the boundary of the compression system and influence the evolution of stall precursors near the rotating stall.29,30And the design of this casing treatment is guided by a 3D compressible stability model rather than trial-and-error methods.30,33Its ability to enhance the stall margin has been confirmed theoretically30,33and experimentally.31,32,34The experimental results from two low-speed compressors and a transonic compressor all showed that the SPS casing treatment could extend the stall margin by 8–18%.More importantly,this technique has obvious advantages over the traditional casing treatment because it can improve the stall margin without obviously changing the variation trend of the pressure-rise characteristics or leading to a recognizable efficiency loss.Further works were also conducted to study its ability to make up the stall margin loss caused by the rotating inletdistortion35and circumferential total pressure inlet distortion.The results showed that the SPS casing treatment,under the co-rotating inlet distortion,could not only make up the stall margin loss caused by inlet distortion but also enhance the stall margin compared to the solid wall casing under the clean inlet condition.35With the circumferential total pressure inlet distortion,the SPS casing treatment could remedy the stall margin loss caused by inlet distortion in small level and partly recover the stall margin loss caused by inlet distortion in large level.At the same time,the additional efficiency variation caused by the SPS casing treatment was less than 1%under both situations,which was within the range of the experimental error.Now,it is natural to wonder whether this technique can also take effects on the stall margin extending with radial inlet distortion,and how it will work in this condition.

    In the present paper,experiments are conducted to study the effectiveness of SPS casing treatment to make up the stall margin loss caused by the radial total pressure inlet distortion on a low speed compressor.The radial distorted inlet flow fields are generated by annular distortion flow generators.The characteristic curves are measured under the conditions with and without SPS casing treatment and radial inlet distor-tion in different levels as well as uniform inlet flow.With these experimental data,the stall margins are accordingly calculated,along with studying the stall margin loss caused by distorted inlet flow and the ability of the SPS casing treatment to recover this stall margin loss.The per-stall behavior of the compressor is also studied to discuss the mechanism behind the stability enhancement.

    2.Experimental facilities and instruments

    The experiments conducted in this paper are based on these facilities and instruments discussed below.

    2.1.Low-speed test compressor

    The experiments conducted in present paper are carried out with a low-speed test compressor called TA36,which is shown in Fig.1.This low-speed compressor is a high load rig with small aspect ratio,consisting of an inlet duct(1 in Fig.1(b)),an axial rotor(8 in Fig.1(b)),an axial stator(9 in Fig.1(b)),an exhaust duct(13 in Fig.1(b))and a throttle composed of a fixed cone(15 in Fig.1(b))and a movable annular sleeve(16 in Fig.1(b))around the cone.

    The rotor is driven by an asynchronous ac-motor(14 in Fig.1(b))with the normal rated power of 15 kW.The rotational speed is controlled by a frequency converter and can vary from 0 to the design value.The load of blades is varied by adjusting the axial position of the annular sleeve to change the area of the outlet.In this way,the operating point can accurately and slowly approach the rotating stall and quickly move away from it.The main design parameters of TA36 are displayed in Table 1.

    2.2.SPS casing treatment

    The SPS casing treatment(10 in Fig.1(b))is installed on the upstream of the rotor blade row and has an overlap with the rotor tip region.A typical configuration of SPS casing treatment is depicted in Fig.2.It contains an annular backchamber and a perforated plate.The chamber’s height is about 50 mm and its length can be adjusted from 100 to 200 mm.The length of the perforated plate is 75 mm and its open area ratio can be adjusted by a controller from 0%to 6.4%so that the casing can switch between the ordinary solid wall casing condition and the SPS casing treatment condition.The SPS casing treatment is designed in terms of the mechanism of the vortex wave and pressure wave interaction(Fig.3)to add some damping into the system,which can influence the stall inception or the evolution of stall precursors.

    Fig.1 Exterior and configuration of TA36.

    2.3.Radial inlet distortion generator

    The radial inlet distortion generator(RIDG)(4 in Fig.1(b))is installed on the casing upstream of the rotor.The distance from the generator to the rotor blade tip(Din Fig.4)is about 450 mm.The RIDGs with different heights(Hin Fig.4)are designed to cover the area of the inlet from 10%to 30%.Using different distortion generators can generate radial distortion inlets with different intensities.

    Table 1 Design parameters of TA36.

    Fig.2 Configuration of SPS casing treatment.

    Fig.3 Mechanism of SPS casing treatment.

    Fig.4 Radial inlet distortion generator.

    2.4.Measurement and data acquisition system

    The measurement system contains the inlet flow field measurement part,the steady-state characteristic measurement part,the dynamic pressure measurement part and data acquisition part.

    2.4.1.Inlet flow field measurement part

    Inlet flow field measurement part is designed to obtain the distortion flow field structure of the selected inlet cross section(plane III in Fig.1(b)).Eight total pressure combs shown in Fig.5(a)are used to measure the inlet total pressure distribution of the whole cross section.There are 40 Kulite highfrequency response pressure sensors installed in these combs(5 for each).As the distortion of the flow field studied in this paper is almost axisymmetric and mainly near the casing wall,a 1D hotwire sensor shown in Fig.5(b)is used to detailedly measure the flow field near the casing wall.

    Fig.5 Inlet flow field measurement part.

    2.4.2.Steady-state characteristic measurement part

    The steady-state characteristics measurement part contains an inlet total pressure comb(2 in Fig.1(b)),4 inlet static pressure sensors(3 in Fig.1(b)),an outlet total pressure comb(11 in Fig.1(b))and 4 outlet static pressure sensors(12 in Fig.1(b)).As shown in Fig.1,the inlet and outlet sensors are mounted in the plane I and plane V to measure the inlet/outlet total and static pressure,respectively.Each total pressure combcontains 5 sensors in the radial direction and 4 inlet/outlet static pressure sensors circumferentially equidistantly mounted on the casing.All characteristic parameters can be calculated from these measures as well as geometrical dimensions of the compressor.

    2.4.3.Dynamic pressure measurement part

    The dynamic pressure measurement part contains 8 highfrequency response dynamic pressure sensors(Fig.1).They are circumferentially equidistantly installed on the casing wall in the plane IV(Fig.6)which is about half a chord length upstream of the rotor blade tip.This part is used to record the static pressure perturbations near the rotor tip region.

    2.4.4.Data acquisition part

    The data acquisition part is mainly a NI-PXI system.This system contains 6 acquisition devices.It can totally provide 48 data acquisition channels.The maximum sampling rate can be up to 100 kHz.

    3.Experiment description

    3.1.Experimental conditions

    In this paper,the radial inlet distortion with 4 kinds of different intensities is studied.Different intensities of distortion are generated by RIDGs with different heights.The area covered by the RIDG in percentage of the total area of the inlet is used to reflect the intensities of the distortion flow field.For example,10%means that the RIDG’s area covers 10%of inlet area.The intensities of the distortion flow field studied in the present work and their corresponding height of RIDG are shown in Table 2.

    Under each intensity of inlet distortion together with the uniform inlet flow,experiments are conducted at 100%and 80%design speeds.Under each rotor speed,the inlet distor-tion flow field is measured at plane III when the downstream throttle is completely opened.Then the throttle is closed slowly to move the work point from large flow condition to rotating stallcondition.During this period,the pressure rise characteristic and efficiency characteristic are recorded.From the pre-stall to rotating stall condition,the near casing wall static pressure near the rotor tip region is recorded to analyze the evolution of the pressure perturbations just before stall occurs.

    Fig.6 Dynamic pressure measurement part.

    Table 2 Intensities of distortion flow field.

    3.2.Parameter calculation

    The parameters used in this paper are calculated based on the following formulas.Thecompressorcharacteristicisin the form of the total-to-static pressure rise coefficient versus the flow coefficient.The total-to-static pressure rise coefficient w is calculated as

    Fig.7 Total inlet pressure field structure under 10%radial inlet distortion.

    Fig.8 Radial distribution of axial velocity under different radial inlet distortions.

    whereNris the rotor speed andrmthe radius at mid-span of the blade.The flow coefficient u is calculated as

    whereVxis the axial inlet velocity,QVthe volume flow,andAthe area of the inlet.The efficiency g is expressed as

    Fig.9 Stall margin enhancement ability of SPS casing treatment with uniform inlet flow.

    Fig.10 Pressure rise curves under 100%design speed.

    To evaluate the stall margin enhancement ability of SPS casing treatment,a parameter SM is used here.Its definition is

    The subscripts s and d stand for the corresponding parameters at stall point and design point,respectively.

    4.Experiment results

    4.1.Intensity of radial inlet distortion

    As mentioned above,different RIDG can generate radial inlet distortion of different intensities.The total inlet pressure field structures of 80%and 100%design speeds under the 10%radial inlet distortion are shown in Fig.7.These figures are recorded by 8 total pressure combs.

    From the figure,the total pressure loss can be observed by the RIDG near the casing wall.The flow field is nearly axisymmetric and changes only in the radial direction.To understand this change in detail,the axial velocity is measured along a radial line(Fig.7(b))using a 1D hotwire sensor.As the static pressure on plan III is nearly uniform,the radial distribution of the axial velocity can re flect the total pressure distribution.

    The axial velocity distribution is depicted in Fig.8.In this figure,80%-UI(uniform inlet flow)means the result of experiment at 80%design speed with uniform inlet flow,while 100%-30 means the result at 100%design speed with 30%inlet distortion.From the figure,it is seen that the axial velocity loss appears near the casing wall.The larger the intensity of inlet distortion is,the more obvious the difference between the velocities of main flow and blade tip regions becomes.For different rotor speeds,the axial velocity distribution is similar with some quantitative differences.Thus,the RIDG can generate radial inlet distortions of different intensities as it is planned to simulate the inlet flow conditions.

    4.2.Stall margin enhancement

    The stall margin enhancement ability of SPS casing treatment with uniform and distortion inlet flow is studied in this part.Different experimental conditions are illustrated as‘Rotating speed-Inlet flow condition-Casing condition’.For example,100%-UI-SC(solid casing)means that the experiment is conducted under 100%design speed and uniform inlet flow with a solid casing(the open area ratio of the perforated plate is 0),while 80%-30-CT(casing treatment)means that the experiment is carried out under 80%design speed and 30%distortion inlet flow with the casing treatment(the open area ratio of the perforated plate is 6.4%).

    The pressure rise curves corresponding to 100%and 80%design speeds with uniform inlet flow are depicted in Fig.9.As calculated by Eq.(5),SM is 45.95%and 45.61%at 100%and 80%design speeds with SC,respectively.After the SPS casing treatment is opened,SM is enhanced to 50.72%and 48.98%.Thus,the SPS casing treatment can improve the stall margin of the compressor with uniform inlet flow.

    Fig.11 Efficiency curves under 100%design speed.

    Fig.12 Pressure rise curve under 80%design speed with inlet distortion.

    As we know,inlet distortion can cause stall margin loss.The pressure rise curves under 100%design speed with 10%intensity of inlet distortion are shown in Fig.10.From the figure,it is seen that the 10%radial inlet distortion causes 3.11%loss of stall margin.Under this condition,the SPS casing treatment can not only make up the stall margin loss caused by the inlet distortion but also improve the stall margin by 2.57%.

    When the intensity of inlet distortion increases to 20%,which causes 12.05%loss of stall margin,the SPS casing treatment cannot totally cover the loss caused by distorted inlet,but it can enable the stall margin loss decrease from 12.05%to 2.99%(Fig.10).

    The efficiency curves corresponding to 100%design speed with 10%and 20%inlet distortion are shown in Fig.11.From these figures,it is known that the radial inlet distortion can cause big efficiency loss.The 20%distortion can cause larger loss of efficiency than the 10%distortion.The use of SPS casing treatment can lead to additional efficiency loss in both cases.However,the corresponding efficiency loss from the use of the SPS casing treatment is much smaller than the loss caused by inlet distortion.The calculation results show that the maximum efficiency losses caused by the SPS casing treatment are only 0.4%and 0.62%under 10%and 20%inlet distortions,respectively.These are basically within the error range(1%)of the measurement system.

    Additional results under 25%and 30%inlet distortion intensities are shown in Figs.10 and 11.With the rise of inlet distortion intensity,both the stall margin loss and efficiency loss are increased.The stall margin decreases by 20.18%under 25%inlet distortion,while this loss rises to 32.52%when intensity of inlet distortion reaches up to 30%.The SPS casing treatment still can partly cover the stall margin loss from distorted inlet flow and does not bring obvious additional eff iciency.Under 25%and 30%inlet distortions,the SPS casing treatment can offset the stall margin losses of 8.34%and 8.79%,meanwhile bringing maximum efficiency losses of 0.82%and 0.3%,respectively.

    Similar studies are also conducted under 80%design speed.The results are shown in Figs.12 and 13.From these figures,the same qualitative law can be observed.

    For a detailed analysis,all the results discussed above are summarized and shown in Table 3.From the table,it is found that radial inlet distortion can cause the obvious stall margin loss under both rotor speeds.The stronger the intensity of inlet distortion is,the more the stall margin loss it brings.The SPS casing treatment can make up or partly make up the stall margin loss about 5%–10%caused by radial inlet distortion.When the intensity of the inlet distortion is relatively weak,the SPS casing treatment can not only recover the lost stall margin but also extend it(minus value in Table 3).The radial inlet distortion can always bring efficiency loss.And the appli-cation of SPS casing treatment may cause additional efficiency loss.However,this additional efficiency loss is less than 1%,which is basically within the error range of the measurement system and much smaller than the loss caused by the radial inlet distortion.

    Fig.13 Efficiency curves under 80%design speed with inlet distortion.

    4.3.Experimental results of pre-stall behavior

    In the previous work,the mechanism of the SPS casing treatment in enhancing the stall margin of the compression system was validated theoretically and experimentally.The relevant investigations29,30show that adding some damping on the boundary of the system can weaken the energy of unsteady flow perturbations and delay the nonlinear ampli fication of the stall precursor waves.In Refs.21,35,time-frequency analysis is conducted and a parameter called power spectral density(PSD)of pressure wave is used to analyze the pre-stall behavior with solid casing and casing treatment in the compressor under uniform,rotating distorted and circumferential pressure distorted inlet flow,respectively.Here,the similar technical method is applied to investigate the pre-stall behavior in the compressor under radial inlet distortion.The definition of PSD is

    Fig.14 PSD shown in time-frequency domain.

    Fig.15 PSD at 100%design speed with 10%radial inlet distortion.

    Fig.16 PSD at 100%design speed with 20%radial inlet distortion.

    Table 3 The enhanced stall margin and max g loss caused by SPS casing treatment

    The analyzed result of one case discussed above is shown in Fig.14.Thexaxis is the relative frequency normalized by the shaft frequency of the compressor.The relative frequency is shown from 0.2 to 0.8 because pressure disturbance within this frequency interval is related to the occurrence of rotating stall.Theyaxis is the time axis.And at the same moment,the two compared cases have the same operating mass flow.Thezaxis is the value of PSD.From the view of time domain,we can see that the SPS casing treatment can suspend the nonlinear increase of the pressure disturbance and postpone the occurrence of rotating stall.

    However,the comparison of the pre-stall behavior of the compression system with solid casing and casing treatment in this kind of figure is not clear.In this 3D view,it is difficult to distinctly distinguish the differences between using the solid casing and SPS casing treatment.In order to make this comparison clearly,the data within the intervalABare extracted and shown only in the frequency domain(projecting the figure on the relative frequency-PSD plane).The pointAis the stall point of the solid casing case and pointBis a near stall work point.PointBis very close to the stall pointAbecause the change of the work point is very slow and the time interval betweenAandBis only 1 s.For every time point in this interval,the flow range is equivalent for the two situations.

    More results of the pre-stall behavior are shown in Fig.17.To show the variation with and without SPS casing treatment clearly,all the results about pre-stall behavior at 100%and 80%design speeds are summarized in Table 4,which shows the comparison of PSD peak values and its decrease due to the application of SPS casing treatment.From this table,it is seen that the value of PSD decreases when the operating speed decreases.In all of the cases,the results show that the SPS casing treatment can clearly weaken the power of the wideband pressure disturbance at pre-stall stage,and especially reduce its peak value.

    Fig.17 PSD at 100%design speed with 25%and 30%radial inlet distortions.

    Table 4 PSD peak value and its decrease due to SPS casing treatment.

    5.Conclusions

    In the present work,experiments are carried out to study the influence of radial inlet distortion on the compression system and the effects of SPS casing treatment under this condition.Some conclusions can be summarized as follows:

    (1)The radial inlet distortion can cause a loss of stall margin and overall efficiency.The stronger the intensity of inlet distortion is,the more stall margin and efficiency loss it brings.

    (2)With uniform inlet flow,the application of SPS casing treatments can improve the stall margin without obviously changing the pressure-rise characteristics.

    (3)With radial inlet distortion,the SPS casing treatment can recover or partly recover the stall margin loss caused by radial inlet distortion.When the intensity of the inlet distortion is relatively weak,the SPS casing treatment can even extend the stall margin.

    (4)The SPS casing treatment may cause slightly additional efficiency loss.This additional efficiency loss is less than 1%,which is basically within the error range of the measurement system and much smaller than the loss caused by the radial inlet distortion.

    (5)The pre-stall behavior with solid casing and casing treatment in the compressor is analyzed by using the timefrequency analysis.The power spectral density is calculated to evaluate the intensity of pressure disturbance at pre-stall stage.The results show that the SPS casing treatment can enhance the stall margin of the compressor by weakening the power of the unsteady pressure perturbations in broadband and suppressing the nonlinear increase of stall precursor.

    Acknowledgements

    This research was supported by National Natural Science Foundation of China(51236001 and 51576008),National Basic Research Program of China(2012CB720201),Aeronautical Science Foundation of China(2014ZB51018)and Fundamental Research Funds for the Central Universities.

    1.Plourde GA,Stenning AH.Attenuation of circumferential inlet distortion in multistage axial compressors.J Aircraft1968;5(3):236–42.

    2.Cook R.Design and test of a highly-loaded three-stage,axial-flow compressor.Reston:AIAA;1976.Report No.:AIAA 1976–6.

    3.Tamaki T,Nagano S.Effects of inlet distortions on a multi-stage compressor.Reston:AIAA;1979.Report No.:AIAA 1979–7003.

    4.Maekawa A,Higashi T,Tanaka S.Performance of rotating cascades under the inlet-distortions:1st report,performance under the radial distortions.Bull JSME1979;22(165):333–9.

    5.Seo N,Nakatani T,Tanaka S.Performance of rotating cascades under the inlet-distortions:2nd report,performance under the circumferential distortions.Bull JSME1981;24(190):628–33.

    6.Hynes TP,Greitzer EM.A method for assessing effects of circumferential flow distortion on compressor stability.J Turbomachinery1987;109(3):371–9.

    7.Hah C,Rabe DC,Sullivan TJ,Wadia AR.Effects of inlet distortion on the flow field in a transonic compressor rotor.J Turbomachinery1998;120(2):233–46.

    8.Lesser A,Niehuis R.Transonic axial compressors with total pressure inlet flow field distortions.New York:ASME;2014.Report No.:ASME GT2014-26627.

    9.Longley JP,Shin H,Plumley RE,Silkowski PD,Day IJ,Greitzer EM,et al.Effects of rotating inlet distortion on multistage compressor stability.J Turbomachinery1996;118(2):181–8.

    10.Peters T,Bu¨rgener T,Fottner L.Effects of rotating inlet distortion on a 5-stage hp-compressor.New York:ASME;2001.Report No.:ASME 2001-GT-0300.

    11.Peters T,Fottner L.Effects of co-and counter-rotating inlet distortions on a 5-stage hp-compressor.New York:ASME;2002.Report No.:ASME 2002-GT-30395.

    12.Pazur W,Fottner L.The influence of inlet swirl distortions on the performance of a jet propulsion two-stage axial compressor.J Turbomachinery1991;113(2):233–40.

    13.Schmid NR,Leinhos DC,Fottner L.Steady performance measurements of a turbofan engine with inlet distortions containing co-and counter rotating swirl from an intake diffuser for hypersonic flight.J Turbomachinery2001;123(2):379–85.

    14.Sheoran Y,Bouldin B,Krishnan PM.Compressor performance and operability in swirl distortion.J Turbomachinery2012;134(4):1–13.

    15.Epstein AH,Ffowcs Williams JE,Greitzer EM.Active suppression of aerodynamic instabilities in turbomachines.J Propulsion Power1989;5(2):204–11.

    16.Ffowcs Williams JE,Huang XY.Active stabilization of compressor surge.J Fluid Mech1989;204(1):245–62.

    17.Day IJ.Active suppression of rotating stall and surge in axial compressors.J Turbomachinery1993;115(1):40–7.

    18.Paduano JD,Epstein AH,Valavani L,Longley JP,Greitzer EM,Guenette GR.Active control of rotating stall in a low speed axial compressor.J Turbomachinery1993;115(1):48–56.

    19.Haynes JM,Hendricks GJ,Epstein AH.Active stabilization of rotating stall in a three-stage axial compressor.J Turbomachinery1994;116(2):226–39.

    20.Van Schalkwyk CM,Paduano JD,Greitzer EM,Epstein AH.Active stabilization of axial compressors with circumferential inlet distortion.J Turbomachinery1998;120(3):431–9.

    21.Spakovszky ZS,Weigl HJ,Paduano JD,Van Schalkwyk CM,Suder KL,Bright MM.Rotating stall control in a high-speed stage with inlet distortion:Part I—Radial distortion.J Turbomachinery1999;121(3):510–6.

    22.Spakovszky ZS,Van Schalkwyk CM,Weigl HJ,Paduano JD,Suder KL,Bright MM.Rotating stall control in a high-speed stage with inlet distortion:Part II—Circumferential distortion.J Turbomachinery1999;121(3):517–24.

    23.Koch CC.Experimental evaluation of outer case blowing or bleeding of single stage axial flow compressor.Washington,D.C.:NASA;1970.Report No.:NASA CR-54592.

    24.Kang CS,Mckenzie AB,Elder RL.Recessed casing treatment effects on fan performance and flow field.New York:ASME;1995.Report No.:ASME 95-GT-197.

    25.Kroeckel T,Hiller SJ,Jeschke P.Application of a multistage casing treatment in a high speed axial compressor test rig.New York:ASME;2011.Report No.:ASME GT2011-46315.

    26.Pixberg CT,Schiffer H,Ross MH,Cameron JD,Morris SC.Stall margin improvement by use of casing treatments.New York:ASME;2013.Report No.:ASME GT2013-95842.

    27.Langston CE.Distortion tolerance:by design instead of by accident.New York:ASME;1969.Report No.:ASME 69-GT-115.

    28.Madden DS,West MA.Effects of inlet distortion on the stability of an advanced military swept fan stage with casing treatment.New York:ASME;2005.Report No.:ASME GT2005-68693.

    29.Sun X,Sun D,Yu W.A model to predict stall inception of transonic axial flow fan/compressors.Chin J Aeronautics2011;24(6):687–700.

    30.Sun X,Sun D,Liu X,Yu W,Wang X.Theory of compressor stability enhancement using novel casing treatment,Part I:Methodology.J Propulsion Power2014;30(5):1224–35.

    31.Sun D,Liu X,Jin D,Gui X,Sun X.Theory of compressor stability enhancement using novel casing treatment,Part II:Experiment.AIAA J Propulsion Power2014;30(5):1236–47.

    32.Dong X,Liu X,Sun D,Sun X.Experimental investigation on sps casing treatment with bias flow.Chin J Aeronautics2014;27(6):1352–62.

    33.Sun D,Liu X,Sun X.An evaluation approach for the stall margin enhancement with stall precursor-suppressed casing treatment.J Fluids Eng2015;137(8):1–16.

    34.Sun D,Nie C,Liu X,Lin F,Sun X.Further investigation on transonic compressor stall margin enhancement with stall precursor-suppressed casing treatment.J Turbomachinery2016;138(2):1–13.

    35.Dong X,Sun D,Li F,Jin D,Gui X,Sun X.Effects of rotating inlet distortion on compressor stability with stall precursorsuppressed casing treatment.J Fluids Eng2015;137(11):1–15.

    19 April 2016;revised 3 September 2016;accepted 4 September 2016

    Available online 22 December 2016

    *Corresponding author.

    E-mail address:sundk@buaa.edu.cn(D.Sun).

    Peer review under responsibility of Editorial Committee of CJA.

    国产精品野战在线观看| 日日爽夜夜爽网站| 久久性视频一级片| 国产一区二区三区在线臀色熟女| avwww免费| 在线观看www视频免费| 丁香六月欧美| 熟女少妇亚洲综合色aaa.| 国产亚洲av嫩草精品影院| 岛国视频午夜一区免费看| 国产三级在线视频| 亚洲三区欧美一区| 免费观看精品视频网站| 在线看三级毛片| 99riav亚洲国产免费| 国产精品av久久久久免费| 久久精品国产清高在天天线| 老司机靠b影院| 免费无遮挡裸体视频| 欧美一级a爱片免费观看看 | 天天添夜夜摸| 宅男免费午夜| 久久久国产欧美日韩av| 国产精品亚洲一级av第二区| 日日爽夜夜爽网站| 国产人伦9x9x在线观看| 日本撒尿小便嘘嘘汇集6| 中国美女看黄片| 国产精品爽爽va在线观看网站 | 国产精品二区激情视频| 精品一区二区三区视频在线观看免费| 久久久久久久久免费视频了| 亚洲aⅴ乱码一区二区在线播放 | 少妇粗大呻吟视频| 丁香六月欧美| 无人区码免费观看不卡| 欧美日韩亚洲综合一区二区三区_| 精品国产一区二区三区四区第35| 色婷婷久久久亚洲欧美| 欧美另类亚洲清纯唯美| 少妇 在线观看| 两性夫妻黄色片| 1024视频免费在线观看| 天天添夜夜摸| 久久久国产精品麻豆| 中文字幕人妻丝袜一区二区| 免费av毛片视频| 久久狼人影院| 亚洲欧美日韩无卡精品| 亚洲真实伦在线观看| 黄色 视频免费看| 成人特级黄色片久久久久久久| 国产精品乱码一区二三区的特点| 久久精品国产综合久久久| 久久 成人 亚洲| 最新美女视频免费是黄的| 亚洲国产精品久久男人天堂| 欧美成人免费av一区二区三区| 18禁黄网站禁片午夜丰满| 亚洲色图 男人天堂 中文字幕| 欧美成人一区二区免费高清观看 | 两个人看的免费小视频| 黄色丝袜av网址大全| www.熟女人妻精品国产| 少妇的丰满在线观看| 女人被狂操c到高潮| 伦理电影免费视频| 亚洲男人的天堂狠狠| 中出人妻视频一区二区| 久久伊人香网站| 免费一级毛片在线播放高清视频| 久久香蕉国产精品| 国产免费av片在线观看野外av| 黄片小视频在线播放| 亚洲av电影在线进入| 久久中文字幕一级| 一个人免费在线观看的高清视频| 亚洲成av片中文字幕在线观看| 国内揄拍国产精品人妻在线 | 岛国视频午夜一区免费看| 美女午夜性视频免费| 午夜影院日韩av| 亚洲av成人不卡在线观看播放网| 久久久久久免费高清国产稀缺| 国产精品免费视频内射| 婷婷亚洲欧美| 久久精品亚洲精品国产色婷小说| 人妻久久中文字幕网| 国产欧美日韩一区二区精品| 18美女黄网站色大片免费观看| 日韩欧美免费精品| 90打野战视频偷拍视频| 色综合站精品国产| 国产午夜精品久久久久久| 欧美激情 高清一区二区三区| 国内揄拍国产精品人妻在线 | 黄色视频不卡| 久久久久久久久久黄片| 香蕉av资源在线| 一区福利在线观看| 亚洲国产看品久久| 日韩精品青青久久久久久| ponron亚洲| 不卡一级毛片| 在线观看www视频免费| 国产亚洲精品久久久久久毛片| 色老头精品视频在线观看| 老司机深夜福利视频在线观看| 热re99久久国产66热| 亚洲av中文字字幕乱码综合 | av欧美777| avwww免费| 日韩欧美一区视频在线观看| 午夜a级毛片| 一a级毛片在线观看| 18禁黄网站禁片免费观看直播| 2021天堂中文幕一二区在线观 | 在线观看日韩欧美| 欧美人与性动交α欧美精品济南到| 久久久久九九精品影院| 18禁裸乳无遮挡免费网站照片 | www国产在线视频色| 首页视频小说图片口味搜索| 十分钟在线观看高清视频www| 国产区一区二久久| 亚洲国产中文字幕在线视频| 一级毛片精品| 久久国产亚洲av麻豆专区| 在线观看舔阴道视频| 亚洲国产日韩欧美精品在线观看 | 男男h啪啪无遮挡| 男男h啪啪无遮挡| 搡老熟女国产l中国老女人| 黄色毛片三级朝国网站| 久久国产亚洲av麻豆专区| 波多野结衣巨乳人妻| 99精品在免费线老司机午夜| 国产精品野战在线观看| 中文字幕高清在线视频| 丁香六月欧美| 看黄色毛片网站| 亚洲国产日韩欧美精品在线观看 | 国产精品1区2区在线观看.| 国产av在哪里看| 动漫黄色视频在线观看| 国产精品98久久久久久宅男小说| 亚洲专区字幕在线| 麻豆一二三区av精品| 我的亚洲天堂| 最近最新中文字幕大全电影3 | 丁香六月欧美| 黄片播放在线免费| √禁漫天堂资源中文www| 精品午夜福利视频在线观看一区| 三级毛片av免费| 黄色a级毛片大全视频| 国产成人精品久久二区二区免费| 国产又黄又爽又无遮挡在线| 成年女人毛片免费观看观看9| 麻豆成人午夜福利视频| 国产成人影院久久av| 国产精品久久久久久人妻精品电影| 国产亚洲精品久久久久5区| 久久久国产精品麻豆| 国产黄a三级三级三级人| 久久久久久久久久黄片| 国产精品免费一区二区三区在线| 亚洲片人在线观看| 少妇被粗大的猛进出69影院| 91成年电影在线观看| 久久亚洲真实| 老司机午夜福利在线观看视频| 黑人欧美特级aaaaaa片| 美国免费a级毛片| 精品不卡国产一区二区三区| 亚洲成人久久性| 久久国产精品男人的天堂亚洲| 亚洲男人天堂网一区| 久久人妻福利社区极品人妻图片| 欧美性猛交黑人性爽| 亚洲精品中文字幕一二三四区| 午夜精品在线福利| 非洲黑人性xxxx精品又粗又长| 国产成人啪精品午夜网站| 成人欧美大片| 美女国产高潮福利片在线看| 国产高清视频在线播放一区| 日韩 欧美 亚洲 中文字幕| 久久精品影院6| 两个人视频免费观看高清| 久久九九热精品免费| 可以免费在线观看a视频的电影网站| 午夜久久久久精精品| 国产精品亚洲美女久久久| 精品一区二区三区视频在线观看免费| 亚洲美女黄片视频| 久久久久久久久久黄片| 欧美激情高清一区二区三区| 男人的好看免费观看在线视频 | 亚洲色图 男人天堂 中文字幕| 婷婷丁香在线五月| 一级毛片女人18水好多| 色老头精品视频在线观看| 国产高清视频在线播放一区| 亚洲国产看品久久| 久久久久久久午夜电影| 久久九九热精品免费| 老司机午夜福利在线观看视频| 丁香六月欧美| 国产成人欧美| 午夜免费观看网址| 香蕉丝袜av| www.www免费av| 亚洲欧美激情综合另类| 亚洲国产中文字幕在线视频| 看免费av毛片| 在线十欧美十亚洲十日本专区| 动漫黄色视频在线观看| 欧美日本亚洲视频在线播放| 少妇被粗大的猛进出69影院| 琪琪午夜伦伦电影理论片6080| 两个人视频免费观看高清| 制服诱惑二区| 久久久精品国产亚洲av高清涩受| 老司机福利观看| 成人一区二区视频在线观看| 色婷婷久久久亚洲欧美| 亚洲熟女毛片儿| 大型黄色视频在线免费观看| 99re在线观看精品视频| 久久久水蜜桃国产精品网| 亚洲一码二码三码区别大吗| 国产真实乱freesex| 成人亚洲精品av一区二区| 50天的宝宝边吃奶边哭怎么回事| 搞女人的毛片| 亚洲三区欧美一区| 2021天堂中文幕一二区在线观 | 亚洲欧美精品综合久久99| 久久人妻av系列| 欧美乱码精品一区二区三区| 少妇 在线观看| 男女午夜视频在线观看| 91麻豆精品激情在线观看国产| 成在线人永久免费视频| 18禁裸乳无遮挡免费网站照片 | 免费观看精品视频网站| 色婷婷久久久亚洲欧美| 久久久久久人人人人人| 亚洲精品久久成人aⅴ小说| 日本免费a在线| 69av精品久久久久久| 大香蕉久久成人网| 欧美日韩一级在线毛片| 性欧美人与动物交配| 亚洲国产精品久久男人天堂| 国产成人av激情在线播放| 亚洲av成人一区二区三| 超碰成人久久| 在线观看免费视频日本深夜| 每晚都被弄得嗷嗷叫到高潮| 无遮挡黄片免费观看| 国产精品av久久久久免费| 在线看三级毛片| 国产亚洲精品av在线| 女性被躁到高潮视频| 中文字幕人成人乱码亚洲影| 叶爱在线成人免费视频播放| 免费在线观看亚洲国产| 可以在线观看毛片的网站| 国产亚洲精品久久久久久毛片| 丁香欧美五月| 成人免费观看视频高清| 久久九九热精品免费| 成人18禁在线播放| 日韩国内少妇激情av| 午夜福利在线在线| 国产精品二区激情视频| 久久精品91无色码中文字幕| 99久久综合精品五月天人人| 久久精品国产综合久久久| 12—13女人毛片做爰片一| 午夜免费观看网址| 亚洲国产精品成人综合色| 两人在一起打扑克的视频| 一级毛片高清免费大全| 一区二区三区激情视频| 久久亚洲精品不卡| 免费女性裸体啪啪无遮挡网站| 亚洲人成伊人成综合网2020| 免费在线观看日本一区| 亚洲中文字幕日韩| 亚洲va日本ⅴa欧美va伊人久久| 激情在线观看视频在线高清| 国产男靠女视频免费网站| 国产欧美日韩精品亚洲av| 国产精品久久视频播放| 精品久久蜜臀av无| 亚洲一区二区三区不卡视频| 免费在线观看影片大全网站| 久久99热这里只有精品18| 精品免费久久久久久久清纯| 十分钟在线观看高清视频www| 99国产精品一区二区三区| 美女免费视频网站| а√天堂www在线а√下载| 亚洲avbb在线观看| 91九色精品人成在线观看| 老司机靠b影院| 国产91精品成人一区二区三区| 欧美在线一区亚洲| 国产亚洲欧美98| 日日摸夜夜添夜夜添小说| 久久久久久免费高清国产稀缺| 午夜福利一区二区在线看| 亚洲欧美一区二区三区黑人| 亚洲成a人片在线一区二区| 真人做人爱边吃奶动态| 亚洲国产欧美日韩在线播放| 久久久久久久久久黄片| 这个男人来自地球电影免费观看| 精品日产1卡2卡| 757午夜福利合集在线观看| 国产精品免费一区二区三区在线| 精品熟女少妇八av免费久了| 国产精品乱码一区二三区的特点| 国产精品一区二区三区四区久久 | 亚洲国产精品999在线| 国产高清有码在线观看视频 | 欧美另类亚洲清纯唯美| 波多野结衣高清无吗| 99riav亚洲国产免费| 香蕉久久夜色| 国产精品一区二区三区四区久久 | 在线十欧美十亚洲十日本专区| 哪里可以看免费的av片| 久久精品国产99精品国产亚洲性色| 欧美亚洲日本最大视频资源| 国产成人精品无人区| 免费观看精品视频网站| 婷婷六月久久综合丁香| 成年免费大片在线观看| 欧美成人午夜精品| 一夜夜www| 一区二区三区国产精品乱码| 久久亚洲精品不卡| 老司机午夜十八禁免费视频| 成人特级黄色片久久久久久久| 国产高清有码在线观看视频 | 午夜成年电影在线免费观看| 日韩成人在线观看一区二区三区| 欧美日韩福利视频一区二区| 精品少妇一区二区三区视频日本电影| 黄色视频,在线免费观看| 黄色丝袜av网址大全| 嫩草影视91久久| 给我免费播放毛片高清在线观看| 美女国产高潮福利片在线看| 日韩视频一区二区在线观看| xxx96com| 亚洲精品久久成人aⅴ小说| 欧美日韩福利视频一区二区| 亚洲专区字幕在线| 久久中文看片网| 久久午夜亚洲精品久久| 亚洲七黄色美女视频| 在线av久久热| 又黄又爽又免费观看的视频| 在线观看www视频免费| 国产精品二区激情视频| 别揉我奶头~嗯~啊~动态视频| 久久国产亚洲av麻豆专区| 久久久久久人人人人人| 91成年电影在线观看| 女人被狂操c到高潮| 欧美日韩瑟瑟在线播放| 男女下面进入的视频免费午夜 | 精品欧美国产一区二区三| 日日爽夜夜爽网站| 国产av在哪里看| 国产高清有码在线观看视频 | 国产精品久久视频播放| 两个人看的免费小视频| 99久久久亚洲精品蜜臀av| 91成年电影在线观看| 少妇的丰满在线观看| 国语自产精品视频在线第100页| 免费av毛片视频| 亚洲色图 男人天堂 中文字幕| 亚洲 欧美一区二区三区| 亚洲av成人不卡在线观看播放网| 听说在线观看完整版免费高清| 精品第一国产精品| 97人妻精品一区二区三区麻豆 | 国产亚洲欧美精品永久| 国产精品1区2区在线观看.| 香蕉av资源在线| 一级作爱视频免费观看| 香蕉久久夜色| 国内揄拍国产精品人妻在线 | 国产精品永久免费网站| 亚洲人成网站在线播放欧美日韩| 亚洲av第一区精品v没综合| 波多野结衣巨乳人妻| 男女下面进入的视频免费午夜 | 国产区一区二久久| 免费在线观看视频国产中文字幕亚洲| 久久久国产欧美日韩av| 亚洲九九香蕉| www国产在线视频色| 欧美国产精品va在线观看不卡| 美女午夜性视频免费| 亚洲av熟女| 国产久久久一区二区三区| 少妇 在线观看| 国产私拍福利视频在线观看| 一区二区三区精品91| 免费观看精品视频网站| 午夜福利高清视频| 久久久国产精品麻豆| 三级毛片av免费| 中文亚洲av片在线观看爽| 91成年电影在线观看| 国产亚洲欧美98| 黄色 视频免费看| 18美女黄网站色大片免费观看| 黄色毛片三级朝国网站| 搡老岳熟女国产| 香蕉国产在线看| 99久久99久久久精品蜜桃| 日韩欧美一区视频在线观看| 麻豆国产av国片精品| 亚洲va日本ⅴa欧美va伊人久久| 亚洲,欧美精品.| av片东京热男人的天堂| 在线观看66精品国产| 亚洲人成77777在线视频| 亚洲中文av在线| 十八禁人妻一区二区| 免费在线观看日本一区| 亚洲欧美日韩无卡精品| 国产亚洲av嫩草精品影院| 搡老熟女国产l中国老女人| 1024香蕉在线观看| av片东京热男人的天堂| 国产精品免费视频内射| 一个人免费在线观看的高清视频| 正在播放国产对白刺激| 18禁黄网站禁片免费观看直播| 热99re8久久精品国产| 国产精品免费视频内射| 日韩视频一区二区在线观看| 我的亚洲天堂| 国产av在哪里看| 日韩高清综合在线| 亚洲中文日韩欧美视频| 欧美黑人精品巨大| 女警被强在线播放| 欧美精品亚洲一区二区| 成人免费观看视频高清| 午夜成年电影在线免费观看| 亚洲 国产 在线| 最近最新免费中文字幕在线| 高清在线国产一区| 男人舔奶头视频| 这个男人来自地球电影免费观看| 精品国产国语对白av| 免费电影在线观看免费观看| 成年免费大片在线观看| 动漫黄色视频在线观看| 黄频高清免费视频| 日韩大尺度精品在线看网址| 美女午夜性视频免费| 精品熟女少妇八av免费久了| 亚洲av五月六月丁香网| 亚洲无线在线观看| 午夜久久久久精精品| 2021天堂中文幕一二区在线观 | 黑人操中国人逼视频| 亚洲男人天堂网一区| 91麻豆av在线| 999久久久国产精品视频| 视频在线观看一区二区三区| 好看av亚洲va欧美ⅴa在| 亚洲久久久国产精品| 国产亚洲欧美精品永久| www.自偷自拍.com| 大香蕉久久成人网| 亚洲一区中文字幕在线| 国产精品自产拍在线观看55亚洲| 久久青草综合色| 巨乳人妻的诱惑在线观看| 免费在线观看黄色视频的| 亚洲 国产 在线| 日日干狠狠操夜夜爽| 欧美一级a爱片免费观看看 | 老司机靠b影院| 亚洲精华国产精华精| 亚洲一区二区三区不卡视频| 国产精品日韩av在线免费观看| 成人18禁高潮啪啪吃奶动态图| 美女午夜性视频免费| x7x7x7水蜜桃| 黄网站色视频无遮挡免费观看| 听说在线观看完整版免费高清| 男女那种视频在线观看| 午夜福利欧美成人| 一进一出抽搐动态| 国产亚洲精品久久久久久毛片| а√天堂www在线а√下载| 亚洲av中文字字幕乱码综合 | 夜夜夜夜夜久久久久| 精品一区二区三区视频在线观看免费| 黄色视频,在线免费观看| 我的亚洲天堂| 国产免费av片在线观看野外av| 午夜久久久在线观看| 丁香欧美五月| 亚洲一区高清亚洲精品| 一进一出抽搐gif免费好疼| 久久久久国内视频| 日本熟妇午夜| 女性生殖器流出的白浆| 久久热在线av| 午夜福利欧美成人| 国产区一区二久久| 日韩免费av在线播放| 在线观看免费日韩欧美大片| 手机成人av网站| 免费看美女性在线毛片视频| 亚洲精品av麻豆狂野| 黄片小视频在线播放| 午夜成年电影在线免费观看| 男女午夜视频在线观看| 最近在线观看免费完整版| 变态另类成人亚洲欧美熟女| 成人午夜高清在线视频 | 欧美激情 高清一区二区三区| 亚洲精品国产一区二区精华液| 午夜免费成人在线视频| 国产真实乱freesex| 亚洲 欧美 日韩 在线 免费| 一级片免费观看大全| 亚洲最大成人中文| 久久中文字幕人妻熟女| 久久 成人 亚洲| 国产av不卡久久| 亚洲欧美日韩高清在线视频| 18禁国产床啪视频网站| 欧美中文综合在线视频| 啦啦啦免费观看视频1| 亚洲一区高清亚洲精品| 亚洲欧美日韩无卡精品| 老司机午夜福利在线观看视频| 亚洲国产精品sss在线观看| 桃色一区二区三区在线观看| 午夜免费成人在线视频| 久久这里只有精品19| 久久精品成人免费网站| 女性被躁到高潮视频| 精品免费久久久久久久清纯| 国产伦一二天堂av在线观看| 91成年电影在线观看| 麻豆一二三区av精品| 国产一区二区在线av高清观看| 免费在线观看日本一区| 天堂动漫精品| 高潮久久久久久久久久久不卡| 亚洲男人天堂网一区| 日本一本二区三区精品| 免费在线观看完整版高清| 最新美女视频免费是黄的| 亚洲成人国产一区在线观看| 久久午夜综合久久蜜桃| 国产91精品成人一区二区三区| 精品久久久久久久久久久久久 | 日韩欧美免费精品| 国产亚洲精品一区二区www| 俄罗斯特黄特色一大片| av电影中文网址| 少妇裸体淫交视频免费看高清 | 欧美黄色淫秽网站| 色在线成人网| 久久久国产成人精品二区| 亚洲欧美日韩无卡精品| 久久草成人影院| 中文字幕人成人乱码亚洲影| 亚洲avbb在线观看| 久久精品国产亚洲av香蕉五月| 亚洲第一电影网av| 午夜免费观看网址| 黑人欧美特级aaaaaa片| 国产一卡二卡三卡精品| 九色国产91popny在线| 中文字幕人妻丝袜一区二区| 中文字幕久久专区| 国产三级黄色录像| 99久久99久久久精品蜜桃| 亚洲九九香蕉| 欧美精品啪啪一区二区三区| 啪啪无遮挡十八禁网站| 国产精品1区2区在线观看.| 久久久久国内视频| 欧美日韩黄片免| 免费在线观看黄色视频的| 波多野结衣高清无吗| 亚洲一区高清亚洲精品| 99re在线观看精品视频| 高清在线国产一区| 久久久久久久久久黄片| 精品国内亚洲2022精品成人| 亚洲无线在线观看| 成人永久免费在线观看视频| 成人午夜高清在线视频 | 黄频高清免费视频|