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

    Experimental study of single-row chevron-jet impingement heat transfer on concave surfaces with different curvatures

    2019-12-19 02:06:42YuanweiLYUJingzhouZHANGXichenLIUYongSHAN
    CHINESE JOURNAL OF AERONAUTICS 2019年10期

    Yuanwei LYU, Jingzhou ZHANG, Xichen LIU, Yong SHAN

    Jiangsu Province Key Laboratory of Aerospace Power System, College of Energy and Power Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China

    KEYWORDS

    Chevron nozzle;

    Concave surface;

    Experimental test;

    Heat transfer;

    Jet impingement

    Abstract To address the curvature effect on single-row chevron-nozzle jet impingement heat transfer on concave surface,a series of experiments are conducted in the present investigation.Four concave surfaces including one semi-cylindrical concave surface and three parabolic concave surfaces with different width-to-depth ratios are tested under three typical Reynolds numbers (Re=5000,10000 and 15000) and several dimensionless nozzle-to-surface distances ranging from 1 to 8. The results show that the concave curvature has a clear impact on chevron-nozzle jet impingement heat transfer, tightly dependent on jet Reynolds number and impinging distance. In general, the semicylindrical concave surface produces the highest longitudinally-averaged Nusselt number at the leading line of concave surface. Under a low jet Reynolds number, the parabolic concave surface with a highly curved curvature produces higher longitudinally-averaged Nusselt number at the leading line and more uniform longitudinally-averaged Nusselt number distribution along the curvilinear direction. However, the longitudinally-averaged Nusselt number at the leading line of concave surface is the lowest for the highly curved surface under a high jet Reynolds number and large impinging distance. In comparison with the round-nozzle, chevron nozzle plays a more significant role on improving jet impingement heat transfer at small impinging distances.

    1. Introduction

    Jet impingement has been widely used in the aeronautical engineering, e.g., intensive cooling of gas turbine hot section components (such as turbine blades, combustor liners, exhaust nozzles),1-3anti-icing of aircraft critical surfaces (such as wings, engine inlet guide struts, nacelle lips).4-6Although vast investigations have been performed on this classical heat transfer scheme,is still a remarkable issue by now,in order to quest more and more efficient means for achieving the maximized heat transfer rate with limited bleed air usage.

    It had been well demonstrated that the jet impingement heat transfer could be effectively enhanced by shaping the nozzle geometry.7Among various shaped nozzles, the tabbed or chevron nozzle received much attention recently. This kind nozzle was initially motivated for the purpose of co-flow mixing enhancement.8-10As the tabs or chevrons could initiate large- and small-scale vortical motion, the near-field jet was effectively excited, resulting in strong momentum transfer between coaxial flows.Owing to this feature,the coherent vortical structure of round impinging jet would be changed by the tabbed or chevron jet. Violato et al.11performed an investigation on 3-D vortex dynamics of chevron jet impingement on a flat plate. Their results indicated that the vortical motion and turbulence intensity of chevron jet approaching to target surface were promoted by the streamwise vortices shedding from chevron notches.As a consequence,the chevron jet resulted in a stronger heat transfer enhancement. The benefit of using tabbed nozzle on improving jet impingement heat transfer was also demonstrated by Yu et al.12,13in a confined straight channel with/without presence of an initial crossflow. It was illustrated that the tabbed jets increased heat transfer in the vicinity of impingement stagnation up to 25% in the situation without initial crossflow but only 10%with initial crossflow,in comparison with the corresponding round jets. Vinze et al.14experimentally investigated the influence of chevron nozzle geometry on jet impingement heat transfer distribution over a flat plate. Ten chevron nozzle configurations were designed with different chevron number and penetration. It was revealed that an increase of 26%-38%in local Nusselt number could be achieved by the chevron nozzle in comparison with the round nozzle.Lyu et al.15performed an experimental study of chevron jet impingement on a flat surface under a series of jet Reynolds numbers ranging from 5000 to 20000 and nozzleto-surface distances ranging from 1 to 8. It was found that 6-chevron nozzle with a chevron length-to-nozzle diameter of 0.6 achieved more favorable heat transfer enhancement.By adopting the chevron nozzle, Guan et al.16,17performed numerical and experimental investigations to evaluate its application in the hot-air anti-icing configuration of engine intake surfaces.In their studies,the overall heating effectiveness was concerned by using conjugate heat transfer model (considering external cold-air forced flow, wall thermal conduction, and internal hot-air impingement). The area-averaged heating effectiveness in the vicinity of conical-wall leading edge could be increased approximately 20%by the chevron nozzle in comparison with the round nozzle.

    To our knowledge, little attention had been paid for illustrating the chevron jet impingement heat transfer on curved surfaces.As the target surface curvature is an important factor affecting the wall jet flow dynamics and consequently the heat transfer process,the conclusion derived on a flat surface is generally distinct from that on a curved surface.18-21For example,Zhou et al.22experimentally studied the local and average heat transfer characteristics of a single round jet impingement on concave surfaces. It was illustrated that the relative surface curvature(d/D)shows different effect on the heat transfer performance. At a fixed surface diameter (D), the increase of relative surface curvature (d/D) increased both stagnation and average heat transfer. However, the case was opposite at a fixed nozzle diameter (d).

    The assessment of curvature effect on single-row chevronjet impingement heat transfer on curved surfaces is the main concern of current research. To address this issue, a series of experiments are conducted in the present investigation. Four concave surfaces including one semi-cylindrical concave surface and three parabolic concave surfaces with different curvatures are taken into considerations. The comparison with the round nozzle is also made.

    2. Experimental setup

    The experimental setup is schematically displayed in Fig.1(a).The impinging air is provided by a compressed air supply. A filter is installed upstream of the pressure regulator. The flow rate is controlled and measured by a pre-calibrated mass flow controller.After routing through the supply passage,the compressed air is fed into a transmission chamber and then distributed into a row of jet tubes. Each jet tube has a length (l)of 120 mm and inner diameter (d) of 10 mm. This length-todiameter ratio (l/d) of 12 ensures a fully developed turbulent velocity profile at the nozzle exit.

    Fig. 1 Schematic test setup of jet impingement.

    The target plate assembly is composed of a supporting substrate and a thin stainless steel sheet, as seen in Fig. 1(b). The supporting substrate is made of a transparent plastic plate with a thickness of 5 mm, by manufacturing a rectangular opening at the center. On the supporting substrate, a 0.05 mm thick stainless steel sheet is tightly mounted, covering the central opening of transparent plastic plate. This stainless steel sheet is heated by DC current. By recording the electric voltage(V) between two copper bars connecting to the heating sheet and the current (I), the joule heating power (qjoule) is determined.Three jet tubes are arranged in a line along the longitudinal direction with a fixed jet-to-jet pitch (P/d) of 5. The length of concave surface (L) is 15 times nozzle diameter(15d). Both longitudinal ends of concave cavity are covered by thermally-adiabatic plate.

    To minimize the heat losses from the rear surface of heating sheet, an insulator box is adopted to enclose the target plate assembly, as seen in Fig. 1(b). On the front plate of insulator box, an infrared glass with a high transmissivity (about 0.97)in the long infrared band is placed for the temperature measurement by infrared camera. Its size is 60 mm (width)×100 mm (length)×5 mm (thickness).

    Fig. 2 Schematic diagram of test model.

    Four concave surfaces are tested in the current study,including one semi-cylindrical concave surface (baseline case)and three parabolic concave surfaces with different curvatures,as shown in Fig.2(a).The parabolic concave surface curvature is changed by varying the base width (Wc=8d, 6d and 4d,respectively) while maintaining the depth of concave cavity(Hc) at 5d. The concave profiles are listed in Table 1. The jet-nozzle centerline is perpendicular to the base plane and orthogonal to the leading line of concave cavity. A specific 6-chevrons nozzle is used in the present study, which is schematically shown in Fig. 2(b). It has six chevrons distributed uniformly in the circumferential direction.Each chevron apex axis is parallel to the centerline of the nozzle. The chevron length(a)is set as 6 mm and the arc width of chevron base (B) was πd/6, giving a length-to-width ratio of 1.15.According to the previous works,11,14the impinging distance or nozzle-to-surface distance (H) is defined as the distance between target surface and chevron orientation section. By using a traverse system, the impinging distance is adjusted in a certain range between H/d=1 and H/d=8. The coordination system is also specified in Fig. 2. The longitudinal direction of the concave surface is defined as x-direction and the curvilinear direction is defined as s-direction.For comparisons,the round-nozzle jet impingement heat transfer is also tested.

    3. Measurements and data processing

    For the multi-jet impingement,the jet Reynolds number(Re)is determined according to the nozzle diameter (d) and bulkaveraged jet velocity (uj)

    where v is the kinematic viscosity of the jet. In the current study, the jet Reynolds number is adjusted by a calibrated mass flow controller, ranging from 5000 to 15000.

    The target temperature, measurement by infrared camera,is same as that adopted by Guan et al16.Both the surface treatment(pre-sprayed with a uniform black paint)and calibration test (by using thermocouples) are conducted in advance to ensure the accurate temperature measurement.

    Considering the heat loss from the backside of heating sheet to ambient,the local convective heat transfer coefficient on the target surface is evaluated as

    where Twis the target temperature, Tjis the impinging jet temperature, qjouleis the joule heating power per unit area or the input heat flux, qsis the heat loss per unit area from the backside of heater surface.

    The heat loss from the backside surface of heating sheet to the ambient is predicted approximately by using the empirical relation of natural convection heat transfer from a vertical flat surface,23according to the measured mean temperature over the outer surface of the insulator box by eight thermocouples.

    where Taand Tbare the ambient temperature and the averaged wall temperature on the outer surface of insulator box exposed to the ambient surrounding, respectively. hnbis the natural convective heat transfer coefficient.

    In the current tests, the joule heating power per unit area(qjoule) is fixed at 4500 W/m2. The heat loss is found to be within 5% of the power input.

    The local Nusselt number (Nu) on the target surface is defined as

    where k is the thermal conductivity of the impinging jet.

    Table 2 Standard uncertainties of independent parameters.

    In the present, longitudinally-averaged Nusselt number on the target surface along curvilinear direction is adopted. As the jet-to-jet pitch is fixed at 5d in the current study, thus the longitudinally-averaged Nusselt number is determined by Total uncertainty,for the local Nu measurement,comprises the uncertainties of all parameters that affect the experiment in different ways. According to the methodology of error transfer,24the standard uncertainty of the measurement is obtained by making a combination of the input standard uncertainties and expressed as following equations:

    where f is a function of Nu in terms of input independent parameters (x1, x2, ..., xn). Δ(xi) is the standard uncertainty of independent parameter xi. In the current estimation, six independent parameters are selected according to Eq. (2) and Eq. (4), such as qjoule, qs, Tw, Tj, k and d.

    The standard uncertainties of independent parameters are summarized in Table 2.By using Eq.(7),the maximum uncertainty in the measurement of Nusselt number is approximately±7.3%.

    4. Results and discussion

    4.1. Local Nusselt number distribution

    Fig. 3 Local Nusselt number distribution on semi-cylindrical concave target by round-jet impingement (Re=10000).

    Fig. 4 Local Nusselt number distribution on semi-cylindrical concave target by chevron-jet impingement (Re=10000).

    Fig. 3 and Fig. 4 present the local Nusselt number distributions on the semi-cylindrical or baseline concave surface under Re=10000, impinged by single round jet and single chevron jet respectively, at different nozzle-to-surface distances. For the round-nozzle jet impingement, as seen in Fig. 3(a), the‘‘double peaks” feature of the local Nusselt number distribution remains vaguely on the concave surface at a small dimensionless nozzle-to-surface distance of H/d=1. As previously illustrated by Livingood and Hrycak,25for the turbulent impinging jets on flat plates, at a small nozzle-to-surface distance, the impinging jet undergoes a sudden turning of flow direction in the stagnation region, resulting in a strong flow acceleration along the radial direction to promote significant flow transition from a laminar to a turbulent boundary layer at a radial position apart from the stagnation point. Corresponding to this significant flow transition, an annular‘‘hump”with higher heat transfer coefficient occurs surrounding the jet stagnation, forming the secondary or outer peak of the local Nusselt number distribution.With regard to the concave target surface, the transition from a laminar to a turbulent boundary layer is conjectured to be faster due to the destabilization effect,and the transitioning region is less effectively cooled because the natural confinement causes the increase of the mean temperature of the ambient air around the jet.26Consequently, in relative to the flat target surface,the alleviation of the secondary peak on the concave surface is caused by the combined effect of destabilization and confinement due to the concave curvature. However, for the chevron jet impingement, the ‘‘double peaks” feature of the local Nusselt number distribution produced by the round jet impingement is disappeared but replaced by the ‘‘lobed pattern”feature that resembles the shape of the 6-chevron nozzle configuration at H/d=1,as seen in Fig.4(a).Of particular interest is that the lobe peaks of local heat transfer contours are shifted in the radial direction. In the inner lobe-shaped contour, the lobe peaks are in correspondence to the chevron notches. Far away from the jet stagnation, the lobe peaks in outer lobe-shaped contour in correspondence to the chevron apexes.At H/d=2,the lobed feature becomes weaker,as seen in Fig.4(b).With the further increase of nozzle-to-surface distance,as seen in Fig.4(c)and Fig.4(d),the characteristic lobeshaped pattern of local Nusselt number distribution observed at H/d=1 produced by the chevron jet impingement vanishes.At a large dimensionless nozzle-to-surface distance, the chevron jet impingement behaves more similar to the round jet impingement.

    Fig. 5 presents the local Nusselt number distributions on the concave surfaces,impinged by a single row of chevron jets under Re=10000 and H/d=4.It is seen clearly that the local convective heat transfer contours in the vicinity of the concave leading edge change gradually from a nearly circular-shape on the semi-cylindrical concave surface (as seen in Fig. 5(a)) to a oval-shape on the highly curved surface(as seen in Fig.5(c)).It is also seen that on the highly curved surface,the local convective heat transfer decreases rapidly along the curvilinear direction far away from the leading line of concave surface.

    Fig.7 and Fig.8 present the influences of nozzle-to-surface distance on mean and Root Mean Square(RMS)velocity contours of chevron jet impingement on a semi-cylindrical concave surface under Re=5000. At a small nozzle-to-surface distance of H/d=2, the approaching jet has a relatively narrow shear layer and big normal mean velocity, resulting in a strong flow acceleration and consequently large RMS velocity at the wall jet region. With the increase of nozzle-to-surface distance, the shear layer width of approaching jet becomes wider but the peak mean velocity decays more seriously. As the wall jet flow along the concave surface is relatively weak at a large nozzle-to-surface distance,it is the recirculation flow inside a concave cavity becomes obviously.Evidently,onto the impact of the recirculation flow inside a concave cavity, the development of impinging jet is conjectured to be destabilized.In addition, due to the confinement effect, the mean temperature of the ambient air around the jet is certainly increased.These two effects due to the concave curvature lead to an alleviation of the jet impingement heat transfer in the vicinity of jet stagnation.From these figures,it is also found that the wall jet along the longitudinal direction is stronger than that along the curvilinear direction.

    Fig. 5 Local Nusselt number distribution of a row of chevronnozzle jet impingement (Re=10000, H/d=4).

    The concave curvature has an important influence on the impinging jet flow inside a concave cavity, as demonstrated in Fig. 9. At a highly curved surface, the confinement effect due to the concave cavity becomes significantly.It is suggested that the concave curvature has two opposite roles on the jet impingement.Firstly,the approaching process of an impinging jet to the target and the development of wall jet along the curvilinear direction suffer the serious confinement effect.Secondly, the lateral diffusion of the wall jet along the longitudinal direction is otherwise promoted. Therefore, the convective heat transfer on a concave surface is tightly dependent on the competition of these two factors having the contrary roles.

    Fig. 6 Computed spatial distribution of coherent structure indicated by iso-contours of Q=500000.

    4.2. Comparison between chevron and round nozzles

    Fig. 7 Computed mean velocity contours of chevron jet impingement on concave surface under Re=5000.

    Fig. 8 Computed RMS velocity contours of chevron jet impingement on concave surface under Re=5000.

    Fig.10 and Fig.11 present the longitudinally-averaged Nusselt number distributions along curvilinear direction, on the concave surfaces with Wc/Hc=0.8 and Wc/Hc=1.6, respectively. It is confirmed that the chevron jet impingement produces stronger convective heat transfer than the corresponding round jet impingement, especially at small nozzleto-surface distances. Approximately 20%-30% increase of longitudinally-averaged Nusselt number in the vicinity of concave leading is obtained by the chevron-jet impingement under small nozzle-to-surface distances. This value is about 10%-15% at large nozzle-to-surface distances. As illustrated in Fig. 6, the chevron nozzle promotes the vortical motion of impinging jet, which is the main cause of chevron-jet heat transfer enhancement. With the increase of nozzle-to-surface distance, the vortical motion of chevron-jet will be more seriously weakened in comparison with that at a smaller nozzleto-surface distance. For this cause, the role of chevron nozzle on improving the jet impingement heat transfer on concave surfaces is alleviated at larger dimensionless nozzle-to-surface distances.

    Fig. 9 Computed flow field inside concave cavity under Re=10000.

    Fig.10 Longitudinally-averaged Nusselt number distribution of jet impingement under concave surface with Wc/Hc=0.8.

    On the highly curved surface, as seen in Fig. 10, the longitudinally-averaged Nusselt number decreases rapidly beyond s/d=2 (s is the distance from the stagnation point along the curve surface), at small nozzle-to-surface distances.As illustrated in Fig. 9, the highly curved surface presents a strong confinement to the jet impingement. The wall jet flow along the curvilinear direction is seriously affected so that a rapid decrease occurs beyond s/d=2. At large nozzle-tosurface distances, the stagnation heat transfer is obviously reduced. However, as the impinging jet has a long distance to develop,the sheer layer is broadened,as seen in Fig.8.Thus wider leading zones of concave surface are subjected to the approaching jet, which makes the longitudinally-averaged Nusselt number distributions in the vicinity of concave leading edge more smoothly when compared to the small nozzle-tosurface distances. Consequently, an intersection point of longitudinally-averaged Nusselt number distributions between small and larger nozzle-to-surface distances appears at about s/d=2.5 under the highly curved concave surface. On the slightly curved surface, as seen in Fig. 11, the confinement effect due to the concave cavity is alleviated in comparison with that on the highly curved surface. The longitudinallyaveraged Nusselt number decreases gradually along the curvilinear direction at small nozzle-to-surface distances, in the vicinity of concave leading edge. The intersection point of longitudinally-averaged Nusselt number distributions between small and larger nozzle-to-surface distances is put far away from the jet stagnation,even disappeared in the vicinity of concave leading edge.

    Fig.11 Longitudinally-averaged Nusselt number distribution of jet impingement under concave surface with Wc/Hc=1.6.

    4.3. Effect of concave curvature

    Fig. 12 and Fig. 13 present longitudinally-averaged Nusselt number distributions along the curvilinear direction,impinged by a row of chevron jets under Re=5000 and Re=15000 respectively.

    Under Re=5000, as seen in Fig. 12, the longitudinallyaveraged Nusselt number at the leading edge of semicylindrical concave surface is higher than that on the other concave surfaces. It is interest to find that the narrowest concave surface (Wc/Hc=0.8) produces nearly the same or even a little higher peak longitudinally-averaged Nusselt number at the leading edge than that on the concave surfaces with Wc/Hc=1.2 and Wc/Hc=1.6. As pointed out early, for the jet impingement on a concave surface,the wall jet flow has less momentum to flow along the curvilinear direction but larger flow capacity along the longitudinal direction of the concave surface. The strong lateral-diffusion of wall jet is conjectured to be helpful for improving the uniformity of convective heat transfer in the longitudinal direction. It is also noticed that the concave surface with Wc/Hc=0.8 produces more smooth longitudinally-averaged Nusselt number distribution along the curvilinear direction at large nozzle-to-surface distances, as seen in Fig. 12(c) and Fig. 12(d). Downstream the location of s/d=1.5, the longitudinally-averaged Nusselt number on the concave surface with Wc/Hc=0.8 is obviously higher that the corresponding one on the other concave surface.

    However, the situations under Re=15000 are not consistence with the corresponding ones under Re=5000, as seen in Fig. 13. As the jet penetration capacity is tremendously enhanced by the increase of jet ejection velocity,the approaching jet at a small nozzle-to-surface distance (such as H/d=2)is not sensitive to the concave curvature, producing nearly the same longitudinally-averaged Nusselt number at the leading edge, as seen in Fig. 13(a). Also, as the flow momentum of the wall jet along the curvilinear direction is also enhanced under higher jet Reynolds number, the longitudinallyaveraged Nusselt number distribution along the curvilinear direction on the semi-cylindrical concave surface changes more smooth than that on the other concave surfaces. As the increase of nozzle-to-surface distance, the difference of longitudinally-averaged Nusselt number at the leading edge between different concave surfaces is appeared. At large nozzle-to-surface distances, such as H/d=6 and H/d=8, as seen in Fig. 13(c) and Fig. 13(d), the longitudinally-averaged Nusselt number at the leading edge is reduced gradually with the decrease of Wc/Hc. The reason is that the approaching jet toward the target at a large nozzle-to-surface distance is more seriously affected by the strong recirculation flow in the vicinity of concave leading edge, in the case of highly curved surface and large jet Reynolds number.

    The effect of concave surface curvature shows complicated influence on the jet impingement. In order to present a global viewing of the average heat transfer for various target structures, the area-averaged Nusselt numbers (Nuav) are adopted here, as displayed in Fig. 14. Two specified areas used for the purpose of average are defined as the zone of 0 ≤s/d ≤2 and 0 ≤s/d ≤5,respectively.In general,at a low jet Reynolds number,as seen in Fig.14(a),the concave surface with a highly curved curvature (Wc/Hc=0.8) is favorable for producing a high longitudinally-averaged Nusselt number at the leading edge and uniform longitudinally-averaged Nusselt number distribution along the curvilinear direction. At a high jet Reynolds number, as seen in Fig. 14(b), the semi-cylindrical concave surface seems to be more favorable for the impingement heat transfer in the vicinity of concave leading.

    5. Conclusions

    (1) For the chevron jet impingement, the ‘‘double peaks”feature of the local Nusselt number distribution occurred in the round jet impingement is not appeared but instead by the‘‘lobe pattern” feature that resembles the shape of the 6-chevron nozzle configuration at a small nozzle-to-surface distance. At a large nozzle-tosurface distance, the chevron jet impingement behaves more similar to the round jet.

    (2) On concave surfaces,the role chevron nozzle on improving jet impingement heat transfer in related to the round nozzle is more pronounced at a small nozzle-to-surface distance. At small nozzle-to-surface distances, the longitudinally-averaged Nusselt number at leading line is increased about 20%-30%. At large nozzle-tosurface distances, this value is only about 10%-15%.

    Fig. 12 Longitudinally-averaged Nusselt number distribution of chevron-nozzle jet impingement under Re=5000.

    (3) The concave curvature has a clear impact on chevronnozzle jet impingement heat transfer, which is tightly dependent on jet Reynolds number and impinging distance. Under a low jet Reynolds number, the concave surface with a highly curved curvature produces higher longitudinally-averaged Nusselt number at the leading line and more uniform longitudinally-averaged Nusselt number distribution along the curvilinear direction.However, the longitudinally-averaged Nusselt number at the leading line of concave surface is the lowest for the highly curved surface under a high jet Reynolds number and large impinging distance.

    Fig. 14 Area-averaged Nusselt number distribution of chevronnozzle jet impingement.

    Acknowledgments

    The authors gratefully acknowledge the financial supports for this project from the National Natural Science Foundation of China (No. 51776097), the Open Fund of Jiangsu Province Key Laboratory of Aerospace Power System (No.APS20A7002), and the Postgraduate Research and Practice Innovation Project of Jiangsu Province (No. KYCX17 0280).

    又粗又硬又长又爽又黄的视频 | 国产精品久久电影中文字幕| 精品久久久久久久久久免费视频| 亚洲中文字幕日韩| 久久久久久伊人网av| 午夜亚洲福利在线播放| 亚洲成a人片在线一区二区| 真实男女啪啪啪动态图| 国产亚洲精品久久久久久毛片| 深夜a级毛片| 97在线视频观看| 99久久无色码亚洲精品果冻| 成年av动漫网址| 欧美日韩国产亚洲二区| 日产精品乱码卡一卡2卡三| 欧美精品国产亚洲| 一本久久精品| 日本三级黄在线观看| 国产精品乱码一区二三区的特点| 人人妻人人澡欧美一区二区| 噜噜噜噜噜久久久久久91| 国内精品宾馆在线| 男女那种视频在线观看| 午夜福利在线观看免费完整高清在 | 美女黄网站色视频| 一本久久精品| 国产色婷婷99| 不卡一级毛片| 国产精品人妻久久久影院| 精品久久久久久久末码| 久久久欧美国产精品| 亚洲中文字幕日韩| 91麻豆精品激情在线观看国产| 亚洲aⅴ乱码一区二区在线播放| 欧美色欧美亚洲另类二区| 国产欧美日韩精品一区二区| av在线蜜桃| 日韩中字成人| 男女视频在线观看网站免费| 亚洲国产精品国产精品| 一级黄色大片毛片| 日韩精品有码人妻一区| av卡一久久| 一卡2卡三卡四卡精品乱码亚洲| 色5月婷婷丁香| 亚洲精品日韩在线中文字幕 | 亚洲,欧美,日韩| 亚洲av中文av极速乱| 国产成人福利小说| 国产视频内射| 日韩国内少妇激情av| 2021天堂中文幕一二区在线观| 久久久精品94久久精品| 麻豆精品久久久久久蜜桃| 麻豆一二三区av精品| 亚洲成人久久性| 亚洲精品色激情综合| 日韩欧美精品免费久久| av在线老鸭窝| 成人毛片60女人毛片免费| 久久亚洲精品不卡| 在线免费观看的www视频| 日韩欧美三级三区| 18禁在线播放成人免费| 国产精品伦人一区二区| 国产精品一区二区三区四区免费观看| 少妇丰满av| 岛国在线免费视频观看| 国内久久婷婷六月综合欲色啪| 在线观看66精品国产| 欧美日韩在线观看h| www日本黄色视频网| 久久精品夜夜夜夜夜久久蜜豆| 少妇被粗大猛烈的视频| 免费看a级黄色片| 日本免费一区二区三区高清不卡| 亚洲一区二区三区色噜噜| av天堂中文字幕网| 久久鲁丝午夜福利片| 欧美激情国产日韩精品一区| 欧美日韩国产亚洲二区| 97超碰精品成人国产| 成人午夜精彩视频在线观看| 日韩欧美在线乱码| 国产黄片美女视频| 日韩制服骚丝袜av| 午夜精品在线福利| 午夜免费男女啪啪视频观看| 日本色播在线视频| 一区二区三区免费毛片| 成人美女网站在线观看视频| 麻豆久久精品国产亚洲av| 国产精品免费一区二区三区在线| 精品一区二区三区视频在线| 赤兔流量卡办理| 日韩av不卡免费在线播放| 91久久精品国产一区二区三区| 久久精品夜色国产| 伊人久久精品亚洲午夜| 国产高清激情床上av| 99国产精品一区二区蜜桃av| 久久人妻av系列| 久久精品91蜜桃| 亚洲18禁久久av| 老司机影院成人| 国产精品电影一区二区三区| 国产精品福利在线免费观看| 国产成人a∨麻豆精品| 国产乱人偷精品视频| 免费观看在线日韩| 国产一区二区亚洲精品在线观看| 久久久国产成人精品二区| 特级一级黄色大片| 亚洲综合色惰| 国产中年淑女户外野战色| kizo精华| 亚洲欧洲国产日韩| 一本久久精品| 国产精品不卡视频一区二区| 日韩精品有码人妻一区| 亚洲四区av| 成人综合一区亚洲| 神马国产精品三级电影在线观看| 少妇裸体淫交视频免费看高清| 免费看日本二区| 特级一级黄色大片| 精品一区二区免费观看| 亚洲欧美中文字幕日韩二区| 熟妇人妻久久中文字幕3abv| 久久久久久久久大av| 午夜福利视频1000在线观看| 午夜视频国产福利| 日韩中字成人| 欧美色视频一区免费| 日韩欧美在线乱码| 亚洲av男天堂| 美女 人体艺术 gogo| 尾随美女入室| 婷婷色综合大香蕉| 老司机影院成人| kizo精华| 亚洲美女视频黄频| 熟女人妻精品中文字幕| 黄色配什么色好看| 国产一区二区在线av高清观看| 中文字幕精品亚洲无线码一区| 久久久久免费精品人妻一区二区| 国产伦精品一区二区三区视频9| 最近手机中文字幕大全| 午夜免费男女啪啪视频观看| 日韩制服骚丝袜av| 久久人人爽人人片av| 黄色配什么色好看| 色哟哟·www| 成年免费大片在线观看| 少妇熟女aⅴ在线视频| 国产伦精品一区二区三区视频9| 国产亚洲av嫩草精品影院| 成人欧美大片| 国产精品av视频在线免费观看| 高清毛片免费观看视频网站| 小蜜桃在线观看免费完整版高清| 日韩欧美一区二区三区在线观看| 国产淫片久久久久久久久| 久久人人爽人人爽人人片va| 久久人人爽人人片av| 黄色一级大片看看| 亚州av有码| 狂野欧美白嫩少妇大欣赏| 熟妇人妻久久中文字幕3abv| 亚洲av成人av| 如何舔出高潮| 亚洲中文字幕日韩| a级毛片免费高清观看在线播放| a级毛色黄片| 国内精品一区二区在线观看| 变态另类丝袜制服| 校园春色视频在线观看| 午夜福利高清视频| 99久久无色码亚洲精品果冻| 免费一级毛片在线播放高清视频| 日韩精品有码人妻一区| 秋霞在线观看毛片| 亚洲人成网站在线播放欧美日韩| 日韩成人av中文字幕在线观看| 别揉我奶头 嗯啊视频| 亚洲一区二区三区色噜噜| 一级二级三级毛片免费看| 欧美极品一区二区三区四区| 亚洲精品自拍成人| 一边摸一边抽搐一进一小说| 午夜视频国产福利| 久久这里只有精品中国| 亚洲一区高清亚洲精品| 亚洲欧美精品自产自拍| 国产精品女同一区二区软件| 成人美女网站在线观看视频| 亚洲欧洲国产日韩| 久久人人爽人人片av| 一个人看的www免费观看视频| 精品午夜福利在线看| 精品无人区乱码1区二区| 国产毛片a区久久久久| 久久精品影院6| 国产69精品久久久久777片| 人人妻人人澡人人爽人人夜夜 | 免费观看在线日韩| 亚洲国产精品成人综合色| 午夜精品一区二区三区免费看| 国产极品天堂在线| 99久国产av精品国产电影| 亚洲国产日韩欧美精品在线观看| 亚洲综合色惰| 男人和女人高潮做爰伦理| 亚洲国产精品国产精品| 日韩国内少妇激情av| 日韩av不卡免费在线播放| 亚洲欧美成人综合另类久久久 | 欧美成人免费av一区二区三区| 日韩大尺度精品在线看网址| 亚洲第一区二区三区不卡| 国产成年人精品一区二区| 好男人视频免费观看在线| 99热这里只有精品一区| 伊人久久精品亚洲午夜| 久久人人精品亚洲av| 亚洲av一区综合| 色视频www国产| 日日摸夜夜添夜夜添av毛片| 日日啪夜夜撸| 久久久精品94久久精品| 国产日本99.免费观看| 久久久欧美国产精品| 2021天堂中文幕一二区在线观| 成人漫画全彩无遮挡| 夜夜夜夜夜久久久久| 国产精品人妻久久久久久| 村上凉子中文字幕在线| 国产精品精品国产色婷婷| 麻豆乱淫一区二区| 看黄色毛片网站| ponron亚洲| 91麻豆精品激情在线观看国产| 一本久久精品| 中文在线观看免费www的网站| 人体艺术视频欧美日本| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品一区二区三区人妻视频| 此物有八面人人有两片| 日本黄大片高清| 国产一区亚洲一区在线观看| 亚洲精品久久国产高清桃花| 美女高潮的动态| 日韩成人伦理影院| 国产精品福利在线免费观看| av在线天堂中文字幕| 国产一级毛片在线| 99久久成人亚洲精品观看| 看黄色毛片网站| 亚洲精品自拍成人| 深夜精品福利| 久久久国产成人免费| 国产v大片淫在线免费观看| 97超碰精品成人国产| 亚洲国产高清在线一区二区三| 精品久久久久久久久亚洲| av在线播放精品| 在线播放国产精品三级| 国产精品福利在线免费观看| 美女 人体艺术 gogo| 成人特级黄色片久久久久久久| 欧美一区二区亚洲| 99久久精品一区二区三区| 国产精品av视频在线免费观看| 国产精品一二三区在线看| 在线天堂最新版资源| 最近中文字幕高清免费大全6| 联通29元200g的流量卡| 国产伦精品一区二区三区视频9| 99热这里只有是精品在线观看| 三级国产精品欧美在线观看| 97热精品久久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 丰满乱子伦码专区| 啦啦啦观看免费观看视频高清| 色噜噜av男人的天堂激情| 国产日本99.免费观看| 色哟哟哟哟哟哟| 嫩草影院入口| 国内精品一区二区在线观看| 一个人免费在线观看电影| 欧美最黄视频在线播放免费| 免费电影在线观看免费观看| 天堂影院成人在线观看| 午夜a级毛片| 久久人人爽人人片av| 日本黄大片高清| 国产精品久久视频播放| 黄色配什么色好看| 国产黄片视频在线免费观看| 亚洲在久久综合| 九九热线精品视视频播放| 亚洲人成网站高清观看| 美女国产视频在线观看| 99热这里只有精品一区| avwww免费| 成人鲁丝片一二三区免费| 变态另类成人亚洲欧美熟女| av视频在线观看入口| 在线播放无遮挡| 高清日韩中文字幕在线| 国产精品永久免费网站| 91av网一区二区| 人妻系列 视频| 成人毛片a级毛片在线播放| 一区二区三区免费毛片| av天堂中文字幕网| 毛片一级片免费看久久久久| 看免费成人av毛片| 欧美性感艳星| 国产成人午夜福利电影在线观看| 国产一区二区在线av高清观看| 亚洲av成人精品一区久久| 亚洲欧洲国产日韩| 亚洲av不卡在线观看| 天天躁夜夜躁狠狠久久av| 国产精华一区二区三区| 免费观看在线日韩| 亚洲四区av| 美女脱内裤让男人舔精品视频 | 亚洲精华国产精华液的使用体验 | 欧美日韩国产亚洲二区| 偷拍熟女少妇极品色| 成人午夜精彩视频在线观看| 一级毛片电影观看 | 国产伦在线观看视频一区| 精品久久久久久久久av| 国产伦精品一区二区三区四那| 亚洲国产欧洲综合997久久,| 国产成人精品一,二区 | 亚洲人成网站在线观看播放| 亚洲国产欧美人成| 美女被艹到高潮喷水动态| 国产精品,欧美在线| 亚洲av不卡在线观看| 深夜a级毛片| 色视频www国产| 国产极品天堂在线| 麻豆乱淫一区二区| 中文欧美无线码| 非洲黑人性xxxx精品又粗又长| 国内精品美女久久久久久| 99国产精品一区二区蜜桃av| 看十八女毛片水多多多| 一区二区三区免费毛片| 亚洲成人av在线免费| 日韩av在线大香蕉| 黄色日韩在线| 久久午夜福利片| h日本视频在线播放| 日韩av在线大香蕉| .国产精品久久| 综合色av麻豆| 免费黄网站久久成人精品| 成年免费大片在线观看| 国产高清三级在线| 日韩欧美 国产精品| 少妇的逼好多水| 少妇人妻精品综合一区二区 | 国产成人精品久久久久久| 久久精品国产清高在天天线| 国产成人aa在线观看| 赤兔流量卡办理| 免费观看a级毛片全部| 亚洲五月天丁香| 99精品在免费线老司机午夜| 亚洲自拍偷在线| 99riav亚洲国产免费| 国内精品久久久久精免费| 男的添女的下面高潮视频| 一本一本综合久久| 12—13女人毛片做爰片一| 亚洲国产色片| ponron亚洲| 亚洲人与动物交配视频| 最好的美女福利视频网| 成人鲁丝片一二三区免费| 国产真实伦视频高清在线观看| 婷婷色av中文字幕| 午夜a级毛片| 亚洲国产精品久久男人天堂| 日日撸夜夜添| 免费观看精品视频网站| 色哟哟·www| 男女那种视频在线观看| eeuss影院久久| 51国产日韩欧美| 亚洲欧美精品自产自拍| 嫩草影院入口| 欧美日韩国产亚洲二区| or卡值多少钱| 国产av麻豆久久久久久久| 男女做爰动态图高潮gif福利片| 亚洲人与动物交配视频| 身体一侧抽搐| 嫩草影院入口| 久久鲁丝午夜福利片| 乱人视频在线观看| 欧美潮喷喷水| 亚洲18禁久久av| 国产 一区 欧美 日韩| 在线观看一区二区三区| 午夜爱爱视频在线播放| 三级男女做爰猛烈吃奶摸视频| 国产一区二区三区在线臀色熟女| 最近最新中文字幕大全电影3| 自拍偷自拍亚洲精品老妇| 插阴视频在线观看视频| 麻豆一二三区av精品| eeuss影院久久| 国产人妻一区二区三区在| 亚洲最大成人中文| 美女黄网站色视频| 国产一级毛片在线| 欧美三级亚洲精品| 美女xxoo啪啪120秒动态图| 国产亚洲av片在线观看秒播厂 | 日本一二三区视频观看| 国产午夜精品论理片| 国产精品人妻久久久影院| 精品日产1卡2卡| 又粗又爽又猛毛片免费看| 少妇熟女欧美另类| 深夜a级毛片| 精品99又大又爽又粗少妇毛片| 精品一区二区免费观看| 亚洲三级黄色毛片| 欧美性猛交╳xxx乱大交人| 日韩,欧美,国产一区二区三区 | 蜜桃久久精品国产亚洲av| 热99在线观看视频| 久久人妻av系列| 国产av不卡久久| 18+在线观看网站| 一夜夜www| 国产精品一区二区在线观看99 | 欧美精品国产亚洲| 日韩在线高清观看一区二区三区| 岛国在线免费视频观看| 夜夜看夜夜爽夜夜摸| 亚洲欧美中文字幕日韩二区| 日本与韩国留学比较| 亚洲精品456在线播放app| www.色视频.com| 男女那种视频在线观看| 欧美一区二区国产精品久久精品| 亚洲性久久影院| av在线播放精品| 精品不卡国产一区二区三区| 乱系列少妇在线播放| av天堂在线播放| 日本撒尿小便嘘嘘汇集6| 我的老师免费观看完整版| 嫩草影院新地址| 欧美激情在线99| 少妇的逼水好多| 亚洲国产欧洲综合997久久,| 国产av不卡久久| 久99久视频精品免费| 午夜免费男女啪啪视频观看| 国产高清激情床上av| 蜜桃亚洲精品一区二区三区| 免费看a级黄色片| 亚洲国产色片| 日本黄色视频三级网站网址| 可以在线观看的亚洲视频| 国产一区二区三区av在线 | 精华霜和精华液先用哪个| 精品人妻一区二区三区麻豆| 变态另类丝袜制服| 日韩国内少妇激情av| 国产成人影院久久av| 男女啪啪激烈高潮av片| 99热这里只有精品一区| 国产黄a三级三级三级人| 欧美+日韩+精品| 欧美xxxx黑人xx丫x性爽| 99久久久亚洲精品蜜臀av| 欧美丝袜亚洲另类| 国内揄拍国产精品人妻在线| 亚洲av成人精品一区久久| 噜噜噜噜噜久久久久久91| 国产精品一区二区三区四区免费观看| 久久久色成人| 精品人妻视频免费看| 级片在线观看| a级毛片a级免费在线| 少妇高潮的动态图| 久久国内精品自在自线图片| 2021天堂中文幕一二区在线观| 别揉我奶头 嗯啊视频| 精品久久久久久久久av| 久久鲁丝午夜福利片| 男人的好看免费观看在线视频| 亚洲精品自拍成人| 99热这里只有精品一区| 一边亲一边摸免费视频| 国产午夜精品论理片| 好男人视频免费观看在线| 国产免费一级a男人的天堂| 亚洲人与动物交配视频| 看非洲黑人一级黄片| 国产 一区 欧美 日韩| 青春草亚洲视频在线观看| 国产白丝娇喘喷水9色精品| 亚洲国产精品成人久久小说 | 老师上课跳d突然被开到最大视频| 黑人高潮一二区| 欧美日韩国产亚洲二区| 午夜福利视频1000在线观看| 国产极品天堂在线| 亚洲天堂国产精品一区在线| 听说在线观看完整版免费高清| 成人亚洲欧美一区二区av| 美女 人体艺术 gogo| 久久精品国产亚洲av香蕉五月| 亚洲av中文字字幕乱码综合| 男女视频在线观看网站免费| 日韩av不卡免费在线播放| 亚洲美女搞黄在线观看| 精品午夜福利在线看| 亚洲性久久影院| 一个人观看的视频www高清免费观看| 免费无遮挡裸体视频| 嫩草影院入口| 人妻久久中文字幕网| .国产精品久久| 亚洲av熟女| 国产午夜精品一二区理论片| 亚洲国产欧美在线一区| 国产精品麻豆人妻色哟哟久久 | 午夜激情福利司机影院| 男的添女的下面高潮视频| 夜夜爽天天搞| 岛国毛片在线播放| 亚洲av熟女| www.av在线官网国产| 欧美性猛交╳xxx乱大交人| 成人综合一区亚洲| 男女做爰动态图高潮gif福利片| 国产v大片淫在线免费观看| 美女大奶头视频| 极品教师在线视频| 草草在线视频免费看| a级一级毛片免费在线观看| 亚洲国产精品成人综合色| 国内少妇人妻偷人精品xxx网站| 欧美成人精品欧美一级黄| 久久久久久久久久黄片| 欧美+日韩+精品| 久久九九热精品免费| 亚洲欧美日韩无卡精品| 国产色爽女视频免费观看| 午夜亚洲福利在线播放| 亚洲中文字幕一区二区三区有码在线看| 久久久久性生活片| 色综合亚洲欧美另类图片| 精品久久久久久久人妻蜜臀av| 夜夜爽天天搞| 日韩精品青青久久久久久| 国产在视频线在精品| 嫩草影院精品99| 全区人妻精品视频| 一级av片app| 国产精品免费一区二区三区在线| 麻豆成人av视频| 日本欧美国产在线视频| 美女大奶头视频| 大又大粗又爽又黄少妇毛片口| av视频在线观看入口| 男的添女的下面高潮视频| 久久草成人影院| 男人舔女人下体高潮全视频| 国产精品永久免费网站| 欧美精品国产亚洲| 人人妻人人看人人澡| 舔av片在线| 岛国毛片在线播放| 能在线免费看毛片的网站| 日本黄色视频三级网站网址| 一进一出抽搐动态| 国产真实伦视频高清在线观看| 欧洲精品卡2卡3卡4卡5卡区| 色综合亚洲欧美另类图片| 黄色欧美视频在线观看| 麻豆国产97在线/欧美| 日韩一区二区视频免费看| 中国国产av一级| 看十八女毛片水多多多| 午夜福利视频1000在线观看| 国产高潮美女av| 非洲黑人性xxxx精品又粗又长| 免费看a级黄色片| 蜜桃亚洲精品一区二区三区| 免费观看a级毛片全部| 欧洲精品卡2卡3卡4卡5卡区| 美女黄网站色视频| 午夜精品一区二区三区免费看| 哪里可以看免费的av片| 日本一本二区三区精品| 有码 亚洲区| 亚洲国产欧洲综合997久久,| 国产精品美女特级片免费视频播放器| 亚洲精品自拍成人| 久久久久久久久久久丰满|