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

    Experimental research in rotating wedge-shaped cooling channel with multiple non-equant holes lateral inlet

    2018-04-19 08:29:00HuiminZHOUHaiwangLIZhiTAOHongwuDENGMeisongYANGShuqingTIAN
    CHINESE JOURNAL OF AERONAUTICS 2018年3期

    Huimin ZHOU,Haiwang LI,Zhi TAO,Hongwu DENG,*,Meisong YANG,Shuqing TIAN

    aNational Key Laboratory of Science and Technology on Aero-Engine Aero-Thermodynamics,Beihang University,Beijing 100083,China

    bCollaborative Innovation Center of Advanced Aero-Engine,Beihang University,Beijing 100083,China

    cResearch and Development Center,AVIC Commercial Aircraft Engine Co.,Ltd.,Shanghai 201108,China

    1.Introduction

    Demands of higher efficiency of modern advanced gas turbines require them operate in a temperature higher than the melting point of the blade,to promote the development of both internal and external cooling techniques to enhance turbine heat transfer.Previous experimental studies on internal cooling technology of advanced gas turbines have been summarized by Han and Huh1–3Internal cooling approaches vary with different parts of the turbine to respond to different needs and achieve the best cooling effect.The trailing edge is one of the most challenging regions due to typically narrow cooling channels.The structural features of the turbine trailing edge can be summarized asfollows:wedge-shaped converging crosssection,pin- fin or other turbulence promoters,lateral coolant ejection and large installation angle.Coriolis force and buoyancy force generated in a rotating state make the flow and heat transfer inside the channel further complicated.

    Under static condition,many studies have been done and confirmed that the flow and heat transfer of turbine trailing edge are mainly influenced by pin- fin cooling configuration and cooling channel structure.Lateral ejection is one of the most effective characteristics.Lau et al.4,5studied a pin- fin cooling channels with lateral ejection structure.The results showed that the decreasing rate of mass flow along the channel caused by the lateral ejection reduced the downstream heat transfer of the channel and furt6her affected the overall heat transfer.Wright and Gohardanifound that the heat transfer of both rectangular channel and wedge-shaped channel is uniform in the cross-section without lateral ejection,and the lateral ejection structure is very effective in increasing the heat transfer coefficient near the out flow holes.Kumaran et al.7confirmed Lau’s conclusion and found that the length of ejection holes had little effect on the channel heat transfer by comparing the heat transfer three different channel configurations pin- fin roughened channel with short lateral ejection holes,with long ones and without lateral ejection holes.McMillin and Lau8studied a pin- fin arrayed rectangular channel with naphthalene sublimation technique and drew the conclusion that the turning of main flow caused by side-wall ejection generated higher mass transfer near the ejection but lower near the opposite inner region.Taslim et al.9studies a wedge-shaped channel with lateral ejection and observed a significant span-wise heat transfer difference due to the trapezoidal converging configuration.It was also found that this span-wise heat transfer variation can be apparently weakened by the side-wall ejection.Hwang et al.10conducted experiments involving eight wedge-shaped channels with different pin- fin and outlet configurations to study the effect of lateral-to-total flow ratio.It was found that the heat transfer was the worst when the ratio approached to 0.3.Hwang and Lu11confirmed that the effect of lateral-tototal was a main factor as well as pin shape after further study via transient liquid crystal technique.The worst heat transfer happened when the ratio approached 0.3–0.4 depending on different pin shape.Kan et al.12did numerical simulation of flow and heat transfer in a high aspect ratio channel with lateral ejection.It was shown that the heat transfer near the end of channel reduced due to a significant flow separation happened there.

    Rotating experimental researches are more challenging than the static ones.Chang et al.13found that the heat transfer in a rib roughened wedge-shaped channel with side-wall ejection holes was lower than that without the side-wall ejection ones,especially in a rotating state.However,the experiments were conducted when the installation angle of channel was 90°rather than the real configuration of a trailing edge cooling channel.Wright et al.14studied the heat transfer characteristics of a rotating smooth wedge-shaped cooling channel with radial out flow under the channel configuration of 135°.The experimental results showed that the rotation number and buoyancy number are ideal parameters for describing the rotational heat transfer characteristics of the channel.The effect of rotation on the leading surface was greater than that on the trailing surface.There was a critical rotation leading to the obvious change of the heat transfer on the leading surface.Liu et al.15did similar investigation but the experiment channel was with lateral ejection.It was found that the lateral ejection weakened the rotation effects and enhanced the overall heat transfer inside the channel for inducing turbulence mixing.The local heat transfer near the out flow slots was especially increased due to a gas film formed near the ejection.Rallabandi et al.16added pin- fins inside the channel which enhanced the heat transfer and weakened the effect of rotation.Liu et al.17also did experiments with turbulators and obtained similar conclusions.Qiu18and Tao19et al.experimentally investigated the effects of channel orientation,outlet boundary condition on flow and heat transfer in a rotating wedge-shaped channel with lateral fluid extraction.The results indicated that the channel orientation was an in fluential factor on channel heat transfer,thus should be chosen carefully to simulate the actual conditions of the turbine.The heat transfer of the channel’s inner top region approached the worst point and was sensitive to the rotation effect under the channel configuration of 90°.Srinivasan et al.20examined Qiu’s experiments18by numerical simulation and confirmed that the heat transfer on the trailing wall was enhanced significantly by rotation.Pardeshi et al.21numerically investigated more complex situation additionally considering pin- fin/rib turbulators.Li et al.22added an additional inlet to wedge-shaped channel at the inner top region and experimental investigated the flow and heat transfer under both rotating and non-rotating conditions.As a result,the second inlet notably improved the heat transfer at the top of the channel and compensated the negative effects induced by the rotation.When the mass flow rate ratio(second inlet mass flow rate/major inlet mass flow rate)came to 0.3,the second inlet benefited overall heat transfer at high rotation numbers.

    From the foregoing,previous studies have confirmed that both lateral outlets and inlets have positive effect on heat transfer of the channel and can weaken the impact of rotation.However,adding a second coolant inlet at the top of the channel can only improve the local heat transfer in the vicinity of the inlet but does not benefit the overall heat transfer enough.From the overall perspective,the negative effect of rotation can lead to local overheating of the channel and low efficiency of coolant utilization.The objective of this paper is to further improve the overall heat transfer of turbine trailing edge cooling channel and decrease adverse impact of rotating.Concretely,the second inlet is expanded from one single hole to multiple lateral holes.The radius of those holes varies depending on their location to enhance channel heat transfer uniformity and improve the utilization efficiency of the cooling air by controlling the mass flow rate through different holes.Experiments have been done to research the effect of this new inlet structure under both static and rotating states and the results are summarized.

    2.Experiment setup

    2.1.Rotating facilities

    Fig.1 Rotating facility.

    As shown in Fig.1,the rotating facility consists of four main modules connected by standard interface,which are electric motor,airing system,rotating arm with support and measurements module.Before being supplied into the rotating facility through a Deublin rotary union,two passes of cooling air cooling air individually pass through a FCI-ST98 thermal flow meter and a GEMU electromagnetic valve.The FCIST98 thermal flow meter can collect the mass flow rate in real time,the measuring range of which is from 0 to 29.5 kg/h and from 0 to 60 kg/h.The GEMU electromagnetic valve cooperating with thermal flow meter achieves automatic control of mass flow rate by regulating valve opening according to the difference between the real-time mass flow rate and set value.The test section is placed at the top of the rotating arm with the maximum rotating diameter of 1.3 m.Cooling air through test section is directly discharged to the atmosphere.In order to ensure a large enough range of the rotation number,a throttle valve is set at the exit of the channel so that the pressure inside the channel is maintained at about 5 times the atmospheric pressure.Temperature measurement is accomplished by thermocouples.The cold ends of the thermocouples are placed in a rotating cavity isolated from hot parts to ensure that the temperature in this cavity could be considered as constant and uniform.The analog signals are converted into digital ones before being transmitted to the non-rotating facilities by slip rings.More details about the rotating facilities can be found in our previous works.22–24

    2.2.Test section

    The test channel is simplified and amplified by five times on the base of actual cooling channel.The dimensions of test model are shown in Fig.2.The channel consists of 230 mm long heated section and a 63 mm long unheated entrance region.The cross-section of the heated channel is a trapezoid with a height of 57.5 mm.The 20 mm wide edge of this wedge narrows to 6 mm and gives the apex an angle of 15.Thus,the Area of the cross section(A)is about 702 mm2and the hydraulic diameter(Dh)is 22.9 mm.The height of the bottom inlet is 2/3 of the heated channel,thus there is a sudden expansion of the entrance,which is similar with actual channel.The distance from the axis of rotation to the major inlet of heated section is 350 mm.The second inlet consists of 8 circular holes located on the side-wall whose diameters vary from 6 mm to 12 mm.Eight slots are drilled on the narrow side-wall,through which the coolant is ejected from the heated channel.The pressure inside the channel is maintained at 5×105Pa by manual regulation through a throttle valve.

    Forty-eight copper plates with a surface area of 22 mm×16 mm are embedded on Nylon made leading and trailing surfaces.A T-type thermocouple is inserted into the back side blind hole of each copper plate to measure the average surface temperature of the copper.One film heater is glued on the back of copper plate to supply the heat energy.In the span-wise direction of the channel,the channel is divided into three parts(inner,mid-span and outer)according to the location of copper plates.And in the stream-wise direction,there are eight locations numbered from P1to P8.Since the resistances of those heaters that connected in series are basically the same,the overall thermal boundary could be considered as iso- flux.Ω in Fig.2 is the rotating angular velocity of the channel.To simulate the real turbine blade trailing edge configuration,the channel is oriented at 135°respect to the rotating plane.1,18The inlet temperature(Tin)is 292 K.The wall temperature ranges from 315 K to 370 K and the average wall temperature is about 335 K.Other details of working conditions are listed in Table 1.

    Fig.2 Dimensions of test channel.

    Table 1 Working conditions of test channel.

    3.Data reduction

    To facilitate the experiment,the mass flow rate of major inlet is confirmed while that of second inlet is adjustable.The Mass flow rate Ratio(MR)is defined as follows based on main inlet to describe the proportion of bottom inletand lateral inlet

    or

    The total inlet mass flow rate ratio is defined as the ratio of lateral flow to total flow so as to research the role of lateral inlet as follows:

    Both entries correspond to a Reynolds number.Reynolds number of working condition is defined with the major inlet mass flow rate for unified comparison as Eqs.(4),(5):

    where ρ is the density of the fluid,μ is the viscosity,and Vinis the velocity of the inlet flow.In this paper,Re refers to Remainwithout special instructions.

    The effective Reynolds number is defined with total inlet mass flow rates to describe conditions while different inlet is major inlet as Eq.(6):

    where Vin,sumis the equivalent velocity of coolant from both inlets.

    Rotation number is defined as follows to represent the magnitude of rotation:

    where n is the rotational speed.

    Nusselt number is defined as follows:

    where λ is the thermal conductivity of the coolant,h is heat transfer coefficient.

    The regionally averaged heat transfer coefficient depends on the net heat flux of the inner wall of the channel(Qnet),heat transfer surface area,regional average wall temperature(Tw)and the overall average bulk temperature(Tb).And the specific definition is shown in Eq.(9):

    The net heat flux can be calculated as Eq.(10)with total electrical energy and the heat loss(Qloss),which can be estimated by heat loss coefficient(α)and the difference between wall temperature and environment temperature(Te):

    where I is the electric current,and Riis the resistance of each temperature measuring point.

    The heat loss coefficient is a function of angular velocity which varies with the change of dimensionless position and can be calibrated by heat loss experiments as follows:

    The local bulk temperature(Tb,x)is difficult to measure or calculate directly because it changes greatly under different operating conditions.So in this experiment,it is described by overall averaged bulk temperature as follows:

    where Toutis the temperature at the outlet,and cpis the heat capacity of the coolant.

    As coolant from two inlets is at the same temperature,Tbcan be further simplified as follows:

    Thermal conductivity is defined as follows:

    where λ0is the thermal conductivity of air at standard atmospheric pressure,and the qualitative temperature Tf,xis selected as Tb,xin this paper.

    The regional buoyancy number is defined as follow:

    in which TR represents the local temperature ratio of the channel and RLrepresents the radius of gyration.The temperature ratio is fixed to 0.14 as the real turbine working conditions on the basis of previous studies.The temperature ratio has an effect on heat transfer characteristics by changing the physical properties and the buoyancy force in the channel.At low rotational speeds,the weakening effect on the heat transfer due to change of the physical property is greater than the enhancement of that due to the increase of buoyancy force.But this effect is not obvious in the actual experimental process as shown in Fig.3,where the ordinate Nu/Nusrepresents the ratio of span-wise average Nusselt number under rotating conditions and that under static conditions.With the increase of the rotation number,the effect of buoyancy increase on heat transfer becomes more and more obvious and the two effects of the temperature ratio cancel each other out.When the rotation number further increases and the influence of buoyancy becomes dominant,especially at the top of the channel,the heat transfer is enhanced with the temperature ratio increasing but the maximum of heat transfer differences between conditions where TR=0.14 and TR=0.08 is less than 8%.So it can be considered that the change of temperature ratio is negligible in this paper.For another,the RL/Dhis a constant in this experiment.So buoyancy number has the same effect as the rotation number in this paper.

    The measured Nusselt numbers are customarily compared with Dittus-Boelter/McAdams correlation or data under static conditions to normalize the data as Eq.(16).In this study,the Prandtl number(Pr)is a constant of 0.4.In the following discussion,Nu0is the average Nusselt number of fully developed turbulence in a static smooth pipe with the same hydraulic diameter as the channel,and Nu0is calculated with the effective Reynolds number without special instructions.It has been proved in previous studies that the Nu/Nu0is a good way to normalize the Nusselt number.22

    Fig.3 Effect of local temperature ratio on heat transfer in a channel with bottom inlet under rotating conditions(Re=20000,MR=0).

    Fanning friction factor(f)is employed as the flow resistance coefficient to describe the frictional losses within the channel.Thus,in Eq.(17),pinand poutrespectively represent static pressure at the inlet and outlet,is the dynamic head at the radial inlet channel cross-section, Δpe,cenis a correction term for pressure changes caused by centrifugal force,and L is the length of the channel.In this experiment,the pressure difference between inlet and outlet is not influenced by centrifugal force because the pressure measuring points are on the same rotating radius.

    To comprehensively analyze heat transfer and flow resistance caused by different channel structure and other factors,the thermal performance(η)is described as Eq.(18)based on the constant pumping power condition25:

    The resistance coefficient and the thermal performance are normalized in the same way as Nusselt numbers.Nus,ηsand fsare data under static conditions;Nu0,η0and f0are data of fully developed turbulence in a static smooth pipe with the same hydraulic diameter as the channel.

    Actually the experiment model and real blade are incompletely similar,which is subject to the rotation radius and rotating speed under laboratory conditions.Two main problems can be caused by incomplete similarity.First,it is difficult to match critical criterion numbers;Second,the small rotation radius can result in a small buoyancy number.As this experiment is conducted at a fixed temperature ratio,the correction of buoyancy number is not discussed in this essay.To match critical criterion numbers,the pressure inside the experiment channel is maintained at about 5 times the atmospheric pressure,which decreases the coolant velocity and increases the coolant density at the inlet.Thus,the rotation number is raised without changing the Reynolds number and similar experimental studies can be carried out.26,27The experimental error is given in Appendix A.

    4.Results and discussion

    4.1.Non-rotating channel with single inlet

    Fig.4 Distribution of Nusselt number ratios at non-rotating conditions with either bottom inlet or lateral inlet(TR=0.14,MR=0,Ro=0).

    The Nu/Nu0distributions along the channel with only bottom inlet or lateral inlet is shown in Fig.4.The results show that under non-rotating conditions,the major inlet Reynolds number barely affect the heat transfer distribution curve,which allows easier analysis of the data under different Reynolds number conditions and wider application of the experimental results.The heat transfer difference between the leading and trailing wall is inconspicuous because of geometric symmetry.Under conditions with only bottom inlet,the heat transfer capacity significantly decreases in stream-wise direction.This phenomenon can be attributed to the coolant continuously spilling from lateral ejection,so the cooling air at the top of the channel is much less than that at the bottom of the channel.In span-wise direction,except P1,the heat transfer at the outer region is higher than that at the mid-span and inner area,for the air flow is accelerated when turning nearly 90°and being ejected through side-wall slots.However,at P1,a flow dead zone formed due to the sudden expansion of the entrance as shown in Fig.5,which slightly weakens the heat transfer near the outer surface.As a comparison,lateral inlet apparently improves the heat transfer at high radius area of the channel.In the ideal case,the coolant should be injected in exact span-wise direction and make the heat transfer uniform throughout the whole channel.In the actual experiment,however,the air flow incidence angle is not exact 90°especially in the vicinity of the root of the channel as the cooling air is injected radially initially and sharply turns to the axial direction,which gives rise to a flow dead zone and reduces the heat transfer near the bottom inlet as shown in Fig.6.When the airflow reaches the middle of the channel,the injection direction is gradually close to axial and the impact effect is strengthened,which enhanced the heat transfer around the middle.On the contrary,the heat transfer at the top of the channel slightly decreases.The reason can be inferred that less cooling air is injected from the top holes.In span-wise direction,the midspan region of the channel present better heat transfer characteristics than the outer and inner areas for being directly impinged by the coolant as shown in Fig.6.Compared with conditions with only bottom inlet,the average heat transfer in the channel with lateral inlet is respectively enhanced by 30.01%and 28.87%when the Reynolds number is 10000 and 20000.And the heat transfer at the high-radius half channel(P5–P8)respectively increases by 168.2%and 128.3%with the Reynolds number 10000 and 20000.Partial data at P8are missing due to thermocouple damage but it has little impact on interpretation of the trend of the whole curve.

    Fig.5 Flow distribution at non-rotating conditions with bottom inlet.

    Fig.6 Flow distribution at non-rotating conditions with lateral inlet.

    4.2.Non-rotating channel with two inlets

    Fig.7 Distribution of the Nusselt number ratios at non-rotating conditions with two inlets(Re=12000,TR=0.14,Ro=0).

    Fig.7 shows the Nu/Nu0distribution in the non-rotating channel with two inlets varying with MR.Reynolds number of the main inlet is kept at 12000 and the MR changes from 0 to 1.0.Considering the good symmetry of the heat transfer on the leading wall and the trailing one,data at corresponding points on both sides are averaged.In Li’s study22,the heat transfer on the top of the channel grows with the ratio of the second coolant injection to total flow.But a decrease in bottom inlet coolant leads to the decrease in Nu/Nu0at the first few locations of the channel.Similarly,in this experiment the lateral inlet expands the influence of the second injection.In short,bottom injection is good for the first few locations while lateral injection benefits the heat transfer from the middle to the top of the channel.It can be inferred that there is a proper MR that is relatively best for the overall heat transfer.

    In order to further study the influence of the mass flow rate ratio on the local heat transfer of the channel,48 temperature measurement points were divided into four areas and numbered 1–4.Averaged data of these four areas are organized by the total inlet mass flow rate ratio(MR′)and the result is shown in Fig.8.It is shown that before MR′reaches 0.2,the heat transfer in Area 1 is almost the same,and that in other 3 areas decreases due to the reduction of the proportion of bottom inlet coolant.When MR′exceeds 0.2,the heat transfer in Area 1 decreases while heat transfer in other 3 areas is enhanced by lateral injection.The standard deviations of Nu/Nu0in these 4 areas are organized by MR′and shown in Fig.9.When MR′approaches 0.56,i.e.when the coolant from two inlets is roughly equal,the heat transfer inside the channel is the most uniform.Considering both intensity and uniformity of the heat transfer,it can be concluded that the lateral mass flow rate of about half of the total mass flow rate is the best option under non-rotating conditions.

    Fig.9 Effect of mass flow rate ratio on standard deviation of averaged Nusselt number ratios at non-rotating conditions with two inlets(Re=12000,TR=0.14,Ro=0).

    Fig.8 Effect of mass flow rate ratio on Nusselt number ratios at non-rotating conditions with two inlets(Re=12000,TR=0.14,Ro=0).

    4.3.Rotating channel with single inlet

    Figs.10 and 11 present the Nu/Nusdistribution of each measuring point with rotation number at different Reynolds numbers(Re=10000 and Re=20000)in rotating channel with either bottom inlet or lateral inlet.The curves at different Reynolds numbers are quite similar but the rotation effect reduces under higher Reynolds number at the top of the channel with bottom inlet so the curve fluctuation amplitude decreases in Fig.10(b),which may be attributed to the increase of coolant.In the channel with bottom inlet,rotation enhances the heat transfer on the trailing wall and weakens that on the leading wall.This phenomenon is especially obvious in the vicinity of the sudden expansion entrance and the top of the channel.Lateral injection reduces negative effect of rotation to a certain extent.In comparison,Nu/Nusof channel with lateral inlet is almost independent with the rotation number and the heat transfer disparity between two walls is apparently reduced.Two counter-current zones are formed at the top of the channel(P8)and the inner region in low radius area(P2and P3),where the average air flow velocity is from the top of the channel to the bottom and the heat transfer on the leading wall is enhanced.

    Fig.12 shows the Nu/Nu0distribution along a rotating channel with lateral inlet at Reynolds number of 10000.Due to the lateral inlet,rotation has hardly affected the distribution except enhancing the heat transfer near P2and P3on the leading wall.The heat transfer on the leading wall and the trailing wall still has a good symmetry and the overall heat transfer is uniform.Similar statement is also applicable to the condition at Reynolds number of 20000(Fig.13).

    4.4.Rotating channel with two inlets

    For a clearer study of the effect of the mass flow rate ratio on the heat transfer characteristics along the channel,the parameters are averaged along the span direction as follows:

    Fig.10 Effect of rotation on Nu/Nusin channel with bottom inlet(TR=0.14,MR=0).

    Fig.11 Effect of rotation on Nu/Nusin channel with lateral inlet(TR=0.14,MR=0).

    Fig.12 Effect of rotation on distribution of Nusselt number ratios in channel with lateral inlet(TR=0.14,Re=10000,MR=0).

    Fig.13 Effect of rotation on distribution of Nusselt number ratios in channel with lateral inlet(TR=0.14,Re=20000,MR=0).

    Fig.14 Effect of rotation on distribution of Nusselt number ratios in channel with two inlets(TR=0.14,Re=12000).

    where the subscript‘T” is for trailing surface,the subscript‘L”is for leading surface,andthe subscript‘outer”,‘mid” and ‘inner”represent different area of the channel as shown in Fig.2.

    It has been discussed that when the total mass flow rate ratio is about 0.5,it is proper for both intensity and uniformity of heat transfer under non-rotating conditions.Experiments have been done to acquire the appropriate MR′considering the effect of rotation and the Nui,ave/Nu0distribution is shown in Fig.14.When MR′is less than 0.5,i.e.coolant from bottom inlet is the main flow,the effect of lateral cooling air on the high radius of the channel is more obvious.When the proportion of lateral coolant is over 50%,lateral air flow’s sphere of influence is expanded to the total channel.When keeping increasing the lateral coolant,heat transfer at the middle and top of the channel is further enhanced with that near the bottom inlet decreasing.

    Fig.15 shows the distribution of normalized resistance coefficient and the thermal performance in a two-inlet channel under the rotating condition.With the increase of the lateral injection,the flow resistance in the channel increases firstly and then decreases obviously.The effect of rotation on the parameters is not significant.Comparing with the condition where MR′=0.5,the f/f0decreases by 8.57%and the η/ηsincreases by 5.34%.

    Considering the overall heat transfer and flow resistance of the channel,the total mass flow rate ratio of 0.67 is an appropriate option.

    Fig.15 Effect of rotation on normalized resistance coefficient and synthesis coefficient in channel with two inlets(TR=0.14,Re=12000).

    5.Conclusions

    The flow and heat transfer in a smooth wedge-shaped channel with two inlets(bottom inlet and lateral inlet formed by multiple non-equant holes)and lateral coolant ejection is studied.48 temperature measurement points distributed on the inner surface of the channel and other measuring equipment provide the regionally averaged heat transfer coefficients.All the experiments are conducted at the inlet Reynolds numbers from 10000 to 20000 and rotation number from 0 to 1.16.The total mass flow rate ratio(lateral mass flow rate/total mass flow rate)varies from 0 to 1.0.The temperature ratio is fixed to 0.14 as the real working conditions of turbine.Conclusions based on experimental facts can be summarized as follows.

    (1)In non-rotating conditions,lateral injection improves the heat transfer at the middle and the top of the channel but reduces at the low-radius areas.With lateral inlet,the average heat transfer in the channel with lateral inlet is respectively enhanced by 30.01%and 28.87%while the heat transferat the high-radius half channel(P5–P8)respectively increases by 168.2%and 128.3%with the Reynolds number 10000 and 20000 compared with conditions with only bottom inlet.

    (2)In non-rotating conditions,when the lateral flow rate is about half of the total flow rate,the intensity and uniformity of overall heat transfer is satisfying.

    (3)In rotating channel,the negative effects induced by rotation are significantly compensated by lateral injection and the heat transfer disparity on the trailing side and leading side of the channel is reduced to an acceptable range.The heat transfer at the middle and top of the channel is notably enhanced with that at the bottom of the channel reducing.The heat transfer distribution is almost free from the impact of rotation.

    (4)Proper proportion of lateral coolant also helps decrease the flow resistance inside the channel.The flow and heat transfer performance is the best when the total flow rate ratio is about 0.67.

    Sincere gratitude goes to AVIC Commercial Aircraft Engine Co.,Ltd.for the help in the research process,including providing information and helping confirm experimental project.

    Appendix A

    The main control variables used in this paper are the Reynolds number and the number of rotations while the main variable is the Nusselt number.The direct measurement quantities in the experiment mainly includes temperature(T),pressure(p),rotating speed,current,mass flow rate(m˙)and geometric dimensions.The direct measurement errors caused by the measurement tools and methods in this experiment are shown in Table A1.In Table A1,M represents the measuring range of flowmeter.

    The errors of the direct measurement quantities further cause the errors of the indirect measurement quantities.Suppose the indirect measurement quantity Y is the general function of the direct measurement quantities X1,X2,...,Xnas Eq.(A1).

    The relative error of the indirect measurement is as follows:

    As mentioned in Table A1,the absolute error of the mass flow rate is relevant to the measuring range of flowmeter.Experiment flowmeter range is 30 kg/h or 60 kg/h,and the relative error of mass flow rate under each condition is shown in Table A2.

    The relative error of the Reynolds number is defined as follows according to Eq.(A3)and the detailed results are shown in Table A3:

    Table A1 Errors of direct measurement quantities.

    Table A2 Relative errors of mass flow rate.

    Table A3 Relative errors of Reynolds number.

    The relative error of the rotation number is described as follows and the detailed results are shown in Table A4:

    The error of the local Nusselt number is as follows:

    Thermal conductivity of air is a function of temperature,

    Then the error of air thermal conductivity can be calculated as Eq.(A7):

    Due to the large number of working conditions and measuring points,the detailed local Nusselt number error is not list here.

    The relative errors of fanning friction factor are calculated as Eqs.(A8)-(A11)and list in Table A5.

    Table A4 Relative errors of rotation number(Remain=12000).

    Table A5 Relative error of lf ow resistance coefficient.

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

    2.Han JC.Turbine blade cooling studies at Texas A&M University:1980–2004.J Thermophys Heat Transfer 2006;20(2):161–87.

    3.Han JC,Huh M.Recent studies in turbine blade internal cooling.Heat Transfer Res 2010;41(8):803–28.

    4.Lau SC,Han JC,Batten T.Heat transfer,pressure drop,and mass flow rate in pin fin channels with long and short trailing edge ejection holes.J Turbomach 1989;111(2):116–23.

    5.Lau SC,Mcmillin RD,Kukreja RT.Segmental heat transfer in a pin fin channel with ejection holes.Int J Heat Mass Transfer 1992;35(6):1407–17.

    6.Wright LM,Gohardani AS.Effect of coolant ejection in rectangular and trapezoidal trailing edge cooling passages.New York:ASME;2008.Report No.:GT2008-50414.

    7.Kumaran TK,Han JC,Lau SC.Augmented heat transfer in a pin if n channel with short or long ejection holes.Int J Heat Mass Transfer 1991;34(10):2617–28.

    8.McMillin RD,Lau SC.Effect of trailing-edge ejection on local heat(mass)transfer in pin fin cooling channels in turbine blades.J Turbomach 1994;116(1):159–68.

    9.Taslim ME,Li T,Spring SD.Experimental study of the effects of bleed holes on heat transfer and pressure drop in trapezoidal passageswith tapered turbulators.J Turbomach 1995;117(2):281–9.

    10.Hwang JJ,Lai DY,Tsia YP.Heat transfer and pressure drop in pin- fin trapezoidal ducts.J Turbomach 1999;121(2):264–71.

    11.Hwang JJ,Lu CC.Lateral- flow effect on end wall heat transfer and pressure drop in a pin- fin trapezoidal duct of various pin shapes.New York:ASME;2000.Report No.:2000-GT-0232.

    12.Kan R,Li Y,Ren J,Jiang HD.Effect of rib configuration and lateral ejection on a high aspect ratio trailing edge channel.New York:ASME;2013.Report No.:GT2013-94549.

    13.Chang SW,Liou TM,Chiou SF,Chang SF.Heat transfer in highspeed rotating trapezoidal duct with rib-roughened surfaces and air bleeds from the wall on the apical side.J Heat Transfer 2008;130(6):142–50.

    14.Wright LM,Liu YH,Han JC,Chopra S.Heat transfer in trailing edge,wedge-shaped cooling channels under high rotation numbers.J Heat Transfer 2008;130(7):135–45.

    15.Liu YH,Huh M,Wright LM,Han JC.Heat transfer in trailingedge channels with slot ejection under high rotation numbers.J Thermophys Heat Transfer 2015;23(2):305–15.

    16.Rallabandi AP,Liu YH,Han JC.Heat transfer in trailing edge wedge-shaped pin- fin channels with slot ejection under high rotation numbers.J Therm Sci Eng Appl 2011;3(2):021007.

    17.Liu YH,Huh M,Han JC.High rotation number effect on heat transfer in a trailing edge channel with tapered ribs.Int J Heat Fluid Flow 2011;33(1):182–92.

    18.Qiu L,Deng H,Tao Z.Effect of channel orientation in a rotating smooth wedge-shaped cooling channel with lateral ejection.New York:ASME;2013.Report No.:GT2013-94758.

    19.Tao Z,Qiu L,Deng H.Heat transfer in a rotating smooth wedgeshaped channel with lateral fluid extraction.Appl Therm Eng 2015;87:47–55.

    20.Srinivasan B,Dhamarla A,Jayamurugan C,Rajan AB.Numerical Studies on effect of channel orientation in a rotating smooth wedge-shaped cooling channel.New York:ASME;2014.Report No.:GT2014-26560.

    21.Pardeshi I,Shih TI,Bryden KM,Ames R,Dennis RA,Ding S,et al.Flow and heat transfer in a rotating and non-rotating wedgeshaped cooling passage with ribs and pin fins.Reston:AIAA;2015.Report No.:AIAA-2015-1444.

    22.Li Y,Xu G,Deng H,Qiu L,Yu X.Effects of coolant mass flow rate ratio on heat transfer in a two-inlet rotating wedge-shaped channel.Int J Heat Mass Transfer 2016;96:353–61.

    23.Deng H,Qiu L,Tao Z,Tian S.Heat transfer study in rotating smooth square U-duct at high rotation numbers.Int J Heat Mass Transfer 2013;66:733–44.

    24.Qiu L,Deng H,Sun J,Tao Z,Tian S.Pressure drop and heat transfer in rotating smooth square U-duct under high rotation numbers.Int J Heat Mass Transfer 2013;66:543–52.

    25.Fu WL,Wright LM,Han JC.Rotational buoyancy effects on heat transfer in five different aspect-ratio rectangular channels with smooth walls and 45 degree ribbed walls.J Heat Transfer 2006;128(11):1130–41.

    26.Cheng J,Ni B,Deng H.Experimental technique and its validation for heat transfer in cooling channels at high rotation numbers.J Aerospace Power 2014;29(8):1817–23[Chinese].

    27.Deng H,Cheng HJ,Chen H,Tian S.Heat transfer in two-pass smooth square channel under large rotation numbers.J Beijing Univ Aeronaut Astronaut 2014;40(5):701–6[Chinese].

    kizo精华| 18在线观看网站| 久久精品人人爽人人爽视色| 久久天躁狠狠躁夜夜2o2o| 色婷婷av一区二区三区视频| √禁漫天堂资源中文www| 国产亚洲午夜精品一区二区久久| 亚洲欧美一区二区三区久久| 精品人妻1区二区| www.自偷自拍.com| 精品少妇黑人巨大在线播放| 女性生殖器流出的白浆| 精品福利永久在线观看| 欧美另类亚洲清纯唯美| 精品少妇久久久久久888优播| 一级a爱视频在线免费观看| 夜夜骑夜夜射夜夜干| 亚洲专区中文字幕在线| 国产午夜精品久久久久久| 老鸭窝网址在线观看| 久久久精品国产亚洲av高清涩受| 精品国内亚洲2022精品成人 | 变态另类成人亚洲欧美熟女 | 免费久久久久久久精品成人欧美视频| 五月天丁香电影| 国产深夜福利视频在线观看| 免费高清在线观看日韩| 日本wwww免费看| 久久国产精品人妻蜜桃| 美女福利国产在线| 国产深夜福利视频在线观看| 大香蕉久久网| 国产av又大| 午夜久久久在线观看| 久久精品亚洲精品国产色婷小说| 国产日韩欧美亚洲二区| 欧美激情高清一区二区三区| 可以免费在线观看a视频的电影网站| 国产免费视频播放在线视频| 999精品在线视频| 他把我摸到了高潮在线观看 | 午夜两性在线视频| 国产精品香港三级国产av潘金莲| 电影成人av| 日韩中文字幕视频在线看片| 国产精品国产高清国产av | 大陆偷拍与自拍| 欧美人与性动交α欧美精品济南到| 欧美精品av麻豆av| 丰满饥渴人妻一区二区三| 国产高清激情床上av| 亚洲七黄色美女视频| 日韩一卡2卡3卡4卡2021年| 50天的宝宝边吃奶边哭怎么回事| 少妇猛男粗大的猛烈进出视频| 成人永久免费在线观看视频 | 色精品久久人妻99蜜桃| 正在播放国产对白刺激| 麻豆乱淫一区二区| 交换朋友夫妻互换小说| 精品福利永久在线观看| 桃花免费在线播放| 国产精品久久久久成人av| av线在线观看网站| 国内毛片毛片毛片毛片毛片| 操美女的视频在线观看| 女人爽到高潮嗷嗷叫在线视频| 乱人伦中国视频| 性高湖久久久久久久久免费观看| 夜夜夜夜夜久久久久| 少妇猛男粗大的猛烈进出视频| 三上悠亚av全集在线观看| 99国产精品免费福利视频| 青草久久国产| 精品国产乱子伦一区二区三区| 麻豆乱淫一区二区| 中文字幕人妻丝袜一区二区| bbb黄色大片| 日日摸夜夜添夜夜添小说| 超碰97精品在线观看| 蜜桃在线观看..| 男女边摸边吃奶| 啪啪无遮挡十八禁网站| 汤姆久久久久久久影院中文字幕| 黄色怎么调成土黄色| 另类精品久久| 露出奶头的视频| 午夜老司机福利片| 不卡av一区二区三区| 波多野结衣av一区二区av| 俄罗斯特黄特色一大片| 成年女人毛片免费观看观看9 | 美国免费a级毛片| 好男人电影高清在线观看| 亚洲精品国产精品久久久不卡| 一区在线观看完整版| 精品一区二区三区av网在线观看 | 少妇裸体淫交视频免费看高清 | 国产精品二区激情视频| 欧美日韩一级在线毛片| 精品国产国语对白av| 精品人妻1区二区| 国内毛片毛片毛片毛片毛片| 成年女人毛片免费观看观看9 | 侵犯人妻中文字幕一二三四区| 欧美黑人精品巨大| 欧美精品一区二区免费开放| xxxhd国产人妻xxx| 三上悠亚av全集在线观看| 国产一区二区激情短视频| 视频区图区小说| 国产成人精品久久二区二区免费| 色在线成人网| av网站免费在线观看视频| √禁漫天堂资源中文www| 亚洲色图av天堂| 久久精品国产亚洲av香蕉五月 | 丰满人妻熟妇乱又伦精品不卡| 亚洲精品一卡2卡三卡4卡5卡| 欧美午夜高清在线| 在线观看免费午夜福利视频| 夜夜骑夜夜射夜夜干| 午夜福利视频在线观看免费| 69av精品久久久久久 | 亚洲熟妇熟女久久| 精品一区二区三区四区五区乱码| 99国产精品一区二区三区| 巨乳人妻的诱惑在线观看| 交换朋友夫妻互换小说| 久久ye,这里只有精品| 麻豆乱淫一区二区| 色94色欧美一区二区| 99久久国产精品久久久| 三上悠亚av全集在线观看| 十八禁网站免费在线| 日本一区二区免费在线视频| 一级毛片女人18水好多| 十八禁高潮呻吟视频| 天天添夜夜摸| 国产成人免费观看mmmm| 亚洲午夜精品一区,二区,三区| 亚洲欧美色中文字幕在线| 久久人妻福利社区极品人妻图片| 日韩大片免费观看网站| 下体分泌物呈黄色| 黄片小视频在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 免费在线观看视频国产中文字幕亚洲| 日韩欧美一区二区三区在线观看 | 国产麻豆69| 欧美乱码精品一区二区三区| 男女无遮挡免费网站观看| 一边摸一边抽搐一进一小说 | 蜜桃国产av成人99| 俄罗斯特黄特色一大片| 久久精品国产99精品国产亚洲性色 | 777久久人妻少妇嫩草av网站| 久久毛片免费看一区二区三区| 国产黄色免费在线视频| 亚洲成人手机| 精品一区二区三卡| 久久午夜亚洲精品久久| 国产成+人综合+亚洲专区| 国产又色又爽无遮挡免费看| 制服诱惑二区| 国产亚洲精品一区二区www | 国产精品久久久久久人妻精品电影 | 国产亚洲一区二区精品| videosex国产| 国产精品 欧美亚洲| 丰满饥渴人妻一区二区三| 日本黄色日本黄色录像| 国产区一区二久久| 动漫黄色视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 91大片在线观看| 中文字幕制服av| 国产有黄有色有爽视频| 在线播放国产精品三级| 亚洲精品国产区一区二| 日韩欧美免费精品| 亚洲欧洲日产国产| 亚洲中文字幕日韩| 亚洲av电影在线进入| 日韩中文字幕视频在线看片| 热re99久久精品国产66热6| 亚洲精品粉嫩美女一区| 久久久国产一区二区| 九色亚洲精品在线播放| 久久精品人人爽人人爽视色| 多毛熟女@视频| 妹子高潮喷水视频| 一夜夜www| www.熟女人妻精品国产| 在线观看免费视频日本深夜| 国产免费现黄频在线看| 男女高潮啪啪啪动态图| 亚洲国产成人一精品久久久| 亚洲av欧美aⅴ国产| 精品少妇久久久久久888优播| 午夜免费成人在线视频| 免费看十八禁软件| 一本久久精品| 天堂动漫精品| 777米奇影视久久| 天堂俺去俺来也www色官网| 五月开心婷婷网| 婷婷丁香在线五月| 国产无遮挡羞羞视频在线观看| 色94色欧美一区二区| 精品欧美一区二区三区在线| 999久久久精品免费观看国产| 日韩视频一区二区在线观看| 搡老乐熟女国产| 日韩制服丝袜自拍偷拍| 777米奇影视久久| 精品一区二区三区视频在线观看免费 | 美国免费a级毛片| 他把我摸到了高潮在线观看 | 一级片免费观看大全| 色视频在线一区二区三区| 国产男靠女视频免费网站| 69av精品久久久久久 | 国产伦人伦偷精品视频| 脱女人内裤的视频| 免费黄频网站在线观看国产| 免费观看a级毛片全部| 国产成+人综合+亚洲专区| 久久毛片免费看一区二区三区| 午夜老司机福利片| 下体分泌物呈黄色| 视频区欧美日本亚洲| 国产亚洲午夜精品一区二区久久| 国产日韩欧美亚洲二区| 午夜福利在线观看吧| 中文字幕av电影在线播放| 日韩 欧美 亚洲 中文字幕| 麻豆国产av国片精品| 成人永久免费在线观看视频 | 亚洲欧美一区二区三区久久| www.999成人在线观看| 日韩一卡2卡3卡4卡2021年| 一级a爱视频在线免费观看| 欧美激情极品国产一区二区三区| svipshipincom国产片| 成人国产一区最新在线观看| 久久99一区二区三区| 青青草视频在线视频观看| 少妇精品久久久久久久| 欧美成狂野欧美在线观看| 黑人操中国人逼视频| 精品人妻在线不人妻| 精品国产乱码久久久久久小说| 久久久久网色| 丝袜美腿诱惑在线| 看免费av毛片| 欧美变态另类bdsm刘玥| 夜夜夜夜夜久久久久| 一区二区三区激情视频| 国产福利在线免费观看视频| 啦啦啦在线免费观看视频4| avwww免费| 欧美日韩福利视频一区二区| 侵犯人妻中文字幕一二三四区| 久久精品人人爽人人爽视色| 50天的宝宝边吃奶边哭怎么回事| 一区福利在线观看| 国产区一区二久久| 自线自在国产av| 丝袜人妻中文字幕| 国产精品国产高清国产av | 日本wwww免费看| 中文亚洲av片在线观看爽 | 最新在线观看一区二区三区| 在线观看www视频免费| 18禁美女被吸乳视频| 怎么达到女性高潮| 99re在线观看精品视频| 免费不卡黄色视频| 午夜视频精品福利| 超色免费av| 日韩免费高清中文字幕av| 精品乱码久久久久久99久播| 国产成人啪精品午夜网站| 久久人人97超碰香蕉20202| 大型av网站在线播放| 亚洲精品av麻豆狂野| 国产免费av片在线观看野外av| 国产国语露脸激情在线看| 国产视频一区二区在线看| 欧美日韩亚洲综合一区二区三区_| 丰满迷人的少妇在线观看| 亚洲欧美色中文字幕在线| 国产色视频综合| 亚洲va日本ⅴa欧美va伊人久久| 国产午夜精品久久久久久| 午夜成年电影在线免费观看| 国产男靠女视频免费网站| 青草久久国产| 一级a爱视频在线免费观看| 国产男女超爽视频在线观看| 一级毛片女人18水好多| a级片在线免费高清观看视频| 少妇的丰满在线观看| 久久久久久人人人人人| 一区二区三区精品91| 成年人午夜在线观看视频| 成人免费观看视频高清| 国产av国产精品国产| 他把我摸到了高潮在线观看 | 亚洲视频免费观看视频| 日韩中文字幕视频在线看片| 精品人妻1区二区| 欧美精品一区二区免费开放| 黄网站色视频无遮挡免费观看| 国产成人精品久久二区二区91| 亚洲av片天天在线观看| 岛国毛片在线播放| 最新美女视频免费是黄的| 免费少妇av软件| 久久免费观看电影| 正在播放国产对白刺激| av在线播放免费不卡| 性少妇av在线| 亚洲国产欧美日韩在线播放| 久久免费观看电影| 亚洲色图av天堂| 中文字幕人妻丝袜制服| 午夜福利视频精品| av有码第一页| 日本av手机在线免费观看| 久久亚洲精品不卡| 一级片'在线观看视频| 亚洲专区中文字幕在线| 亚洲人成电影观看| 亚洲一区中文字幕在线| 黄频高清免费视频| 午夜福利影视在线免费观看| 欧美日韩中文字幕国产精品一区二区三区 | av不卡在线播放| 成人免费观看视频高清| 窝窝影院91人妻| h视频一区二区三区| 岛国在线观看网站| 久久久国产精品麻豆| 久久影院123| 国产又色又爽无遮挡免费看| 精品久久蜜臀av无| 天天躁狠狠躁夜夜躁狠狠躁| 精品久久久久久电影网| 亚洲精品在线观看二区| 女警被强在线播放| 午夜精品国产一区二区电影| 成年人黄色毛片网站| 超碰成人久久| 极品教师在线免费播放| 亚洲av片天天在线观看| 日韩视频一区二区在线观看| 欧美日韩国产mv在线观看视频| 五月开心婷婷网| 久久久久久久久免费视频了| 亚洲,欧美精品.| 国产男女内射视频| 亚洲第一av免费看| 老司机福利观看| 精品亚洲成a人片在线观看| 人人妻人人澡人人看| 久久久久久久大尺度免费视频| 久久久久国内视频| 亚洲午夜精品一区,二区,三区| 亚洲av第一区精品v没综合| 国产主播在线观看一区二区| 亚洲avbb在线观看| 精品国产一区二区三区四区第35| 2018国产大陆天天弄谢| 亚洲国产看品久久| 久久久国产一区二区| 免费黄频网站在线观看国产| 成年女人毛片免费观看观看9 | 天天影视国产精品| av网站在线播放免费| 无遮挡黄片免费观看| 欧美日韩福利视频一区二区| 飞空精品影院首页| 麻豆国产av国片精品| 亚洲精品国产区一区二| 香蕉国产在线看| 国产精品一区二区精品视频观看| 久久人妻福利社区极品人妻图片| 色综合婷婷激情| 一二三四社区在线视频社区8| 国产精品一区二区在线不卡| 国产不卡av网站在线观看| 国产激情久久老熟女| 天天躁狠狠躁夜夜躁狠狠躁| 国产单亲对白刺激| 国产在线观看jvid| 老司机福利观看| 美女午夜性视频免费| 日韩大片免费观看网站| 国产又爽黄色视频| 亚洲精品国产一区二区精华液| 国产深夜福利视频在线观看| 91精品国产国语对白视频| 他把我摸到了高潮在线观看 | 老司机午夜福利在线观看视频 | 精品国产一区二区久久| 极品少妇高潮喷水抽搐| av网站在线播放免费| 99热国产这里只有精品6| 婷婷成人精品国产| 一进一出抽搐动态| 精品人妻1区二区| 亚洲成国产人片在线观看| 国产精品久久久av美女十八| 免费久久久久久久精品成人欧美视频| 国产91精品成人一区二区三区 | 国产av精品麻豆| 一区二区三区乱码不卡18| 精品国产乱码久久久久久男人| 青青草视频在线视频观看| 欧美黑人欧美精品刺激| 男男h啪啪无遮挡| 午夜福利免费观看在线| 美国免费a级毛片| 在线看a的网站| a级毛片黄视频| 国产黄色免费在线视频| 少妇粗大呻吟视频| 亚洲国产av影院在线观看| 免费观看a级毛片全部| 成人18禁在线播放| 中文字幕制服av| 麻豆av在线久日| 9191精品国产免费久久| 欧美黑人精品巨大| 国产欧美日韩综合在线一区二区| 亚洲专区字幕在线| 他把我摸到了高潮在线观看 | 日本黄色视频三级网站网址 | 亚洲精华国产精华精| 香蕉久久夜色| 在线观看免费视频网站a站| 国产精品美女特级片免费视频播放器 | 麻豆国产av国片精品| 亚洲 欧美一区二区三区| 91成年电影在线观看| 一本大道久久a久久精品| 亚洲黑人精品在线| 亚洲欧美一区二区三区久久| 中文字幕另类日韩欧美亚洲嫩草| 99国产精品一区二区三区| 老鸭窝网址在线观看| 黄网站色视频无遮挡免费观看| 一区在线观看完整版| 亚洲全国av大片| 久久精品91无色码中文字幕| 欧美日韩黄片免| 黄色片一级片一级黄色片| 少妇精品久久久久久久| 99久久精品国产亚洲精品| 亚洲欧美色中文字幕在线| 欧美日韩一级在线毛片| 女人高潮潮喷娇喘18禁视频| 午夜成年电影在线免费观看| 亚洲国产看品久久| 国产在线精品亚洲第一网站| 国产精品.久久久| 日日夜夜操网爽| 国内毛片毛片毛片毛片毛片| 黄色成人免费大全| av网站免费在线观看视频| 精品国产超薄肉色丝袜足j| 视频在线观看一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 黄片播放在线免费| 亚洲成人国产一区在线观看| 啦啦啦视频在线资源免费观看| 国产在线视频一区二区| 王馨瑶露胸无遮挡在线观看| 美女高潮到喷水免费观看| netflix在线观看网站| 国产伦人伦偷精品视频| 在线观看舔阴道视频| 1024香蕉在线观看| 亚洲欧洲精品一区二区精品久久久| 啦啦啦免费观看视频1| 777米奇影视久久| 国产又爽黄色视频| 黄片播放在线免费| 99精品欧美一区二区三区四区| 亚洲精品av麻豆狂野| 久久中文字幕人妻熟女| 91九色精品人成在线观看| 在线看a的网站| 一边摸一边做爽爽视频免费| 亚洲av国产av综合av卡| www.自偷自拍.com| 欧美中文综合在线视频| 午夜福利乱码中文字幕| 国产成人av激情在线播放| 99国产综合亚洲精品| 精品少妇内射三级| 肉色欧美久久久久久久蜜桃| 一级黄色大片毛片| 在线看a的网站| 成年版毛片免费区| 国产欧美日韩一区二区精品| 国产精品亚洲一级av第二区| 亚洲国产av影院在线观看| 亚洲欧美色中文字幕在线| 狂野欧美激情性xxxx| 人妻 亚洲 视频| 日韩大片免费观看网站| 久久亚洲精品不卡| 黄色毛片三级朝国网站| 99在线人妻在线中文字幕 | 在线观看免费高清a一片| 国产精品九九99| 人人澡人人妻人| 丝袜美腿诱惑在线| 亚洲成人免费电影在线观看| 这个男人来自地球电影免费观看| 亚洲成人免费电影在线观看| 精品福利永久在线观看| 三上悠亚av全集在线观看| 嫁个100分男人电影在线观看| av电影中文网址| 啦啦啦 在线观看视频| 岛国在线观看网站| 麻豆国产av国片精品| 狠狠精品人妻久久久久久综合| 2018国产大陆天天弄谢| 午夜福利乱码中文字幕| 最新美女视频免费是黄的| 亚洲国产看品久久| 中文字幕最新亚洲高清| 精品视频人人做人人爽| 国产高清国产精品国产三级| 日韩精品免费视频一区二区三区| 亚洲国产中文字幕在线视频| 免费少妇av软件| av网站在线播放免费| 他把我摸到了高潮在线观看 | 亚洲九九香蕉| 中文亚洲av片在线观看爽 | 50天的宝宝边吃奶边哭怎么回事| 日本精品一区二区三区蜜桃| 欧美人与性动交α欧美软件| 香蕉国产在线看| 欧美 日韩 精品 国产| 成人三级做爰电影| 欧美在线黄色| 国产亚洲午夜精品一区二区久久| 欧美激情 高清一区二区三区| 91九色精品人成在线观看| 中国美女看黄片| 午夜两性在线视频| 成人国语在线视频| 国产一区有黄有色的免费视频| 人人妻人人添人人爽欧美一区卜| 狠狠狠狠99中文字幕| 色婷婷久久久亚洲欧美| 亚洲天堂av无毛| 三上悠亚av全集在线观看| 午夜成年电影在线免费观看| 国产成人精品在线电影| 成人永久免费在线观看视频 | 美女国产高潮福利片在线看| 最黄视频免费看| 国产淫语在线视频| 正在播放国产对白刺激| 国产有黄有色有爽视频| 中文字幕精品免费在线观看视频| 香蕉国产在线看| 国产xxxxx性猛交| 欧美另类亚洲清纯唯美| 色94色欧美一区二区| 日韩一区二区三区影片| 日本wwww免费看| 精品亚洲成国产av| 怎么达到女性高潮| 欧美精品人与动牲交sv欧美| 一区二区三区激情视频| 欧美成人免费av一区二区三区 | 一本一本久久a久久精品综合妖精| 一本综合久久免费| 亚洲七黄色美女视频| 老司机靠b影院| av福利片在线| 国产精品国产av在线观看| 99久久99久久久精品蜜桃| 久久人妻熟女aⅴ| 中文字幕制服av| 国产高清videossex| 美女高潮喷水抽搐中文字幕| 久久久久网色| 久久久精品国产亚洲av高清涩受| 搡老岳熟女国产| 婷婷丁香在线五月| 久久国产精品男人的天堂亚洲| 久久午夜综合久久蜜桃| 午夜福利在线免费观看网站| 变态另类成人亚洲欧美熟女 | 咕卡用的链子| 久久久久久久精品吃奶| av线在线观看网站| 国产精品一区二区在线不卡| 亚洲精品国产精品久久久不卡| 亚洲欧美日韩另类电影网站| 中文字幕人妻熟女乱码| 亚洲第一av免费看| 中文字幕av电影在线播放| 一级,二级,三级黄色视频| 91字幕亚洲|