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

    Effect of polymer and fiber additives on pressure drop in a rectangular channel *

    2017-11-02 09:09:15AmirEshghinejadfardKashyapaSharmaDominiqueThvenin
    水動力學研究與進展 B輯 2017年5期

    Amir Eshghinejadfard, Kashyapa Sharma, Dominique Thévenin

    Laboratory of Fluid Dynamics and Technical Flows, University of Magdeburg “Otto von Guericke”,Universit?tsplatz 2, 39106 Magdeburg, Germany, E-mail: amir.eshghinejadfard@ovgu.de

    Effect of polymer and fiber additives on pressure drop in a rectangular channel*

    Amir Eshghinejadfard, Kashyapa Sharma, Dominique Thévenin

    Laboratory of Fluid Dynamics and Technical Flows, University of Magdeburg “Otto von Guericke”,Universit?tsplatz 2, 39106 Magdeburg, Germany, E-mail: amir.eshghinejadfard@ovgu.de

    The influence of minute amounts of additives on pressure drop is an interesting fundamental phenomenon, potentially with important practical applications. Change of the pressure drop in a quasi-two-dimensional channel flow using various additives is experimentally investigated. Tests were conducted for a wide range of concentrations (100 ppm-500 ppm) and Reynolds numbers(16 000-36 000) with two polymers and four rigid fibers used as additive. Maximum drag reduction of 22% was observed for xanthan gum. However, xanthan gum loses its drag-reducing property rapidly. It was also seen that drag reduction percentage of xanthan gum remains almost constant for different Reynolds numbers. Guar flour demonstrated good drag reduction property at high Reynolds numbers. Drag reduction of 17.5% at R e = 33200 using 300 ppm solution was observed. However, at low Reynolds numbers guar flour will cause an increase in pressure drop. Fiber fillers (aspect ratio=21) have been tested as well. In contrast to polymers, they increased the drag for the range of examined concentrations and Reynolds numbers. Polyacrylonitrile fiber with three different aspect ratios (106, 200, 400) was also used, which showed an increase in pressure drop at low aspect ratios. Polyacrylonitrile fibers of larger lengths (6 mm) demonstrated minor drag-reducing effects (up to 3%).

    Pressure drop, polymer, fiber, channel flow

    Introduction

    Effect of suspended particles and additives on flow properties is always an important topic for both experimental and numerical studies. In this regard,some types of additives lead to lower pressure drop in the flow[1-3]while some others increase it[4,5]. Particle and fluid type in addition to flow regime and geometry play the main roles in this regard. Therefore, recognizing the conditions that can result in reduction or enhancement of pressure drop is of great importance.

    Reduction of pressure drop and consequently pumping cost is an interesting topic for many liquid transportation systems. This can be achieved by addition of small amounts of drag-reducing agents(DRA)[6]. The concept of drag reduction (DR) can be traced back to the work of Toms[7]. He discovered that introduction of a few parts per million (ppm) of longchain polymers can result in reduction of wall friction.DR leads to an increase of pipeline flow rate for the same pressure head. Most studies focus on polymers as drag reduction agent[8-11]. In this regard, long-chain and flexible polymers such as polyacrylamide and polyethylene oxide (PEO) are recognized as DRA.Onset of DR occurs when polymers are unraveled from the coiled to an extended conformation. Therefore, DR due to polymer solutions is expected to happen beyond a certain Reynolds number for a unique concentration. Below this Reynolds number usually no DR is observed.

    Two types of explanation have been given for the onset of DR when polymers are utilized. Some researchers[12-14]attribute DR to viscous effects. They reason that stretching of polymers in the regions of strong deformations of flow field leads to a higher effective viscosity in this turbulent region and consequently the thickness of buffer region is increased and velocity gradient is reduced. On the other hand, the elastic theory of Tabor and De Gennes[15]postulates that DR occurs when the elastic energy of polymers is large enough to terminate the Kolmogorov-type energy cascade. A major practical problem associated with flexible polymers is their fast thermal, biological, chemical, or mechanical degradation. This would limit their application to one-passsystems. Nowadays, rigid polymers[16], surfactants[17],bubbles[18], nanoparticles[19]and fibers[20,21]have been proven as possible alternatives for flexible polymers.

    In most cases from dilute to dense regimes,however, pressure drop increases after addition of rigid particles[22-27]. This effect is highly dependent on particle shape, Stokes number ()St and solid phase volume fraction ()φ. There are many practical cases,in which, fibers or polymers are transported by the piping system. It would be thus important to recognize the conditions that can lead to each of above mentioned pathways and to understand the effect of different additives on decrease or increase of pressure drop.

    The current work aims at experimentally investigating the influence of different additives on pressure drop in a straight rectangular channel. The experimental set-up includes a long channel of rectangular cross section to mimic a quasi-two-dimensional flow.Pressure drop measurements are performed and results are compared with those of pure water to calculate the pressure drop change. The first part of the results deals with behavior of semi-rigid polymers (xanthan gum, guar flour) and investigates their effect for different flow rates. Since less attention has been paid to rigid particles so far, the second part considers the effect of rigid particle shape, concentration and size on pressure drop in semi-dilute regime. In the end, we aim to arrive at a final method on the selection of rigid fibers to lower the pressure drop in particle-laden systems.

    1. Turbulent channel flow characteristics

    For theoretical calculations, a fully developed turbulent flow through a channel is assumed in this study. Pressure drop and flow characteristics are thus related through

    where is the fluid density,Df is Darcy friction factor,Ff is Fanning friction factor, L is distance,Dhis the hydraulic diameter and Ubis the mean(bulk) velocity. Due to large width to height ratio of the channel, Dhis considered to be equal to 4H with H being the half-channel height. Consequently,the friction velocity (uτ) is given by

    withwτ being the mean shear stress at the wall and Δ p a constant pressure gradient. The Reynolds number based on the friction velocity reads

    where is the fluid kinematic viscosity. The bulk Reynolds number is given based on the channel height by

    The Fanning friction factor for turbulent fluid flow in smooth pipes is usually described by the Prandtl-von Karman (PK) expression[28]

    In the case of water as carrier fluid, the Blasius approximation can be used as well

    Based on the Virk's maximum drag reduction(MDR) asymptote[29], the lowest friction factor that can be attained in a turbulent flow is defined by

    Therefore, the friction factor of a drag-reducing system must fall between the PK law and Virk's MDR asymptote for a given Reynolds number. In multiphase flows, DR percent is defined as the ratio of reduction in the pressure drop to the pressure drop without DRA

    Here, subscripts 0 and DR denote the states in the absence and presence of DRA, respectively. Both pressure drops are measured at the same Reynolds number based on water properties.

    Another important parameter in particle-laden turbulent flows (mainly for rigid particles) is the particle Stokes number. It is defined as the ratio of the particle response time (τp) to the characteristic time scale of the turbulent flow. A particle with low Stokes number responds quickly to the flow. Characteristic time scale of the fluid flow can be either theone is denoted by S tKand the latter one by S t+.Elghobashi[30,31]presented a classification map for particle-laden turbulent flows based on the volume fraction and Stokes number. This map shows whether the turbulent kinetic energy (TKE) in the presence of particles in dilute regime will increase or decrease.For dilutesuspension of rigid particleshe recognized two main zones based on the Kolmogorov Stokes number. Normally, for very low Stokes numbers(Microparticles) TKE increases.Forparticles do not respond immediately to velocity fluctuations of surrounding fluid. Thus,reduction of both TKE and dissipation relative to the single-phase flow is expected. For10≥ , due to vortex shedding behind the particle, TKE increases.However, in dense regimes with finite-size particles reduction of TKE and Reynolds shear stress was observed[4,32]. In contrast to TKE whose changes can be distinguished based on the solid phase volume fraction and Stokes number, changes in the pressure drop show strong dependency on particle shape and Reynolds number as well. It means that reduction of turbulence intensity does not necessarily lead to drag reduction and many other parameters must be considered. For example, in numerical studies of Zhao and van Wachem[25], the turbulence kinetic energy decreased by 6.0% for the flow with spherical particles and by 4.8% for the ellipsoidal particles, but the mean flow velocity is reduced as compared with single-phase flow. These effects will be discussed in the following sections.

    Fig.1(a) Sketch of the experimental set-up

    Fig.1(b) Channel cross section

    2. Experimental set-up and procedure

    Figure 1(a) depicts a sketch of the experimental set-up. A rectangular channel is used as the main test section. Channel length is 6 m and its cross section is rectangular with aspect ratio of 200 mm to 10 mm(Fig.1(b)). A stainless steel tank, pressure transmitter,temperature sensor (thermometer), flow meter and flow regulating pump were used in the set-up as well.The set-up works as a closed loop where the fluid runs from an open stainless steel tank through aΦ50mm diameter pipe section followed by the channel, and finally returns to the tank. Transition from the circular cross section of the pipe to the rectangular section of the channel has very gentle slope to prevent any sudden change in the flow. The screw pump regulates the fluid flow rate.

    The pressure drop measurements were performed between two pressure taps that are placed 1 m apart using a Hottinger-Baldwin differential pressure transmitter. To obtain a hydrodynamically fully developed flow region, the pressure sensors were located far (2.8 m) from the transition zone between circular and rectangular cross-section. The zero error of the sensors was noted before each measurement and subsequently used in the final results calculation. The thermometer installed on the discharge line of the tank records the temperature; the flow rate of the fluid is measured using a magnetic flow meter (Endress and Hauser). Tap water was used as working fluid and the influence of the following additives on the pressure drop was investigated:

    (1) Xanthan gum, which is a high molecular weight polysaccharide with a semi-rigid structure.Its primary structure consists of repeated pentasaccharide units formed by two glucose units, two mannose units, and one glucuronic[33]. Xanthan gum increases the viscosity of the liquid when it is added in very low concentrations at moderate temperatures.Owing to its distinct rheological properties, it is known as an important industrial biopolymer and has broad applications in food and cosmetic products.The xanthan gum E-415 from Carl Roth Co. (size 80 mesh) was used in the experiments.

    (2) Guar flour, from Carl Roth Co. (CAS 9000-30-0) was used. Guar flour mainly consists of high molecular weight polysaccharides with an average size of up to some hundred microns obtained from the guar plant seed and composed of watersoluble galactomannans which are formed from a linear chain of (14→)-linked β-D-mannopyranosyl units with (16→)-linked α-D-galactopyranosyl residues as side chains. Its main applications are in the food industry.

    (3) Carbon filler, SFC-B10 (from Schwarzw?lder Textil-Werke group) which is a short length carbon fiber consisting of very thin strands of the element carbon with density of 1 780 kg/m3was also tested. Its average diameter and length are 7 μm and 150 μm,respectively, which results in aspect ratio equal to 21.The elastic modulus of this fiber is 190 GPa.

    (4) Polyacrylonitrile fiber, which is a high strength and high modulus of elasticity white-yellowish fiber from Schwarzw?lder Textil-Werke group known as PAC251 with a specific gravity of 1.18 was used. It has various applications in textile, composite materials and aircraft industry. The fiber diameter is approximately 15 μm and available in various lengths.For the experiment, 1.6, 3 and 6 mm fibers at different concentrations were used. The elastic modulus of this type of additive is 10 GPa-12 GPa.

    Different preparation procedures were used for individual additives. For xanthan gum, the required quantity of the gum was added directly into the tank and mixed for 30 min using a hydrofoil blade mixer(tank stirrer). The system was left alone for at least 24 h to let the gum dissolve gradually in the water.The drag-reducing capacity of the xanthan gum was not instantly retrieved. Therefore, the experiments were carried out with batches of 8 h continuous operation each. Above this duration, the DR effect of the xanthan solution decreased noticeably. For each concentration, the entire set-up was cleaned and a fresh batch of sample and water was used.

    To obtain a homogeneous mixture for the fibers the samples were prepared outside the system. Small quantities of the sample (50 g-70 g) were dispersed in 3 to 5 liters of water and stirred using a small stirrer in a separate set-up. Small quantity of soap solution was used in order to obtain good dispersion of the sample.Batches of prepared sample solutions were later transferred into the main tank. Along with the above mentioned procedure, the tank stirrer was finally used for at least 2 h to obtain a homogeneous dispersed phase. For higher concentration, new batches of fiber solutions were added to the existing mixture. Each experiment took at least 3 h of continuous run.

    Fig.2 Friction factor vs. Re for water and xanthan gum

    Fig.3 DR percentage vs. Re for xanthan gum and carbon fillers

    3. Experimental results

    3.1Part I: Polymers

    3.1.1 Xanthan gum

    In this section the measurement results for water and 100 ppm and 200 ppm solutions of xanthan gum are presented. Friction factor versus Reynolds number is plotted in Fig.2. It is observed that, for the examined range of Reynolds number (Re=18000-36000), the pressure drop measured from water circulation is in close consistency with the Prandtl-von Karman (PK) correlation (Eq.(5)). Brown squares represent Virk's MDR theory (Eq.(7)). DR capability of xanthan gum can be clearly observed for both 100 and 200 ppm solutions. As expected, a higher DR percentage is achievable at higher concentrations of polymer.However,for current concentrations(100 ppm and 200 ppm) DR is lower than maximum DR asymptote.

    On the other hand, it is seen that xanthan gum loses its effectiveness rapidly. In this regard, degraded solution of 200 ppm has a lower efficiency than 100 ppm fresh solution. This occurs after only 15 h of running the apparatus. DR percentage for various Reynolds numbers is reported in Fig.3. DR percentage is computed from Eq.(8) at the same temperature of water and xanthan solution to exclude the effect of temperature on water density and viscosity. It is seen that DR initially shows some dependency on Reynolds number but tends to reach a constant value. This transition to constant DR percentage occurs at higher Reynolds number when the polymer concentration is higher. Maximum DR percentage in the range of examined flow rates was observed atRe=33890 standing at 22% for a solution containing 200 ppm xanthan gum.

    It must be noted that the current channel is rather narrow. In larger systems, use of higher concentrations of polymers is easier without formation of clumps. The reduced efficiency of polymers can be obviously observed from Fig.3. This reduction of efficiency is mainly attributed to de-aggregation of rigid polymers[34].

    3.1.2 Guar flour

    The next polymer that is studied in this part is guar flour. Two concentrations of 100 ppm and 300 ppm are considered. The results of changes in pressure drop are shown in Fig.4. This figure illustrates that guar flour does not demonstrate any drag reduction at Reynolds numbers less than 21 000.In the range of low Reynolds numbers, both concentrations led to pressure drop increase. A maximum increase of 9.4% was observed for 100 ppm guar gum solution. Furthermore, the amount of this increase is almost similar for both 100 ppm and 300 ppm solutions. AtRe? 2 2000 guar flour neither increases, nor decreases the pressure drop. At higher Reynolds numbers (Re>2 2000),a noticeable drag reduction is seen from the graph. Higher concentration (300 ppm) yields higher DR percentage.Its corresponding graph has a higher slope as well.This trend is approximated by a third order polynomial fit for both concentrations. One can get 17.5% drag reduction atRe= 33200 when 300 ppm guar flour is used. This percentage remained constant after 3 h of running the system. Maximum DR percentage attained by using 100 ppm solution was 10.8% atRe= 31300.

    Fig.4 DR percentage vs. Re for guar flour

    Fig.5 DR percentage vs. Re for Polyacrylonitrile fiber

    3.2Part II: Fibers

    3.2.1 Carbon fillers

    The influence of carbon fillers (SFC-B10) (with concentrations of 300 ppm and 500 ppm) on pressure drop is shown in Fig.3 as a function of Reynolds number and compared with those from xanthan gum.It is interesting to see that carbon fillers produce no DR effect. In contrast, this type of rigid fibers leads to extra pressure drop in the system. The pressure drop increase was up to 5%. This effect can be attributed to their low length to diameter ratio and relatively high specific gravity and will be discussed in more details in Section 3.3. Carbon fillers are not able to modify the near-wall turbulence structure even at higher concentration of 500 ppm. Although the increase in pressure drop is slightly higher at low Re values, the slope of the trend line is very low. It can be seen also that results of both 300 ppm and 500 ppm are very close to each other.

    3.2.2 Polyacrylonitrile fiber

    The results corresponding to this type of fiber(PAC251) are given in Fig.5. Three different particle lengths namely 1.6 mm, 3 mm and 6 mm have been examined. It is seen from the figure that, the shortest fiber increases the pressure drop in the pipeline for the whole range of flow rates examined. Maximum increase of 5% in the pressure drop was observed for 300 ppm of 1.6 mm length fiber. The results for both 300 and 500 ppm concentrations are in the same range.This effect is due to low Stokes number of particles.By using 3 mm long fibers, pressure drop remains close to that of pure water but minor drag reduction effect was observed. Interestingly, 6 mm length fibers atc=500ppm showed drag reduction up to 3%.This part of the experiment proves that, in the case of rigid fibers, the ratio of fiber length to diameter plays the main role in DR. At low aspect ratios, particles cannot change the flow field structures. They can lead to a lower pressure drop only when their aspect ratio is sufficiently high.

    It must be pointed out that, we should not expect high levels of drag reduction in dilute regimes of fiber-laden flows. This statement is supported by numerical studies of Zhao et al.[35]. They examined elongated ellipsoids with major and minor radii of 480 μm and 96 μm at volume fraction ofφ= 3.6×10-4andSt+=τp/τV=30. Drag reduction of less than 1% was observed in their simulations.

    3.3Physical effect of additives

    We are interested in the effect of additives on flow field based on the results given. In particular,finding the conditions that fibers decrease the pressure drop is important.

    In DR by polymers the change of polymer structure from a coiled to a stretched one is an important factor. The strain rate and vorticity fields in the buffer layer make it possible that this stretching leads to a larger thickness of the buffer layer[11]. Therefore,higher effective viscosity in this region leads to the reduction of momentum flux to the wall more than viscous drag enhancement and consequently dissipative eddies are damped. This effect is also influenced by Reynolds number because this transition to stretched formation is triggered easier at higher Reynolds numbers. For xanthan gum the transition occurs at relatively low Reynolds numbers.

    Now the reason behind the change in pressure drop by rigid fibers will be discussed. Turbulence intensity and pressure drop of particle-laden flows are influenced by particle shape, density, diameter (Stokes number) and volume fraction.

    The particle response time for ellipsoid particles can be evaluated by[36]

    where and a are the aspect ratio and length of minor-axis of the particle, respectively.

    For carbon fillers examined in this paper(SFC-B10, Section 3.2.1) particle size and density lead to very low Stokes number of particle based on friction velocity+0.04-0.13 St ≈ . These particles are referred to as “microparticles”[30]. Due to relatively small size of particles they can follow the streamlines which leads to higher flow turbulence.

    In general, total change of stresses in the system can be understood by the streamwise momentum budget which can be expressed by three terms[26]

    where is the total stress, τTis the turbulent Reynolds shear stress, τVis the viscous stress and τSis the stress induced by solid phase. In dilute regimes of carbon fillers, the increase of pressure drop can be justified by particle-induced stresses along with the increase of turbulence intensity of the fluid itself.The effect of particle-induced stresses on pressure drop enhancement was also reported in numerical studies of Ref.[26].

    Polyacrylonitrile fibers' behavior depends strongly on particle size and aspect ratio. By increasing the particle length from lp=1.6 mm to 6 mm particle Stokes number ( S t+) increases (from 0.58 to 0.75).Consequently, particle effect on the flow is expected to change from “ghost” particles (for lp=1.6mm) to“critical” particles (for lp=6mm). Ghost particle does not influence the TKE. However, viscous stressesVτ and particle-induced stressesSτ lead to higher pressure drop in the system. By further increasing the particle length (and aspect ratio) critical regime is reached. In the scope of critical regime,TKE of the fluid phase is reduced. These two effects of reduction of fluid flow turbulent intensity and particle-induced stress, counteract each other at some point and further increase of particle Stokes number results in drag reduction for 6 mm long fibers.

    Larger particleStokes numbers+( 1)St≥which can occur as a result of big particle selection(not present in this study) induce extra effect on turbulence intensity augmentation in dilute regime.Vortex shedding behind these particles gives rise to turbulence intensity even in the dilute regime.

    Effect of volume fraction in this study can be explained by tendency of the system to move from a dilute to dense regime. Approaching dense regime means that particle can change the structure of streamwise vortices, weaken their strength, thus reducing turbulence intensity. It must be noted also that fibers are more efficient in systems with large enough pipes and less elbows so that high volume fraction without any worry on system blockage can be reached. Furthermore, fibers are more resistant,particularly at high Reynolds numbers, where polymers are subject to high scission rates. It would also be helpful to select particles of better elasticity.

    4. Conclusions

    The influence of polymers and fibers on pressure drop in a channel flow of rectangular cross section has been experimentally investigated. Xanthan gum, guar flour, carbon filler and polyacrylonitrile fibers have been considered as additive. Based on our observations, following conclusions can be drawn:

    (1) A solution of 200 ppm xanthan gum revealed a noticeable DR of 22% at R e = 33000. However,this DRA loses its efficiency rapidly. Guar flour is another polymer examined in this experiment. At low Reynolds numbers, guar flour increased pressure drop for both 100 ppm and 300 ppm solutions. However, a drag reduction of up to 18% can be achieved with 300 ppm of guar flour.

    (2) Carbon fillers at 300 ppm and 500 ppm were tested as well. In both concentrations, an increase in pressure drop was seen, which is justified by their low ratio of length to diameter and low Stokes number.

    (3) Polyacrylonitrile fiber is another additive which is examined and introduced as a potential DRA with good rigidity. It demonstrates up to 3% drag reduction at high particle lengths with only 500 ppm concentration. Therefore, rigid particles decrease the pressure drop in the channel flow only when they have sufficiently high aspect ratios. However, due to their better resistance, mechanical degradation does not occur. This fiber is expected to yield more drag reduction if used in higher concentrations. In this study, due to small diameter of piping we tested concentrations up to 500 ppm.

    Based on the results shown, for future studies with rigid particles it would be helpful to use particles with following characteristics:

    (1) It is recommended to use particles of high aspect ratios (>100) to get more DR. For the same aspect ratio, moderate particle Stokes numbers are preferred. For Stokes numbers in the upper range of criteria mentioned, lower particle aspect ratios can be used. However, it must be noted that very large particles (with higher Stokes numbers) experience wake shedding, which increases turbulence and pressure drop consequently. Then, very large volume fractions and Reynolds numbers are probably required to get small levels of DR.

    (2) The next point is the proper selection of solid-phase volume fraction. Higher volume fractions of moderate Stokes number particles usually lead to more attenuation of turbulence.

    (3) Fibers can be preferably used as DRA in systems of relatively large size with few elbows to make high solid-phase volume fractions possible without any pipe blockage.

    (4) Better elasticity of particle is useful to reduce the pressure drop more but it is not a necessary parameter.

    For future studies, using higher volume fractions of rigid particles and high Stokes number particles are recommended.

    Acknowledgements

    The financial support of the German research foundation (DFG) within the graduate college for micro-macro-interactions in structured media and particle systems (GRK 1554) is gratefully acknowledged.

    [1] Shao X., Lin J., Wu T. et al. Experimental research on drag reduction by polymer additives in a turbulent pipe flow [J]. Canadian Journal of Chemical Engineering,2002, 80(2): 293-298.

    [2] Edomwonyi-Otu L. C., Chinaud M., Angeli P. Effect of drag reducing polymer on horizontal liquid–liquid flows[J]. Experimental Thermal and Fluid Science, 2015, 64:164-174.

    [3] Steele A., Bayer I. S., Loth E. Pipe flow drag reduction effects from carbon nanotube additives [J]. Carbon, 2014,77(2): 1183-1186.

    [4] Shao X., Wu T., Yu Z. Fully resolved numerical simulation of particle-laden turbulent flow in a horizontal channel at a low Reynolds number [J]. Journal of Fluid Mechanics, 2012, 693: 319-344.

    [5] Uhlmann M. Interface-resolved direct numerical simulation of vertical particulate channel flow in the turbulent regime [J]. Physics of Fluids, 2008, 20(5): 53305.

    [6] Gyr A., Bewersdorff H.-W. Drag reduction of turbulent flows by additives [M]. Dordrecht, The Netherlands:Kluwer Academic, 1995.

    [7] Toms B. A. Some observations on the flow of linear polymer solutions through straight tubes at large Reynolds numbers [C]. First International Congress on Rheology.Scheveningen, The Netherlands, 1948, 135-141.

    [8] Warholic M. D., Massah H., Hanratty T. J. Influence of drag-reducing polymerson turbulence:Effectsof Reynolds number, concentration and mixing [J]. Experiments in Fluids, 1999, 27(5): 461-472.

    [9] Ptasinski P. K., Nieuwstadt F. T. M., Van Den Brule B. et al. Experiments in turbulent pipe flow with polymer additives at maximum drag reduction [J]. Flow, Turbulence and Combustion, 2001, 66(2): 159-182.

    [10] Mowla D., Naderi A. Experimental study of drag reduction by a polymeric additive in slug two-phase flow of crude oil and air in horizontal pipes [J]. Chemical Engineering Science, 2006, 61(5): 1549-1554.

    [11] White C. M., Mungal M. G. Mechanics and prediction of turbulent drag reduction with polymer additives [J].Annual Review of Fluid Mechanics, 2008, 40(1): 235-256.

    [12] Lumley J. L. Drag reduction by additives [J]. Annual Review of Fluid Mechanics, 1969, 1(1): 367-384.

    [13] Ryskin G. Turbulent drag reduction by polymers: A quantitative theory [J]. Physical Review Letters, 1987,59(18): 2059.

    [14] Gillissen J. J. J. Polymer flexibility and turbulent drag reduction [J]. Physical Review E Statistical Nonlinear and Soft Matter Physics, 2008, 78(4): 046311.

    [15] Tabor M., De Gennes P. G. A cascade theory of drag reduction [J]. Europhysics Letters, 1986, 2(7): 519-522.

    [16] Japper-Jaafar A., Escudier M. P., Poole R. J. Turbulent pipe flow of a drag-reducing rigid “rod-like” [J]. Journal of Non-Newtonian Fluid Mechanics, 2009, 161(1-3):86-93.

    [17] Cai S. P., Higuchi Y. Drag-reduction behavior of an unusual nonionic surfactant in a circular pipe turbulent flow[J]. Journal of Hydrodynamics, 2014, 26(3): 400-405.

    [18] Elbing B. R., Winkel E. S., Lay K. A. et al. Bubble-induced skin-friction drag reduction and the abrupt transition to air-layer drag reduction [J].Journal of Fluid Mechanics, 2008, 612: 201-236.

    [19] Pouranfard A. R., Mowla D., Esmaeilzadeh F. An experimental study of drag reduction by nanofluids through horizontal pipe turbulent flow of a Newtonian liquid [J].Journal of Industrial and Engineering Chemistry, 2014,20(2): 633-637.

    [20] Radin I.,. Zakin J. L., Patterson G. K. Drag reduction in solid-fluid systems [J]. AIChE Journal, 1975, 21(2):358-371.

    [21] You Z. J., Lin J. Z., Shao X. M. et al. Stability and drag reduction in transient channel flow of fibre suspension [J].Chinese Journal of Chemical Engineering, 2004, 12(3):319-323.

    [22] Ko G. H., Heo K., Lee K. et al. An experimental study on the pressure drop of nanofluids containing carbon nanotubes in a horizontal tube [J]. International Journal of Heat and Mass Transfer, 2007, 50(23-24): 4749-4753.

    [23] Liu Z. H., Liao L. Forced convective flow and heat transfer characteristics of aqueous drag-reducing fluid with carbon nanotubes added [J]. International Journal of Thermal Sciences, 2010, 49(12): 2331-2338.

    [24] Lin J. Z., Xia Y., Ku X. K. Flow and heat transfer characteristics of nanofluids containing rod-like particles in a turbulent pipe flow [J]. International Journal of Heat and Mass Transfer, 2016, 93: 57-66.

    [25] Zhao F., van Wachem B. G. M. Direct numerical simulation of ellipsoidal particles in turbulent channel flow [J].Acta Mechanica, 2013, 224(10): 2331-2358.

    [26] Picano F., Breugem W.-P., Brandt L. Turbulent channel flow of dense suspensions of neutrally buoyant spheres [J].Journal of Fluid Mechanics, 2015, 764: 463-487.

    [27] Lin J. Z., Xia Y., Ku X. K. Pressure drop and heat transfer of nanofluid in turbulent pipe flow considering particle coagulation and breakage [J]. Journal of Heat Transfer,2014 , 136(11): 111701.

    [28] von Kármán T. Mechanische ?hnlichkeit und turbulenz[C]. Proceedings of the 3rd International Congress for Applied Mechanics. Stockholm, Sweden, 1930, 85-105.

    [29] Virk P. S. Drag reduction fundamentals [J]. AIChE Journal, 1975, 21(4): 625-656.

    [30] Elghobashi S. An updated classification map of particleladen turbulent flows [C]. IUTAM Symposium on Computational Approaches to Multiphase Flow. Berlin,Germany: Springer, 2006, 3-10.

    [31] Elghobashi S. On predicting particle-laden turbulent flows[J]. Applied Scientific Research, 1994, 52(4): 309-329.

    [32] Fornari W., Formenti A., Picano F. et al. The effect of particle density in turbulent channel flow laden with finite size particles in semi-dilute conditions [J]. Physics of Fluids, 2016, 28(3): 033301.

    [33] Garc?a-Ochoa F., Santos V. E., Casas J. A. et al. Xanthan gum: Production, recovery, and properties [J]. Biotechnology Advances, 2000, 18(7): 549-579.

    [34] Soares E. J., Sandoval G. A. B., Silveira L. et al. Loss of efficiency of polymeric drag reducers induced by high Reynolds number flows in tubes with imposed pressure [J].Physics of Fluids, 2015, 27(12): 125105.

    [35] Zhao F., George W. K., van Wachem B. G. M. Four-way coupled simulations of small particles in turbulent channel flow: The effects of particle shape and Stokes number [J].Physics of Fluids, 2015, 27(8): 083301.

    [36] Shapiro M., Goldenberg M. Deposition of glass fiber particles from turbulent air flow in a pipe [J]. Journal of Aerosol Science, 1993, 24(1): 65-87.

    March 16, 2017, Revised April 10, 2017)

    * Biography: Amir Eshghinejadfard (1983-), Male,Ph. D. Candidate

    av国产久精品久网站免费入址| 亚洲欧美清纯卡通| 亚洲av免费高清在线观看| 色视频www国产| 国产色婷婷99| 我的老师免费观看完整版| 免费播放大片免费观看视频在线观看| 成人漫画全彩无遮挡| av不卡在线播放| 另类精品久久| 亚洲精华国产精华液的使用体验| 亚洲国产精品一区二区三区在线| 日日撸夜夜添| av卡一久久| 熟女人妻精品中文字幕| av在线老鸭窝| 大陆偷拍与自拍| av免费观看日本| 久久久精品94久久精品| 国产黄色视频一区二区在线观看| 晚上一个人看的免费电影| 简卡轻食公司| 久久热精品热| 欧美日韩一区二区视频在线观看视频在线| 女人久久www免费人成看片| 亚洲经典国产精华液单| 九九在线视频观看精品| 日韩人妻高清精品专区| 国产男人的电影天堂91| 五月天丁香电影| 老女人水多毛片| 性高湖久久久久久久久免费观看| 日本wwww免费看| 少妇裸体淫交视频免费看高清| 一级黄片播放器| 青春草视频在线免费观看| 美女脱内裤让男人舔精品视频| 久久精品久久久久久噜噜老黄| 亚洲国产精品专区欧美| 嫩草影院新地址| 久久精品夜色国产| 国产在线男女| 乱人伦中国视频| 又爽又黄a免费视频| 亚洲精品成人av观看孕妇| 免费黄色在线免费观看| 国产老妇伦熟女老妇高清| 在线免费观看不下载黄p国产| 国产69精品久久久久777片| 亚洲精品国产av蜜桃| 在线观看www视频免费| 青春草国产在线视频| 中文字幕制服av| 精品亚洲乱码少妇综合久久| 高清在线视频一区二区三区| 精品人妻一区二区三区麻豆| 成人美女网站在线观看视频| 国产伦理片在线播放av一区| 国产永久视频网站| av一本久久久久| 在线 av 中文字幕| 久久久精品免费免费高清| 国产成人精品久久久久久| 亚洲av成人精品一二三区| 九色成人免费人妻av| 黄色毛片三级朝国网站 | a级毛色黄片| 欧美 日韩 精品 国产| 成年人免费黄色播放视频 | 精品国产一区二区三区久久久樱花| 99视频精品全部免费 在线| 久久热精品热| 精品久久久久久久久亚洲| 亚洲自偷自拍三级| 成年人午夜在线观看视频| 丁香六月天网| 亚洲第一区二区三区不卡| 色吧在线观看| 下体分泌物呈黄色| 久久久欧美国产精品| 精品人妻偷拍中文字幕| 亚洲欧洲精品一区二区精品久久久 | 中文天堂在线官网| 高清av免费在线| 精品一区二区三区视频在线| 国产淫语在线视频| 午夜日本视频在线| 亚洲精品久久午夜乱码| 嘟嘟电影网在线观看| 国产精品99久久久久久久久| 久久免费观看电影| 成人国产av品久久久| 国产精品欧美亚洲77777| 久久午夜综合久久蜜桃| videossex国产| 久久精品国产亚洲网站| 久久久久久久久久久丰满| 国产在线男女| 精品国产露脸久久av麻豆| 亚洲精品日韩在线中文字幕| 精品国产乱码久久久久久小说| 我的女老师完整版在线观看| 亚洲精品国产成人久久av| 少妇猛男粗大的猛烈进出视频| 老司机亚洲免费影院| 色视频www国产| 在线播放无遮挡| 寂寞人妻少妇视频99o| 一区在线观看完整版| 国产伦精品一区二区三区视频9| 日韩三级伦理在线观看| 欧美少妇被猛烈插入视频| 久久97久久精品| 精品酒店卫生间| 免费看不卡的av| 岛国毛片在线播放| 成人毛片60女人毛片免费| 只有这里有精品99| 日韩三级伦理在线观看| 国产成人免费无遮挡视频| 99热这里只有精品一区| 美女脱内裤让男人舔精品视频| 18禁裸乳无遮挡动漫免费视频| 一本一本综合久久| 精品久久久噜噜| 午夜福利,免费看| 国产成人免费无遮挡视频| 国产黄片视频在线免费观看| 国产免费一区二区三区四区乱码| 秋霞在线观看毛片| 老女人水多毛片| 91精品一卡2卡3卡4卡| 精品国产一区二区三区久久久樱花| 亚洲真实伦在线观看| 丰满人妻一区二区三区视频av| 18禁在线播放成人免费| 色婷婷久久久亚洲欧美| 日韩电影二区| 3wmmmm亚洲av在线观看| 日韩亚洲欧美综合| 国产毛片在线视频| 亚洲欧美成人综合另类久久久| 日韩欧美精品免费久久| 少妇被粗大的猛进出69影院 | 嘟嘟电影网在线观看| 久久久国产精品麻豆| 99热网站在线观看| 久久久久国产精品人妻一区二区| 日本黄色日本黄色录像| 中国三级夫妇交换| 老司机亚洲免费影院| 少妇的逼水好多| 一区二区三区精品91| 黄色怎么调成土黄色| 亚洲不卡免费看| 边亲边吃奶的免费视频| 亚洲欧美精品自产自拍| 99热国产这里只有精品6| 黑人巨大精品欧美一区二区蜜桃 | 人人澡人人妻人| 日本爱情动作片www.在线观看| 成人无遮挡网站| 免费看av在线观看网站| 亚洲三级黄色毛片| 国产男女超爽视频在线观看| 亚洲精品久久久久久婷婷小说| 中国三级夫妇交换| 一区二区三区四区激情视频| 精品人妻熟女av久视频| 特大巨黑吊av在线直播| 十八禁高潮呻吟视频 | 亚洲天堂av无毛| 欧美精品国产亚洲| 最近手机中文字幕大全| 三级国产精品片| 亚洲三级黄色毛片| 国产熟女欧美一区二区| 观看av在线不卡| 国产精品一区二区在线不卡| 精品一区二区三区视频在线| 国产亚洲av片在线观看秒播厂| www.av在线官网国产| 寂寞人妻少妇视频99o| 爱豆传媒免费全集在线观看| 国产黄片视频在线免费观看| 麻豆成人午夜福利视频| 99热这里只有精品一区| 尾随美女入室| 欧美精品人与动牲交sv欧美| www.色视频.com| 黄色怎么调成土黄色| 亚洲欧美日韩另类电影网站| 国产色婷婷99| 日本与韩国留学比较| 美女福利国产在线| 两个人免费观看高清视频 | 极品少妇高潮喷水抽搐| 2022亚洲国产成人精品| 六月丁香七月| 国产又色又爽无遮挡免| 亚洲精品日韩在线中文字幕| 亚洲欧洲精品一区二区精品久久久 | 国产精品蜜桃在线观看| 91精品一卡2卡3卡4卡| av.在线天堂| 国产精品福利在线免费观看| av免费在线看不卡| 91精品国产国语对白视频| 特大巨黑吊av在线直播| 亚洲av成人精品一区久久| 黄色毛片三级朝国网站 | 国产午夜精品一二区理论片| 大话2 男鬼变身卡| 国产成人免费观看mmmm| 久久97久久精品| 久久久久久久亚洲中文字幕| 国产一区二区在线观看日韩| 在线免费观看不下载黄p国产| 国产黄片美女视频| 在线观看美女被高潮喷水网站| 乱人伦中国视频| 一级,二级,三级黄色视频| 免费观看性生交大片5| 99热这里只有是精品50| 日韩精品有码人妻一区| 久久人人爽人人片av| 黑丝袜美女国产一区| 国产精品欧美亚洲77777| 亚洲人成网站在线播| 精品人妻偷拍中文字幕| 欧美 亚洲 国产 日韩一| 国产精品99久久久久久久久| 九色成人免费人妻av| 精品99又大又爽又粗少妇毛片| 两个人的视频大全免费| 亚洲高清免费不卡视频| 在线精品无人区一区二区三| 男女国产视频网站| 欧美精品高潮呻吟av久久| 欧美区成人在线视频| 日韩中文字幕视频在线看片| 国产精品久久久久久久久免| 日韩人妻高清精品专区| av卡一久久| 18禁动态无遮挡网站| 亚洲经典国产精华液单| 亚洲精品国产成人久久av| 亚洲欧美一区二区三区国产| 欧美激情极品国产一区二区三区 | 欧美bdsm另类| 国产成人aa在线观看| 波野结衣二区三区在线| 一级毛片久久久久久久久女| 亚洲三级黄色毛片| 十分钟在线观看高清视频www | 中文天堂在线官网| 成人无遮挡网站| 午夜福利,免费看| 亚洲一级一片aⅴ在线观看| 成人二区视频| 成人毛片a级毛片在线播放| 中文字幕av电影在线播放| 中文字幕久久专区| 丰满人妻一区二区三区视频av| 婷婷色av中文字幕| 亚洲精品日韩在线中文字幕| 午夜福利影视在线免费观看| 免费大片黄手机在线观看| 精品久久久久久久久av| tube8黄色片| 搡老乐熟女国产| 99热全是精品| 久热这里只有精品99| 亚洲精品亚洲一区二区| 丰满乱子伦码专区| 性色av一级| 久久这里有精品视频免费| 国产黄色免费在线视频| 我要看日韩黄色一级片| 国产爽快片一区二区三区| 亚洲欧美精品自产自拍| 一级毛片电影观看| 国产无遮挡羞羞视频在线观看| 亚洲在久久综合| 婷婷色综合大香蕉| 国产精品久久久久久av不卡| 乱人伦中国视频| 大话2 男鬼变身卡| 日本猛色少妇xxxxx猛交久久| a 毛片基地| 97超碰精品成人国产| 丁香六月天网| 又爽又黄a免费视频| 亚洲av二区三区四区| 国产成人91sexporn| 精品久久久精品久久久| 91成人精品电影| a级毛片在线看网站| 久久影院123| 老女人水多毛片| 国产一区二区在线观看日韩| 麻豆精品久久久久久蜜桃| 成人毛片60女人毛片免费| 亚洲国产最新在线播放| 精品少妇黑人巨大在线播放| 国产高清国产精品国产三级| 日韩一区二区视频免费看| 国产一区二区在线观看av| 国产无遮挡羞羞视频在线观看| 天天操日日干夜夜撸| 欧美日韩国产mv在线观看视频| 国产在线视频一区二区| 色视频在线一区二区三区| av视频免费观看在线观看| 黑人高潮一二区| 免费看av在线观看网站| 最近中文字幕2019免费版| 亚州av有码| 亚洲综合精品二区| 久久精品国产a三级三级三级| 99视频精品全部免费 在线| 大片电影免费在线观看免费| 狂野欧美白嫩少妇大欣赏| 国产伦在线观看视频一区| 97超视频在线观看视频| 精品亚洲成a人片在线观看| 少妇丰满av| 亚洲av综合色区一区| 日日啪夜夜撸| 国产中年淑女户外野战色| 永久免费av网站大全| 婷婷色综合www| 伦理电影大哥的女人| 亚洲精品国产色婷婷电影| 国产欧美日韩一区二区三区在线 | 男人舔奶头视频| 亚洲精品日本国产第一区| 人体艺术视频欧美日本| 高清不卡的av网站| 久久精品国产亚洲av涩爱| 成人特级av手机在线观看| 夜夜爽夜夜爽视频| 高清av免费在线| 99热这里只有是精品50| 亚洲国产欧美在线一区| 中文字幕久久专区| 国产精品免费大片| 精品一区二区三卡| 亚洲内射少妇av| 黄色欧美视频在线观看| 国产亚洲91精品色在线| 男女无遮挡免费网站观看| 日韩中文字幕视频在线看片| 欧美一级a爱片免费观看看| 亚洲一级一片aⅴ在线观看| 日日啪夜夜爽| 亚洲国产精品专区欧美| 欧美成人午夜免费资源| 日本黄大片高清| 国产一区二区在线观看日韩| 日韩人妻高清精品专区| 激情五月婷婷亚洲| 男人狂女人下面高潮的视频| 黄色日韩在线| 国产男女内射视频| 丝袜在线中文字幕| 中文字幕人妻丝袜制服| 成人特级av手机在线观看| 自线自在国产av| 22中文网久久字幕| 91久久精品国产一区二区成人| 免费看日本二区| 成年av动漫网址| 男女免费视频国产| av在线播放精品| 久久毛片免费看一区二区三区| 久久99精品国语久久久| 99久久精品国产国产毛片| 99国产精品免费福利视频| 日韩中字成人| 亚洲欧美一区二区三区国产| 国产高清不卡午夜福利| 有码 亚洲区| 黑人巨大精品欧美一区二区蜜桃 | 精品一区在线观看国产| 日本91视频免费播放| 国产熟女午夜一区二区三区 | 国产视频内射| 日本免费在线观看一区| 日韩强制内射视频| 日本91视频免费播放| 午夜av观看不卡| av女优亚洲男人天堂| 五月开心婷婷网| 99九九线精品视频在线观看视频| 欧美一级a爱片免费观看看| 日韩制服骚丝袜av| 人妻少妇偷人精品九色| 亚洲国产精品999| 极品人妻少妇av视频| 国产黄色视频一区二区在线观看| 亚洲国产精品999| 爱豆传媒免费全集在线观看| 亚洲成人一二三区av| 中文字幕人妻丝袜制服| 一级二级三级毛片免费看| 日本欧美视频一区| 日韩不卡一区二区三区视频在线| 在线观看免费高清a一片| 精品少妇黑人巨大在线播放| 美女福利国产在线| 欧美日韩一区二区视频在线观看视频在线| 老司机亚洲免费影院| 中文天堂在线官网| 青春草国产在线视频| 国产精品无大码| 女的被弄到高潮叫床怎么办| 精品一区二区免费观看| 亚洲第一av免费看| 激情五月婷婷亚洲| 在线精品无人区一区二区三| 男男h啪啪无遮挡| 欧美区成人在线视频| 99久久精品国产国产毛片| 国产免费一级a男人的天堂| 午夜福利视频精品| 欧美激情国产日韩精品一区| 日本欧美视频一区| 亚洲精品亚洲一区二区| 亚洲色图综合在线观看| 成人午夜精彩视频在线观看| 视频区图区小说| 精品人妻熟女av久视频| 波野结衣二区三区在线| 久久人人爽av亚洲精品天堂| 日本欧美视频一区| 日韩精品有码人妻一区| 99热这里只有是精品在线观看| 免费观看在线日韩| 一级毛片我不卡| 欧美xxⅹ黑人| 久久人人爽人人片av| 少妇人妻一区二区三区视频| 九九爱精品视频在线观看| 国产极品天堂在线| a 毛片基地| 在线天堂最新版资源| 欧美性感艳星| 久久人人爽人人片av| 97超碰精品成人国产| 午夜日本视频在线| 男女边吃奶边做爰视频| 久久久久国产精品人妻一区二区| a级片在线免费高清观看视频| 22中文网久久字幕| √禁漫天堂资源中文www| 最近的中文字幕免费完整| 另类亚洲欧美激情| 日韩一本色道免费dvd| 国产成人免费无遮挡视频| 国产毛片在线视频| 久久午夜综合久久蜜桃| 精品一品国产午夜福利视频| 七月丁香在线播放| 国产精品一二三区在线看| 欧美激情极品国产一区二区三区 | 久久国产乱子免费精品| 国产片特级美女逼逼视频| 久久久久久久久久久免费av| 国产熟女午夜一区二区三区 | 菩萨蛮人人尽说江南好唐韦庄| 天堂俺去俺来也www色官网| 精品人妻偷拍中文字幕| 免费观看的影片在线观看| 久久久国产精品麻豆| 中文字幕av电影在线播放| 国产成人一区二区在线| 国产精品一区二区在线不卡| 国产欧美日韩精品一区二区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品一区蜜桃| 一级a做视频免费观看| 久久亚洲国产成人精品v| 久久6这里有精品| 免费av不卡在线播放| 亚洲国产精品成人久久小说| 亚洲欧美精品专区久久| 久久精品国产a三级三级三级| 久久精品久久久久久噜噜老黄| 国产欧美日韩一区二区三区在线 | 日韩 亚洲 欧美在线| 色94色欧美一区二区| 男女国产视频网站| 久久精品熟女亚洲av麻豆精品| 免费不卡的大黄色大毛片视频在线观看| 三上悠亚av全集在线观看 | 三级国产精品片| 午夜福利网站1000一区二区三区| 国产免费视频播放在线视频| 在线精品无人区一区二区三| 一级片'在线观看视频| 亚洲国产成人一精品久久久| 亚洲电影在线观看av| 永久免费av网站大全| 亚洲国产最新在线播放| 人人妻人人爽人人添夜夜欢视频 | 我要看黄色一级片免费的| 久久久a久久爽久久v久久| 人人妻人人澡人人看| 亚洲精品中文字幕在线视频 | 亚洲av电影在线观看一区二区三区| 久久国产精品男人的天堂亚洲 | 亚洲欧美日韩东京热| 中文字幕av电影在线播放| 你懂的网址亚洲精品在线观看| 97精品久久久久久久久久精品| 日韩精品有码人妻一区| 精品一区二区三卡| 91在线精品国自产拍蜜月| 高清午夜精品一区二区三区| 午夜激情福利司机影院| 久久精品久久久久久噜噜老黄| 成人美女网站在线观看视频| 2018国产大陆天天弄谢| 精品久久久噜噜| 亚洲精品色激情综合| 久久久久久久亚洲中文字幕| 国产极品粉嫩免费观看在线 | 久久国内精品自在自线图片| 国产av精品麻豆| 午夜老司机福利剧场| 国产av精品麻豆| 高清欧美精品videossex| 欧美三级亚洲精品| 欧美日韩av久久| 久久免费观看电影| 一级黄片播放器| 91aial.com中文字幕在线观看| a级毛色黄片| 如日韩欧美国产精品一区二区三区 | 少妇丰满av| 免费黄网站久久成人精品| 亚洲欧洲国产日韩| 国产深夜福利视频在线观看| 精品国产一区二区久久| 丝瓜视频免费看黄片| 99久久精品热视频| 十八禁高潮呻吟视频 | 亚洲国产精品国产精品| 草草在线视频免费看| av线在线观看网站| 国产成人91sexporn| 丁香六月天网| 一本色道久久久久久精品综合| 又黄又爽又刺激的免费视频.| 亚洲av男天堂| 丝袜喷水一区| 九色成人免费人妻av| 久久这里有精品视频免费| 亚洲精品国产成人久久av| 午夜免费男女啪啪视频观看| 色视频在线一区二区三区| 久久鲁丝午夜福利片| 好男人视频免费观看在线| 一级毛片黄色毛片免费观看视频| 日韩亚洲欧美综合| 日韩人妻高清精品专区| 欧美成人午夜免费资源| 一级毛片电影观看| 欧美区成人在线视频| 国产亚洲一区二区精品| 中国三级夫妇交换| 国产探花极品一区二区| 国产精品久久久久久精品电影小说| 精华霜和精华液先用哪个| 极品人妻少妇av视频| 日韩一区二区视频免费看| 热re99久久国产66热| 亚洲激情五月婷婷啪啪| 精品亚洲成a人片在线观看| 国产精品女同一区二区软件| 免费观看av网站的网址| 91精品国产九色| 国产精品免费大片| h日本视频在线播放| 最近2019中文字幕mv第一页| 人人妻人人澡人人看| 欧美变态另类bdsm刘玥| 一级毛片我不卡| 日本与韩国留学比较| 久久国产精品男人的天堂亚洲 | 狂野欧美激情性xxxx在线观看| 午夜视频国产福利| 色视频在线一区二区三区| 国产av一区二区精品久久| 欧美三级亚洲精品| 国产片特级美女逼逼视频| 22中文网久久字幕| 国产69精品久久久久777片| 午夜av观看不卡| 亚洲精品乱码久久久v下载方式| 日本vs欧美在线观看视频 | 免费人妻精品一区二区三区视频| 亚洲av不卡在线观看| 久久久久精品久久久久真实原创| 国产在线一区二区三区精| 国产欧美另类精品又又久久亚洲欧美| 国产色婷婷99| av天堂久久9| 人妻制服诱惑在线中文字幕| 亚洲国产最新在线播放| 在线看a的网站| 少妇 在线观看| 黄片无遮挡物在线观看| 黑人巨大精品欧美一区二区蜜桃 |