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

    Effect of inlet rotating swirl on endwall film cooling for two representative hole arrangements

    2018-05-17 10:07:36YangZHOUYangZHANGXinrongSUXinYUAN
    CHINESE JOURNAL OF AERONAUTICS 2018年5期

    Yang ZHOU,Yang ZHANG,Xinrong SU,Xin YUAN

    Key Laboratory for Thermal Science and Power Engineering of Ministry of Education,Tsinghua University,Beijing 100084,China

    1.Introduction

    In modern turbomachinery,the constantly increasing turbine inlet temperature proposes stringent challenge to the design,and advanced cooling technology is heavily employed to protect the high temperature turbine components.1To reduce the NOxemission and preserve the combustion stability,swirl generator is widely used in lean-burn combustor.2Inside the combustor,the swirling flow is not fully smeared and contributes to the non-uniform flow field at the combustorturbine interface.As a result,the real inlet flow condition of turbine is non-uniform and the non-uniformity would have negative impact on the aerodynamics and heat transfer features of the turbine.Consequently,the investigation on inlet swirl and Inlet Guide Vane(IGV)endwall film-cooling3is exceedingly indispensable.

    In the last decade,several related studies concerning the combustor-turbine interactions were made.Shih and Lin4studied the interaction mechanism between inlet swirl and leading-edge airfoil fillet.They controlled secondary flow using these two methods to reduce aerodynamic loss and vane surface heat transfer and found that the swirls and fillets increased the size of the stagnation region on the endwall around the airfoil’s leading edge.Fillets accomplished this by geometry

    Nomenclature

    contouring and swirl flow achieved this by more kinetic energy.But they did not discuss the effect of inlet swirl on the film covering on the endwall.Researchers in Oxford University proposed the Integrated Combustor Vane(ICV)concept,Rosic et al.5investigated numerically the interaction between the combustor and the first vane in an industrial turbine inside which 16 can-com bus tors and 32 vanes were installed.It was claimed that the combustor walls had dreadful effects on the first vane film cooling.Promising combustor-vane integration minimized the axial distance between the combustor and vane.In addition,the cooling and aerodynamics performance of the turbo system was not detrimentally affected.Beneficial in reducing secondary flows in the ICV vane,the integrated design also saved up to 25%of the total coolant.Jacobi and Rosic6compared a ero thermal experimental data and their inhouse Computational Fluid Dynamics(CFD)code TBLOCK simulation to analyse the decrease in total pressure loss coefficient and heat transfer coefficient levels of integrated vanes.Their method halved the number of vanes,eliminated the coolant need of leading edge and considerably reduced manufacturing and development costs.Furthermore,Khanal et al.studied the effect of combustor hot streak and swirl on the unsteady aerothermal performance of high pressure turbine stage based on MT1,7then proposed a combined performance parameter which enabled aerodynamics and heat transfer behaviour to be quanti fied and marked in their solver HYDRA.According to He et al.8,the combined interplay of hot-streak and swirl resulted in crucial variations in heat transfer and aerodynamic performance.Qureshi et al.9investigated the effect of residual swirl from combustor exit on thermodynamics by comparing with and without inlet swirl.Inside their research,experimental part presented time-averaged heat transfer and static pressure measurement in an unshrouded turbine,and numerical simulation part based on Rolls-Royce in-house code HYDRA also were proceeded.In view of Oxford colleagues’research,Beard et al.10dig deeper to study the relationship between inlet swirl and stage efficiency.They conducted the experiment about the effect of swirl on the turbine stage efficiency with an efficiency measurement system capable of resolving efficiency changes within±0.16%.However,their current research hardly focused on endwall film cooling which is very significant in modern gas turbine system.Yang et al.11predicted film cooling and heat transfer characteristic by numeric al simulations with consideration of stator-rotor purge flow and discrete hole flows.The calculation was carried out based on Reynolds stress turbulence model with a non-equilibrium wall function.Las kowski et al.12conducted several research cases about turbine vanes with leading edge film cooling in the conjugate calculation method which can be in analogy with endw all film cooling.Dong et al.13constructed a CFD model to study the impact of vane wakes and leading edge bow wakes on ingestion,meanwhile,predicted turbine wheels pace cooling flow interactions with transonic hot gas path.These simulation research was meaningful in the study of industrial complicated calculations,especially for actual turbine working conditions,such as inlet swirl.

    As for experimental investigation,Cha et al.14measured temperature and flow fields of full-annual,rich-burn com bus tors and high pressure turbines using passive scalar tracing methods which were isothermal and non-reacting.The complementary results generally provided a datasets benchmark for conventional turbines design in the future.Koupper et al.15developed a kind of engine representative combustor simulator to generate hot streak pro file.This trisector rig was equipped with an automatic control system which allowed the facility to operate at stable situations for long time runs.The research results mentioned above were conducted on the industrial sized rigs,it is more difficult to monitor the exact mechanism about the flow phenomenon,than the same experiments on simple cascades.Especially for catalytic and dry low NOxcombustor,Ames et al.16measured endwall heat transfer distributions under turbulence from the combustor exit in a large-scale,low speed linear cascade facility.They made the assumption that the individual vortex had less evident influences.With consideration of a large range of Re numbers and turbulence conditions,the database of the exercise could be an assessment of endwall thermal simulation capabilities.Hedlund et al.17observed heat transfer and flow phenomena in a swirl combustor and relative turbine blade passage,and found that surface Nusselt numbers and time-averaged flow characteristics had important variations because of arrays of passage vortex pairs.Papa et al.18used naphthalene-saturated air and oil-dot visual-ization to measure the mass/heat transfer on the endwall surface in a linear cascade.The experimental results by Barringer et al.19indicated that the inlet total pressure pro files affected the aerodynamics loading by as much as 10%,revealing that the combination of different total pressure and total temperature pro files had significant effect.The bold exercise was conducted with the endwall cooling hole sets located along is ovelocity lines,this idea was further analogized to holes sets decorated along is o-pressure lines in this paper.

    Some novel experimental methodologies and new test facilities were developed to measure more information of a erothermal fields in detail.Luque et al.20,21proposed a new technique to accurately evaluate cooling system performance.Taking engine-representative geometry,and internal and external cooling flow into consideration,the facility based on thermochromics liquid crystal provided a totally integrated approach to assessing the cooling system.Besides the above method,Luque et al.also made use of five-hole probes and transient infrared thermography technique to monitor flow information such as:turbulence,swirl,leakage flow,shed wake distributions,and surface coolant flows.Krichbaum et al.22introduced an extensively redesigned large scale turbine rig,which was a full annular 1.5-stage axial low speed turbine and equipped with combustor simulator modules including swirl generators.

    As is known,endwall is always the emphasis of turbine cooling.Due to the difficulty of extremely accurate CFD prediction of endwall film cooling phenomenon,most of the investigation on endwall film cooling were conducted based on experimental measurements.Zhang and Jaiswal23compared the cooling results of double-row hole injection and singlerow discrete slot injection using Pressure Sensitive Paint(PSP).They argued the secondary flow which dominated the near endwall flow is at low mass flow ratio.On the contrary,the cooling film flow with a higher momentum dominated the near wall flow field and suppressed the secondary flow at higher mass flow ratios.Still based on PSP technique,Li et al.24investigated experimental impacts of swirl purge flows on platform cooling and phantom cooling of suction surface.The main parameters on cooling effectiveness were mass flow ratio and swirl ratio,which could change coolant distribution on platform close to suction side and phantom cooling result.On a multi-stage rotating test facility,Suryanarayanan et al.imitated real rotating conditions to measure cooling effectiveness on blade platform25,and suggested that film cooling effectiveness increased with an increase in the coolant-to mainstream blowing ratio(M)for different all rotating speeds,respectively.26

    So far,related research literatures27mainly focused on the effect analysis of interaction between inlet non-uniform flow and the whole turbine aerothermal performance in the downstream.But few investigations were about the most important heat transfer feature,turbine component cooling.Moreover,scant studies about endwall film-cooling hole arrangement criterions were searchable in this field.These literatures did not explain clearly the mechanism of swirl effect on film cooling.There were almost no related research about the effect of inlet swirl on endwall cooling and how to suppress this turbulence influence with the method of changing holes arrangement strategy.Besides,limited by calculation ability,it must be admitted that current numerical simulation had less accuracy and less precision than experimental measurements.Thus,using experimental methods to clarify the effect of inlet swirl on endwall cooling is significant.A more in-depth arrangement plan comparing two representative cooling-hole to expound the reason of the difference is quite meaningful.Based on our view that specially designed holes arrangement could be beneficial to cooling effectiveness and the development of cooling film,it was in novatively inferred that new holes arrangement under inlet swirl conditions,are supposed to contribute better cooling effectiveness.In this paper,to resist the film lift-off effect of inlet flow non-uniformities and protect the endwall from being ablated by high temperature gas,a new hole arrangement was optimized and chosen by CFD simulations according to Design Guidelines that film-cooling holes lay along the isobars on the endwall.As for this experimental research,we used a clockwise swirl generator model embedded in a removable plexiglas plate to imitate swirl flow induced from the combustor in the real turbine set.The cooling effectiveness were measured under variable blowing ratios on two typical endwalls,then the mechanism of result physical phenomenon were analysed.The novel discoveries are that inlet swirl has hugely devastating effect on endwall film coverage,and the endwall on which holes lie on isobars has strong resistance against swirl effect.

    As for the structure of this paper,experimental theory and test rig facility will be introduced firstly,then measurement plan in different groups will be clari fied;the final part is conclusion and discussion about the results.

    2.Experimental methodology

    In this experiment,PSP technique based on oxygen-quenched photoluminescence was used as the measurement method.PSP is a kind of contactless measurement paint,with the reacting theory resolved by Zhang and Jaiwal.23In current experiments,Light-Emitting Diode(LED)lights whose wavelength are 450 nm are used to activate PSP so that PSP could emit reflected light,then the light will be captured by a Charge-Coupled Device(CCD)camera with high spectral sensitivity.The CCD camera has a band pass filter whose traversing wavelength is 600 nm,to catch the 600 nm reflected light by PSP in the energized state.In this study,the main stream is air which contains about 21%oxygen.The coolant is industrial nitrogen for which the partial pressure of oxygen is almost 0.With the consideration of the analogy between mass transfer and heat transfer,oxygen concentration and partial pressure could express the film effectiveness in terms of Eq.(1)as below:

    where η means the film cooling effectiveness,Cairrepresents the oxygen concentration of the main stream,Cmixrepresents the oxygen concentration of the mixture of air/nitrogen,CN2represents the oxygen concentration of the nitrogen,which equals to 0.Meanwhile,(PO2)airrepresents the oxygen partial pressure of the main stream,(PO2)mixrepresents the oxygen partial pressure of the mixture of air/nitrogen.With this technique,the measured film cooling effectiveness ranges between 0%(at far upstream and far downstream locations)and 100%(inside the film cooling hole).

    In each PSP film-cooling test,CCD camera takes 4 pictures at the same main stream temperature value to measure the film-cooling effectiveness accurately.28The No.1 picture is a‘dark’image that is photographed without LED light and is the main stream flow condition.The No.2 picture is the‘reference’image photographed without the main stream flow,but with the LED light on.The No.3 picture is ‘a(chǎn)ir’image is photographed with both the main stream flow and LED light on,and the coolant is totally air injection as mainstream.The No.4 picture is ‘N2’image photographed in the same condition as No.3,except that the coolant is pure nitrogen injection instead of air.By visualization program,the effectiveness distribution could be obtained from the four figures for every test sequence.The information on the change in oxygen concentration due to the change in oxygen partial pressure,which cannot affect the magnitude of local film-cooling effectiveness.This could be provided by the distribution of reference ratio data in the‘a(chǎn)ir’image.Compared with the reference ratio data,the other one from the ‘nitrogen’image could make contributions to effectiveness by absolute oxygen concentrations.Referring to the analogy between mass transfer and heat transfer,the post algorithm obtains the local adiabatic film-cooling effectiveness distributions by these two groups of reference ratio distributions.

    It is necessary to calibrate the PSP employed along the test surface to obtain the relationships between light intensities and local partial pressures of oxygen,before every experimental measurement.For this goal,a calibration system is established which is shown in Fig.1.During the calibration of the PSP system,a PSP-coated copper coupon is employed as the model of the experimental surface,underneath which three thermocouples are mounted to monitor the unsteady change of surface temperature.This PSP-coated copper coupon is placed in a partial or total vacuum room,sealed by machinery.A heater located at the back of the coupon increases both the temperature of the copper coupon and the heater to keep the desired temperature of the sample with accuracy better than±0.5 K,continuously.Through a transparent plexiglas window,the sample coupon is in the vision of CCD camera.Since the surface temperatures in this system varies from 298 K to 308 K,the calibrations are conducted for surface temperatures of 276.5 K,298 K,and 308 K,in addition,for pressures from partial vacuum(almost 1/3 of an atmosphere)to one standard atmosphere,respectively.These three temperature points of calibration undertaken make sure the completeness of influence information of temperature variation.The interdependence between the light intensity factor and pressure factor is shown in Fig.1,which shows the calibration result curve.In this figure,the abscissa is dimensionless relative pressure ratio,where P means local nitrogen partial pressure,and Prefmeans the reference nitrogen partial pressure.The ordinate is dimensionless relative light intensity,where I means local light intensity,and Irefmeans reference light intensity.It is obvious that there is little dependence of calibration data on surface temperature,and the quantitative change of the three data sets of different temperature values are tiny.

    Based on the calibration results,the conventional dimensionless temperature around the cooling areas also decided by light intensity,could be expressed in terms of Eq.(2)as below:

    where θ means dimensionless relative temperature,T∞means temperature of freestream,Tcmeans temperature of injected coolant,T means local temperature.

    For constant property and low speed flow,after importing the adiabatic wall temperature parameter,the equation yields dimensionless film-cooling effectiveness.Furthermore,the final adiabatic film effectiveness could be represented in terms of Eq.(3)as below:

    where Tawmeans the temperature of adiabatic wall,η means the same cooling effectiveness as in Eq.(1).

    As for the uncertainty and tolerance,the confidence of the calibration is on 95%extent.The measurement relative tolerance of the adiabatic film-cooling effectiveness is almost 3%,and the magnitude of uncertainty value is less than 0.5.However,if the effectiveness is close to 0,the non-confidence will increase,even to a magnitude of almost 20%when the value is about 0.05.The cumulative uncertainty is the superposition result from calibration error(4%)and deviation of capturing pictures(1%).Thus,the absolute error of effectiveness measurement changes during 0.01–0.02.As for extreme assumption,maybe the relative uncertainty for extremely low effectiveness(0.01)would be exceedingly high,about 100%.

    3.Experimental facility

    The test facility is composed of an inlet segment,a linear turbine cascade and an exit section.In the experiments,the free-stream air velocity is 24.3 m/s,which corresponds to a Mach number of 0.075.In the measurement section,four GE-E3guide vanes are mounted.The tested cascade29has a scaling ratio of 2.2,with its chord length equal to 79 mm and height 129 mm.In the current study,the effect of the non-uniform flow field at the combustor-turbine interface is generated with a swirl generator,as demonstrated in Fig.2.

    Also,for comparative study,the swirl generator module is removable and in this manner the experiment with uniform freestream can be conducted very easily.The swirl intensity is defined as Swirl Number in the terms of Eq.(4)as below,which is a widely used parameter in both research field and engineering:

    where S represents Swirl Number,Rinrepresents the inside radius of swirl generator,Routrepresents the outside radius,and the φ represents the mounting angles of blades in the swirl generator.As for the current con figuration as given Fig.2,the inside radius is 14 mm and the outside radius is 48 mm,with which the swirl number is about 1 and this swirl intensity is representative in modern lean-burn combustor.

    During the entire experimental process,the swirl generator was installed at the same position,the accurately middle of the span wise test-section width,to keep the swirl flow conditions constant.Almost all components of the test section can be replaced including endwall parts,swirl generator,cooled vanes and coolant supply cavities which is installed below endwall plane.Since the holes arrangement on the endwall is the primary concern during this research,no other components except two representative endwall parts is changed.

    The test rig30is shown in Fig.3(the real system diagram)and Fig.4(the schematic diagram).

    To compare the film-cooling effectiveness distribution accurately,the relative position between the CCD camera and the main endwall module is held stable constantly during the experiments.A high power compressor supplies enough coolant air for the secondary air system,and the compressor delivers coolant into cavities located under the wind tunnel test section.To realize the practical conditions inside industry turbine as real as possible,the structure of coolant cavities was designed as below:all secondary air supply share one communal passage which is connected with Air or Nitrogen gas cylinder,the total mass flow could be controlled by flow meter.After the shared general coolant passage,several cavities are linked between coolant supply outlet and cooling holes inlet,respectively.With the consideration of processing convenience,each row of cooling holes on endwall is allocated with one divided cavity,by which the coolant could go spouted outside the fan-shaped holes.Since all cavities share one coolant supply which means the same total pressure,the distribution of coolant would be dominated by pressure difference of each cavity,instead of active interference.During the entire experiments to all tested cases,the total coolant mass flow was changed to simulate mutative M,which is a macro variable,and no more active control on coolant distribution was adopted.Above is the working mechanism of secondary air system on this test rig,which is shown in Fig.4(a).Consequently,the coolant supply structure and working principle are permanent,no matter which tested part is installed.

    As for the hole shapes,it is widely known that the jet de flection of fan-shaped holes could result in a wider coolant diffusion area.31Therefore,typical fan-shaped holes with the geometry features shown in Fig.5 were chosen as the holes on the experimental endwalls.

    There are two kinds of endwalls used in the test.The first one is called ‘Type A’endwall.Taking the vortex sweeping coolant from endwall to the mid-span region into consideration,all holes are arranged in straight lines perpendicular to the axial to re flect the blowing-off phenomenon more clearly.To cover the gill region fully,the cooling structure consists of four rows of fan-shaped film-cooling holes.This endwall is shown in Fig.6(a).As shown in the figure,the first row of holes is fixed at upstream of the vane leading-edge,and the other three rows of holes are located in the passage.The angle between the axial of fan-shaped holes and platform surface is 30°,and the compound angles for the main axial direction are 0°,30°,45°and 60°for the four rows of holes,respectively.Both the lateral expansion and forward expansion into the platform are 10°.The hole diameter of the cylindrical part below the expansion is 1 mm,as suggested in open literatures and used widely in industry.

    The second one is called ‘Type B’endwall,for the purpose of resisting the film blow-off effect.With the consideration of the same outlet back pressure making sure the coolant flow evenly,all holes are arranged along the isobars except the first row.Under the permission of structural strength,coolant flows out from five rows of holes to protect the whole gill region,as shown in Fig.6(b).The first row is located upstream of the leading-edge,and the row is separated to be two parts,the part around the pressure-side of channel is moved toward downstream to obtain better cooling effectiveness.The other holes is fixed up with both the consideration of isobars and the structural strength,so this endwall could be manufactured for industry.As shown in the figure,the second row includes ten holes,with three discrete holes located as a triangle near the pressure-side among the ten.The third row of holes starts from the head of suction-side in the passage and ends at the mid-Caxof pressure-side in the passage.The fourth and fifth rows are arranged downstream,nearly parallel to the main axial direction of the turbine stage.Overall,the cooling holes arrangement could be seen in the figure.For every single hole,the geometry size and expansion angles of the ‘Type A’endwall are the same as the ‘Type B’endwall.Under the platform with the holes,every row is installed correspondingly with a coolant cavity to simulate the secondary-air system,and all cavities share one coolant source,like real turbine cooling condition.The coolant flow could be changed by rotameters connected with the cavities and coolant source.

    The information about the geometric and position data of film holes on ‘Type A’and ‘Type B’endwall is shown in Tables 1 and 2,respectively.In these tables,X/Caxmeans relative position of each hole,Caxmeans actual chord length of scaled-up blade pro file at axial chord direction.D means diameter of each cooling hole.Meanwhile,the coolant flow in the experiments is listed in Table 3,the detail information about the geometric and flow conditions of test rig is shown in Table 4.In Table 3,M represents blowing ratio in the experiments plan.

    4.Results and discussion

    In this study,the two types of endwalls with different film cooling hole arrangements were tested with both uniform and swirling inlet conditions to study their performance within the practical turbine environment.A broad range of M value from 0.2 to 1.3 were investigated in the experiments,with emphasis on the off-design non-uniform performance of the endwall film cooling.

    4.1.Cooling effect difference between two endwalls without swirl

    The contours of film-cooling effectiveness on two gill regions with uniform freestream are shown in Figs.7–10.For the series of cooling effectiveness contours,abscissa is dimensionless chord length,X/Caxmeans the ratio of axial position and actual chord length of scaled-up blade pro file at axial chord direction.Ordinate is dimensionless pitchwise length,ZPmeans length of one pitch at the pitchwise direction.

    It is pretty obvious that the ‘Type B’endwall could obtain a more satis fied cooling effectiveness distribution under the same conditions.For the designed working M situation,Fig.7 shows film-cooling effectiveness distribution on ‘Type B’endwall without inlet swirl when the M is 0.7.It is obvious in the image that almost the whole gill region around film-cooling holes could be protected by coolant coverage,except the single holes located closest to leading edge of pressure side,where the outlet backpressure on the platform surface is too high toallow the coolant flow out heavily.Thus,this area of gill region near leading-edge of pressure-side is crucial to the cooling blind area.With consideration of the blind area,the total outcome of film-cooling effectiveness is pretty satisfying.There is almost noplace without coolant coverage in the downstream of the first row.Fig.8 shows film-cooling effectiveness distribution on ‘Type A’endwall without inlet swirl when M is 0.7,and the film-cooling effectiveness distribution is also satisfying and there forms a successive coverage.However,besides the blind area around leading-edge of pressure-side,some area around holes in the middle of first and second rows cannot supply enough coolant to protect the endwall surface,especially for the 3rd,4th,6th,7th,9th hole in the second row.During the second row and the third row,the coolant film is thin and the cooling effectiveness is lower than downstream behind the third row.Figs.9 and 10 show film-cooling effectiveness distribution on ‘Type B’and ‘Type A’endwall without inlet swirl when the M is 0.3,respectively.This is for further off-designed work conditions,with decreasing coolant out flow.Compared with Figs.7 and 8,an intensive difference of filmcooling effectiveness could be observed.For ‘Type B’endwall,there still is a successive film of coolant,although one single hole in the blind area cannot allow coolant flow out heavily and the effectiveness of some areas near the pressure-side is somewhat poor.In general,the effect of film-cooling is pretty good,even M is so low as 0.3.In contrast,for ‘Type A’endwall,the film of coolant is not continuous when M is low,as shown in Fig.10.The pitchwise lines at 20%Caxand 75%Caxare chosen to reflect quantitatively the cooling effectiveness of upstream and downstream,and the line locations are shown in Figs.7–10,as 1–1,2–2,3–3,and 4–4.The cooling effectiveness distributions along pitchwise lines are shown in

    Figs.11 and 12.For these series of curve graphs,abscissa is pitchwise location,where Span means one pitch size,Z means local pitchwise coordinate.Ordinate is cooling effectiveness.At the upstream,the film-cooling effectiveness of the‘Type B’endwall is overwhelmingly more favourable than the ‘Type A’endwall.At the downstream,the averaged film-cooling effectiveness of‘Type A’endwall is slightly more positive than‘Type B’endwall.However,the effectiveness value of the ‘Type B’endwall is close to 0.3.This value is near that of the ‘Type A’endwall and is satis fied.Compared with Figs.11 and 12,it is found that the upstream effectiveness value is close to the downstream for the ‘Type B’endwall,the value is around 0.35,meaning the coolant flow out is generally very uniform.On the contrary,the situation for ‘Type A’is unsatisfied.The upstream effectiveness value is around 0.2 and the downstream effectiveness value is around 0.4,with consideration of effectiveness contours,obtaining the result that the coolant flow out is not uniform at all.In Fig.12,there are two peaks on the curve of‘Type B’endwall,because the 75%Caxpitchwise line passes through film-cooling holes on the platform.Thus,the maximum effectiveness appears without significant

    meaning,the averaged effectiveness is more persuasive,therefore,‘Type B’endwall has more homogeneous cooling effect.

    Table 2 Film holes location and orientation for ‘Type B’endwall.

    Table 3 Injection coolant flow in experiments for two endwalls(irrelevant to swirl).

    Table 4 Geometric and flow conditions.

    To sum the performance comparison between two endwalls without inlet swirl, ‘Type B’generally overwhelms ‘Type A’,according to the detailed experimental results.It is mainly because with almost all holes lay on the isobars,the pressure gradients between secondary-air pressure and outlet back pressure on the isobars are close,so the flow out of cooling holes are more uniform,this phenomenon will lead to less coolant momentum difference,which means less interaction between coolant secondary flow and free stream,retaining film in longer distance and bigger lateral expansion area for shaped holes.More uniform coolant flow allocation also avoids the problem that some holes cannot appear unobstructed coolant flow.On the other side,for the ‘Type B’endwall,due to smaller distance between adjacent rows than ‘Type A’,the coolant from upstream rows is able to survive until covering the near downstream rows,this overlay consequent also makes contribution to better cooling effect for ‘Type B’endwall.

    4.2.Effect of swirl on endwall cooling for each endwall

    The effects of inlet non-uniform swirl to endwall film cooling are key points of this research.Two different endwall cooling structure design forms were studied respectively under representative swirl conditions whose intensity is equal to 1.

    With the inlet swirl,generally on the endwall,there appears a massive blow away instead of continuous film.Fig.13 is the film-cooling effectiveness distribution on ‘Type A’endwall with inlet swirl of the M equal to 0.7,corresponding to the result without swirl shown in Fig.8.The comparison indicates the huge adverse influence of swirl disturbance when M is high.Under this situation,the effect of film cooling is unsatis fied,and there in no clear connection between adjacent rows,meaning the coolant film is not successive at all.Fig.14 is filmcooling effectiveness distribution on ‘Type A’endwall with inlet swirl of M equal to 0.3,corresponding to Fig.10.The comparison shows that inlet swirl has a devastating influence on the cooling effectiveness.For instance,half of the rows of holes cannot allow coolant to flow out,which means the mainstream with high temperatures would flow backward into endwall to melt the components.Laterally averaged film-cooling effectiveness distribution along the axial direction on the‘Type A’endwall with different M values without inlet swirl is shown in Fig.15.In the figure,the result of the averaged cooling effectiveness gets decreasing as M decreases,monotonically.The infomation with inlet swirl is shown in Fig.15,too.It is evident that the swirl has a huge influence on end wall film cooling effectiveness,especially at the upstream of the gill region,the effectiveness is sensitive to swirl,and the sensitivity at downstream is slight,relatively.Also,the comparation clarifies there exists over-cooled accumulation coolant in downstream area,on which inlet swirl has just slight effect.

    Fig.16 is the film-cooling effectiveness distribution on‘Type B’endwall with the inlet swirl and M equals to 0.7,corresponding to Fig.7.The comparison indicates that the inlet swirl also has a damage impact on the cooling effectiveness on ‘Type B’endwall.Although the whole pattern of coolant film is not changed,the coolant traces behind the cooling holes almost disappears,instead of long wakes in the Fig.7.Fig.17 is the film-cooling effectiveness distribution on ‘Type B’endwall with the inlet swirl of M equal to 0.3,corresponding to Fig.9.The difference between Figs.9 and 17 suggests the swirl impact on cooling effectiveness with low M.It is observed that four holes in the first row inside the blind area cannot flow coolant,and the film near pressure-side is nothing.Under strong inlet swirl,the general film cooling effectiveness distribution is somewhat terrible.Laterally averaged film-cooling effectiveness distribution along the axial direction on the‘Type B’endwall with different M values without inlet swirl is shown in Fig.18.In the figure the trend of averaged cooling effectiveness is as low as the M,monotonically.The same information with inlet swirl is shown in Fig.18,too.

    Aaccording to Fig.18,it is claimed that the swirl has considerable effects on the endwall film cooling effectiveness.However,due to uniform coolant allocation of‘Type B’endwall,the changing trends under one specific M are approximately similar.It is found that for the ‘Type B’endwalls on which the holes are along isobars,only under the low M and with swirl disturbed,the averaged cooling effectiveness is poor(lower than 0.2).The film effectiveness of other cases is still satisfying.

    According to the experimental results and film lift-off theory,the mechanism could be modelled as below.Fig.19 shows the lift-off physical process that normally happens several D distance far from the ejection hole.The mainstream with momentum entrainment always conducts reciprocal momentum transfer with coolant kidney vortex pair, finally the vortex pair grow up and get lift-off,then no coolant exists between endwall and freestream hot gas,which means cooling film blows away.

    Taking swirl perturbation into consideration,it is clear that the mutual effect between coolant and hot gas will be more intense.Due to powerful swirl flow,the secondary flow intensity and momentum entrainment of freestream will be more.Freestream has an extremely strong trend to merge vortex pair of coolant,and mixes with coolant vortices to conduct momentum and energy transfer, finally,one side of vortex pair of coolant will get enough energy and momentum to detach from endwall earlier,thus the balance of coolant vortex pair is broken.As the result, film coverage also cannot be maintained immediately.The final outcome is freestream and coolant vortices get blended together,like swirl flow sweeps coolant film away from endwall.It is the sweeping and blending effect of swirl perturbation that destroy the film coverage on the endwall.The sweeping increases the mixing process between mainstream and coolant ejection.During the interaction,swirl transfers more momentum to cooling ejection,these movements have a negative effect on film attachment and could damage cooling effectiveness.In the current paper,it could be found that ‘Type B’endwall has more resistance against swirl flow than traditional‘Type A’endwall,less sensitivities to inlet swirl,because of the uniform coolant allocation and holes laying on isobars.

    4.3.Cooling effect difference between two endwalls with inlet swirl

    Based on the above measurement data,the experiment results with consideration of both swirl and different endwall showed that under the same intensity of disturbance,‘Type B’endwall has stronger resistance against disturbance and the film coverage effectiveness is relatively high even after swirl disturbance.For both two endwalls,it is difficult for the coolant to flow outside the upstream area when M is low.However,there also is severe excessive cooling downstream on ‘Type B’.This is because the ‘Type B’endwall has less holes downstream where pressure is fairly low.The cost of coolant near downstream is less,so the upstream endwall gets more coolant allocation,the distribution of coolant can be more averaged overall.In upstream part,the coolant trail from the first row can cover the next row due to the fact that the holes are distributed in a high density.As a result,this superposition effect makes the better cooling effect.

    Figs.13 and 16 show the film-cooling effectiveness distribution on ‘Type A’and ‘Type B’endwall with inlet swirl of M equal to 0.7,respectively.The comparison suggests that both the endwalls could be protected by successive coolant film.For ‘Type A’endwall the situation is more adverse.For instance,the cooling blind area around pressure-side and the area between adjacent rows of holes is hardly covered with film due to short of coolant.The momentum of coolant injection will join the convection and exchange momentum with swirl flow and mainstream.As a result,the cooling film would be swept by secondary flow in the passage.Figs.20 and 21 compare the film-cooling effectiveness distribution on ‘Type B’and‘Type A’endwall with inlet swirl of a mild M value,which is 0.5.

    Furthermore,Figs.14 and 17 compare the film-cooling effectiveness distribution on ‘Type B’and ‘Type A’endwall respectively with inlet swirl of a low M value,which is 0.3.The tendency is clear that the film cooling effectiveness is getting lower as M gets lower,and the effectiveness distribution of‘Type B’endwall is more favourable than ‘Type A’endwall.When the M is extremely low at the value of 0.2,as shown in Figs.22 and 23,both distributions of cooling effectiveness are terrible under inlet swirl interference.For ‘Type B’endwall,the large area around the cooling blind area(dimensionless Pitch from 0.6 to 0.9,Caxfrom-0.2 to 0.6)almost has nothing in the coverage,and the situation is more dreadful for ‘Type A’endwall.On the ‘Type A’endwall,only three holes have the ability of supplying coolant out flow normally in the first row,and half of the holes in the second row cannot work as designed.Fig.24 shows pitchwise distribution of filmcooling effectiveness at 20%Caxwith inlet swirl when M is 0.5.Fig.25 is the same distribution of film-cooling effectiveness at 75%Caxwith inlet swirl when M is 0.5.Except the two peaks in the curves of‘Type B’endwall shown in Figs.24 and 25,for‘TypeB’endwalltheaveraged cooling effectiveness at upstream is approximately 0.3,and the value at downstream is nearly 0.2.In contrast,for the ‘Type A’endwall the averaged effectivenss at upstream is lower than 0.1,and the value at downstream is still about 0.25.Although the downstream of‘Type A’endwall is higher than ‘Type B’endwall,with the consideration of the accumulation of the coolant downstream,the effectiveness(0.25 level)is not necessary,and the effectiveness of‘Type B’endwall(equal to 0.2)is enough to protect the endwall.As for the paramount of cooling,the upstream of‘Type A’endwall has poor cooling effectiveness,but the situation of‘Type B’endwall is fairly promising.

    To sum up the charcteristic of‘Type B’endwall and the reasons why the cooling is better,the keypoints are the uniform coolant allocation and overlay priciple.Inlet swirl does have a huge impact on film-cooling effectiveness,however,under high M values,the film coverage is not fragile and on endwall surfaces still form continuous coolant film.It means changes of cooling effect is not very sensitive to swirl.But when M is extremely low,the film coverage will be susceptible to swirl effect,coolant will blow away and be lift-off,the changes of cooling effect is sensitive to swirl disturbance at that time.

    5.Conclusions

    This research studied the effect of inlet swirl on the film effectiveness of‘Type B’endwall and ‘Type A’endwall under different M values by a test cascade tunnel rig.The difference of test conditions with varying M values and inlet uniformities indicates the reasons why film-cooling effectiveness is low somewhere.Film effectiveness is measured by pressure sensitive paint technology.In the experiments,the complete endwallsurfaceismainly compared with contentsunder different cooling structure protection,with the M varying from 1.3 to 0.2.In general,the inlet swirl has a huge in fluence on the endwall film-cooling.The coolant ejection and mainstream with strong secondary flow system get mixing together near the platform surface,therefore,sweeping the coolant around the mixture fluid vortex.

    From this work,the following conclusions can be drawn:

    (1)Hole arrangements of the endwall with holes along isobars have the ability to resist the inlet swirl disturbance and to protect the endwall surface under strong swirl conditions.

    (2)The sensitivity of distribution of film cooling effectiveness on the endwall is low when the M is high.In contrast,the sensitivity is high when M is fairly low.

    (3)The film cooling effectiveness of endwall with cooling holes along isobars is better than endwall with holes along straight lines vertical to axial,with other conditions remaining the same.

    In the further step,more non-uniform features in the actual gas turbine will be investigated in experimental method,and the affected objects will be included but not limited to blade surface and tip plane.

    Acknowledgements

    This study was supported by National Natural Science Foundation of China(Nos.51476082 and 51506107).The language writing assistance from Mr.Ross Boyd and Mr.Zach Kossow was very helpful to this paper.

    References

    1.Han J.Fundamental gas turbine heat transfer.J Therm Sci Eng Applicat 2013;5(2):021007–21015.

    2.Zhang M,Fu ZB,Lin YZ,Li JB.CFD study of NOxemissions in a model commercial aircraft engine combustor.Chin J Aeronaut 2012;25(6):854–63.

    3.Yang CF,Zhang JZ.In fluence of multi-hole arrangement on cooling film development.Chin J Aeronaut 2012;25(2):182–8.

    4.Shih TI,Lin YL.Controlling secondary- flow structure by leadingedge airfoil fillet and inlet swirl to reduce aerodynamic loss and surface heat transfer.J Turbomach 2003;125(1):48–56.

    5.Rosic B,Denton JD,Horlock JH,Uchida S.Integrated combustor and vane concept in gas turbines.J Turbomach 2011;134(3):031005–31010.

    6.Jacobi S,Rosic B.Development and aerothermal investigation of integrated combustor vane concept.J Turbomachin 2015;138(1):011002–11010.

    7.Rahim A,Khanal B,He L,Romero E.Effect of nozzle guide vane lean under in fluence of inlet temperature traverse.J Turbomach 2013;136(7):071001–71012.

    8.Khanal B,He L,Northall J,Adami P.Analysis of radial migration of hot-streak in swirling flow through high-pressure turbine stage.J Turbomach 2013;135(4):041005–41011.

    9.Qureshi I,Beretta A,Chana K,Povey T.Effect of aggressive inlet swirl on heat transfer and aerodynamics in an unshrouded transonic HP turbine.J Turbomach 2012;134(6):061023-11.

    10.Beard PF,Smith AD,Povey T.Effect of combustor swirl on transonic high pressure turbine efficiency.J Turbomach 2013;136(1):011002–12.

    11.Yang H,Gao Z,Chen HC,Han J,Schobeiri MT.Prediction of film cooling and heat transfer on a rotating blade platform with stator-rotor purge and discrete film-hole flows in a 1–1/2 turbine stage.J Turbomach 2009;131(4):041003–12.

    12.Dong P,Wang Q,Guo Z,Huang H,Feng G.Conjugate calculation of gas turbine vanes cooled with leading edge films.Chin J Aeronaut 2009;22(2):145–52.

    13.Laskowski GM,Bunker RS,Bailey JC,Ledezma G,Kapetanovic GM,Itzel GM,et al.An investigation of turbine wheelspace cooling flow interactions with a transonic hot gas path—part II:CFD simulations.J Turbomach 2011;133(4):041020,041020-12.

    14.Cha CM,Hong S,Ireland PT,Denman P,Savarianandam V.Experimental and numerical investigation of combustor-turbine interaction using an isothermal,nonreacting tracer.J Eng Gas Turb Power 2012;134(8):081501–18.

    15.Koupper C,Caciolli G,Gicquel L,Duchaine F,Bonneau G,Tarchi L,et al.Development of an engine representative combustor simulator dedicated to hot streak generation.J Turbomach 2014;136(11):111007–10.

    16.Ames FE,Barbot PA,Wang C.Effects of catalytic and dry low NOxcombustor turbulence on endwall heat transfer distributions.J Heat Transf 2005;127(4):414–24.

    17.Hedlund CR,Ligrani PM,Moon HK,Glezer B.Heat transfer and flow phenomena in a swirl chamber simulating turbine blade internal cooling.J Turbomach 1999;121(4):804–13.

    18.Papa M,Srinivasan V,Goldstein RJ.Film cooling effect of rotorstator purge flow on endwall heat/mass transfer.J Turbomach 2011;134(4):041014–8.

    19.Barringer MD,Thole KA,Polanka MD.Effects of combustor exit pro files on vane aerodynamic loading and heat transfer in a high pressure turbine.J Turbomach 2009;131(2):021008–10.

    20.Luque S,Povey T.A novel technique for assessing turbine cooling system performance.J Turbomach 2010;133(3):031013–31019.

    21.Luque S,Kanjirakkad V,Aslanidou I,Lubbock R,Rosic B,Uchida S.A new experimental facility to investigate combustorturbine interactions in gas turbines with multiple can combustors.J Eng Gas Turb Power 2015;137(5):051503–9.

    22.Krichbaum A,Werschnik H,Wilhelm M,Schiffer HP,Lehmann.A large scale turbine test rig for the investigation of high pressure turbine aerodynamics and heat transfer with variable in flow conditions.ASME turbo eExpo 2015:Turbine technical conference and exposition;2015 June 15–19;Montreal,Quebec,Canada.New York:ASME;2015.

    23.Zhang LJ,Jaiswal RS.Turbine nozzle endwall film cooling study using pressure-sensitive paint.J Turbomach 2001;123(4):730–8.

    24.Li SJ,Lee J,Han JC,Zhang LZ,Moon HK.Turbine platform cooling and blade suction surface phantom cooling from simulated swirl purge flow.J Turbomach 2016;138(8):081004–11.

    25.Suryanarayanan A,Ozturk B,Schobeiri MT,Han JC.Filmcooling effectiveness on a rotating turbine platform using pressure sensitive paint technique.J Turbomach 2010;132(4):041001–13.

    26.Suryanarayanan A,Mhetras SP,Schobeiri MT,Han JC.Filmcooling effectiveness on a rotating blade platform.J Turbomach 2008;131(1):011014-12.

    27.Chung H,Hong CW,Kim SH,Cho HH,Moon HK.Heat transfer measurement near endwall region of first stage gas turbine nozzle having platform misalignment at combustor-turbine interface.Int Commun Heat Mass Transf 2016;78(11):101–11.

    28.Zhang Y,Yuan X.Experimental investigation of nonperiodic endwall film cooling in neighboring passages with upstream rotating flow.ASME turbo expo 2014:Turbine technical conference and exposition;2014 June 16–20;Düsseldorf,Germany.New York:ASME;2014.

    29.Zhang Y,Yuan X.Endwall film cooling using the staggered combustor-turbine gap leakage flow.ASME turbo expo 2014:Turbine technical conference and exposition;2014 June 16–20;Düsseldorf,Germany.New York:ASME;2014.

    30.Zhang Y,Li Y,Yuan X.Effects of inlet swirl on suction side phantom cooling.ASME turbo expo 2016:Turbomachinery technical conference and exposition;2016 June 13–17;Seoul,South Korea.New York:ASME;2016.

    31.Shi HH,Kiriyama K,Itoh M.Effect of diffusively shaped holes on the turbulent flow field of a film cooling plate.Chin J Aeronaut 2002;15(1):6–11.

    国产精品99久久99久久久不卡| 精品欧美一区二区三区在线| 久久精品亚洲av国产电影网| 日本一区二区免费在线视频| 大香蕉久久成人网| 天天操日日干夜夜撸| 久久国产精品大桥未久av| 人妻 亚洲 视频| 人妻 亚洲 视频| 亚洲人成网站在线观看播放| 亚洲av在线观看美女高潮| 好男人电影高清在线观看| 丝袜在线中文字幕| 少妇 在线观看| 国产片特级美女逼逼视频| 久久国产亚洲av麻豆专区| 王馨瑶露胸无遮挡在线观看| 婷婷色综合www| 老司机亚洲免费影院| 老司机在亚洲福利影院| 一边摸一边做爽爽视频免费| 国产精品国产三级国产专区5o| 黄色视频在线播放观看不卡| 性色av一级| 婷婷色综合www| 热99国产精品久久久久久7| 精品少妇一区二区三区视频日本电影| 妹子高潮喷水视频| 久久精品国产a三级三级三级| 亚洲中文字幕日韩| 亚洲精品国产区一区二| 亚洲熟女毛片儿| 麻豆乱淫一区二区| 精品少妇一区二区三区视频日本电影| 老司机在亚洲福利影院| 欧美 亚洲 国产 日韩一| 国产免费一区二区三区四区乱码| 我要看黄色一级片免费的| 九色亚洲精品在线播放| avwww免费| 欧美日韩亚洲综合一区二区三区_| 黄色片一级片一级黄色片| 啦啦啦中文免费视频观看日本| 建设人人有责人人尽责人人享有的| 成人午夜精彩视频在线观看| 中国美女看黄片| 亚洲欧洲日产国产| 国产亚洲av高清不卡| 日韩欧美一区视频在线观看| 嫁个100分男人电影在线观看 | 亚洲国产最新在线播放| 国产成人系列免费观看| 亚洲欧美精品自产自拍| 久久ye,这里只有精品| 亚洲精品一卡2卡三卡4卡5卡 | 人妻一区二区av| 色精品久久人妻99蜜桃| av一本久久久久| 97人妻天天添夜夜摸| 青春草视频在线免费观看| 国产成人精品久久久久久| 一本—道久久a久久精品蜜桃钙片| 咕卡用的链子| 日本五十路高清| 国产1区2区3区精品| 欧美 日韩 精品 国产| 在线精品无人区一区二区三| bbb黄色大片| 国产成人精品久久二区二区免费| 久久久国产欧美日韩av| 老鸭窝网址在线观看| 黑人猛操日本美女一级片| 狂野欧美激情性bbbbbb| 一本一本久久a久久精品综合妖精| 免费高清在线观看视频在线观看| 久久国产亚洲av麻豆专区| 成人国产一区最新在线观看 | 欧美性长视频在线观看| 赤兔流量卡办理| 亚洲激情五月婷婷啪啪| 一级毛片女人18水好多 | 在线av久久热| 日韩 亚洲 欧美在线| 80岁老熟妇乱子伦牲交| 亚洲天堂av无毛| 日日夜夜操网爽| 中文字幕亚洲精品专区| 成人亚洲精品一区在线观看| 久久久国产一区二区| 亚洲一区中文字幕在线| 成人午夜精彩视频在线观看| 欧美日韩一级在线毛片| 午夜91福利影院| 国产野战对白在线观看| 亚洲国产欧美一区二区综合| 欧美大码av| 超碰97精品在线观看| 国产欧美日韩精品亚洲av| 69精品国产乱码久久久| 一区二区三区乱码不卡18| 丝袜喷水一区| 亚洲av国产av综合av卡| 看免费av毛片| tube8黄色片| 亚洲av日韩在线播放| av网站在线播放免费| 亚洲国产精品成人久久小说| 中文字幕人妻丝袜一区二区| 一本—道久久a久久精品蜜桃钙片| 成年人午夜在线观看视频| av欧美777| 中国美女看黄片| 精品一区二区三区av网在线观看 | 亚洲精品久久午夜乱码| 精品亚洲成国产av| 久久精品成人免费网站| 丝袜美足系列| 午夜免费观看性视频| 国产黄色免费在线视频| av国产久精品久网站免费入址| 国产日韩欧美亚洲二区| 亚洲av欧美aⅴ国产| 精品一区在线观看国产| 天天躁夜夜躁狠狠躁躁| 亚洲av日韩在线播放| 国产人伦9x9x在线观看| 黑人巨大精品欧美一区二区蜜桃| 国产伦人伦偷精品视频| 久久久久久久久免费视频了| videos熟女内射| 满18在线观看网站| 两个人看的免费小视频| 十八禁高潮呻吟视频| 肉色欧美久久久久久久蜜桃| 在线观看国产h片| 黄色 视频免费看| 久久久国产精品麻豆| 十分钟在线观看高清视频www| 国产片特级美女逼逼视频| kizo精华| 亚洲伊人久久精品综合| 丝袜脚勾引网站| 亚洲国产精品999| 亚洲欧洲国产日韩| 老司机亚洲免费影院| 91麻豆精品激情在线观看国产 | 久久青草综合色| 十八禁高潮呻吟视频| 精品福利永久在线观看| 国产成人av激情在线播放| 午夜精品国产一区二区电影| 天天操日日干夜夜撸| 亚洲情色 制服丝袜| 三上悠亚av全集在线观看| 黄频高清免费视频| 亚洲情色 制服丝袜| 老熟女久久久| 国产女主播在线喷水免费视频网站| 中国国产av一级| 韩国精品一区二区三区| 又大又爽又粗| 国产成人免费无遮挡视频| 精品亚洲成国产av| 欧美日韩视频精品一区| 免费日韩欧美在线观看| 国产精品.久久久| 久久久久视频综合| 国产精品久久久久久精品古装| 久久久国产精品麻豆| 欧美日韩亚洲国产一区二区在线观看 | 久久久久久人人人人人| 男女下面插进去视频免费观看| 国产精品一区二区在线不卡| 亚洲精品一区蜜桃| 电影成人av| 日本av手机在线免费观看| 亚洲少妇的诱惑av| 久久九九热精品免费| 热re99久久国产66热| 一区福利在线观看| 亚洲国产av影院在线观看| 中文乱码字字幕精品一区二区三区| 黄片小视频在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 校园人妻丝袜中文字幕| 另类精品久久| 在线观看人妻少妇| 男女床上黄色一级片免费看| 免费看十八禁软件| 一级片'在线观看视频| 手机成人av网站| 国产xxxxx性猛交| 免费观看人在逋| 高清av免费在线| 亚洲欧美激情在线| 亚洲成av片中文字幕在线观看| 免费看十八禁软件| 午夜免费鲁丝| 欧美+亚洲+日韩+国产| 久久精品人人爽人人爽视色| 青草久久国产| 一边亲一边摸免费视频| 又黄又粗又硬又大视频| 国产免费又黄又爽又色| 少妇裸体淫交视频免费看高清 | 亚洲av片天天在线观看| 在线亚洲精品国产二区图片欧美| 国产欧美日韩一区二区三 | 国产xxxxx性猛交| 亚洲一码二码三码区别大吗| 亚洲av成人精品一二三区| 丁香六月欧美| 国产高清视频在线播放一区 | 在线天堂中文资源库| 18在线观看网站| 高清av免费在线| 亚洲午夜精品一区,二区,三区| 亚洲国产日韩一区二区| 国产精品一区二区在线不卡| 亚洲国产看品久久| 久久这里只有精品19| 水蜜桃什么品种好| 国产亚洲精品久久久久5区| 久久影院123| 免费黄频网站在线观看国产| 亚洲欧美成人综合另类久久久| 欧美日韩精品网址| 亚洲九九香蕉| 亚洲精品av麻豆狂野| 日韩精品免费视频一区二区三区| 七月丁香在线播放| 亚洲专区中文字幕在线| 男女边摸边吃奶| 午夜日韩欧美国产| 国产在线视频一区二区| 黑人猛操日本美女一级片| 亚洲精品一区蜜桃| 精品国产乱码久久久久久小说| 国产精品99久久99久久久不卡| 亚洲成色77777| 久久影院123| 欧美av亚洲av综合av国产av| 一区在线观看完整版| www日本在线高清视频| 老汉色av国产亚洲站长工具| 欧美大码av| 天天躁狠狠躁夜夜躁狠狠躁| 日韩制服骚丝袜av| 黄网站色视频无遮挡免费观看| 精品国产乱码久久久久久男人| 亚洲av电影在线观看一区二区三区| 国产亚洲精品第一综合不卡| 婷婷色综合大香蕉| 97精品久久久久久久久久精品| 丝袜美腿诱惑在线| 好男人视频免费观看在线| 国产又色又爽无遮挡免| 国产精品人妻久久久影院| 妹子高潮喷水视频| av有码第一页| 欧美人与性动交α欧美精品济南到| 国产在线一区二区三区精| 19禁男女啪啪无遮挡网站| 日本午夜av视频| 欧美性长视频在线观看| 久久亚洲精品不卡| 欧美精品一区二区免费开放| 国语对白做爰xxxⅹ性视频网站| 日韩视频在线欧美| 国产精品成人在线| 99久久99久久久精品蜜桃| 精品人妻一区二区三区麻豆| 纵有疾风起免费观看全集完整版| av在线播放精品| 亚洲精品久久成人aⅴ小说| 亚洲国产日韩一区二区| 一本色道久久久久久精品综合| 国产老妇伦熟女老妇高清| 婷婷色av中文字幕| 90打野战视频偷拍视频| 国产极品粉嫩免费观看在线| 久久热在线av| 另类精品久久| 成人亚洲精品一区在线观看| 精品一区二区三区四区五区乱码 | 日韩中文字幕欧美一区二区 | 国产高清videossex| av不卡在线播放| 成年动漫av网址| 久久久亚洲精品成人影院| 熟女少妇亚洲综合色aaa.| 久久久久国产精品人妻一区二区| 午夜福利乱码中文字幕| 亚洲色图 男人天堂 中文字幕| 免费一级毛片在线播放高清视频 | 丁香六月欧美| 日韩,欧美,国产一区二区三区| 91麻豆av在线| 少妇粗大呻吟视频| 亚洲av日韩在线播放| 欧美成人精品欧美一级黄| 人人妻,人人澡人人爽秒播 | 男女免费视频国产| 国产精品亚洲av一区麻豆| 亚洲一区中文字幕在线| 一二三四社区在线视频社区8| 在线观看一区二区三区激情| 男女之事视频高清在线观看 | 国产精品久久久久久精品古装| 你懂的网址亚洲精品在线观看| 欧美久久黑人一区二区| 亚洲国产欧美日韩在线播放| 黑丝袜美女国产一区| svipshipincom国产片| 久久久久久免费高清国产稀缺| 亚洲精品美女久久av网站| 日本vs欧美在线观看视频| 日本欧美国产在线视频| 精品国产一区二区三区四区第35| 两人在一起打扑克的视频| 国产精品免费大片| 老司机在亚洲福利影院| 另类亚洲欧美激情| 国产精品成人在线| videosex国产| 欧美成人精品欧美一级黄| av线在线观看网站| 一边摸一边做爽爽视频免费| av在线播放精品| 精品一区二区三区av网在线观看 | 欧美日韩视频高清一区二区三区二| 香蕉国产在线看| 欧美97在线视频| 国产淫语在线视频| 九草在线视频观看| 色网站视频免费| 叶爱在线成人免费视频播放| 叶爱在线成人免费视频播放| 99国产精品一区二区蜜桃av | 啦啦啦啦在线视频资源| 黑人巨大精品欧美一区二区蜜桃| 午夜视频精品福利| 男人舔女人的私密视频| 建设人人有责人人尽责人人享有的| 黄网站色视频无遮挡免费观看| 99久久人妻综合| av在线app专区| 黑人巨大精品欧美一区二区蜜桃| 久久人人97超碰香蕉20202| 国产主播在线观看一区二区 | 国产视频首页在线观看| 啦啦啦在线免费观看视频4| 日日夜夜操网爽| 一边摸一边抽搐一进一出视频| 一级黄色大片毛片| 中文字幕制服av| 国产真人三级小视频在线观看| 亚洲国产看品久久| 久久久久精品人妻al黑| 男女免费视频国产| 国产黄频视频在线观看| 国产精品偷伦视频观看了| 亚洲激情五月婷婷啪啪| 日本av免费视频播放| 可以免费在线观看a视频的电影网站| 午夜久久久在线观看| 一区二区三区四区激情视频| 超碰成人久久| 亚洲熟女毛片儿| 久久久精品免费免费高清| 男女床上黄色一级片免费看| 免费女性裸体啪啪无遮挡网站| 日韩制服骚丝袜av| 欧美日韩国产mv在线观看视频| 国产欧美亚洲国产| 亚洲一码二码三码区别大吗| 国产人伦9x9x在线观看| 老汉色∧v一级毛片| 国产欧美亚洲国产| 日本vs欧美在线观看视频| 亚洲国产欧美网| 在线观看一区二区三区激情| 日韩人妻精品一区2区三区| 精品一区在线观看国产| 蜜桃国产av成人99| 十八禁网站网址无遮挡| 亚洲av成人精品一二三区| 少妇猛男粗大的猛烈进出视频| 国产亚洲av片在线观看秒播厂| 精品免费久久久久久久清纯 | 69精品国产乱码久久久| 精品一区二区三区av网在线观看 | 国产精品香港三级国产av潘金莲 | 丰满少妇做爰视频| 国产xxxxx性猛交| 亚洲av欧美aⅴ国产| 青春草亚洲视频在线观看| 欧美亚洲日本最大视频资源| 国产免费现黄频在线看| 欧美黄色淫秽网站| 欧美精品亚洲一区二区| 亚洲,欧美,日韩| 久久久精品国产亚洲av高清涩受| 天天躁日日躁夜夜躁夜夜| 少妇精品久久久久久久| 一级片免费观看大全| 亚洲欧洲精品一区二区精品久久久| 欧美亚洲日本最大视频资源| 午夜福利在线免费观看网站| 国产成人欧美在线观看 | 99久久99久久久精品蜜桃| 欧美在线黄色| 1024视频免费在线观看| xxxhd国产人妻xxx| 考比视频在线观看| 精品一区二区三区四区五区乱码 | 国产精品国产av在线观看| 亚洲熟女精品中文字幕| 国产99久久九九免费精品| 亚洲专区中文字幕在线| 午夜福利在线免费观看网站| 久久久久国产一级毛片高清牌| 久热爱精品视频在线9| 人妻一区二区av| 国产一区二区在线观看av| 国产99久久九九免费精品| 午夜福利视频在线观看免费| 最近最新中文字幕大全免费视频 | 青春草视频在线免费观看| 国产精品偷伦视频观看了| 午夜福利视频精品| 国产在线观看jvid| 这个男人来自地球电影免费观看| 91精品国产国语对白视频| 男女午夜视频在线观看| 一级片免费观看大全| 亚洲欧美日韩高清在线视频 | 国产成人免费观看mmmm| 国产精品久久久久成人av| 久久毛片免费看一区二区三区| 免费观看a级毛片全部| 我要看黄色一级片免费的| 18禁国产床啪视频网站| 欧美黄色淫秽网站| 男人爽女人下面视频在线观看| 一本一本久久a久久精品综合妖精| av不卡在线播放| 亚洲图色成人| 久久久亚洲精品成人影院| 啦啦啦在线免费观看视频4| 在线观看一区二区三区激情| 中文欧美无线码| 欧美精品啪啪一区二区三区 | 黄频高清免费视频| 欧美日韩视频高清一区二区三区二| 一区二区日韩欧美中文字幕| 久久天堂一区二区三区四区| 爱豆传媒免费全集在线观看| 国产高清视频在线播放一区 | 男女免费视频国产| 啦啦啦视频在线资源免费观看| 人妻 亚洲 视频| 免费人妻精品一区二区三区视频| 99九九在线精品视频| 赤兔流量卡办理| bbb黄色大片| 狂野欧美激情性bbbbbb| 精品少妇黑人巨大在线播放| 亚洲综合色网址| 亚洲精品国产av蜜桃| 国产精品av久久久久免费| 欧美大码av| 国产日韩欧美在线精品| 女警被强在线播放| 秋霞在线观看毛片| 少妇裸体淫交视频免费看高清 | 久久久国产一区二区| 欧美亚洲日本最大视频资源| 国产女主播在线喷水免费视频网站| 一级毛片我不卡| 老汉色av国产亚洲站长工具| √禁漫天堂资源中文www| 一区二区三区乱码不卡18| 国产av一区二区精品久久| 国产极品粉嫩免费观看在线| 亚洲欧美日韩另类电影网站| 亚洲欧美一区二区三区黑人| 精品高清国产在线一区| 欧美人与善性xxx| av一本久久久久| 国产精品久久久人人做人人爽| 国产一区亚洲一区在线观看| 欧美+亚洲+日韩+国产| 制服诱惑二区| 97精品久久久久久久久久精品| 欧美在线一区亚洲| 午夜福利影视在线免费观看| videosex国产| 一区二区三区四区激情视频| 人妻一区二区av| 中文字幕最新亚洲高清| 欧美黑人精品巨大| 中文精品一卡2卡3卡4更新| 久久久国产欧美日韩av| 日本一区二区免费在线视频| 国产欧美日韩一区二区三区在线| 青草久久国产| 九色亚洲精品在线播放| 在线看a的网站| 国产精品熟女久久久久浪| 你懂的网址亚洲精品在线观看| 好男人视频免费观看在线| 黄色视频不卡| 久久久久久久久免费视频了| 久久毛片免费看一区二区三区| 国产精品免费大片| 日韩 欧美 亚洲 中文字幕| 色婷婷av一区二区三区视频| 国产精品久久久久久精品古装| 中国美女看黄片| 久热爱精品视频在线9| 色精品久久人妻99蜜桃| 国产精品熟女久久久久浪| 狠狠婷婷综合久久久久久88av| 亚洲成人手机| 亚洲国产日韩一区二区| av国产精品久久久久影院| www.999成人在线观看| 亚洲av电影在线观看一区二区三区| 国产午夜精品一二区理论片| av网站在线播放免费| av欧美777| 精品欧美一区二区三区在线| 欧美av亚洲av综合av国产av| 女警被强在线播放| 伊人久久大香线蕉亚洲五| 人人妻人人澡人人看| 狠狠精品人妻久久久久久综合| 日本一区二区免费在线视频| 在线亚洲精品国产二区图片欧美| 国产精品成人在线| 一边摸一边做爽爽视频免费| 日韩大码丰满熟妇| 一二三四在线观看免费中文在| 一级片免费观看大全| 狂野欧美激情性bbbbbb| 视频在线观看一区二区三区| 水蜜桃什么品种好| 国产成人影院久久av| 777久久人妻少妇嫩草av网站| 51午夜福利影视在线观看| 久久精品亚洲熟妇少妇任你| 中文字幕另类日韩欧美亚洲嫩草| 老司机深夜福利视频在线观看 | 电影成人av| 欧美日韩av久久| 大话2 男鬼变身卡| 91九色精品人成在线观看| 色网站视频免费| 国产精品熟女久久久久浪| cao死你这个sao货| av不卡在线播放| 亚洲黑人精品在线| 日韩一区二区三区影片| 欧美日韩视频精品一区| 美女视频免费永久观看网站| 欧美成狂野欧美在线观看| 欧美 日韩 精品 国产| 少妇人妻 视频| 亚洲成人免费av在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 精品一品国产午夜福利视频| 亚洲欧洲精品一区二区精品久久久| 尾随美女入室| 一区二区三区四区激情视频| 一边亲一边摸免费视频| 999精品在线视频| 日韩av不卡免费在线播放| 国产亚洲一区二区精品| 男女边摸边吃奶| 国产精品免费视频内射| 免费人妻精品一区二区三区视频| 欧美日韩亚洲高清精品| av片东京热男人的天堂| 国产淫语在线视频| 搡老乐熟女国产| 亚洲国产日韩一区二区| 亚洲国产欧美网| 欧美xxⅹ黑人| 一级a爱视频在线免费观看| 亚洲av国产av综合av卡| 韩国高清视频一区二区三区| 香蕉国产在线看| 自线自在国产av| 中文字幕色久视频| 80岁老熟妇乱子伦牲交| 人人澡人人妻人| 久久女婷五月综合色啪小说| 2021少妇久久久久久久久久久| 满18在线观看网站| 少妇 在线观看| 亚洲美女黄色视频免费看| 三上悠亚av全集在线观看| 精品国产一区二区久久| 丝袜喷水一区| 免费av中文字幕在线| 欧美黑人欧美精品刺激| 免费一级毛片在线播放高清视频 | 精品一区在线观看国产| 亚洲成av片中文字幕在线观看| 丁香六月欧美| 人体艺术视频欧美日本| 亚洲欧美一区二区三区久久| 一级毛片黄色毛片免费观看视频| 欧美精品啪啪一区二区三区 |