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

    Investigation on subsurface vortices within a closed pump intake under different pressure conditions *

    2020-12-16 02:19:02MiaoGuoXuelinTangXiaoqinLiFujunWangXiaoyanShi

    Miao Guo , Xue-lin Tang , Xiao-qin Li , Fu-jun Wang , Xiao-yan Shi

    1. College of Water Resources and Civil Engineering, China Agricultural University, Beijing 100083, China

    2. Beijing Engineering Research Centre of Safety and Energy Saving Technology for Water Supply Network System, China Agricultural University, Beijing 100083, China

    Abstract: An experimental study to investigate various features of subsurface vortices and their evolutions inside a closed pump intake under different pressure conditions at free surface, and to analyze the influence of a pipe bell with different mean axial velocities at the same water level on flow patterns is presented. For all cases, three different flow regimes are closely observed: various side-wall attached, back-wall attached and submerged vortices structures, and their evolutions are subject to the free-surface pressure for a critical flow rate, the number and strength of all these vortex rise as the absolute value of the pressure increases for a higher flow rate,for a lower flow rate, whether the side-wall vortices occur or not is related closely to the pressure magnitude, and there are no vortices for lower absolute pressure values, the number and the strength of the side-wall vortices increase as the absolute pressure values increase. Besides, based on statistical analyses, the side-wall vortex probability on the right side along the stream-wise direction is much higher than that on the opposite side for all operation cases, particularly single vortex probability, and the back-wall vortex probabilities are quite higher with values of more than 90%, meanwhile, the multi-vortices probability is much higher than that of the single-vortex. The captured vortices have been qualitatively classified into different “vortex zones”: no vortex zones and vortex zones including transition, single- and multi-vortex. The Reynolds number and the absolute value pressure are the most important parameters which influence the existence and the number of these vortices, and the increase of the two parameters reduces the no vortex transition zone and causes bigger difference in the back-wall and side-wall vortex zones.

    Key words: Evolution of subsurface vortices, particle-image velocimetry (PIV), vortex probability, closed pump intake

    Introduction

    Closed pump intakes are widely used in various industrial fields like the cooling water systems of nuclear power stations, sodium-cooled fast reactors,pumping stations with great variations in water level,and so on. The very complex flows inside pump intakes commonly contain many coherent structures[1-2], especially near the surrounding area of the pump bell. The subsurface vortices[3-5]influence safe and steady operation of pump stations, which causes noise and vibrations, impeller damage due to cavitation[5], and uneven impeller loadings[6-7]. This not only results in severe deterioration of flow patterns inside closed pump intakes, and even influences pump performances, but it also leads to a significant increase in operational and maintenance costs[8].

    It is often very useful and important to understand the various complex vortices inside pump intakes by means of observations and measurements of physical models. Generally, visual observation technologies,such as the acoustic Doppler velocimeter (ADV)technique, particle-image velocimetry (PIV), high-speed video, and laser Doppler velocimetry (LDV) techniques,are applied in test models to study subsurface vortices and the flow in pump intakes. Sch?fer et al. observed the flow and subsurface vortices in the front of open and covered intake systems by model testing, and the optimization schemes were used for the devices guiding the cost-effective flows in open and covered intake systems and to improve pump performance[9]. Choi et al.[10]examined flow uniformity in a pump model with a single intake and seven intakes to find out the detailed cause of vortex occurrence, and found that the submerged vortex intensity in the vicinity of side walls was higher than those both near back wall and near the upstream region of the bell mouth, and also that the vorticity around the bell-mouth was higher value in the multi-intake model than that in the single-intake. Ahmad et al. obtained the velocity distribution at nine points with the help of an ADV probe in a scaled pump-intake model. Type-2 (dye core) subsurface vortices were noticed and the swirl angle was high in most of the tests without vortex breakers, and the axial velocities at the throat of the bell mouth of four test pumps were almost uniform[11]. Ansar obtained the mean velocity distribution in a pump intake model with and without a cross-flow device using an ADV probe. The results showed that there were nearly symmetric flows with respect to the longitudinal axis of the intake and no reverse flow for the no cross-flow case. However, the reverse flow concentrated at lower depths and close to the right intake sidewall for the cross-flow case.Simultaneously, a large recirculation region existed in the upstream of the suction pipe at higher depths and near the left intake side wall[12]. Rajendran and Patel[13]used PIV to obtain the vortices in a laboratory simple intake model under conditions of small surface tension and negligible free-surface elevation effects. They found that some vortices had relatively-fixed mean positions while others were intermittent or meandering and the subsurface vortex strength depended on the distance between the wall or floor to which the vortex was attached and the intake pipe, and smaller distance promoted stronger vortices. Nagahara et al.[14]studied the effect of submerged vortices on hydraulic forces of a mixed flow pump impeller in a model pump intake, and obtained the strength of the submerged vortex, the tangential velocities, and the diameter of the cavitation area by using a high-speed video camera and a digital image processing system. It was found that the strength of submerged vortices was in the control of the upper stream velocity distortion and the diameter of the cavitation area of the submerged vortex correlated with the circulation of the vortex. Besides, Tahershamsi et al.[15]used a submerged water jet to alleviate or even eliminate the vortices, and verified it by experiments.

    All the above experimental investigations of subsurface vortices focus on the open channel pump intakes,and pay less attention to the flow pattern and subsurface vortices in closed pump intakes. Besides, for the same cases, the evolution of the vortex probabilities and geometric scale of these vortices are not clear, which are significantly affect safe and steady operation in pump stations. It is not conducive to design the appropriate geometry of the closed pump intake, nor can it provide good inflow conditions for the inlet of pump. Therefore,it is worthwhile to investigate such subsurface vortices in closed pump intake including the related kinematic and dynamic mechanism.

    The stereo PIV technique is used to measure the complicated flow patterns for more accurately analyzing in the closed pump intake model. The specific impacts of both vortex evolution and pressure variation on velocity distribution under the pipe bell were investigated in a closed pump intake. This paper also systematically analyzed the various subsurface vortices and the relative vorticity distribution on different selected sections, and the time-dependent evolution inside the closed pump intake under different free-surface pressures. The locations and the strength of different subsurface vortices, the occurrence probability, and the statistical characteristics of various subsurface vortices, as well as vortex probability distribution of various subsurface vortices in the closed pump intake, were also comprehensively analyzed.

    The remaining parts of the present paper are organized as follows. In section of experimental system, the apparatus of the closed pump intake model and the experimental procedure adopted in this study are described. The experimental study focuses three different kinds of flow rates and five pressure conditions in the closed pump intake. In section of results and discussion, the various subsurface vortices and the vorticity distribution on two side walls and back wall, and the related evolution with different pressures of the closed pump intake are analyzed in detail. In section of experimental statistics and vortex probability, the occurrence probability and statistics of various subsurface vortices on selected sections were also carried out. Finally, in section of conclusions, the summary and conclusions of this experimental study of closed pump intake are provided.

    1. Experimental system

    The experimental tests have been carried out in a pump intake facility[16]. The main component of the experimental facility is a parallelepiped shaped intake and pipe bell in polymethyl methacrylate (PMMA).Other components of the experimental facility (i.e.,pump, regulation valve, pressure/vacuum pump, surge tank, flow meter, data processor, synchronizer, data acquisition system, charge-coupled device (CCD)cameras, Lasers and Diversion barriers et al.)complete the loop in a closed circuit.

    1.1 Experimental facility

    Geometrical and physical similarity are basic principles to allow a definite transfer of model test results to a prototype structure. By these means, the Froude number (Fr) of the physical and prototype model is equal to 1 in order to achieve similar flow conditions in the model and the prototype structure[17].

    In Fig. 1 the closed pump intake experimental facility is shown. The test rig is characterized by a closed loop and pump intake with different pressure conditions to obtain a steady state test flow pattern. The four main operational functions, which can be controlled in the experimental loop are: (1) draining and refilling the closed pump intake, (2) controlling the water flow rate, (3) controlling and steadying the pressure inside of the closed pump intake, (4) PIV operations and taking pictures.

    The pressure inside of the closed pump intake can be obtained by vacuum pump (in this experiment,controlled by the motor control system) connected with the closed pump intake and the valve installed under the surge tank. The state of the valve (open/closed position)and the motor control system allows the investigation of several pressure conditions in the closed pump intake. In the present work, five different pressure conditions are considered. The fluid from the pipe bell discharge flows through the flow meter, then is pumped by the pump system, and is controlled by the regulation valve, and finally into the test section. The measurement devices used in this experiment and its accuracy are given in Table 1.

    The closed PMMA intake channel is 800 mm×450 mm×450 mm (Lx×Ly×H) with walls 20 mm thick. The parameters of the single closed pump intake and various subsurface vortices sketch are exhibited in Fig. 2. The details of the geometry of the closed pump intake are list in Table 2.

    1.2 Experimental procedure

    The initial conditions for each case strongly affect the vortices formation. For this reason, the experimentally captured procedure took an initial waiting time to reach a better condition within the closed circuit and pump intake. According to the first experimental campaign, 20 min has been recommended as the minimum waiting time. This particular procedure has been applied by the need to compare consistent experimental data, avoiding initial conditions affecting the results. Flow patterns and various subsurface vortices formation are certainly affected by non-uniform boundary conditions or by geometry shapes; however, these instabilities are unavoidable in any experimental campaign and are not related to the particular geometry. It is need to highlight that the repeatability of the different experimental cases was verified concerning the general validity of the test,replicating each observation with the aim of eliminating statistical errors due to these instabilities as much as possible.

    Fig. 1 (Color online) A schematic description of the present experimental system. The numbers marked in the figure represent(1-Pump system, 2-Regulation valve, 3-Vacuum pump system, 4-Surge tank, 5-Magnetic flow meter, 6-Data processor,7-Synchronizer, 8-Data acquisition system, 9-CCD cameras, 10-Lasers, 11-Test pump intake, 12-Diversion barriers,13-Controlling valve, 14-Pressure sensor, 15-Motor control system)

    Table 1 A list of details of the experimental instrument employed in the present experiment

    Fig. 2 (Color online) The details of closed pump intake

    Table 2 Geometry parameters of closed pump intake

    Particles in the fluid are illuminated by a sheet of light that is pulsed, which are the 8 μm-12 μm diameter glass balloons with good optical capabilities.The size of the interrogation area is 64×64 pixels in the experiment, with a pixel being approximately 0.15 mm. The parameters of PIV such as exposure timetδ, the time step of system ΔT, the number of framesNpic[18]and the number of the total picturesNtalare listed in Table 3. In order to capture the evolutionary process of various vortices in a closed pump intake accurately and comprehensively, 90 images in each successful experimental step were considered sufficient to obtain the properties of the subsurface vortices. The PIV system is combined in the experimental facility shown as Fig. 1.

    Table 3 The paramount parameters of PIV. The given parameter refers to the case 1

    1.3 Experimental conditions

    Three different flow regimes have been investigated: a “l(fā)ower flow rate” (case 1), a “critical flow rate” (case 2), and a “higher flow rate” (case 3).Referencing to Denny’s definition[19], the critical submergence is 260 mm in this experiment, within the flow rate in case 2. Three different conditions with different flow rates,Q, and same submergence water depth,Hf, of the pump bell listed in Table 4 are considered. The detailed pump intake non-dimensional parameters are listed in Table 5. Additionally,the pressure conditions performed in the present experiment are given in Table 6. Two parameters have been changed within the following ranges: the flow rate between 19.78 m3/h to 25.43 m3/h and the pressure between -0.07 MPa to 0 MPa.

    Table 4 The definition of the operating conditions in terms of three cases

    Table 5 Pump intake dimensionless parameters

    Table 6 The pressure conditions of the present experiments

    2. Results and discussion

    2.1 Velocity and vorticity analysis and discussion

    When the local pressure in the closed pump intake was lower than the vapor pressure of the fluid,cavitation occurred and the vortices attached to side walls or back wall of the closed pump intake would be suctioned in the pump and decreased the efficiency of the system. Hence, 3 selected sections (two side walls and the back wall) around the pipe bell were chosen to reproduce the real flow in relative planes through the PIV technique. Figure 3 shows the measurement section of the side walls 0.35Dfrom the physical side wall 1 and the same distance from the real side wall 2.The domain of two side sections presented by the coordinates inxandydirections are [5.45,7.45]and [0,2].

    Fig. 3 (Color online) PIV measurement section of side walls

    The combination of streamline and vorticity is used to identify the vortices in the present experiment.The evolution of the velocity and vorticity on the side walls section from 0 MPa to -0.07 MPa for case 1 is displayed in Fig. 4. Obviously, the number and scale of the vortices on the side wall 1 were almost increasing with the absolute value of the pressure.There were only one or two small scale vortices formed on side wall 1 at the first three pressure conditions, but two bigger scale vortices at least were captured on the same section at higher absolute value of the pressure conditions (p=-0.06 MPa, -0.07 MPa).Inverse, the velocity and vorticity on the side wall 2 were different to those on the side wall 1. No vortex was found on this section at the lower absolute value of pressure conditions (p=0 MPa, -0.02 MPa,-0.04 MPa and -0.06 MPa), and a vortex with weak strength was formed on the side wall 2 at highest absolute value of pressure condition (p=-0.07 Μ Pa) of this experiment.

    Fig. 4 (Color online) The velocity and vorticity in y-direction (ω y) on the side walls section for case 1. (a1)-(e1), (a2)-(e2) refer to the velocity distributions under the conditions with p=0 MPa, -0.02 MPa, -0.04 MPa, -0.06 MPa and -0.07 MPa on the section of side walls 1 and 2, respectively. (f1)-(j1), (f2)-(j2) refer to the ωy distributions under the conditions with p =0 MPa, -0.02 MPa, -0.04 MPa, -0.06 MPa and -0.07 MPa on the section of side wall 1 and 2, respectively

    The evolution of the velocity and vorticity on the side walls section for case 2 at all pressure conditions is shown in Fig. 5. With the constant tendency of case 1, the number and scale of the vortices on side wall 1 were slightly increasing with the increase of the absolute value of the pressure. Compared with case 1,only one distinct vortex was captured on side wall 1 at the first three pressure conditions, and at least two vortices were captured on the same section at higher absolute pressure value conditions (p=-0.06 MPa,-0.07 MPa), which was the common point of velocity evolution between cases 1 and 2 on side wall 1.However, the velocity and vorticity on side wall 2 have different tendencies with those on side wall 1.There was not any distinct vortex that could be found until the absolute value of the pressure was elevated to 0.06 MPa or higher. Between the evolution of velocity and vorticity in case 1 and 2, the remarkable progress is the absolute pressure value of a distinct vortex formed on side wall 2, which dropped from 0.07 MPa to 0.06 MPa. This development was attributed to the higher mean velocity for case 2.

    When the mean axial velocity in the intake pipe rose to 0.9 m/s, the Reynold numberRe=9× 105, the more turbulent flow was performed in the closed pump intake. For case 3, velocity and vorticity on the side walls section evolved from 0 MPa to -0.07 MPa is indicated in Fig. 6. A distinct vortex with increasing strength located in the same domain was captured on side wall 1 at the first three pressure conditions, then,more vortices meandered on this section with fluctuating vorticities at higher absolute pressure values. It is noted that an obvious vortex located at the same domain on side wall 2 was found at all of the pressure conditions.

    Fig. 6 (Color online) The velocity and vorticity in y-direction (ω y) on the side walls section for case 3. (a1)-(e1), (a2)-(e2) refer to the velocity distributions under the conditions with p=0 MPa, -0.02 MPa, -0.04 MPa, -0.06 MPa and -0.07 MPa on the section of side walls 1 and 2, respectively. (f1)-(j1), (f2)-(j2) refer to the ωy distributions under the conditions with p =0 MPa, -0.02 MPa, -0.04 MPa, -0.06 MPa and -0.07 MPa on the section of side wall 1 and 2, respectively

    After analyzing and comparing the velocity and vorticity evolution on the side walls sections for cases 1, 2 and 3, it was concluded that the flow on section of side wall 1 was more turbulent than that of side wall 2,and the vortex/vortices were progressing with the rising of the absolute value of the pressure conditions on section of side wall 2, which became more and more unsteady and meandering.

    The back-wall section was chosen to study the real flow on the back wall of the closed pump intake.Figure 7 shows the measurement section of the back wall 0.35Dfrom the physical one. The domain of selected sections located by theyandzcoordinates directions are [1.25,3.25] and [0,2].

    Fig. 7 (Color online) PIV Measurement section of back wall

    The velocity and vorticity evolution on the back-wall section from 0 MPa to -0.07 MPa for Case 1 is displayed in Fig. 8. Clearly, the number and scale of the vortices on the back wall were almost increasing along with the increase of the absolute pressure value. Specifically, only one distinct vortex on the back wall could be captured at the first two pressure conditions, then, more vortices were found on the same section with the rising of the absolute value of the pressure. It should be noted that the core of the vortex on this section was unsteady and meandered.

    Fig. 8 (Color online) The velocity and vorticity in x-direction(ω x) on the back-wall section for case 1. (a)-(e) refer to the velocity distributions under the conditions with p=0 MPa,-0.02MPa,-0.04MPa,-0.06MPa and-0.07 MPa respectively. (f)-(j) refer to the ωx distributions under the conditions with p=0 MPa, -0.02 MPa,-0.04 MPa, -0.06 MPa and -0.07 MPa respectively

    The velocity and vorticity evolution on the section of the back wall from 0 MPa to -0.07 MPa for case 2 is illustrated in Fig. 9. Similarly, the number of the vortices on the back wall was increasing with the increase of the absolute pressure value for case 2.Specifically, only one distinct vortex on the back wall formed at the first pressure conditions (i.e.,p= 0 MPa ),then more vortices were found on the same section with the rising absolute pressure value. Additionally,the core of the vortex was unsteady and meandered on this section. Remarkably, comparing the velocity and vorticity evolution between cases 1 and 2, the absolute pressure value conditions of multiple vortices formed on the back wall dropped from 0.04 MPa to 0.02 MPa.This change was contributed by the higher mean velocity for case 2.

    Fig. 9 (Color online) The velocity and vorticity in x-direction(ω x) on the back-wall section for case 2. (a)-(e) refer to the velocity distributions under the conditions with p=0MPa,-0.02MPa,-0.04MPa,-0.06MPaand-0.07 MPa respectively. (f)-(j) refer to the ωx distributions under the conditions with p=0 MPa, -0.02 MPa,-0.04 MPa, -0.06 MPa and -0.07 MPa respectively

    When the mean axial velocity in the intake pipe rose to 0.9 m/s, the velocity and vorticity evolution on the section of the back wall from 0 MPa to -0.07 MPa is indicated in Fig. 10 for case 3. The dramatic progress is that multiple vortices were captured on the back wall for case 3 at all pressure conditions, which was distinctly different from the previous cases.Moreover, the core of the vortex was unsteady and meandered on this section.

    Fig. 10 (Color online) The velocity and vorticity in x-direction (ω x) on the back-wall section for case 3. (a)-(e)refer to the velocity distributions under the conditions with p=0 MPa,-0.02 MPa,-0.04 MPa,-0.06 MPa and -0.07 MPa respectively. (f)-( j) refer to the ωx distributions under the conditions with p=0 MPa,-0.02MPa,-0.04MPa,-0.06MPa and-0.07 MPa respectively

    Based on the results provided by this paper on the side walls and back wall sections for all cases under all pressure conditions, it can conclude that the flow patterns on the back wall were more turbulent than those on the sections of the side wall, and the core of the vortex was more unsteady and flexible on the section of back wall and the number of these vortices was more than that on side walls, which was difficult to predict.

    2.2 Experimental statistics and vortex probability

    The vortices were unsteady and meandered in pump intake[20-22]. After examining the PIV experimental results attentively and comprehensively, it was revealed that the various subsurface vortices were developing with the condition in the closed pump intake and they showed some trends in probability and statistics. Hence, it is necessary and meaningful to statistically analyze and correlate the probability of the subsurface vortices in the closed pump intake through experiments, which are contributing to a better understanding and explaining the formation and development of the various subsurface vortices in the closed pump intake.

    Four vortex probability bar graphs are shown in Figs. 11-14 which include experimental results obtained with 90 groups of continuous experimental images. The vortex probability on different measured sections (side-wall 1 and 2, back-wall) was computed and displayed for each case and different pressure conditions.

    Fig. 11 (Color online) Vortex probability on the section of two side walls for case 1 at different pressure condition. (a),(b) refer to the vortex probability on the section of side walls 1 and 2 respectively

    Fig. 12 (Color online) Vortex probability on the section of two side walls for case 2 at different pressure condition. (a),(b) refer to the vortex probability on the section of side walls 1 and 2 respectively

    Fig. 13 (Color online) Vortex probability on the section of two side walls for case 3 at different pressure condition. (a),(b) refer to the vortex probability on the section of side walls 1 and 2 respectively

    Fig. 14 (Color online) Vortex probability on the section of the back wall for all cases at different pressure condition.(a), (b) and (c) refer to the vortex probability on the section of the back wall for case 1, 2 and 3 respectively

    During the tests, all the four vortex probabilities have been computed, merging PIV experimental data obtained with the different sections and cases. The letterη(light color bars) stands for that there was vortex captured on the selected sections for each cases at different pressure conditions, on the other side, the value of 1-ηmeans that there was not any vortex existing on the measured sections (not displayed in Figs. 11-14). The dark color bars ()θare used to represent the probability of a single vortex forming on the PIV captured domain, similarly, at least two different vortices are observed on the selected sections and can be calculated throughη-θ(not displayed in Figs. 11-14). The relative definitions are

    whereNTis the total number of the successful PIV images of each test,NA,NBare the number of the PIV images of the existing vortices in the PIV results,and only one vortex on the selected section at least in each test, respectively.

    The vortex probabilities on the two side wall sections for case 1 and different pressure conditions are illustrated in Fig. 11. On the section of side wall 1,the value of throughηis rising with the increase of the absolute value of the pressure, and reached the maximum value of 88.33% at -0.07 MPa. Inversely,the probability of only one vortex is down to the minimum value of 15% at -0.07 MPa. On the section of side wall 2, the vortex probabilities are lower than the corresponding value on the side wall 1. The value ofηis increasing with the increase of the absolute pressure value, and reached the maximum value of 80%at -0.07MPa. The variation ofθis small on this section, whose value is maintained around 15%.

    The probabilities on the two side wall sections for case 2 and different pressure conditions are illustrated in Fig. 12. On the section of side wall 1, the value ofηis increasing with the increase of the absolute pressure value, reaching the maximum value 98.33% at -0.07 MPa. In the other side, the value ofθis lowest at 5% at -0.07 MPa. It is evident that the probabilities ofη,η-θand 1-η,θhave a clearly inversed tendency in side wall 1. On the section of side wall 2, the vortex probabilities are lower than the corresponding value on the section of side wall 1. The value ofηis increasing with the increase of the absolute pressure value, and reached the maximum value of 68.33% at -0.07MPa.Compared to case 1, the variation ofθis small too,whose value is maintained around 30%.

    The vortex probabilities on the two side wall sections for case 3 and different pressure conditions are shown in Fig. 13. On the section of side wall 1, the value ofηis maintained a high level with small variation. The maximum and minimum ofηare 90%and 61.67% at 0 MPa and -0.04 MPa or 0.07 MPa,respectively. The probability of only one vortexθis down to the minimum value at -0.07 MPa, which is 16.67%. On the section of side wall 2, compared with the previous conditions, the value ofηshowed an obvious trend of rising. The maximum and minimum values ofηwere 93.33% and 28.33%, respectively.Compared to the previous cases, the variation ofθis small as well, whose value is maintained around 28.33%.

    The vortex probabilities on the back wall for all cases and different pressure conditions are shown in Fig. 14. Clearly, the value ofηis quite high compared with its value on the side walls, reaching the maximum value of 100%. These flow patterns indicated that vortices are easily formed and developed on the back wall for all cases. Contrarily,the probability of only one vortexθis relatively small for all cases, especially at pressured conditions,and its minimum value wasθ=5% for case 2 at-0.07 MPa.

    2.3 Experimental statistics and vortex probability distribution

    The vortex probability distribution and merging experimental data obtained with the different pressure conditions have been charted dramatically as well.According to the statistical results of the above relative parts, the core of the vortices is unsteady and meandering, which is driven by the instantaneous conditions of the flow. Moreover, the vorticities of the vortices changed during the test running. It is not clear if this instability is a peculiar characteristic of the vortex phenomena, or if it is only caused by the closed pump intake geometry and flow perturbations introduced by the pump or other elements. Either way,each experimental condition has been classified considering the worst vortex phenomena observed during the experiment.

    A dimensional analysis was performed in order to identify the dimensionless parameters that influence the various subsurface vortices formation and the number of various subsurface vortices in the closed pump intake.

    Taking the mean axial velocity in the pipeup,the streamwise mean velocity in the channeluc,density of fluidρ, the interior diameter of the pipeD, the dynamic viscosity of waterμ(at 20°C), the angular velocity of rotation of the systemΩ(the angular velocity of rotation of the earth in the present paper) and the width of the intakeLyas the repeating variables:

    Different selected sections in the closed pump intake of experiments have been grouped according to the observed subsurface vortex phenomena, different pressure conditions, and the dimensionless parameters.The subsurface vortex probability distribution are built according to the relative experimental results. A preliminary evaluation of the different subsurface vortex phenomena zone boundaries has been done.However, this must be considered only a preliminary and rough estimation of the transition boundaries,resulting from the qualitative classification criteria adopted in the present paper.

    In detail, the relative zones are partitioned according to the following: No vortex zone means the Vortex probabilityα=[0%,30%), transition zone present the vortex probabilityα=[30%,60%)α=[30%, 60%), and the vortex probability of vortex zone is between 60% and 100%. Furthermore, the vortex zone is separated into single and multi-vortices zones, and the vortex probability of the single vortex zone isβ=(50%,100%], the multi-vortices zone equals that of the vortex probabilityβ=[0%,50%).

    On the basis of a large number of experimental data, combined with dimensionless analysis, the vortex probability distribution is obtained, furthermore, which is also based on rich experimental experience accumulated during the experiment. The vortex probability distribution of the side wall 1 section is illustrated in Fig. 15. The transition zone includes only one case: The Reynolds number is 7×104and the corresponding pressure is 0 MPa. The single vortex zone contains three cases: The Reynolds number is 7×104and the corresponding pressure is-0.02 MPa and -0.04 MPa, the Reynolds number is 9×104and the corresponding pressure is 0 MPa. The rest of the cases in the present experiments are all classified in the multi-vortices zone. Moreover, in this map, a general agreement of the experimental results considering the influence of the variation of the pressure and dimensionless number can still be noticed. The no vortex transition zones are located the area with low pressure, along with the values of the dimensionless number. At first, with the raising of theReand pressure values, the no vortex zone turns into the transition zone, then, the vortex zone is divided into single and multi-vortices zones, and they even appeared in the high pressure and the dimensionless number domains. For the single vortex zone on the side wall 1 for largeRecase, the possible explanation is that the flow near the side wall 1 to gather and then turn back to form a single vortex with clear flow pattern after the impact on the back wall.The same mechanism can also explain the related flow patterns on the section of back wall.

    The vortex probability distribution of the side wall 2 section is demonstrated in Fig. 16. The no vortex zone includes four cases: The Reynolds number is 7×104and the corresponding pressure is 0 MPa, -0.02 MPa and -0.04 MPa, the Reynolds number is 9×104and the corresponding pressure is 0 MPa. The transition zone contains four cases else:the Reynolds number is 7×104and the corresponding pressure is -0.06 MPa, the Reynolds number is 8×104and the corresponding pressure is 0 MPa, -0.02 MPa and -0.04 MPa. The rest of the cases in this section are all classified in the multi-vortices zone. The relative high area of no vortex zone and transition zone indicate that vortex is not a common phenomenon in this section. Moreover, it is need to high-light that there is no unambiguous single vortex zone on the measurement section, which means that there is multi- or no distinct vortex on this measurement section, and the flow patterns on the same section is quite complicated.

    Fig. 15 (Color online) P- Re vortex probability distribution on the section of side wall 1

    Fig. 16 (Color online) P- Re vortex probability distribution on the section of side wall 2

    The vortex probability distribution of the back-wall section is demonstrated in Fig. 17. The no vortex zone and transition zone on this section are not included in any case in the present experiments, which are charted qualitatively by experience and the related fluid mechanics principle. Then, the single vortex zone includes only one case: The Reynolds number is 9×104, the corresponding pressure is 0 MPa, and this vortex pattern may result from the strong circulation and the flow impacting on this wall. The rest of the cases in the present section are all classified in the multi-vortices zone, which intimates that the multi-vortices structures on the back-wall section are very common.

    Fig. 17 (Color online) P- Re vortex probability distribution on the section of back wall

    3. Conclusions

    An experimental study was undertaken in order to investigate the working conditions of a closed pump intake, with the same water level, made up of a vertical intake pipe with variable flow rates and pressure conditions.

    For different cases and pressure conditions, the number of vortices on the side wall sections were almost increasing with the increase of the absolute pressure value and the strength of the vortex fluctuated from 0 MPa to -0.07 MPa. Additionally, the scale of the vortices on the side wall 1 were almost increasing with the increase of the absolute pressure value. The remarkable progress in the evolution of vortices on the side wall sections is the pressure condition of a distinct vortex formed on side wall 2, which dropped with the rising of the absolute pressure value. This progress is attributed to the higher mean velocity for case 2 and indicated that the vortices were formed more easily on side wall 1 than 2 at the same time. Moreover, the evolution of the velocity and vorticity on the back-wall section from 0 MPa to -0.07 MPa for all cases had a similar and stronger tendency compared to the selected vertical planes sections under the pipe bell and side walls.

    The vortex probability on different measured sections were presented for all cases and pressure conditions. The variation ofθwas small in all cases,whose value maintained around 15%, 30% and 28.33%for cases 1, 2 and 3, respectively. On the back-wall section, the value ofηwas quite high compared with that of the side wall sections, and reached the maximum value of 100%. This indicated that vortices are formed and developed on the back wall easily for all cases. The probability of only one vortexθwas relatively small for all cases.

    The captured vortices have been qualitatively classified in different “vortex zones”, and the most important parameters that influenced the physical phenomenon were identified through a dimensionless analysis. The results showed that the no vortex transition zone was expanding with the rising of theReand the absolute pressure value, following the order of the back wall, and side walls 1 and 2, which were quite well distinguished. Moreover, the single vortex transition zone was expanding with the rising ofReand the absolute pressure value, following the order of the back wall, and side walls 1 and 2, but was not recognized as easily as it was in the previous zones. Finally, multiple vortices zones were divided after being analyzed and counted in the experimental results. It was found that the no vortex transition zone is located in the area with low pressure and values of the dimensionless number. Then, with the rising of theReand the absolute pressure value, the single vortex transition zones were separated, and they even appeared in the high pressure and the dimensionless number domains.

    Experiments showed that the vortical features and their evolution in a closed pump intake are quite complicated, but follow some regularity. Therefore,the present paper suggested that the key sections and regions for the suppression of subsurface vortices should be carried out according to the vortex probability distribution. The results of the paper are of guiding significance for optimizing geometry of the closed pumping stations and designing inlet flow conditions of the next installations. The rest of the results of closed pump intakes at pressured conditions will be shown in future work.

    少妇猛男粗大的猛烈进出视频| 久久毛片免费看一区二区三区| 91成人精品电影| 极品少妇高潮喷水抽搐| videos熟女内射| av线在线观看网站| 晚上一个人看的免费电影| 免费看av在线观看网站| 亚洲精品456在线播放app| 中文字幕亚洲精品专区| 国产爽快片一区二区三区| 国产在线一区二区三区精| 欧美97在线视频| 亚洲欧美一区二区三区国产| 三级国产精品片| 国产成人精品一,二区| 黄色配什么色好看| 91午夜精品亚洲一区二区三区| 日韩精品有码人妻一区| 天美传媒精品一区二区| 国产69精品久久久久777片| 一区二区日韩欧美中文字幕 | 久久久久久久精品精品| av电影中文网址| 麻豆精品久久久久久蜜桃| 久久精品国产亚洲av天美| 久热这里只有精品99| 久久国产精品大桥未久av| 久久99一区二区三区| 观看av在线不卡| 精品国产一区二区久久| 99热网站在线观看| 青春草亚洲视频在线观看| 亚洲婷婷狠狠爱综合网| 亚洲欧美一区二区三区国产| 亚洲精品一区蜜桃| 一区二区日韩欧美中文字幕 | 亚洲精品中文字幕在线视频| 性色av一级| 在线观看www视频免费| 三上悠亚av全集在线观看| 蜜桃国产av成人99| 中文字幕精品免费在线观看视频 | 色哟哟·www| 男女边吃奶边做爰视频| 国产视频首页在线观看| 青青草视频在线视频观看| 黄色毛片三级朝国网站| 美女cb高潮喷水在线观看| 丝袜脚勾引网站| 性色av一级| 十八禁高潮呻吟视频| 国产免费又黄又爽又色| 国产精品不卡视频一区二区| 伦理电影大哥的女人| 成人毛片a级毛片在线播放| av免费观看日本| 亚洲综合色网址| 婷婷色综合www| 国产日韩欧美亚洲二区| 亚洲经典国产精华液单| 日本vs欧美在线观看视频| 美女福利国产在线| av黄色大香蕉| 两个人的视频大全免费| 男女高潮啪啪啪动态图| 另类精品久久| 国产白丝娇喘喷水9色精品| 3wmmmm亚洲av在线观看| 国产黄片视频在线免费观看| 嘟嘟电影网在线观看| 欧美人与善性xxx| 国产 一区精品| 人人澡人人妻人| 一边摸一边做爽爽视频免费| 精品少妇黑人巨大在线播放| 99热国产这里只有精品6| 国产高清有码在线观看视频| 99久久人妻综合| 久久女婷五月综合色啪小说| 嘟嘟电影网在线观看| 18禁观看日本| 伊人久久精品亚洲午夜| 草草在线视频免费看| 少妇猛男粗大的猛烈进出视频| 纯流量卡能插随身wifi吗| 欧美激情极品国产一区二区三区 | 成人18禁高潮啪啪吃奶动态图 | 三级国产精品欧美在线观看| 一级毛片电影观看| 91精品伊人久久大香线蕉| 国产一区二区在线观看av| 欧美97在线视频| 亚洲,欧美,日韩| 久久精品国产亚洲av涩爱| 亚洲丝袜综合中文字幕| 日本91视频免费播放| 少妇精品久久久久久久| 欧美激情 高清一区二区三区| 99re6热这里在线精品视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲精品视频女| 街头女战士在线观看网站| 一级毛片 在线播放| 国产亚洲午夜精品一区二区久久| 人人澡人人妻人| 最近手机中文字幕大全| 亚洲五月色婷婷综合| 丁香六月天网| 亚洲精品一二三| 99视频精品全部免费 在线| 一区二区三区精品91| 各种免费的搞黄视频| 天堂8中文在线网| 久久久久国产精品人妻一区二区| 看免费成人av毛片| 18禁裸乳无遮挡动漫免费视频| 午夜久久久在线观看| 丝袜脚勾引网站| 久久久欧美国产精品| 婷婷色麻豆天堂久久| 天堂8中文在线网| 亚洲av免费高清在线观看| 欧美丝袜亚洲另类| 美女cb高潮喷水在线观看| 午夜精品国产一区二区电影| 大又大粗又爽又黄少妇毛片口| 久热这里只有精品99| 在线播放无遮挡| 女人久久www免费人成看片| 日本免费在线观看一区| av播播在线观看一区| 爱豆传媒免费全集在线观看| 亚洲色图综合在线观看| 亚洲国产精品成人久久小说| 熟妇人妻不卡中文字幕| 一区二区日韩欧美中文字幕 | 国产免费又黄又爽又色| 少妇高潮的动态图| 亚洲激情五月婷婷啪啪| 十分钟在线观看高清视频www| av电影中文网址| 亚洲不卡免费看| 夫妻性生交免费视频一级片| 精品一区在线观看国产| 亚洲人成网站在线播| 欧美最新免费一区二区三区| 男女边吃奶边做爰视频| 亚洲美女搞黄在线观看| 欧美精品人与动牲交sv欧美| 日日摸夜夜添夜夜爱| 新久久久久国产一级毛片| 三级国产精品欧美在线观看| av在线播放精品| 久久精品国产鲁丝片午夜精品| 亚洲第一区二区三区不卡| 男女高潮啪啪啪动态图| 国产乱人偷精品视频| 男女免费视频国产| 亚洲欧洲精品一区二区精品久久久 | 五月天丁香电影| 久久这里有精品视频免费| 欧美日韩视频精品一区| 我要看黄色一级片免费的| 成人综合一区亚洲| 美女主播在线视频| a 毛片基地| 亚洲性久久影院| 久久久久精品久久久久真实原创| 国产精品一区www在线观看| 国产毛片在线视频| 亚州av有码| 国产黄色视频一区二区在线观看| 人妻一区二区av| 国产午夜精品久久久久久一区二区三区| 国产成人免费无遮挡视频| 18禁观看日本| 91午夜精品亚洲一区二区三区| 黄色怎么调成土黄色| 伊人久久国产一区二区| 亚洲精品,欧美精品| 免费观看无遮挡的男女| 自线自在国产av| 寂寞人妻少妇视频99o| 亚洲精品国产色婷婷电影| 三上悠亚av全集在线观看| 夜夜看夜夜爽夜夜摸| 男女高潮啪啪啪动态图| 晚上一个人看的免费电影| 精品国产露脸久久av麻豆| 欧美国产精品一级二级三级| 欧美成人午夜免费资源| 最新的欧美精品一区二区| 日日摸夜夜添夜夜爱| 免费少妇av软件| av在线观看视频网站免费| 国产精品久久久久久av不卡| 蜜桃在线观看..| 免费播放大片免费观看视频在线观看| 午夜激情福利司机影院| 国产熟女欧美一区二区| 天美传媒精品一区二区| 国产日韩欧美在线精品| 26uuu在线亚洲综合色| 亚洲av免费高清在线观看| 日韩成人伦理影院| 亚洲伊人久久精品综合| 日韩 亚洲 欧美在线| 国产精品一区二区三区四区免费观看| 色哟哟·www| 麻豆精品久久久久久蜜桃| 熟女av电影| 插逼视频在线观看| 国产精品 国内视频| 日本av免费视频播放| 在线观看国产h片| 99久国产av精品国产电影| 中文字幕亚洲精品专区| 岛国毛片在线播放| 母亲3免费完整高清在线观看 | 欧美激情极品国产一区二区三区 | 亚洲欧美中文字幕日韩二区| 亚洲av二区三区四区| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品久久久久久久性| 精品国产露脸久久av麻豆| 美女xxoo啪啪120秒动态图| 精品亚洲成国产av| 国产一区二区在线观看日韩| 欧美 亚洲 国产 日韩一| 尾随美女入室| 一本—道久久a久久精品蜜桃钙片| 国产亚洲一区二区精品| 22中文网久久字幕| 在现免费观看毛片| 欧美激情极品国产一区二区三区 | 午夜激情av网站| 久久精品人人爽人人爽视色| 欧美日韩视频高清一区二区三区二| 国产成人精品一,二区| 久久99蜜桃精品久久| 国产精品熟女久久久久浪| 青春草视频在线免费观看| 91久久精品电影网| 日韩欧美精品免费久久| 天天影视国产精品| 国产熟女欧美一区二区| av有码第一页| 丝袜脚勾引网站| 少妇的逼水好多| 青春草视频在线免费观看| 不卡视频在线观看欧美| 国产精品无大码| 一级片'在线观看视频| 曰老女人黄片| 国产精品久久久久久精品古装| 午夜免费男女啪啪视频观看| 久久人妻熟女aⅴ| 国产国语露脸激情在线看| xxx大片免费视频| 亚洲精品久久午夜乱码| 久久人人爽人人片av| 日日啪夜夜爽| 蜜桃在线观看..| 交换朋友夫妻互换小说| 亚洲欧洲日产国产| 最近的中文字幕免费完整| 精品久久久噜噜| 久久热精品热| 欧美性感艳星| 中文字幕制服av| 人体艺术视频欧美日本| 我的女老师完整版在线观看| 午夜激情av网站| 少妇被粗大的猛进出69影院 | 99热这里只有精品一区| 夫妻性生交免费视频一级片| 免费不卡的大黄色大毛片视频在线观看| 成人二区视频| 中文欧美无线码| 2021少妇久久久久久久久久久| 亚洲精品美女久久av网站| 又粗又硬又长又爽又黄的视频| 欧美精品亚洲一区二区| 亚洲高清免费不卡视频| 美女主播在线视频| 亚洲色图 男人天堂 中文字幕 | 在线精品无人区一区二区三| av免费观看日本| 制服诱惑二区| 蜜桃在线观看..| 99在线人妻在线中文字幕 | 国产精品久久久久久精品古装| 国产深夜福利视频在线观看| 精品国产乱码久久久久久男人| 别揉我奶头~嗯~啊~动态视频| 法律面前人人平等表现在哪些方面| 建设人人有责人人尽责人人享有的| 久久精品国产亚洲av香蕉五月 | 国产男女超爽视频在线观看| 亚洲av国产av综合av卡| 亚洲人成伊人成综合网2020| 亚洲一区中文字幕在线| 国产淫语在线视频| 99riav亚洲国产免费| 高潮久久久久久久久久久不卡| 他把我摸到了高潮在线观看 | 国产深夜福利视频在线观看| 免费女性裸体啪啪无遮挡网站| 亚洲av成人一区二区三| 1024香蕉在线观看| www.熟女人妻精品国产| www.自偷自拍.com| 国产极品粉嫩免费观看在线| 国产精品久久久久成人av| 精品久久久久久久毛片微露脸| 国产亚洲欧美精品永久| 久久国产精品男人的天堂亚洲| 亚洲一区二区三区欧美精品| 国产精品1区2区在线观看. | 国产精品九九99| 成人国语在线视频| 男人舔女人的私密视频| 午夜福利在线观看吧| 亚洲一区中文字幕在线| 成人免费观看视频高清| 免费在线观看日本一区| 亚洲精品久久成人aⅴ小说| 欧美日韩福利视频一区二区| 夫妻午夜视频| 露出奶头的视频| 老汉色∧v一级毛片| 精品亚洲成国产av| 亚洲精品国产精品久久久不卡| 成在线人永久免费视频| 免费观看av网站的网址| a在线观看视频网站| 老司机午夜福利在线观看视频 | 免费久久久久久久精品成人欧美视频| 9色porny在线观看| 一级黄色大片毛片| 捣出白浆h1v1| 黄片大片在线免费观看| 亚洲情色 制服丝袜| 女人被躁到高潮嗷嗷叫费观| 十八禁高潮呻吟视频| 国产成人精品在线电影| 老司机午夜福利在线观看视频 | 国产在视频线精品| 亚洲国产欧美日韩在线播放| 国产精品国产av在线观看| 69av精品久久久久久 | 狠狠婷婷综合久久久久久88av| 亚洲视频免费观看视频| 中文字幕另类日韩欧美亚洲嫩草| 欧美亚洲 丝袜 人妻 在线| 可以免费在线观看a视频的电影网站| 人人澡人人妻人| 又黄又粗又硬又大视频| 又紧又爽又黄一区二区| 国产亚洲av高清不卡| 精品一品国产午夜福利视频| 国产精品 国内视频| 免费在线观看黄色视频的| 久9热在线精品视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲欧洲日产国产| 色综合婷婷激情| www.999成人在线观看| 男人操女人黄网站| 黑人操中国人逼视频| 免费黄频网站在线观看国产| 99久久国产精品久久久| 亚洲三区欧美一区| 久久性视频一级片| 欧美日韩av久久| 精品福利永久在线观看| 夜夜骑夜夜射夜夜干| 免费观看av网站的网址| 色尼玛亚洲综合影院| 男女无遮挡免费网站观看| tocl精华| 午夜91福利影院| av在线播放免费不卡| 午夜久久久在线观看| 丰满饥渴人妻一区二区三| 欧美成狂野欧美在线观看| 怎么达到女性高潮| 国产成人精品在线电影| 搡老熟女国产l中国老女人| 人妻一区二区av| 精品少妇久久久久久888优播| 亚洲国产看品久久| av一本久久久久| 99久久人妻综合| e午夜精品久久久久久久| 麻豆乱淫一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 777久久人妻少妇嫩草av网站| 国产精品久久久久久人妻精品电影 | 欧美精品亚洲一区二区| 国产精品二区激情视频| 后天国语完整版免费观看| 99国产精品99久久久久| 一级黄色大片毛片| 高清欧美精品videossex| 国产精品久久久av美女十八| 亚洲av美国av| 欧美黄色淫秽网站| 午夜福利影视在线免费观看| 精品亚洲成a人片在线观看| 免费人妻精品一区二区三区视频| 成人18禁高潮啪啪吃奶动态图| 伊人久久大香线蕉亚洲五| 国产精品自产拍在线观看55亚洲 | 少妇猛男粗大的猛烈进出视频| 老司机深夜福利视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 他把我摸到了高潮在线观看 | 久久影院123| 日韩一区二区三区影片| 国产精品久久久人人做人人爽| 午夜久久久在线观看| 久热这里只有精品99| 老司机福利观看| 一区福利在线观看| 久久久久久亚洲精品国产蜜桃av| 午夜福利乱码中文字幕| 久久久精品国产亚洲av高清涩受| 亚洲第一av免费看| 国产高清videossex| 一个人免费看片子| 精品久久蜜臀av无| av天堂在线播放| 亚洲国产欧美在线一区| 嫩草影视91久久| 久久精品国产99精品国产亚洲性色 | 女警被强在线播放| 一本久久精品| h视频一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 久久中文字幕人妻熟女| 亚洲色图综合在线观看| 啦啦啦 在线观看视频| 制服人妻中文乱码| 国产一区二区激情短视频| 亚洲欧洲精品一区二区精品久久久| 精品少妇黑人巨大在线播放| 99热国产这里只有精品6| 午夜免费成人在线视频| 国产福利在线免费观看视频| 日本黄色视频三级网站网址 | av超薄肉色丝袜交足视频| 国产av国产精品国产| 国产精品av久久久久免费| 久久热在线av| 精品人妻在线不人妻| 老司机午夜福利在线观看视频 | 视频区欧美日本亚洲| av网站在线播放免费| 午夜福利一区二区在线看| tocl精华| 精品亚洲成国产av| 视频在线观看一区二区三区| 激情视频va一区二区三区| bbb黄色大片| 久久婷婷成人综合色麻豆| 日本撒尿小便嘘嘘汇集6| 黄色 视频免费看| 精品国产一区二区久久| 亚洲一区中文字幕在线| 制服诱惑二区| 国产在线一区二区三区精| 久9热在线精品视频| 91成人精品电影| 久久久精品国产亚洲av高清涩受| 国产又爽黄色视频| 国产黄频视频在线观看| 韩国精品一区二区三区| a级毛片在线看网站| 超碰97精品在线观看| 精品免费久久久久久久清纯 | 超碰成人久久| 视频在线观看一区二区三区| 9热在线视频观看99| 后天国语完整版免费观看| 在线观看66精品国产| 欧美日韩福利视频一区二区| 国产熟女午夜一区二区三区| 亚洲av成人不卡在线观看播放网| 久久九九热精品免费| 国产av国产精品国产| 一夜夜www| 一本综合久久免费| 日韩视频在线欧美| 国产精品成人在线| 多毛熟女@视频| 高清在线国产一区| 国产视频一区二区在线看| 久久天堂一区二区三区四区| 久久精品国产99精品国产亚洲性色 | 亚洲免费av在线视频| 在线播放国产精品三级| 免费少妇av软件| 18禁黄网站禁片午夜丰满| 伊人久久大香线蕉亚洲五| 高清毛片免费观看视频网站 | 男女下面插进去视频免费观看| 国产片内射在线| 日本黄色视频三级网站网址 | 精品国内亚洲2022精品成人 | 亚洲国产中文字幕在线视频| av天堂久久9| 久久毛片免费看一区二区三区| 国产精品偷伦视频观看了| 国产男女内射视频| 国产精品香港三级国产av潘金莲| 免费在线观看黄色视频的| 欧美日韩亚洲高清精品| 18禁美女被吸乳视频| 777米奇影视久久| 国产亚洲一区二区精品| 久久热在线av| 欧美在线黄色| 国产成人影院久久av| 国产精品欧美亚洲77777| 久久精品熟女亚洲av麻豆精品| 国产成人av教育| 伊人久久大香线蕉亚洲五| 免费人妻精品一区二区三区视频| 欧美成人午夜精品| 黄色视频,在线免费观看| 久久精品国产综合久久久| 女人高潮潮喷娇喘18禁视频| 操美女的视频在线观看| 欧美精品亚洲一区二区| 午夜福利乱码中文字幕| 夫妻午夜视频| 久久精品国产a三级三级三级| 日韩一卡2卡3卡4卡2021年| 女同久久另类99精品国产91| 久久精品91无色码中文字幕| 亚洲av片天天在线观看| 成在线人永久免费视频| 黑人欧美特级aaaaaa片| 亚洲精品美女久久久久99蜜臀| 午夜激情久久久久久久| 久久亚洲精品不卡| 日韩一卡2卡3卡4卡2021年| 久久久精品免费免费高清| 欧美成人免费av一区二区三区 | 亚洲三区欧美一区| 最近最新中文字幕大全免费视频| 成年版毛片免费区| 精品人妻在线不人妻| 99国产极品粉嫩在线观看| 一二三四社区在线视频社区8| 久久久久久久久免费视频了| 日韩免费av在线播放| 亚洲人成伊人成综合网2020| 欧美日韩av久久| av超薄肉色丝袜交足视频| 亚洲精品国产一区二区精华液| 搡老熟女国产l中国老女人| 老汉色∧v一级毛片| 亚洲av日韩在线播放| 国产成人啪精品午夜网站| 日韩免费av在线播放| 国产欧美日韩精品亚洲av| 一级毛片女人18水好多| 在线观看人妻少妇| 99热网站在线观看| 人人妻人人添人人爽欧美一区卜| 国产日韩欧美视频二区| 国产午夜精品久久久久久| 亚洲人成伊人成综合网2020| 国产精品国产高清国产av | 久久久国产欧美日韩av| 女人精品久久久久毛片| 超碰97精品在线观看| 女性生殖器流出的白浆| 美女福利国产在线| 欧美黑人精品巨大| 日韩中文字幕视频在线看片| 熟女少妇亚洲综合色aaa.| 12—13女人毛片做爰片一| 亚洲一卡2卡3卡4卡5卡精品中文| 超碰成人久久| 日韩人妻精品一区2区三区| 日日摸夜夜添夜夜添小说| 国产精品自产拍在线观看55亚洲 | 国产精品98久久久久久宅男小说| 无遮挡黄片免费观看| 高清毛片免费观看视频网站 | 国产成人啪精品午夜网站| 男人舔女人的私密视频| 肉色欧美久久久久久久蜜桃| 在线播放国产精品三级| 美女国产高潮福利片在线看| 日韩 欧美 亚洲 中文字幕| 肉色欧美久久久久久久蜜桃| 50天的宝宝边吃奶边哭怎么回事| 亚洲 国产 在线| 亚洲国产精品一区二区三区在线| 黑人欧美特级aaaaaa片| 日韩免费高清中文字幕av| 热99久久久久精品小说推荐| 波多野结衣av一区二区av| 国产区一区二久久| 亚洲精品一卡2卡三卡4卡5卡| 久久九九热精品免费| 老司机在亚洲福利影院| 狠狠狠狠99中文字幕| 精品国产国语对白av| 18禁美女被吸乳视频|