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

    Thermohydraulics of Turbulent Flow Through Heat Exchanger Tubes Fitted with Circular-rings and Twisted Tapes*

    2013-06-07 11:21:31SmithEiamsaardVichanKongkaitpaiboonandKwanchaiNanan

    Smith Eiamsa-ard, Vichan Kongkaitpaiboon and Kwanchai Nanan

    Department of Mechanical Engineering, Faculty of Engineering, Mahanakorn University of Technology, Bangkok 10530, Thailand

    Thermohydraulics of Turbulent Flow Through Heat Exchanger Tubes Fitted with Circular-rings and Twisted Tapes*

    Smith Eiamsa-ard**, Vichan Kongkaitpaiboon and Kwanchai Nanan

    Department of Mechanical Engineering, Faculty of Engineering, Mahanakorn University of Technology, Bangkok 10530, Thailand

    The influences of circular-ring turbulators (CRT) and twisted tape (TT) swirl generators on the heat transfer enhancement, pressure drop and thermal performance factor characteristics in a round tube are reported. The circular-ring turbulators were individually employed and together with the twisted tape swirl generators in the heated section of the tube. Three different pitch ratios (l/D=1.0, 1.5, and 2.0) of the CRT and three different twist ratios (y/W=3, 4, and 5) of the TT were introduced. The experiments were conducted using air as the working fluid under a uniform wall heat flux condition, for the Reynolds number between 6000 and 20000. The experimental results reveal that the heat transfer rate, friction factor and thermal performance factor of the combined CRT and TT are considerably higher than those of CRT alone. For the range examined, the increases of mean Nusselt number, friction factor and thermal performance, in the tube equipped with combined devices, respectively, are 25.8%, 82.8% and 6.3% over those in the tube with the CRT alone. The highest thermal performance factor of 1.42 is found for the combined device consisting of the CRT with l/D=1.0 and TT with y/W=3. The correlations of the Nusselt number, friction factor and thermal performance factor of the tubes with combined devices are also developed in terms of Reynolds number, Prandtl number, twist ratio and pitch ratio.

    heat transfer enhancement, circular-ring, twisted tape, swirl generator, turbulator, heat exchanger

    1 INTRODUCTION

    Heat exchangers play a significant role in heat transfer processes for numerous applications, for example aviation and spacecraft engineering, power engineering, refrigerating and cryogenic engineering, chemical, petroleum refining and food industries, etc. To reduce size as well as operating cost of heat exchangers, heat transfer enhancement devices are commonly applied. Heat transfer enhancement techniques are classified into two main groups. One is the active method, which requires an external power, for example, surface-fluid vibration, injection and suction of the fluid, jet impingement, and electrostatic fields. The other is the passive method, which does not involve with an external power, such as surface modifications, swirl generators and turbulence promoters. Many passive techniques have been developed for various heat exchangers. The turbulator and swirl generator are recognized as key devices of the passive techniques due to their effectiveness for heat transfer enhancement.

    Turbulators or reverse flow devices are widely employed in thermal engineering applications such as heat exchangers, combustion chambers, gas turbine blades, and electronic devices. The reverse flow or separation flow with reattachment which commonly induces high heat flux and momentum transfer as a result of the large effective driving potential force. Numerous previous investigations on improving convective heat transfer coefficient by various turbulator promoters were reported, for example truncated hollow cone, V-nozzle [1], conical-nozzle [2], conical ring [3-5], diamond-shaped [6], circular-ring [7], perforated conical-ring [8], etc. Yakut et al. [3, 4] examined the effect of pitch ratio of the conical-ring turbulator on the heat transfer and pressure drop and found that the turbulators with smallest pitch ratio provided the highest heat transfer and thermal performance. In addition, vortex-shedding frequencies and amplitude were also determined. Durmus [9] analyzed exergy loss in cut out conical turbulators to evaluate energy saving. Promvonge [5] investigated conical ring turbulator arrangements [conical-ring, diverging conical-ring (DR) array, and converging-diverging conical-ring (CDR) array] on heat transfer rate, friction factor and thermal performance factor. For the thermal performance concern, the DR array was found to be superior to CDR array. Kiml et al. [10, 11] presented the effect of the elliptic ring (angled/transverse ribs) on the flow structure and circumferential heat transfer distribution characteristics. Eventually, the angled ribs offered greater heat transfer than the transverse ribs which is responsible by a development of the rib-induced secondary flow in a form of a pair of vortices, resulting in more efficient fluid transfer between core region and wall region. Kongkaitpaiboon et al. [7] studied the effect of the circular-ring turbulator (CRT) on the heat transfer and fluid friction characteristics with various diameter ratios (d/D=0.5, 0.6 and 0.7) and pitch ratios (l/D=6, 8 and 12). The results also reveal the CRT with the smallest pitch and diameter ratios provided the highest heat transfer rate which was accompanied with the largest pressure loss. Again, Kongkaitpaiboon et al. [8] examined the influences of the perforated conical-ring (PCR) on the turbulentconvective heat transfer, friction factor ( f ) and thermal performance factor characteristics in a circular tube with three pitch ratios (l/D=4, 6 and 12) and three numbers of perforated holes (N=4, 6 and 8 holes). Evidently, the PCRs gave higher thermal performance factor than the typical conical-ring (CR), at the same pumping power. Recently, Eiamsa-ard and Promvonge [6] investigated the heat transfer and friction factor characteristics of the fully developed turbulent airflow through a uniform heat flux tube fitted with diamond-shaped turbulators in tandem arrangements at different cone angles (θ=15°, 30° and 45°), and tail length ratios (ltail-length/lhead-length=1.0, 1.5 and 2.0). At the given cone angle, heat transfer enhancement was considerably increased with decreasing tail length ratio.

    Another key device for heat transfer enhancement is swirl generator. Twisted tape insert as the common swirl generator is widely employed to generate swirl into the bulk flow for thinning the thermal/velocity boundary layer. The twisted tape also plays a vital role as the blockage of tube flow cross section which leads to a higher flow velocity and hence momentum transfers. It was reported that the performance of heat transfer enhancement of twisted tapes is strongly dependent on their geometries [12].

    Apart from the utilization of each enhancing heat transfer device, efforts have also been made by combining different types of enhancement devices known as compound devices, to further improve heat transfer. Twisted tapes were applied together with several modified tubes for improving the heat transfer rate such as tube with internal extended surfaces [13], single/three-start spirally corrugated tubes [14, 15], converging-diverging tube [16], square ducts with internal transverse rib turbulator [17], spirally grooved tube [18], micro-fin tube [19], dimpled tube [20]. In addition twisted tapes were also combined with other types of turbulators such as conical-ring [21], wire coils [22, 23], interrupted ribs in square channels [24], duct with periodic transverse rib [25], and wire nails [26].

    In general, the use of compound devices results in the increase in both heat transfer coefficient over that of individual devices. However, other factors such as the overall performance indicated by thermal performance factor above unity, the ease of installation/operation, the durability as well as the strength of devices should be taken into consideration for selecting compound devices.

    Although twisted tapes have been extensively combined with several enhancing heat transfer devices, to our best knowledge, a combination of twisted tape and circular ring has rarely been reported. Basically, both twisted tape and circular ring are acceptable for their ease of installation/operation, durability and strength. The superior performance factor of the compound device over those of the individual devices is expected owing to the synergetic enhancing effect of the reverse/separation flow induced by circular-rings and the swirling flow generated by twisted tape. Therefore, the compound device consisting of the two devices is proposed in the present work. The effects of pitch ratio (l/D=1.0, 1.5, and 2.0) of the circular-ring turbulator and twist ratio (y/W=3, 4, and 5) of the twisted tape swirl generator are also examined.

    2 EX PERIMENTAL SETUP AND OPERATING PROCEDURE

    The experimental work was conducted to reach a more general and practical look at the effect of the circular-ring and twisted-tape for both heat transfer and friction factor and thermal performance characteristics. The test section was made of copper tube with length (L) of 1500 mm, inner diameter (D) of 64 mm, and outer diameter of 67 mm. The twisted tape is made of a 0.8 mm thick and 50 mm wide aluminum strip. The tapes were prepared at three different twist ratios of y/W=3.0, 4.0 and 5.0 (y=150, 200 and 250 mm) with constant thickness (δ) of 0.8 mm for generating different swirl intensities. The optimal twist ratios of twisted tapes have been found in the range 2.5<y/W<3, indicated by reasonable heat transfer rate, pressure loss and thermal performance factor, when they being used individually. However, in the present work, the tapes at larger twist ratios are chosen to compromise the friction effect as combined with the circular-ring. A circular-ring was made of 3 mm thick aluminum sheet with its inner diameter (d) of 0.8D. The ratio d/D=0.8 of circular ring is chosen based on optimization on both heat transfer and pressure loss results. The use of a ring at smaller d/D causes significantly higher pressure drop. On the other hand, the ring with larger d/D is inefficient for fluid disturbance, resulting in poor heat transfer. The circular-rings were attached on the inner wall of the test tube with three different pitch lengths (l/D=1, 1.5 and 2). For the compound device setup, each twisted tape was subsequently put into the circular-rings. The schematic drawings of both circular-ring and twisted-tape setup are shown in Fig. 1.

    The schematic diagram of the experimental facility is demonstrated in Fig. 2. The experiments were carried out in an open-loop experimental facility consisting of a blower, orifice flow meter, test section, variac transformer, U-tube manometer, inverter, thermocouple or resistance temperature detector (RTD) and data logger. The hydrodynamic development section of a length of 2.0 m (31D) was provided to eliminate any entrance effect. The mass flow rate of the test air was measured by means of the orifice flow meter built according to Japan society of mechanical engineers (JSME) industrial standard [27]. A mixing chamber located behind the test section was used to provide a uniform temperature. The test tube was heated on outer tube surface using electrical wires wrapped around the tube to provide uniform wall heat flux conditions. The experiments were conducted at a constant nominal power of 660 W by keeping the current to be less than 3 Amps. Fifteen thermocouples were tapped on the local wall surfaces of the tube andthe thermocouples were placed round the tube in line with the same spacing to measure the circumferential temperature variation, which was found to be negligible while the whole system is insulated by wrapping insulation tape around all the exposed portions of the system. The mean local wall temperature was determined by means of calculations based on the reading of chromel-constantan thermocouples. Besides, one RTD was used at the inlet and the other two RTD at the outlet of the mixing chamber. Both the inlet and outlet temperatures of the bulk air from the tube were measured by multi-channel chromel-constantan thermocouples, calibrated within ±0.2 °C deviation by thermostat before being used.

    During the test, the inlet bulk air at 25 °C from a 2.2 kW high speed blower (TECO: model 4EVF) was directed through the orifice meter and passed to the heat transfer test section. The volumetric air flow rates from the blower were adjusted by varying motor speed through the inverter. The data obtained from the thermocouples and RTDs were recorded using a computer via a data logger. All data were taken at the steady-statecondition which was usually reached at about 1 hour after the start-up. The Reynolds number of the inlet air was varied from 6000 to 20000 based on constant volumetric air flow rate. The various characteristics of the flow, the Nusselts numbers were based on the average of tube wall temperature and outlet air temperature. The pressure drop across the test section was measured without turning on the heating unit.

    Figure 1 Diagram and photograph of circular-rings combined with a twisted tape

    Figure 2 Schematic heat transfer experiment

    The uncertainties of the reduced data were experimentally evaluated, adopting method of Ref. [28]. The maximum uncertainties of non-dimensional parameters were found to be ±4% for Reynolds number, ±5% for Nusselt number and ±7% for friction factor. The individual uncertainties of volumetric flow rate, pressure and temperature measurements were estimated to be within ±7%, ±2% and ±0.5%, respectively. All the present tests were conducted at the same inlet conditions of the test tube and the tube flow was expected to be fully developed because of the developing length of 31D before the tested section.

    3 DA TA REDUCTION

    During the test, air in the test section receives heat (Qair) from the electrical energy input to the system principally via the convective heat transfer mechanism. Thereby, Qairis assumed to be equal to the convective heat transfer within the test section.

    The heat supplied by electrical winding in the test tube is found to be 4 to 8% higher than the heat absorbed by the fluid for thermal equilibrium test due to convection and radiation heat losses from the test section to surroundings. Thus, only the heat transfer rate absorbed by the fluid is taken for internal convective heat transfer coefficient calculation. The convection heat transfer rate from the test section can be expressed as

    Substituting Eq. (2) into Eq. (3), yields

    The average Nusselt number, Nu is estimated as follows:

    The Reynolds number based on the average flow inlet velocity and the tube inlet diameter, is given by

    where U is mean air velocity of the tube and ν is a kinematic viscosity of the working fluid. In fully developed flow, the friction factor (f) can be determined by measuring the pressure drop across the test tube length (L).

    where ΔP is the pressure drop across the test tube measured by U-tube manometers, L is the test tube length and U is the mean air velocity at the entrance of the test section which was calculated from volumetric flow rate divided by the cross-section area of the tube. D is the inner diameter of the test tube at the inlet. The values of the thermo-physical properties of air were evaluated at the bulk fluid temperature from equation (4).

    For a constant pumping power,

    Eq. (10) can be rewritten in terms of friction and Reynolds number as:

    Regarding to the suggestion by Webb and Kim [29], the thermal performance factor η, is defined as the ratio of the htof an augmented tube to that of a plain tube and hpat the same pumping power. The thermal performance factor can be expressed in terms of Nusselt number and friction factor as

    where pp means pumping power.

    4 RESUL TS AND DISCUSSION

    4.1 Validation of the Experimental Setup

    In this section, the results of heat transfer (in the form of Nusselt number), pressure drop (in the form of friction factor) and thermal performance factor associated by the tube fitted with ring-turbulators and twisted tapes, are described and discussed. To validate the present plain tube and other experimental facilities, the experiment was firstly conducted without any enhancement device. The present experimental results were then compared with the predictions obtained from theprevious correlations (correlation of Dittus-Boelter for Nusselt number and correlation of Blasius for friction factor), in order to evaluate their agreements.

    Nusselt number correlations:

    Correlation of Dittus-Boelter,

    Friction factor correlations:

    Correlation of Blasius,

    Figure 3 shows that the Nusselt numbers of the present plain tube agree well with those obtained from the Dittus-Boelter correlation within ±4.9%. The results also show that Nusselt number almost increases linearly with increasing Reynolds number. Comparison in Fig. 4 reveals that the friction factors between experimental results and results predicted by Blasius correlation are in good agreement within ±5.7%.

    Figure 3 Confirmation of Nusselt number for plain tubeplain tube; Dittus-Boelter

    Figure 4 Confirmation of friction factor for plain tubeplain tube; Blasius

    The Nusselt number and friction factor for the present plain tube are correlated as follows:

    The comparisons between the experimental and the predicted data of Eqs. (16) and (17) reveal that the mean deviations of the predicted Nusselt number (Nu) and friction factor ( f ) for the present plain tube are 1.5% and 0.6%, and the corresponding maximum deviations for them are 2.3% and 1%, respectively.

    4.2 Heat transfer

    Figure 5 displays the variation of Nusselt number with Reynolds number of the plain tube, the tubes with circular-ring turbulators alone and the tubes with compound devices (circular-ring turbulator together with twisted tapes). The figure also shows the effect of pitch ratio of circular-ring turbulator and twist ratio of twisted tape. For all tubes, Nusselt number increases with increasing Reynolds number, the maximum Nusselt number was obtained at the highest Reynolds number (Re=20000). It is obvious that applying circular-ring turbulator results in a significant increase in Nusselt number above the plain tube. The simultaneous use of circular-ring turbulator and twisted tape as a compound device gives a further increase of Nusselt number, attributed to the combined effect of reverse/ separation flow by the circular-ring turbulator and swirl flow by the twisted tape. At similar operating conditions, the mean Nusselt number of using the compound devices is higher than that of employing the circular-ring alone and the plain tube around 25.8% and 20.6%, respectively. At the same Reynolds number, the Nusselt number increases with decreasing pitch ratio of circular-ring and/or twist ratio of twistedtape. Over the tested range, the maximum Nusselt number is obtained with the use of the compound device at the smallest pitch ratio (l/D=1.0) and twist ratio (y/W=3.0). This is attributed to the increase of heat transfer area as well as the increase of turbulent intensity imparted to the flow, leading to a better mixing of the fluid between the core and the tube wall region. It is also observed that the effect of the CRT pitch ratio is more prominent than that of the twist ratio. The mean Nusselt number associated by the circular ring with the lowest pitch ratio (l/D=1.0) is higher than those induced by the circular rings with l/D=1.5 and 2.0 around 16.9% and 32.6%, respectively. Meanwhile the mean Nusselt number of the tape with y/W=3.0 is higher than those with y/W=4.0 and 5.0 at about 11.4% and 23.3%, respectively.

    Figure 5 Variation of Nusselt number with Reynolds number for the tubes with heat transfer enhancement devices

    Figure 6 shows the Nusselt number ratio, Nut/Nup, against the Reynolds number value. It is clearly seen that the Nusselt number ratio value tends to decrease with the rise of Reynolds number. This signifies that the heat transfer improvement is more efficient at low Reynolds number. This can be simply explained that at lower Reynolds number, the thermal boundary is thicker, thus the disruption of the boundary by the circular-ring and twisted tape insert is more recognizable. The Nusselt number ratio values of all cases are above unity. This indicates an advantageous gain of using the tubes with circular ring and twisted tape insert over the plain tube.

    Figure 6 Variation of Nusselt number ratio with Reynolds number for the tube s with heat tr ansfer enhancement devices

    4.3 Pressure drop

    The variation of the pressure drop in terms of friction factor across the test section with Reynolds number is portrayed in Fig. 7. For all cases, the friction factor tends to reduce with increasing Reynolds number. At a given Reynolds number, the tubes with compound devices possess higher friction factors than the tubes with circular-rings alone and also the plain tube. This is due the dissipation of dynamic pressure of the fluid caused by the combined effects of (1) the flow obstruction of circular-ring and twisted tape (2) the reverse, separation and swirl flows induced by the devices and (3) the increase of flow contact area. As shown in the figure, mean friction factor in the tube with compound devices is higher than those in the tube with circular-ring alone around 82.8% and 17.6 times of the plain tube, respectively. Similar to Nusselt number, friction factor increases with decreasing pitch ratio of circular-ring and decreasing twist ratio of twisted tape. The maximum friction factor is achieved with the compound device consisting of circular-ring turbulator with the smallest pitch ratio (l/D=1.0) and twisted tape with the smallest twist ratio (y/W=3.0). The mean friction factor induced by the circular ring with the lowest pitch ratio (l/D=1.0) is higher than those given by the circular rings with l/D=1.5 and 2.0 at around 42.3% and 75.5%, respectively. The mean friction factor generated by the twisted tape with smallest twist ratio (y/W=3.0) is higher than those induced by the twisted tapes with y/W=4.0 and 5.0 by around 17.2% and 37.6%, respectively.

    Figure 8 shows the variation of friction factor ratio, ft/fp, with Reynolds number. It is seen that the friction factor ratio tends to increase with the rise of Reynolds number and the decreases of pitch ratio (l/D) and twist ratio and (y/W). According to the figure, the average/mean friction factor generated by the compound devices consisting of the twisted tape at the smallest twist ratio (y/W=3.0) and the circular rings with l/D=1.0, 1.5 and 2.0 are respectively, found to be about 33.8, 22.7 and 17 times of that for the plain tube.

    Figure 7 Variation of friction factor with Reynolds number for the tubes with heat transfer enhancement devices

    Figure 8 Variation of frictio n factor ratio with Reynolds number for the tube s with heat tr ansfer enhancement devices

    4.4 Thermal performance factor

    The effect of using the circular-ring turbulator in conjunction with twisted tape on thermal performance factor based on the same pumping power is depicted in Fig. 9. This ultimate factor beyond unity can be used as the criteria for the proper enhancement devices as well as operating condition. As can be seen from Fig. 9, the performance factors are all greater than unity (η≥1) for all enhancement devices. In addition, the performance factor tends to decrease as Reynolds number increases. The results also reveal that the performance factor increases with the reduction of the pitch ratio (l/D) and twist ratio (y/W). The performance factors of the compound devices are consistently higher than those of the circular-rings alone. This signifies the beneficial gain in view point of energy saving by using both circular-ring and twisted tape over the circular-ring alone. Roughly, the mean thermal performance factor of the compound device is higher than that of the circular-ring alone by 6.3%. For the present Re range studied, the maximum thermal performance factor of 1.42, given by the compound enhancement devices is found at the lowest Reynolds number, l/D=1.0 and y/W=3.0.

    Figure 9 Variation of therm al performance factor with Reynolds number for the tu bes w ith he at tr ansfer e nhancement devices

    4.5 Empirical correlations and their accuracies

    The Nusselt number, friction factor and thermal performance factor values of the circular-ring alone as well as the compound device (circular-ring in conjunction with twisted tape) at various pitch ratios (l/D=1.0, 1.5 and 2.0) and twist ratios (y/W=3.0, 4.0 and 5.0) are correlated (with the average relative errors in brackets) as follows:

    For circular-ring alone

    For compound device

    Evidently, the experimental data of Nusselt number, friction factor, thermal performance factor, are well represented by the above correlations.

    4.6 Comparison with the previous work

    Figures 10 to 12 show the comparisons of the present results for Nusselt number, friction factor and thermal performance factor with those in a previous work [7]. As seen in Fig. 10, under similar conditions, the Nusselt numbers in the present work (and the combined device consisting of CRT and TT as well) are consistently higher than those in [7] (only CRT used). The superior heat transfer in the present work is attributed to the combined effects of turbulent flow (by CRT) and swirl flow (by TT) resulting in more effective thermal eruption and better temperature mixing. From Fig. 11, in spite of their higher Nusselt numbers, friction factors for generated by all combined devices in the present work are significantly lower than those obtained from the use of only CRT at l/D=6.0, d/D=0.5 and comparable with those achieved from the utilization of only CRT at l/D=6.0, d/D=0.6in [7]. Regarding to the comparisons of Nusselt number and friction factor mentioned above, the consequent thermal performance factors in the present work are higher (Fig. 12). This is attributed to the prominent effect of an advantage by the use of CRT and TT (an increase of Nusselt number relative to that of the plain tube) over a disadvantage one (an increase of friction factor relative to that of the plain tube) as compared to the results by the use of only CRT [7].

    Figure 10 Comparison of Nusselt number of the presentCRT and TT with those in the previous work [7]

    Figure 11 Comparison of friction factor of the present CRT and TT with those in the previous work [7]

    Figure 12 Comparison of thermal performance factor of the present CRT and TT with those in [7]

    5 CONCLU SIONS

    In the present study, detailed heat transfer measurements were conducted in order to understand the heat transfer and friction factor inside tube inserts. The behaviors of the compound devices consisting of (1) circular-ring turbulators (CRT) with three different pitch ratios of l/D=1.0, 1.5 and 2.0 and (2) twisted tape swirl-generators (TT) with three different twist ratios of y/W=3.0, 4.0 and 5.0, were determined. The circular-ring alone and the plain tube were also tested for comparison. The order of heat transfer enhancement associated by the compound devices are about 2.36 to 4.47 times the plain tube which are accompanied by pressure loss ranging from 10.7 to 31.2 times, depending on the pitch ratio (l/D) and twist ratio (y/W) and Reynolds number (Re) values. It is found that the heat transfer rate increases with increasing Reynolds number and decreasing pitch/twist ratio. Comparatively, the compound enhancement devices provide higher heat transfer rate than the circular-ring alone around 25.8%. At a constant pumping power, the thermal performance factor for all enhancement devices is found to be higher than unity. The maximum thermal performance factor of 1.42 is reached at the lowest Reynolds number, l/D=1.0 and y/W=3.0. In addition, the correlations for Nusselt number, friction factor and thermal performance factor for the compound devices were developed as a function of the pitch ratio (l/D) and twist ratio (y/W), Prandtl number (Pr) and Reynolds number (Re) in this article.

    ACKNOWLEDGEMENTS

    The authors gratefully acknowledge the office of higher education commission CHE and mahanakorn university of technology (MUT).

    NOMENCLATURE

    A heat transfer surface area, m?2

    Cp,airspecific heat capacity of air, kJ·kg?1·K?1

    D inner diameter of test tube or inner diameter of CRT, m

    d outer diameter of CRT, m

    f friction factor

    h mean heat transfer coefficient, W·m?2·K?1

    I current, A

    k thermal conductivity of air, W·m?1·K?1

    L length of the test tube, m

    l pitch length of CRT, m

    Nu Nusselt number

    Pr Prandtl number

    ΔP pressure drop, Pa

    Q heat transfer rate, W

    Re Reynolds number

    T temperature, K

    T? mean temperature, K

    t thickness of the test tube, m

    U mean axial velocity, m·s?1

    V voltage, V

    ν kinematic viscosity, Ns·m?2

    W tape width, m

    y twist length, m

    η thermal performance factor

    ρ fluid density, kg·m?3

    Subscripts

    b bulk

    i inlet

    m mean

    o outlet

    p plain tube

    t

    turbulator

    w wall

    REFERENCES

    1 Eiamsa-ard, S., Promvonge, P., “Experimental investigation of heat transfer and friction characteristics in a circular tube fitted with V-nozzle turbulators”, Int. Comm. Heat Mass Transfer, 33, 591-600 (2006).

    2 Promvonge, P., Eiamsa-ard, S., “Heat transfer and turbulent flow friction in a circular tube fitted with conical-nozzle turbulators”, Int. Comm. Heat Mass Transfer, 34, 72-82 (2007).

    3 Yakut, K., Sahin, B., Canbazoglu, S., “Performance and flow-induced vibration characteristics for conical-ring turbulators”, Appl. Energy, 79, 65-76 (2004).

    4 Yakut, K., Sahin, B., “Flow-induced vibration analysis of conical rings used of heat transfer enhancement in heat exchanger”, Appl. Energy, 78, 273-288 (2004).

    5 Promvonge, P. “Heat transfer behaviors in round tube with conical ring inserts”, Energy Conver. Manag., 49, 8-15 (2008).

    6 Eiamsa-ard, S., Promvonge, P., “Thermal characterization of turbulent tube flows over diamond-shaped elements in tandem”, Int. J. Therm. Sci., 49, 1051-1062 (2010).

    7 Kongkaitpaiboon, V., Nanan, K., Eiamsa-ard, S., “Experimental investigation of convective heat transfer and pressure loss in a round tube fitted with circular-ring turbulators”, Int. Comm. Heat Mass Transfer, 37, 568-574 (2010).

    8 Kongkaitpaiboon, V., Nanan, K., Eiamsa-ard, S., “Experimental investigation of heat transfer and turbulent flow friction in a tube fitted with perforated conical-rings”, Int. Comm. Heat Mass Transfer, 37, 560-567 (2010).

    9 Durmus, A., “Heat transfer and exergy loss in cut out conical turbulators”, Energy Conver. Manag., 45, 785-796 (2004).

    10 Kiml, R., Mochizuki, S., Murata, A., Stoica, V., “Effects of rib-induced secondary flow on heat transfer augmentation inside a circular tube”, J. Enhanced Heat Transfer, 10, 9-20 (2003).

    11 Kiml, R., Magda, A., Mochizuki, S., Murata, A., “Rib-induced secondary flow effects on local circumferential heat transfer distribution inside a circular rib-roughened tube”, Int. J. Heat Mass Transfer, 47, 1403-1412 (2004).

    12 Jaisankar, S., Radhakrishnan, T.K., Sheeba, K.N., Suresh, S., “Experimental investigation of heat transfer and friction factor characteristics of thermosyphon solar water heater system fitted with spacer at the trailing edge of left-right twisted tapes”, Energy Conver. Manag., 50, 2638-2649 (2009).

    13 Liao, Q., Xin, M.D., “Augmentation of convective heat transfer inside tubes with three-dimensional internal extended surfaces and twisted-tape inserts”, Chem. Eng. J., 78, 95-105 (2000).

    14 Zimparov, V., “Enhancement of heat transfer by a combination of three-start spirally corrugated tubes with a twisted tape”, Int. J. Heat Mass Transfer, 44, 551-574 (2001).

    15 Zimparov, V., “Enhancement of heat transfer by a combination of a single-start spirally corrugated tubes with a twisted tape”, Exp. Therm. Fluid Sci., 25, 535-546 (2002).

    16 Hong, M.N., Deng, X.H., Huang, K., Li, Z.W., “Compound heat transfer enhancement of a converging-diverging tube with evenly spaced twisted-tapes”, Chin. J. Chem. Eng., 15, 814-820 (2007).

    17 Pramanik, D., Saha, S.K., “Thermohydraulics of laminar flow through rectangular and square ducts with transverse ribs and twisted tapes”, Int. J. Heat Transfer Trans., 128, 1070-1080 (2006).

    18 Bharadwaj, P., Khondge, A.D., Date, A.W., “Heat transfer and pressure drop in a spirally grooved tube with twisted tape insert”, Int. J. Heat Mass Transfer, 52, 1938-1944 (2009).

    19 Al-Fahed, S., Chamra, L.M., Chakroun, W., “Pressure drop and heat transfer comparison for both microfin tube and twisted-tape inserts in laminar flow”, Exp. Therm. Fluid Sci., 18, 323-333 (1998).

    20 Thianpong, C., Eiamsa-ard, P., Wongcharee, K., Eiamsa-ard, S.,“Compound heat transfer enhancement of a dimpled tube with a twisted tape swirl generator”, Int. Comm. Heat Mass Transfer, 36, 698-704 (2009).

    21 Promvonge, P., Eiamsa-ard, S., “Heat transfer behaviors in a tube with combined conical-ring and twisted-tape insert”, Int. Comm. Heat and Mass Transfer, 34, 849-859 (2007).

    22 Promvonge, P., “Thermal augmentation in circular tube with twisted tape and wire coil turbulators”, Energy Conver. Manag., 49, 2949-2955 (2008).

    23 Eiamsa-ard, S., Nivesrangsan,P., Chokphoemphun S., Promvonge, P., “Influence of combined non-uniform wire coil and twisted tape inserts on thermal performance characteristics”, Int. Comm, Heat Mass Transfer, 34, 849-859 (2010).

    24 Pal, P.K., Saha, S.K., “Thermal and friction characteristics of laminar flow through square and rectangular ducts with transverse ribs and twisted tapes with and without oblique teeth”, J. Enhanced Heat Transfer, 17, 1-21 (2010).

    25 Zhang, Y.M., Azad, G.M., Han, J.C., Lee, C.P., “Turbulent heat transfer enhancement and surface heating effect in square channels with wavy, and twisted tape inserts with interrupted ribs”, J. Enhanced Heat Transfer, 7, 5-49 (2000).

    26 Murugesan, P., Mayilsamy, K., Suresh, S., “Heat transfer and friction factor studies in a circular tube fitted with twisted tape consisting of wire-nails”, Chin. J. Chem. Eng., 18, 1038-1042 (2010).

    27 Japanese Industrial Standard, “Measurement of Fluid Flow by Means of Orifice Plates, Nozzles and Venturi Tubes”, Japanese Standards Association, Japan, JIS Z 8762-1988 (1988).

    28 ANSI/ASME, “Measurement Uncertainty”, PTC 19, 1-1985. Part I (1986).

    29 Webb, R.L., Kim N.H., Principles of Enhanced Heat Transfer, 2nd edition, Taylor& Francis Group, New York (2005).

    2011-09-22, accepted 2011-12-29.

    * Supported by the Thailand Research Fund (TRF).

    ** To whom correspondence should be addressed. E-mail: smith@mut.ac.th

    丁香六月欧美| 国产亚洲av高清不卡| 少妇熟女aⅴ在线视频| 欧美 亚洲 国产 日韩一| 激情视频va一区二区三区| 亚洲精品中文字幕在线视频| 一区二区日韩欧美中文字幕| 国产精品一区二区三区四区久久 | 窝窝影院91人妻| 级片在线观看| 一级a爱视频在线免费观看| 日本黄色视频三级网站网址| 国产精品影院久久| 国产av一区二区精品久久| 老汉色av国产亚洲站长工具| 一夜夜www| 俄罗斯特黄特色一大片| 久久九九热精品免费| 国产成人一区二区三区免费视频网站| 一卡2卡三卡四卡精品乱码亚洲| 黑人操中国人逼视频| 色老头精品视频在线观看| 黄网站色视频无遮挡免费观看| 国产精品免费视频内射| 精品人妻在线不人妻| 十八禁网站免费在线| 99国产综合亚洲精品| 99精品久久久久人妻精品| 十分钟在线观看高清视频www| 欧美中文综合在线视频| 久久久久久久久中文| 欧美亚洲日本最大视频资源| 一级a爱片免费观看的视频| 亚洲第一青青草原| 极品教师在线免费播放| 国产亚洲欧美98| 国产av一区在线观看免费| 免费在线观看日本一区| 久久精品91蜜桃| 亚洲少妇的诱惑av| 岛国视频午夜一区免费看| 欧美日本中文国产一区发布| 母亲3免费完整高清在线观看| 在线观看免费午夜福利视频| 正在播放国产对白刺激| 国产亚洲精品一区二区www| 久久国产精品影院| 欧美老熟妇乱子伦牲交| 精品免费久久久久久久清纯| 午夜精品在线福利| 国产精品一区二区免费欧美| 岛国在线观看网站| 不卡一级毛片| 欧美激情极品国产一区二区三区| 麻豆av在线久日| 最新在线观看一区二区三区| 国产精品电影一区二区三区| 国产精品1区2区在线观看.| 久久国产精品男人的天堂亚洲| 亚洲精品中文字幕一二三四区| 青草久久国产| √禁漫天堂资源中文www| 欧美最黄视频在线播放免费| 久久精品亚洲熟妇少妇任你| 在线观看一区二区三区| 欧美激情极品国产一区二区三区| a在线观看视频网站| 精品日产1卡2卡| 日韩大尺度精品在线看网址 | 久久欧美精品欧美久久欧美| 午夜免费成人在线视频| 老司机靠b影院| 亚洲久久久国产精品| 欧美日韩乱码在线| 麻豆av在线久日| 久9热在线精品视频| 黄色a级毛片大全视频| 欧美 亚洲 国产 日韩一| 欧美亚洲日本最大视频资源| 老熟妇仑乱视频hdxx| 在线免费观看的www视频| 国产精品av久久久久免费| 欧美丝袜亚洲另类 | 国产在线观看jvid| 香蕉丝袜av| 午夜精品在线福利| 日日爽夜夜爽网站| av免费在线观看网站| 免费在线观看亚洲国产| 亚洲第一青青草原| 亚洲av日韩精品久久久久久密| 90打野战视频偷拍视频| 丝袜美腿诱惑在线| 国产熟女午夜一区二区三区| 亚洲精华国产精华精| 一区在线观看完整版| 国产国语露脸激情在线看| 丝袜人妻中文字幕| 国产精品香港三级国产av潘金莲| 亚洲 欧美 日韩 在线 免费| 国产又色又爽无遮挡免费看| 国产真人三级小视频在线观看| 精品人妻在线不人妻| 一个人观看的视频www高清免费观看 | 午夜成年电影在线免费观看| 高清毛片免费观看视频网站| 亚洲自偷自拍图片 自拍| 18禁裸乳无遮挡免费网站照片 | 中文字幕人妻丝袜一区二区| 黄片小视频在线播放| 两个人视频免费观看高清| 亚洲精华国产精华精| 国产成人影院久久av| 国产一区二区激情短视频| 国产真人三级小视频在线观看| 亚洲中文字幕日韩| 精品久久久精品久久久| 日韩欧美在线二视频| 国产精品秋霞免费鲁丝片| 国产亚洲av嫩草精品影院| 香蕉丝袜av| 露出奶头的视频| 波多野结衣高清无吗| 亚洲五月色婷婷综合| 国产激情欧美一区二区| 色综合欧美亚洲国产小说| 亚洲精华国产精华精| 中文字幕人妻丝袜一区二区| 国产乱人伦免费视频| 熟妇人妻久久中文字幕3abv| 看免费av毛片| 正在播放国产对白刺激| 无人区码免费观看不卡| 69av精品久久久久久| 天堂√8在线中文| 午夜福利一区二区在线看| 神马国产精品三级电影在线观看 | 亚洲,欧美精品.| 黄色a级毛片大全视频| a在线观看视频网站| 久久天躁狠狠躁夜夜2o2o| 午夜福利免费观看在线| 在线永久观看黄色视频| 国产男靠女视频免费网站| 制服人妻中文乱码| 两人在一起打扑克的视频| 国产1区2区3区精品| 中文字幕久久专区| 正在播放国产对白刺激| 国产精品永久免费网站| 久久人人爽av亚洲精品天堂| 成人欧美大片| 一级毛片高清免费大全| 日本三级黄在线观看| 校园春色视频在线观看| av片东京热男人的天堂| 操出白浆在线播放| 在线永久观看黄色视频| 亚洲精品av麻豆狂野| 级片在线观看| 亚洲欧美日韩无卡精品| 大型黄色视频在线免费观看| 99在线人妻在线中文字幕| 国产精品av久久久久免费| 51午夜福利影视在线观看| 一区二区三区精品91| 在线观看午夜福利视频| tocl精华| 丝袜人妻中文字幕| 精品熟女少妇八av免费久了| 国产一区二区三区在线臀色熟女| 国产黄a三级三级三级人| 男女床上黄色一级片免费看| av欧美777| 婷婷精品国产亚洲av在线| 999精品在线视频| 亚洲精品国产区一区二| 欧美一区二区精品小视频在线| 夜夜看夜夜爽夜夜摸| 精品国产亚洲在线| 999久久久精品免费观看国产| 一级黄色大片毛片| 亚洲 欧美一区二区三区| 日韩中文字幕欧美一区二区| 欧美一区二区精品小视频在线| 午夜影院日韩av| 99国产精品免费福利视频| 99国产精品免费福利视频| 国产午夜精品久久久久久| 欧美日韩中文字幕国产精品一区二区三区 | 国产伦一二天堂av在线观看| aaaaa片日本免费| tocl精华| 一区二区三区精品91| 国产成人影院久久av| 18美女黄网站色大片免费观看| 亚洲黑人精品在线| 大型av网站在线播放| 高清黄色对白视频在线免费看| 亚洲九九香蕉| 国产成人精品久久二区二区91| 精品免费久久久久久久清纯| 久久久精品国产亚洲av高清涩受| 国产成人av教育| 十分钟在线观看高清视频www| 亚洲精品国产色婷婷电影| 国产在线精品亚洲第一网站| 亚洲精品av麻豆狂野| 欧美性长视频在线观看| 侵犯人妻中文字幕一二三四区| 亚洲熟女毛片儿| ponron亚洲| 精品人妻1区二区| 国产又色又爽无遮挡免费看| 一区二区日韩欧美中文字幕| 国产精品免费一区二区三区在线| 一二三四社区在线视频社区8| 国产精品自产拍在线观看55亚洲| 亚洲人成77777在线视频| 99国产精品99久久久久| 国内毛片毛片毛片毛片毛片| 欧美乱码精品一区二区三区| 日本a在线网址| 国产私拍福利视频在线观看| 亚洲欧美日韩高清在线视频| 视频在线观看一区二区三区| 亚洲男人天堂网一区| 日本免费一区二区三区高清不卡 | 国产极品粉嫩免费观看在线| 色综合站精品国产| av视频免费观看在线观看| 久久精品亚洲熟妇少妇任你| 成人亚洲精品av一区二区| 露出奶头的视频| 男人舔女人下体高潮全视频| 12—13女人毛片做爰片一| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品国产区一区二| 麻豆一二三区av精品| 国产精品一区二区在线不卡| 人妻久久中文字幕网| 操美女的视频在线观看| 久久久久久久久久久久大奶| 久久久精品欧美日韩精品| 午夜两性在线视频| 亚洲欧美日韩另类电影网站| 九色亚洲精品在线播放| 热re99久久国产66热| www.熟女人妻精品国产| tocl精华| 成年版毛片免费区| 亚洲五月天丁香| 免费高清视频大片| 日韩成人在线观看一区二区三区| 伊人久久大香线蕉亚洲五| 亚洲天堂国产精品一区在线| 久久精品国产99精品国产亚洲性色 | 亚洲第一欧美日韩一区二区三区| 咕卡用的链子| 成人18禁在线播放| 欧美 亚洲 国产 日韩一| 最新美女视频免费是黄的| av天堂久久9| 国产精品 欧美亚洲| 99在线视频只有这里精品首页| 久久精品91无色码中文字幕| 久久九九热精品免费| 精品欧美国产一区二区三| 可以在线观看的亚洲视频| 视频在线观看一区二区三区| 日本免费a在线| 亚洲av成人不卡在线观看播放网| 侵犯人妻中文字幕一二三四区| 精品国内亚洲2022精品成人| 人人澡人人妻人| 亚洲国产精品久久男人天堂| 欧美中文综合在线视频| 久久午夜综合久久蜜桃| 久久伊人香网站| 亚洲精品av麻豆狂野| 啦啦啦观看免费观看视频高清 | 国产成人精品久久二区二区免费| 咕卡用的链子| 国产单亲对白刺激| 免费看十八禁软件| 别揉我奶头~嗯~啊~动态视频| 欧美成人一区二区免费高清观看 | 人人妻人人澡欧美一区二区 | 老熟妇仑乱视频hdxx| 亚洲自拍偷在线| 亚洲国产看品久久| 一a级毛片在线观看| 国产私拍福利视频在线观看| АⅤ资源中文在线天堂| 天天躁狠狠躁夜夜躁狠狠躁| 国产蜜桃级精品一区二区三区| 琪琪午夜伦伦电影理论片6080| 精品久久久久久久毛片微露脸| 岛国视频午夜一区免费看| 亚洲精品国产精品久久久不卡| 国产精品九九99| 中亚洲国语对白在线视频| 在线观看日韩欧美| 亚洲成a人片在线一区二区| 丰满人妻熟妇乱又伦精品不卡| 高潮久久久久久久久久久不卡| 99香蕉大伊视频| 别揉我奶头~嗯~啊~动态视频| 欧美日韩一级在线毛片| 好男人在线观看高清免费视频 | 色综合婷婷激情| 黑人操中国人逼视频| 中文字幕精品免费在线观看视频| 精品不卡国产一区二区三区| 亚洲成人精品中文字幕电影| 国产亚洲精品第一综合不卡| 色综合亚洲欧美另类图片| 欧美日韩乱码在线| 麻豆国产av国片精品| 婷婷丁香在线五月| 久久国产精品人妻蜜桃| 精品国产国语对白av| 免费在线观看视频国产中文字幕亚洲| 成人特级黄色片久久久久久久| 亚洲全国av大片| 1024香蕉在线观看| 亚洲精品粉嫩美女一区| videosex国产| 亚洲无线在线观看| 亚洲人成77777在线视频| 亚洲九九香蕉| 最近最新免费中文字幕在线| 热99re8久久精品国产| 母亲3免费完整高清在线观看| 啪啪无遮挡十八禁网站| 国产激情久久老熟女| 黄色成人免费大全| 国产精品免费一区二区三区在线| 男人的好看免费观看在线视频 | 亚洲av成人不卡在线观看播放网| 欧美激情高清一区二区三区| 成人亚洲精品av一区二区| 亚洲欧美精品综合一区二区三区| 两个人免费观看高清视频| 国产亚洲av嫩草精品影院| 亚洲精品av麻豆狂野| 韩国精品一区二区三区| 国产成人精品无人区| 久久久久精品国产欧美久久久| 亚洲熟女毛片儿| 亚洲美女黄片视频| 国产一区二区三区综合在线观看| 国产一区二区三区在线臀色熟女| 露出奶头的视频| 大型av网站在线播放| 免费在线观看亚洲国产| 校园春色视频在线观看| 日日夜夜操网爽| 国产成人精品久久二区二区免费| 咕卡用的链子| 亚洲在线自拍视频| 成人精品一区二区免费| 亚洲五月色婷婷综合| 97碰自拍视频| 亚洲欧美精品综合一区二区三区| 亚洲伊人色综图| 1024香蕉在线观看| 一二三四社区在线视频社区8| 日韩av在线大香蕉| 一级片免费观看大全| 亚洲成av人片免费观看| 女人精品久久久久毛片| 欧美一级a爱片免费观看看 | 色综合欧美亚洲国产小说| 99精品在免费线老司机午夜| 亚洲一区二区三区不卡视频| 日韩精品中文字幕看吧| 制服人妻中文乱码| 老司机深夜福利视频在线观看| 国产视频一区二区在线看| 久久久国产欧美日韩av| 亚洲精品美女久久久久99蜜臀| 99热只有精品国产| 51午夜福利影视在线观看| 国产欧美日韩一区二区精品| 欧美中文日本在线观看视频| 黑人欧美特级aaaaaa片| 欧美日韩精品网址| 又黄又粗又硬又大视频| 欧美乱色亚洲激情| 成人特级黄色片久久久久久久| 人人妻,人人澡人人爽秒播| 亚洲成av片中文字幕在线观看| or卡值多少钱| 欧美日韩福利视频一区二区| svipshipincom国产片| 老司机深夜福利视频在线观看| 嫩草影视91久久| 中文亚洲av片在线观看爽| 国产一区在线观看成人免费| 一区二区日韩欧美中文字幕| 国产精品精品国产色婷婷| bbb黄色大片| 午夜久久久久精精品| 在线十欧美十亚洲十日本专区| 岛国在线观看网站| 日韩欧美国产一区二区入口| 精品国产美女av久久久久小说| av视频在线观看入口| 亚洲五月婷婷丁香| 91九色精品人成在线观看| 欧美性长视频在线观看| 天天躁夜夜躁狠狠躁躁| 男人舔女人下体高潮全视频| 一个人免费在线观看的高清视频| 欧美乱妇无乱码| 午夜久久久久精精品| 超碰成人久久| 欧美日韩黄片免| 国产真人三级小视频在线观看| 黄色 视频免费看| avwww免费| 波多野结衣av一区二区av| 国产黄a三级三级三级人| 色在线成人网| 男人操女人黄网站| 久久精品国产99精品国产亚洲性色 | 女人高潮潮喷娇喘18禁视频| 妹子高潮喷水视频| 无遮挡黄片免费观看| 人人妻人人爽人人添夜夜欢视频| av超薄肉色丝袜交足视频| 国产激情久久老熟女| 亚洲自偷自拍图片 自拍| 黑人巨大精品欧美一区二区蜜桃| 99精品欧美一区二区三区四区| 亚洲人成电影观看| 高清在线国产一区| 日韩一卡2卡3卡4卡2021年| 中亚洲国语对白在线视频| 成人永久免费在线观看视频| 88av欧美| 校园春色视频在线观看| 国产男靠女视频免费网站| 国产成人影院久久av| 欧美性长视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 黄片大片在线免费观看| 女性被躁到高潮视频| 嫩草影院精品99| 黄色a级毛片大全视频| 精品一品国产午夜福利视频| bbb黄色大片| 动漫黄色视频在线观看| 免费搜索国产男女视频| 亚洲自偷自拍图片 自拍| 一个人免费在线观看的高清视频| 50天的宝宝边吃奶边哭怎么回事| 日韩精品免费视频一区二区三区| 99久久国产精品久久久| 亚洲欧美激情在线| 亚洲aⅴ乱码一区二区在线播放 | 久久精品国产亚洲av香蕉五月| 高清毛片免费观看视频网站| www.熟女人妻精品国产| 欧美在线黄色| 夜夜夜夜夜久久久久| 天天躁夜夜躁狠狠躁躁| 91国产中文字幕| 9热在线视频观看99| 国产精品自产拍在线观看55亚洲| 免费观看精品视频网站| 他把我摸到了高潮在线观看| 91麻豆av在线| 村上凉子中文字幕在线| 日韩一卡2卡3卡4卡2021年| 成人av一区二区三区在线看| 视频在线观看一区二区三区| 免费在线观看视频国产中文字幕亚洲| 琪琪午夜伦伦电影理论片6080| 久久青草综合色| 丰满的人妻完整版| 99久久综合精品五月天人人| 色在线成人网| 免费在线观看亚洲国产| 国产欧美日韩一区二区三区在线| 亚洲人成电影观看| 亚洲人成网站在线播放欧美日韩| 国产成人精品无人区| 欧美精品啪啪一区二区三区| 国产一区二区三区在线臀色熟女| 91精品国产国语对白视频| 一二三四社区在线视频社区8| 制服人妻中文乱码| 在线永久观看黄色视频| 少妇的丰满在线观看| 韩国精品一区二区三区| 久99久视频精品免费| 动漫黄色视频在线观看| 亚洲人成电影观看| 欧美精品亚洲一区二区| 精品免费久久久久久久清纯| 国产成+人综合+亚洲专区| 久久久国产精品麻豆| 日日爽夜夜爽网站| 国产精品二区激情视频| 天堂影院成人在线观看| 国产欧美日韩一区二区三| 日本在线视频免费播放| 日本 av在线| 日本三级黄在线观看| 在线观看免费视频网站a站| 国产亚洲欧美精品永久| 波多野结衣av一区二区av| 国产成人欧美在线观看| 99国产极品粉嫩在线观看| 热99re8久久精品国产| 欧美+亚洲+日韩+国产| 69精品国产乱码久久久| 国产午夜福利久久久久久| 在线观看午夜福利视频| 丁香六月欧美| 变态另类丝袜制服| 我的亚洲天堂| 久久久水蜜桃国产精品网| 亚洲情色 制服丝袜| 他把我摸到了高潮在线观看| www.精华液| av在线播放免费不卡| 国产国语露脸激情在线看| 男人舔女人的私密视频| av福利片在线| 欧美一区二区精品小视频在线| 午夜福利在线观看吧| 女同久久另类99精品国产91| 黄色女人牲交| 法律面前人人平等表现在哪些方面| 男人舔女人下体高潮全视频| 国产激情久久老熟女| 国产91精品成人一区二区三区| 国产精品一区二区在线不卡| 亚洲第一青青草原| 国产高清videossex| 一区福利在线观看| 国内久久婷婷六月综合欲色啪| 亚洲片人在线观看| 黄色丝袜av网址大全| 久久中文看片网| 午夜免费鲁丝| 亚洲中文av在线| 9191精品国产免费久久| 91av网站免费观看| 亚洲最大成人中文| 99国产精品一区二区蜜桃av| 777久久人妻少妇嫩草av网站| 亚洲五月天丁香| 九色亚洲精品在线播放| 黄片大片在线免费观看| 色哟哟哟哟哟哟| 一进一出抽搐动态| 亚洲av成人不卡在线观看播放网| 国内精品久久久久久久电影| 村上凉子中文字幕在线| 91麻豆av在线| 久久精品国产亚洲av香蕉五月| 正在播放国产对白刺激| 亚洲av美国av| 丁香六月欧美| 国产精品自产拍在线观看55亚洲| 久久精品亚洲精品国产色婷小说| 久久久水蜜桃国产精品网| 丝袜美足系列| 亚洲九九香蕉| 老司机在亚洲福利影院| 黄片大片在线免费观看| 成人av一区二区三区在线看| 很黄的视频免费| 乱人伦中国视频| www国产在线视频色| 精品国产一区二区三区四区第35| 久久精品国产99精品国产亚洲性色 | 色哟哟哟哟哟哟| 精品欧美国产一区二区三| 别揉我奶头~嗯~啊~动态视频| 国产1区2区3区精品| 国产精品av久久久久免费| 国产亚洲精品综合一区在线观看 | 精品福利观看| 国产成人精品久久二区二区91| 老司机在亚洲福利影院| 久久性视频一级片| 18禁观看日本| 99久久精品国产亚洲精品| 午夜福利欧美成人| av天堂久久9| 国产一区二区激情短视频| 国产成人精品久久二区二区免费| 国产亚洲精品久久久久5区| 两性午夜刺激爽爽歪歪视频在线观看 | 精品午夜福利视频在线观看一区| 久久精品国产亚洲av香蕉五月| 9色porny在线观看| 久久香蕉国产精品| 久久精品国产综合久久久| 满18在线观看网站| 老司机在亚洲福利影院| 视频在线观看一区二区三区| 在线观看免费视频日本深夜| 麻豆国产av国片精品| 日韩免费av在线播放| 老司机午夜十八禁免费视频| 91成人精品电影| svipshipincom国产片| 亚洲成av片中文字幕在线观看| 亚洲国产高清在线一区二区三 |