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

    Flow structures and hydrodynamics of unsteady cavitating flows around hydrofoil at various angles of attack *

    2018-05-14 01:42:58DongmeiJu劇冬梅ChangleXiang項(xiàng)昌樂(lè)ZhiyingWang王志英JunLi李軍NanxiXiao肖南溪
    關(guān)鍵詞:昌樂(lè)南溪李軍

    Dong-mei Ju (劇冬梅), Chang-le Xiang (項(xiàng)昌樂(lè)) , Zhi-ying Wang (王志英) , Jun Li (李軍) ,Nan-xi Xiao (肖南溪)

    1. State Key Laboratory of Vehicle Transmission, Beijing Institute of Technology, Beijing 100081, China

    2. Ordnance Science and Research Academy of China, Beijing 100089, China

    Introduction

    It is well known that unsteady cavitating flow in the turbomachineries and the controlled surface would lead to many problems such as sudden changes of loads, pressure pulsations, vibrations and noises[1].Based on the physical characteristics, four primary forms of cavitation can be described: the inception cavitation[2], as well as the sheet cavitation, the cloud cavitation, and the supercavitation with sustained,transient cavities[3]. The unsteady structure of the cloud cavitation can occur not only in a spatially or temporally varying flow, but also in cases with a stationary hydrofoil and a steady and uniform inlet flow[4].

    Two main classes of instabilities of unsteady cavitating flows were distinguished[5], namely, the intrinsic instabilities and the system instabilities,according to the origin of the unsteadiness. An intrinsic instability might be originated in the cavity itself, e.g., the collapse of individual bubbles through the formation of microjets, or the unsteady breakup and shedding of cloud cavities through the formation and the development of the re-entrant flow. On the other hand, the system instability is caused by dynamic interplay between the cavity and the rest of the system, e.g., the interaction with the flow inlet,outlet, or the pump casing wall, the interaction between rotor and stator.

    Among various cavitation phenomena, special attentions are paid to the unsteady characteristics of the cloud cavitation, with vibrations, noises and erosions[6-9]. The cloud cavitation typically starts with the formation of the re-entrant flow at the closure region of the attached cavity. When the re-entrant flow reaches the cavity leading edge, the cavity breaks,lifts off, and then shed downstream, we have the cloud cavitation[10]. Many experiments, conducted in various testing facilities around the world, with very different hydraulic impedances, found cavity shedding fre-quencies with similar Strouhal numbers[11,12].

    To improve the understanding of the complex structures of cavitating flows, various experimental studies were conducted. Gopalan and Katz[13]used a high-speed camera to measure the details of flow structures of the sheet cavitation. They found that the formation and the collapse of the vapor bubbles in the cavity closure region is an important source of the vorticity production. They also found that the size of the cavity had significant impacts on the turbulent intensity and the momentum thickness of the downstream boundary layer. Arndt et al.[14]investigated the instability of the partial cavitation around the NACA 0015 hydrofoil by experimental and numerical methods. They showed that the transition from the sheet cavitation to the cloud cavitation could induce significant fluctuations in lift, thrust and torque. Ausoni et al.[15]applied the PIV to study the effects of the cavitation and the fluid-structure interaction on the vortex generation mechanism at the trailing edge of a hydrofoil in a moderate Reynolds number (Re=2.5× 104-6.5× 104) flow. They demonstrated that the vortex strength increased linearly with the transverse velocity at the foil trailing edge, hence, the cavitation in the vortex street could not be considered as a pas- sive agent for the turbulent wake flow.Aeschlimann et al.[16]used the PIV-LIF to conduct the velocity measurement in a 2-D cavitating shear layer flow, the experimental results showed that the cavitation played an important role in the interaction between the large and small scale flows.

    Although many studies explored the effects of the cavitation through experimental studies, the complex structure and the hydrodynamic characteristics of the unsteady cloud cavitation are still not well understood. The present study is to obtain a comprehensive set of flow field data and video recordings to reveal the structures and the hydrodynamics of the unsteady cavitation with various combinations of angle of attack α and cavitation number σ, to offer a basis for future numerical validation studies. A high-speed video camera is used to observe the time-evolution of the cavitation structure on a Clark-Y hydrofoil, and the PIV system is applied to obtain the velocity and vorticity fields. The temporal and spatial variations of the flow velocity and vorticity profiles at selected chord wise locations along and aft of the foil are shown to better understand the unsteady cloud cavitating flow structures. Meanwhile, the drag and lift coefficients in different cavitation regimes are measured with a dynamic measurement device.

    1. Experimental apparatus and instrumentation

    In order to better investigate the cavity structure,experimental studies are conducted in a closed-loop cavitation tunnel. An axial flow pump, to drive the flow in the tunnel, is located about 5 m below the test section. A tank with a volume of 5 m3is placed upstream the test section to filter out free stream bubbles. The top of the tank is connected to a vacuum pump for controlling the pressure in the cavitation tunnel. Between the test section and the tank, a corner vane and a straightening vane are applied to reduce the turbulence level of the flow. A schematic diagram of the experimental setup is illustrated in Fig.1. The cavitation phenomena are registered by a high-speed digital camera (HG-LE, by Redlake). In order to maintain desirable spatial resolutions, depending on the focus of the investigation, 1000 fps and 4000 fps are used in this study. The flow field at the mid-span of the foil is illuminated by a continuous laser sheet from the bottom wall.length of the hydrofoil isc=0.07 m. The position of the hydrofoil inside the test section is shown in Fig.2.The test section is 10cin length, 2.7cin height,and 1.0cin width. The flow structures could be

    Fig.1 Schematic diagram of the experimental setup

    Fig.2 Sketch of the foil’s position in the test section

    A two-component high-speed PIV system developed by Dantec is used, as shown in Fig.1. The light source is a dual head, Nd: YAG laser, where the beam is expanded into a 1 mm wide sheet. For the actual data acquisition, a 2000 Hz recording rate at 1024×512 pixels is chosen. The commercial PIV-soft ware,the Dynamic Studio, is applied to process the velocity vector field with an interrogation area of 32×32 pixels and 50% overlap in general. In the cavitating flows,the velocity information inside the cavity is captured using the vapor bubbles as “tracer particles”[17].

    A Clark-Y hydrofoil with different angles of attack is adopted in the present study, and the chord observed at the top, bottom and sides. The cavitation tunnel is capable of generating a free stream with velocities ranging from2 m/s-15 m/s, the minimum cavitation number could be reduced to 0.3. The hydrofoil is mounted horizontally in the tunnel test section at α=0°, 5° and 8° , respectively. The foil is clamped to the back wall and the suction side of the foil is mounted toward the bottom of the test section for the convenience of viewing the cavity structures.

    The data correspond to an inflow speed of 10 m/s and a cavitation number of σ=2.0. To avoid the wall boundary layer effects, the data are collected along the mid-span of the hydrofoil where the axial velocity is nearly uniform and the turbulent intensity is less than 2%. In this experimental research, the reference velocityU∞is fixed at 10 m/s, the pressure is adjusted to vary the cavitation number σ=(p∞-pv)/(0.5ρlU∞2), where p∞ is the tunnel pressure measured at the test section inlet and pv is the saturated vapor pressure, the Reynolds number isRe=U c/ν=7× 105, where ν is the kinamatic vis-∞cosity. And the Froude number iswheregis the gravitation acceleration. In general,the experimental conditions are maintained to within an uncertainty of 1% for the hydrofoil angle of incidence, and an uncertainty of 2% for both the flow velocity and the upstream pressure. The uncertainty of the electromagnetic flow meter is 0.5%, and the uncertainty of the pressure transducer is 0.25%. The cavitation number can be controlled to within an uncertainty of 5%[18,19].

    2. Results and discussions

    2.1 Global multiphase structures in different cavitation regimes

    Fig.3 (Color online) Typical cavity shapes in different cavitation regimes when Re =7× 105, α=8°, viewed from the bottom of the test section

    Depending on the cavitation number σ and the incidence angle α, four cavitation regimes could be identified: the inception cavitation, the sheet cavitation, the cloud cavitation, and the supercavitation.Figure 3 shows the typical cavity shapes in different cavitation regimes at σ=1.8, σ=1.4, σ=0.8 and σ=0.4 as observed in the experimental visualization, when the angle of the hydrofoil is fixed at α= 8°. Figure 4 shows the means and the fluctuations of the measured attached cavity lengths for α= 0°, α=5° and α=8° at various cavitation numbers, respectively. Cavitation inception is observed at about σ=1.8 for α=8°, which is advanced to about σ=1.6 for α=5°, and σ=1.2 for α=0°, respectively. In the regime of the inception cavitation, the cavitating flow shows in white color at the leading edge of the hydrofoil, indicating that it contains a number of micro-sized vapor bubbles.Lowering the cavitation number, the cavitation changes from the inception cavitation to a sheet cavity at about σ=1.4 for α=8°, as shown in Fig.3(b).Although the cavity is attached at the LE, the rear portion of the sheet is unsteady and it rolls up into a series of bubble eddies that shed imtermittently. The sheet cavitation begins to appear at about σ=1.3 for α= 5° and σ=1.0 for α=0°, respectively. Further decreasing the cavitation number to about σ =0.8 for α=8°, the sheet cavitation grows and the trailing edge becomes increasingly unsteady.which is accompanied with massive bubble shedding in the rear portion of the cavity, to form the cloud cavitation. In the cloud cavitation stage, the attached cavity length shows large fluctuations. It should be noted that due to the unsteady shedding of the cavity,the fluctuation amplitudes in the attached cavity are much higher for α=8° than for α=5° and α= 0°. The supercavitation is the final stage of the cavitation, At this stage, the cavitating area covers the entire hydrofoil, extending to the downstream region,the cavity boundary seems to be quite steady, as shown in Fig.3(d), and the attached cavity lengths see less fluctuations, as compared with the cases of the cloud cavitation.

    Fig.4 The measured normalized cavity lengths (L c /c) at various cavitation numbers and attack angles

    2.2 The unsteady multiphase structures of cloud cavitating flows

    In the cloud cavitation regime, the unsteady shedding and the breakdown occur violently, with highly unsteady hydrodynamic forces on the hydrofoil,which may lead to strong dynamic instabilities[20,21].

    In the present experimental study, when the cavitation number is reduced to σ=0.4 for α=0°,σ =0.6 for α=5° and σ=0.8 for α=8°, the cavity grows and becomes increasingly unsteady until the transient cloud cavitation develops. In all above cases, it is observed that the development of the cloud cavitation assumes a distinctly quasi-periodic pattern.Figure 5 shows the side views of the periodic cavity structures, within a single flow cycle in all cases.Based on the experimental visualization, the periodsTcycleof the cloud cavitation in the case of σ=0.6 and α=5° and in the case of σ=0.8 and α=8°are almost the same, taking a value of about 40 ms,which corresponds to 5.7c/U∞, and the period in the case of σ=0.4 and α=0° is about 0.75Tcycle.

    An example of a typical cavity shape observed in the experimental visualization is presented in Fig.6(a),which includes both the original flow visualization captured by the high speed camera and the schematic interpretations drawn by an in-house feature-recognition software package[22]. Figure 6(b) presents the measured time evolution of the nomarlized cavity area over several flow cycles in cloud cavitation regimes(α=0°, σ=0.4, α=5°, σ=0.6, α=8°, σ=0.8) for different attack angles. The curves see a periodic feature. The measured mean cavity areas are 0.49Sc, 0.59Scand 0.72Sc, respectively, whereScis the cross-section area of the hydrofoil.

    Fig.6(b) comparison of the measured time evolution of the normalized cavity area in the typical cloud cavitation regimes (α=0°, σ=0.4, α=5°, σ=0.6,α= 8°, σ=0.8), Re=7× 105, V∞ =10m/s

    The void fraction distribution is of interest, and it should be related to the grayscale contours captured by the video. The mean and the standard deviation of the gray level distributions around the hydrofoil in the cases of α=0°, σ=0.4 and α=8°, σ=0.8 are shown in Fig.7, respectively. The mean cavity area is much larger in the case of σ=0.8 and α=8° than in the case of σ=0.4 and α=0°. It should be noted that the value of the standard deviation in the case of σ=0.8 and α=8° is larger than that in the case of σ=0.4 and α=0°. This is an indication of highly unsteady cavitating behavior and is typical of large-scale cloud fluctuations, as shown in Fig.5(c).

    Figures 8, 9 show the evolution of the cloud shedding process in the case of σ=0.4 and α=0°and in the case of σ=0.8 and α=8°, respectively,viewed from the bottom (foil suction-side) of the test section. In the case of σ=0.8 and α=8°, the attached sheet cavity grows to its maximum length, as shown in Figs. 9(a)-9(d), and a re-entrant flow forms and pushes the flow towards the leading edge. Meanwhile, the interface becomes increasingly unsteady with the presence of bubbles due to the reverse motion of the two-phase mixture. As the re-entrant flow reaches the vicinity of the leading edge, the cavity is lifted away from the wall and shed downstream in the form of the cloud cavitation, as shown in Figs. 9(e),9(f), where the cavity shedding process is very similar to the schematic diagram of the typical cloud cavity transformation shown in Ref. [10]. The complete breakup of the cavity from the leading edge and the unsteady transformation to the cloud cavitation in the case of σ=0.8 and α=8° are different from the experimental visualization in the case of σ=0.4 and α=0°. At a smaller angle of attack, the attached cavity seems much more stable, and the cavity partially breaks at the trailing edge of the hydrofoil,followed by the successive shedding of smaller vapor structures, as shown in Fig.8.

    Fig.7 (Color online) Mean gray level distribution around the hydrofoil in the typical cloud cavitation regimes (α=,0°,σ=0.4, α=8°, σ=0.8), Re=7× 105, V∞=10m/s. A value of 0 corresponds to black and a value of 255 corresponds to pixel saturation (white)

    Fig.6(a) (Color online) Typical flow visualization and schematic interpretation

    Fig.8 (Color online) Photographs of the time sequence of the cloud shedding obtained via high-speed video for σ=0.4, Re=7× 105, α=0°, viewed from the bottom (foil suction-side) of the test section

    Fig.9 (Color online) Photographs of the time sequence of the cloud shedding obtained via high-speed video for σ=0.8, Re=7× 105, α=8°, viewed from the bottom (foil suction-side) of the test section

    Fig.10 (Color online) Numerically predicted time-averaged cavity shapes and re-entrant flow affected regions,Re=7× 105, V=10 m/s∞

    Two parameters are shown to be of the most importance in the analysis of the re-entrant flow instability: the adverse pressure gradient and the cavity thickness as compared with the re-entrant flow thickness[23]. In order to address the physical origin of the instability of the cavity at different angles of attack,Fig.10 shows the predicted time-averaged cavity shapes and the re-entrant flow affected regions in the case of σ=0.8 and α=8° and in the case of σ =0.4 and α=0° , as compared with the numerical simulation, validated with the experimental visualization and data for the same Clark-Y hydrofoil[18].Figure 11 shows the corresponding predicted contours of the negative (upstream) axial (u) velocity as a function of position and time (withTrefandLrefdefined asTref=c/U∞andLref=c). It could be concluded that the cavity thickness is directly controlled by the angle of attack, and the ratio of the thickness of the cavity to that of the re-entrant flow is substantially important to estimate the unsteady behavior of the partial cavity. In the case of σ=0.8 and α= 8°, the cavity is thick enough to limit the interaction between the re-entrant flow and the cavity interface, which allows the re-entrant flow to reach the cavity leading edge and the formation of a large scale vapor structures. However, in the case of σ=0.4 and α=0°, because of the smaller cavity thickness,one sees a strong interaction between the cavity interface and the re-entrant flow. The re-entrant flow can only reach the rear part of the hydrofoil, as shown in Fig.12, contrary to the large scale cloud shedding in the case of σ=0.8 and α=8°, where small vapor structures are formed.

    Fig.11 (Color online) Numerically predicted time evolution of the predicted reverse u-velocity in various sections (the position is in ordinate, from the foil leading edge x/ L0=0 to the trailing edge x/ L0=1), Re=7× 105,V∞=10 m/s

    Fig.12 Schematic representation of the re-entrant flows in the region of an attached cavity in the typical cloud cavitation regimes (α=0°,σ=0.4,α=8°, σ=0.8)

    2.3 The cavitating flow structures

    Fig.13 (Color online) Normalized ensemble averaged axial velocity contours on the suction side of the hydrofoil,Re=7×105 , V∞= 10 m/ s

    Fig.14 Means of the ensemble averaged axial velocity profiles at selected chord wise locations along the foil for cloud cavitating flow, Re=7×105 , V∞= 10 m/ s

    The cavitating flow fields are measured by the PIV technology, and the images are acquired at 2000 Hz,corresponding to a mean flow displacement of 2.3% of the field of view between captures. 500 pairs of PIV images are sufficient to record more than 6 full cloud cavitation shedding cycles. Based on the measured instantaneous velocity fields, the 500 instantaneous velocity fields are averaged in the post- processing to obtain the ensemble averaged in-plane velocity and vorticity fields. Figure 13 compares the contours of the normalized ensemble averaged axial velocity which runs parallel to the inflow direction. As shown in the figure, the axial velocity near the foil trailing edge in the cavitation region is much lower. Figure 14 shows the mean values of the normalized ensemble averaged axial (u) velocity profiles at different chord wise locations in the case of σ=0.8 and α= 8° and in the case of σ=0.4 and α=0°,respectively. TheYaxis represents the vertical distance from the suction side of the hydrofoil.Substantial differences can be observed between the two cloud cavitating cases, especially near the cavity closure region. In fact, in the case of σ=0.8 and α= 8°, the recirculation breaks up and lifts the cavity upward, and the large scale shedding is formed, which significantly increases the thickness of the boundary layer and a more extended low velocity region is observed as compared to the case of σ=0.4 and α= 0°. Consequently, the gradients of the mean axial velocity profiles in the case of σ=0.4 and α=0°tend to be smaller (more evenly distributed) than those in the case of σ=0.8 and α=8°.

    Fig.15 (Color online) Normalized amplitudes of the turbulent velocity fluctuations, Re=7× 105, V∞=10 m/s

    To further investigate the turbulent velocity fields, Figure 15 shows the normalized amplitude of the averages of the fluctuating velocities,I=, in the case of σ=0.8 and α=8°and in the case of σ=0.4 and α=0°. Here,u′andv′ are the horizontal and vertical components of the turbulent velocity fluctuations. In both cloud cavitating cases shown in Fig.16, large velocity fluctuations are seen near the aft half of the hydrofoil and the fluctuation area is extended into the wake,showing the unsteady large-scale fluctuations due to the shedding of the cloud cavity, especially in the case of σ=0.8 and α=8°, meanwhile, the turbulent fluctuations are confined to a thicker layer on the suction side of the hydrofoil. Figure 17 shows the normalized amplitude of the averages of the fluctuating velocities along the selected monitor locations in the case of σ=0.8 and α=8° and in the case of σ=0.4 and α=0°. The fluctuating velocities are much higher in the case of σ=0.8 and α=8°, with a much thicker turbulent boundary layer as compared to the case of σ=0.4 and α= 0°.

    Fig.16 Normalized amplitudes of the turbulent velocity fluctuating velocity profiles at the selected monitoring locations along the foil, Re=7×105 , V∞= 10 m/ s

    Figure 17 shows the side-by-side comparisons of the normalized out-of-plane (z-component) vorticity field, ωz/ωo, in the case of σ=0.8 and α=8°and in the case of σ=0.4 and α=0°. The detailed comparisons of the normalized vorticity profiles in the case of σ=0.8 and α=0° and in the case of σ =0.4 and α=0° at the selected monitoring locations are shown in Fig.18. Here, ωzis defined as:

    Fig.17 (Color online) Normalized averaged z -vorticity contours, ωz /ωo, on the suction side of the hydrofoil,Re=7× 105, V=10 m/s∞

    ( δ≈ 0.3c=0.021m is the approximated turbulent boundary layer thickness at the foil trailing edge based on the experimental measurements of the flow velocity in the cloud cavitating case of σ=0.8 and α= 8°). In both cases, large-scale counter-clockwise and clockwise vortical structures can be observed in the cavity region and at the trailing edge of the hydrofoil, respectively. Comparisons between the vorticity fields in the case of σ=0.8 and α=8° and in the case of σ=0.4 and α=0° suggest that the cavity structures play an important role in the vorticity production, in addition to the turbulent level and the turbulent boundary layer thickness. The detailed normalized vorticity profiles shown in Fig.18 in the case of σ=0.4 and α=0° suggest the existence of a thin vortical layer near the foil surface, with the thickness of the vortical zone growing toward the trailing edge, up to a maximum thickness of about 0.15c. In the case of σ=0.8 and α=8°, on the other hand, the vortical zone is much thicker (near 0.3cat the foil trailing edge), the vorticity profile is drastically different, and the magnitude of the vorticity is much larger than that in the case of σ =0.4 and α=0°.

    Figure 19 shows the measured lift coefficient(CL) and drag coefficient (CD) over a range of cavitation numbers and angles of attack. When the flow is subcavitating, both lift and drag coefficients remain largely unchanged as the cavitation number varies. In the cavitation inception stage, the traveling cavity bubbles appear as shown in Fig.3, and the effect of the cavitation on the lift and drag coefficients is very small. In the subcavitating and inception cavitation stages, the hydrodynamic coefficients show little fluctuations due to the relatively stable cavity dynamics, as indicated by the error bars. Further decreasing the cavitation number, the development of the cavitation increases the drag while reduces the lift coefficient. In the cloud cavitation stage, the larger vortex shedding and the related unsteady movement strongly affect the flow structure around the hydrofoil,leading to a relatively higher magnitude of drag and a lower magnitude of lift. The combination of decreases in lift and increases in drag leads to a decrease in the efficiency, which is a classical issue for the cavitation.As indicated by the error bars of the hydrodynamic coefficients in the cloud cavitation stage, the measured values show large fluctuations, especially in the case of α=8°, which is induced by the unstable cavity dynamics. Finally, while transitioning from the cloud cavitation to the supercavitation, the cavities show more steady characteristics and both drag and lift coefficients are reduced.

    2.4 Hydrodynamic coefficient

    Fig.19 The measured lift coefficients and drag coefficie nts for various cavitation numbers at different angles ofattack,Re=7× 105, V=10 m/s∞

    3. Conclusions

    This paper investigates the unsteady structures and the hydrodynamic characteristics of the cavitating flows, to provide experimental data for future numerical validating studies. Experimental studies are for a Clark-Y hydrofoil at fixed angles of attack α=0°, 5°and 8° atRe=7× 105for various cavitation numbers,from the subcavitating flow to the supercavitation.High-speed videos of the evolution of the cloud cavitation dynamics and PIV measurements of the velocity and vorticity fields reveal the structure of unsteady cavitating flows. Statistics of the cavity lengths, the velocity distributions, the turbulent intensities, as well as the hydrodynamics are presented to better quantify the unsteady process. Detailed analysis of the cavity behavior is applied to identify the re-entrant flow instability which leads to the cloud cavitation at different angles of attack. The primary findings include:

    (1) Cloud cavity is highly unsteady. High-speed photographs show a self-oscillatory behavior of the transformation from sheet to cloud cavitations. In the case of σ=0.4 and α=0°, the thickness of the attached cavity is comparable to the re-entrant flow, a relatively strong interaction is observed between the cavity interface and the re-entrant flow throughout its upstream movement. The cavity in the rear part of the hydrofoil splits into small vapor structures. In the case of σ=0.8 and α=8°, the attached cavity is thick enough to limit the interaction between the re-entrant flow and the cavity interface during its movement to the leading edge and the interface is cut, contrary to the small scale cloud shed in the case of σ=0.4 and α= 0°, the typical large vapor structures are formed.

    (2) Compared to the case of σ=0.4 and α =0°, the velocities in the cavitating region are much lower, and the fluctuating velocities as well as the turbulent intensities are much higher in the cloud cavitating case of σ=0.8 and α=8°, which leads to a much thicker turbulent boundary layer. The results suggest that the periodic shedding and collapse of the vapor cavities are the important mechanisms for vorticity production, as evidenced by the production of large-scale vortical structures observed in the cavitating region of the foil and in the cavitating wake in the case of σ=0.8 and α=8°, the affected zones and the magnitude of the vorticity are much larger than those in the case of σ=0.4 and α=0°.

    (3) The dynamic characteristics of the cavitation vary considerably with the angles of attack and the cavitation numbers. The effect of the inception cavity on the lift is very small, little fluctuation of the dynamic characteristics is obseved. With the decrease of the cavitation numbers, the lift is affected by the cavitating flow structure. At the cloud cavitation, the large scale shedding leads to a relatively higher magnitude of drag and a lower magnitude of lift. As the supercavitation appears, the cavity shapes and the dynamics seem more stable.

    [1] Wang G. Y., Wu Q., Huang B. Dynamics of cavitation–structure interaction [J].Acta Mechanica Sinica, 2017,33(4): 685-708.

    [2] Luo X. W., Ji B., Tsujimoto Y. A review of cavitation in hydraulic machinery [J].Journal of Hydrodynamics, 2016,28(3): 335-358.

    [3] Wang G., Senocak I., Shyy W. et al. Dynamics of attached turbulent cavitating flows [J].Progress in Aerospace sciences, 2001, 37(6): 551-581.

    [4] Peng X. X., Ji B., Cao Y. et al. Combined experimental observation and numerical simulation of the cloud cavitation with U-type flow structures on hydrofoils [J].International Journal of Multiphase Flow, 2016, 79:10-22.

    [5] Franc J. P. Partial cavity instabilities and re-entrant jet [C].Fourth International Symposium on Cavitation, Paris,France, 2001.

    [6] Huang B. Physical and numerical investigation of unsteady cavitating flows [D]. Doctoral Thesis, Beijing,China: Beijing Institute of Technology, 2012(in Chinese).

    [7] Huang B., Ducoin A., Young Y. L. Physical and numerical investigation of cavitating flows around a pitching hydrofoil [J].Physics of Fluids, 2013, 25(10): 102109.

    [8] Ji B., Luo X., Wu Y. et al. Numerical analysis of unsteady cavitating turbulent flow and shedding horse-shoe vortex structure around a twisted hydrofoil [J].International Journal of Multiphase Flow, 2013, 51: 33-43.

    [9] Ji B., Luo X. W., Peng X. X. et al. Three-dimensional large eddy simulation and vorticity analysis of unsteady cavitating flow around a twisted hydrofoil [J].Journal of Hydrodynamics, 2013, 25(4): 510-519.

    [10] Kawanami Y., Kato H., Yamauchi H. et al. Mechanism and control of cloud cavitation [J].Journal of Fluids Engineering, 1997, 119(4): 788-794.

    [11] Fujii A., Kawakami D. T., Tsujimoto Y. et al. Effect of hydrofoil shapes on partial and transitional cavity oscillations [J].Journal of Fluids Engineering, 2007, 129(6):669-673.

    [12] Kjeldsen M., Arndt R. E. A., Effertz M. Spectral characteristics of sheet/cloud cavitation [J].Journal of Fluids Engineering, 2000, 122(3): 481-487.

    [13] Gopalan S., Katz J. Flow structure and modeling issues in the closure region of attached cavitation [J].Physics of fluids, 2000, 12(4): 895-911.

    [14] Arndt R. E. A., Song C. C. S., Kjeldsen M. et al. Instability of partial cavitation: A numerical/experimental approach [C].Proceedings Twenty-Third Symposium on Naval Hydrodynamics, Valde Reuil, France, 2000.

    [15] Ausoni P., Farhat M., Escaler X. et al. Cavitation influence on von Kármán vortex shedding and induced hydrofoil vibrations [J].Journal of Fluids Engineering, 2007,129(8): 966-973.

    [16] Aeschlimann V., Barre S., Djeridi H. Velocity field analysis in an experimental cavitating mixing layer [J].Physics of Fluids, 2011, 23(5): 055105.

    [17] Wosnik M., Fontecha L. G., Arndt R. E. A. Measurements in high void-fraction bubbly wakes created by ventilated supercavitation [C].ASME 2005 Fluids Engineering Division Summer Meeting, Houston, USA, 2005, 531-538.

    [18] Huang B., Yong Y. L., Wang G. et al. Combined experimental and computational investigation of unsteady structure of sheet/cloud cavitation [J].Journal of Fluids Engineering, 2013, 135(7): 071301.

    [19] Li X., Wang G., Zhang M. et al. Structures of super cavitating multiphase flows [J].International Journal of Thermal Sciences, 2008, 47(10): 1263-1275.

    [20] Ji B., Long Y., Long X. P. et al. Large eddy simulation of turbulent attached cavitating flow with special emphasis on large scale structures of the hydrofoil wake and turbulence-cavitation interactions [J].Journal of Hydrodynamics, 2017, 29(1): 27-39.

    [21] Luo X., Ji B., Peng X. X. et al. Numerical simulation of cavity shedding from a three-dimensional twisted hydrofoil and induced pressure fluctuation by large-eddy simulation [J].Journal of Fluids Engineering, 2012,134(4): 041202.

    [22] Zhang M., Song X., Wang G. Design and application of cavitation flow image programs. Design and application of cavitation flow image program [C].International Symposium on Photoelectronic Detection and Imaging 2007: Image Processing, Beijing, China, 2008.

    [23] Callenaere M., Franc J. P., Michel J. M. et al. The cavitation instability induced by the development of a reentrant jet [J].Journal of Fluid Mechanics, 2001, 444:223-256.

    猜你喜歡
    昌樂(lè)南溪李軍
    關(guān)于“借錢”
    新少年(2023年3期)2023-03-20 10:16:30
    唐 南溪詩(shī)序并南溪玄巖銘拓片
    木棉花開
    人民之聲(2022年3期)2022-04-12 12:00:14
    Superconductivity in octagraphene
    A physics-constrained deep residual network for solving the sine-Gordon equation
    干杯!
    新少年(2021年3期)2021-03-28 02:30:27
    Long A Sounds
    擰瓶蓋的男生
    昌樂(lè)傳統(tǒng)石雕工藝的傳承發(fā)展思路探析
    科技視界(2017年10期)2017-08-21 08:46:16
    昌樂(lè)方言詞匯研究
    无遮挡黄片免费观看| 男女那种视频在线观看| 午夜亚洲福利在线播放| 亚洲18禁久久av| 国产精品久久久人人做人人爽| www日本黄色视频网| 日韩成人在线观看一区二区三区| www.自偷自拍.com| 亚洲欧美日韩无卡精品| 草草在线视频免费看| 久久婷婷人人爽人人干人人爱| 日本黄色视频三级网站网址| 亚洲国产看品久久| а√天堂www在线а√下载| 很黄的视频免费| 精品一区二区三区av网在线观看| 欧美日韩国产亚洲二区| 亚洲专区字幕在线| 在线国产一区二区在线| 国产综合懂色| 欧美丝袜亚洲另类 | 中国美女看黄片| 综合色av麻豆| 精品国产三级普通话版| 亚洲成人中文字幕在线播放| 悠悠久久av| 老司机深夜福利视频在线观看| 亚洲精品456在线播放app | 国产精品99久久久久久久久| 国产精品一区二区精品视频观看| 一进一出抽搐gif免费好疼| 91av网一区二区| 十八禁网站免费在线| 国产精品久久久久久久电影 | 真人做人爱边吃奶动态| 变态另类成人亚洲欧美熟女| 天天躁狠狠躁夜夜躁狠狠躁| 岛国视频午夜一区免费看| 久久久精品欧美日韩精品| 国模一区二区三区四区视频 | 国产精品影院久久| 日本精品一区二区三区蜜桃| 我要搜黄色片| 99riav亚洲国产免费| 亚洲国产色片| 国产麻豆成人av免费视频| 久久性视频一级片| 一本一本综合久久| 日韩有码中文字幕| 国产成人一区二区三区免费视频网站| 日本撒尿小便嘘嘘汇集6| 黄色片一级片一级黄色片| 欧美色欧美亚洲另类二区| 欧美3d第一页| 老司机福利观看| 欧美大码av| 黑人操中国人逼视频| 12—13女人毛片做爰片一| 国产精品香港三级国产av潘金莲| 亚洲自拍偷在线| 99久国产av精品| 狂野欧美激情性xxxx| 中出人妻视频一区二区| 色av中文字幕| 久久伊人香网站| 中文字幕最新亚洲高清| 91老司机精品| 黑人欧美特级aaaaaa片| 天堂网av新在线| 麻豆一二三区av精品| 欧美日本视频| 岛国在线观看网站| 国内精品美女久久久久久| 美女午夜性视频免费| 亚洲精品中文字幕一二三四区| 久久草成人影院| 欧美精品啪啪一区二区三区| 国产一区二区在线观看日韩 | a级毛片a级免费在线| 亚洲人与动物交配视频| 一个人看视频在线观看www免费 | 嫩草影院入口| www.精华液| 亚洲色图av天堂| 国内揄拍国产精品人妻在线| 成人无遮挡网站| aaaaa片日本免费| 免费搜索国产男女视频| 男人舔女人的私密视频| 午夜影院日韩av| 国产私拍福利视频在线观看| 噜噜噜噜噜久久久久久91| 免费无遮挡裸体视频| 久久人妻av系列| 久久久久久久久免费视频了| 老司机午夜福利在线观看视频| 欧美又色又爽又黄视频| 好男人电影高清在线观看| 1000部很黄的大片| 久久久久久人人人人人| 国产男靠女视频免费网站| 久久精品aⅴ一区二区三区四区| 国产精品九九99| 亚洲精品中文字幕一二三四区| 欧美黑人巨大hd| 在线观看免费午夜福利视频| 日韩有码中文字幕| 国内精品久久久久久久电影| 在线免费观看不下载黄p国产 | 国产精品综合久久久久久久免费| av国产免费在线观看| 亚洲av美国av| 18禁美女被吸乳视频| 一本久久中文字幕| www.精华液| 国产精品亚洲一级av第二区| 久久久久久久午夜电影| 欧美丝袜亚洲另类 | 免费搜索国产男女视频| 少妇的丰满在线观看| 国产熟女xx| 99re在线观看精品视频| 男人和女人高潮做爰伦理| 国产激情偷乱视频一区二区| 日韩 欧美 亚洲 中文字幕| 中亚洲国语对白在线视频| 久久久久精品国产欧美久久久| 欧美三级亚洲精品| 久久中文字幕人妻熟女| 亚洲男人的天堂狠狠| 国产精品99久久99久久久不卡| 色吧在线观看| 国产淫片久久久久久久久 | 亚洲欧美激情综合另类| 又紧又爽又黄一区二区| 99国产综合亚洲精品| 国产精品亚洲av一区麻豆| www日本在线高清视频| 久久99热这里只有精品18| 久久久久国内视频| 一级毛片高清免费大全| 97人妻精品一区二区三区麻豆| 99久久无色码亚洲精品果冻| 少妇的逼水好多| 亚洲欧美日韩无卡精品| 免费看十八禁软件| 国产精品野战在线观看| 精华霜和精华液先用哪个| 欧美xxxx黑人xx丫x性爽| 看免费av毛片| 日韩 欧美 亚洲 中文字幕| 午夜免费激情av| 波多野结衣巨乳人妻| 成人三级黄色视频| 好男人在线观看高清免费视频| 欧美日本视频| 日本精品一区二区三区蜜桃| 国产亚洲av嫩草精品影院| 欧美日韩亚洲国产一区二区在线观看| 亚洲第一欧美日韩一区二区三区| 欧美另类亚洲清纯唯美| 日日摸夜夜添夜夜添小说| 亚洲在线自拍视频| 国产激情偷乱视频一区二区| 网址你懂的国产日韩在线| 人妻夜夜爽99麻豆av| 国产主播在线观看一区二区| 精品国产乱码久久久久久男人| 亚洲精品久久国产高清桃花| 色综合欧美亚洲国产小说| 国产精品野战在线观看| 国产免费男女视频| ponron亚洲| 精品国产三级普通话版| 欧美日韩乱码在线| 日韩有码中文字幕| 久久伊人香网站| 亚洲一区二区三区不卡视频| 日本黄大片高清| 色吧在线观看| 国产成人影院久久av| 两个人看的免费小视频| 色综合婷婷激情| 又大又爽又粗| 欧美一级毛片孕妇| 免费一级毛片在线播放高清视频| 久久久久国产精品人妻aⅴ院| 亚洲精品在线观看二区| 黄片小视频在线播放| 久久中文字幕一级| 日韩有码中文字幕| 欧美一级a爱片免费观看看| 成年女人看的毛片在线观看| 嫩草影视91久久| 老司机福利观看| av福利片在线观看| 国产免费av片在线观看野外av| 午夜福利免费观看在线| 在线看三级毛片| 午夜日韩欧美国产| 国产乱人视频| 波多野结衣高清作品| 色播亚洲综合网| 首页视频小说图片口味搜索| 级片在线观看| 国产熟女xx| 两个人的视频大全免费| 宅男免费午夜| 国产成年人精品一区二区| 亚洲成av人片免费观看| 国产激情欧美一区二区| 99视频精品全部免费 在线 | 亚洲五月婷婷丁香| 欧美3d第一页| 国产乱人伦免费视频| 男女下面进入的视频免费午夜| 国产黄色小视频在线观看| 精品无人区乱码1区二区| 欧美日本视频| 精华霜和精华液先用哪个| 国产午夜精品久久久久久| 国产成人影院久久av| 精品久久蜜臀av无| 不卡av一区二区三区| 国产不卡一卡二| 99国产精品一区二区蜜桃av| 99热这里只有精品一区 | 69av精品久久久久久| 99久久综合精品五月天人人| 亚洲va日本ⅴa欧美va伊人久久| 精品一区二区三区视频在线观看免费| 丰满人妻一区二区三区视频av | 一个人免费在线观看电影 | 国产男靠女视频免费网站| 人人妻人人澡欧美一区二区| 99re在线观看精品视频| 免费看十八禁软件| 蜜桃久久精品国产亚洲av| 真实男女啪啪啪动态图| 亚洲自偷自拍图片 自拍| 亚洲最大成人中文| 动漫黄色视频在线观看| 婷婷六月久久综合丁香| a在线观看视频网站| 丰满人妻熟妇乱又伦精品不卡| 国产精品一区二区三区四区免费观看 | 亚洲美女视频黄频| 日韩中文字幕欧美一区二区| 国产精品一及| 青草久久国产| 中文字幕av在线有码专区| 午夜两性在线视频| 一二三四社区在线视频社区8| 18禁黄网站禁片午夜丰满| 国产免费男女视频| 国产精品久久久人人做人人爽| 一本久久中文字幕| 亚洲人成网站高清观看| 亚洲一区二区三区色噜噜| 日韩 欧美 亚洲 中文字幕| a在线观看视频网站| 久久午夜综合久久蜜桃| 亚洲成人久久性| 国产成人福利小说| 舔av片在线| 亚洲欧美精品综合久久99| 亚洲欧美精品综合一区二区三区| 国产精品,欧美在线| 国产精华一区二区三区| 日韩欧美在线二视频| 成人高潮视频无遮挡免费网站| 欧美av亚洲av综合av国产av| 色播亚洲综合网| 欧美黑人欧美精品刺激| 男女做爰动态图高潮gif福利片| 国产午夜精品论理片| 身体一侧抽搐| 脱女人内裤的视频| 九九在线视频观看精品| 99久久精品热视频| 一本一本综合久久| 成人欧美大片| 国内精品美女久久久久久| 他把我摸到了高潮在线观看| 国模一区二区三区四区视频 | 女人高潮潮喷娇喘18禁视频| 99久久无色码亚洲精品果冻| 亚洲欧洲精品一区二区精品久久久| 日韩欧美在线乱码| 国产高清videossex| 欧美日韩精品网址| 国产主播在线观看一区二区| 国产精华一区二区三区| 亚洲成a人片在线一区二区| 麻豆成人午夜福利视频| 99热这里只有是精品50| 成人国产综合亚洲| 亚洲欧美日韩无卡精品| 免费大片18禁| 欧美日韩亚洲国产一区二区在线观看| av福利片在线观看| 国产欧美日韩一区二区三| 久久精品91无色码中文字幕| 精品久久蜜臀av无| 久久国产精品影院| 美女黄网站色视频| 99riav亚洲国产免费| 国产综合懂色| 国内精品美女久久久久久| 欧美一区二区国产精品久久精品| 成人三级做爰电影| 人妻夜夜爽99麻豆av| 免费观看人在逋| svipshipincom国产片| 两个人的视频大全免费| 深夜精品福利| 国产主播在线观看一区二区| 国产精品久久视频播放| 欧美国产日韩亚洲一区| 成人特级黄色片久久久久久久| 精品乱码久久久久久99久播| 国产淫片久久久久久久久 | 亚洲国产精品sss在线观看| 午夜免费激情av| 国产探花在线观看一区二区| 精品人妻1区二区| 岛国在线免费视频观看| 天堂网av新在线| 亚洲自拍偷在线| 最近在线观看免费完整版| 夜夜看夜夜爽夜夜摸| 麻豆av在线久日| 网址你懂的国产日韩在线| 最近最新中文字幕大全免费视频| 天天一区二区日本电影三级| 51午夜福利影视在线观看| 99精品在免费线老司机午夜| 91av网一区二区| 亚洲avbb在线观看| 美女扒开内裤让男人捅视频| 亚洲熟女毛片儿| 亚洲精品粉嫩美女一区| 亚洲一区高清亚洲精品| 一个人看视频在线观看www免费 | 国产不卡一卡二| 两人在一起打扑克的视频| 亚洲 欧美 日韩 在线 免费| 免费看美女性在线毛片视频| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美国产一区二区入口| 这个男人来自地球电影免费观看| 日韩欧美在线乱码| 国产一区在线观看成人免费| 真人一进一出gif抽搐免费| 久久精品国产综合久久久| 在线观看一区二区三区| 午夜福利视频1000在线观看| 9191精品国产免费久久| 国产精品久久久av美女十八| 亚洲精华国产精华精| 精品国内亚洲2022精品成人| 91av网一区二区| 欧美色欧美亚洲另类二区| 国产精品影院久久| 又爽又黄无遮挡网站| 亚洲国产精品sss在线观看| 国产精品亚洲av一区麻豆| 久久人人精品亚洲av| e午夜精品久久久久久久| 国产三级中文精品| 欧美色欧美亚洲另类二区| 亚洲真实伦在线观看| 97超级碰碰碰精品色视频在线观看| 久久欧美精品欧美久久欧美| 亚洲一区高清亚洲精品| 国产一区二区在线观看日韩 | 热99re8久久精品国产| 黄色日韩在线| 女生性感内裤真人,穿戴方法视频| www.精华液| av中文乱码字幕在线| 免费观看的影片在线观看| 日本熟妇午夜| 亚洲精品粉嫩美女一区| 高清在线国产一区| 亚洲美女黄片视频| 成人av在线播放网站| 一级黄色大片毛片| 91av网站免费观看| 99精品在免费线老司机午夜| 国产精品久久久久久精品电影| 首页视频小说图片口味搜索| 欧美一区二区精品小视频在线| 最近最新中文字幕大全免费视频| 免费一级毛片在线播放高清视频| 精品国产三级普通话版| 久久久久国内视频| 国产精品99久久久久久久久| 欧美一区二区国产精品久久精品| 国产精品一区二区免费欧美| 国产一区在线观看成人免费| 又大又爽又粗| 欧美精品啪啪一区二区三区| av在线天堂中文字幕| 人人妻人人看人人澡| 男人舔女人下体高潮全视频| 麻豆av在线久日| 欧美绝顶高潮抽搐喷水| 九九热线精品视视频播放| 啦啦啦观看免费观看视频高清| 精华霜和精华液先用哪个| 免费在线观看成人毛片| 两个人的视频大全免费| 最近最新免费中文字幕在线| 欧美乱色亚洲激情| 国产午夜精品久久久久久| 午夜福利成人在线免费观看| 伊人久久大香线蕉亚洲五| 可以在线观看的亚洲视频| 亚洲国产欧美网| 久久这里只有精品中国| 国产真人三级小视频在线观看| 国产av不卡久久| 国产激情久久老熟女| 少妇人妻一区二区三区视频| 亚洲国产色片| 亚洲欧美精品综合一区二区三区| 国产成人精品无人区| 国产成人精品久久二区二区免费| 日韩三级视频一区二区三区| 国产精品一区二区免费欧美| 久久精品aⅴ一区二区三区四区| 一进一出抽搐动态| 在线国产一区二区在线| 美女黄网站色视频| 色综合欧美亚洲国产小说| 午夜免费成人在线视频| av女优亚洲男人天堂 | 成在线人永久免费视频| 怎么达到女性高潮| 国产精品免费一区二区三区在线| 日本三级黄在线观看| 97人妻精品一区二区三区麻豆| 国产成人欧美在线观看| 精品国产三级普通话版| 亚洲欧美日韩高清在线视频| a级毛片a级免费在线| 免费在线观看视频国产中文字幕亚洲| 欧美成人一区二区免费高清观看 | 国产日本99.免费观看| 精华霜和精华液先用哪个| 高清在线国产一区| 亚洲精品国产精品久久久不卡| 日韩欧美在线乱码| 国产高潮美女av| 免费搜索国产男女视频| 欧美中文综合在线视频| 淫妇啪啪啪对白视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产精华一区二区三区| 欧美乱色亚洲激情| 亚洲无线在线观看| xxx96com| 久久这里只有精品中国| 99久国产av精品| 午夜影院日韩av| 免费在线观看成人毛片| 久久久久亚洲av毛片大全| 97超级碰碰碰精品色视频在线观看| 久久精品人妻少妇| 久99久视频精品免费| 日韩精品青青久久久久久| 中出人妻视频一区二区| 欧美黄色片欧美黄色片| 黑人操中国人逼视频| 欧美成人性av电影在线观看| 九九在线视频观看精品| 美女 人体艺术 gogo| 日韩 欧美 亚洲 中文字幕| 最新美女视频免费是黄的| 国产成人福利小说| www日本在线高清视频| 黄片大片在线免费观看| av天堂在线播放| 法律面前人人平等表现在哪些方面| 在线免费观看的www视频| 丁香六月欧美| 99国产精品一区二区蜜桃av| av在线蜜桃| 色视频www国产| 久久亚洲精品不卡| 免费在线观看日本一区| 国产黄色小视频在线观看| 九九在线视频观看精品| 国内精品久久久久久久电影| 国内揄拍国产精品人妻在线| 亚洲av第一区精品v没综合| 在线国产一区二区在线| 90打野战视频偷拍视频| 国产精品野战在线观看| 久久香蕉国产精品| 看免费av毛片| 亚洲精品粉嫩美女一区| 小蜜桃在线观看免费完整版高清| 露出奶头的视频| 黑人操中国人逼视频| 国产伦人伦偷精品视频| 亚洲专区国产一区二区| 欧美黑人巨大hd| 网址你懂的国产日韩在线| 嫩草影院入口| 三级男女做爰猛烈吃奶摸视频| 母亲3免费完整高清在线观看| 成年版毛片免费区| 一边摸一边抽搐一进一小说| 亚洲中文字幕日韩| 成年人黄色毛片网站| 人妻夜夜爽99麻豆av| 性欧美人与动物交配| 国产91精品成人一区二区三区| 黄色丝袜av网址大全| 久久天堂一区二区三区四区| 欧美大码av| 欧美黄色片欧美黄色片| 天堂网av新在线| 九色成人免费人妻av| 99国产精品99久久久久| 国产精品av久久久久免费| 国产精品亚洲一级av第二区| 很黄的视频免费| 一进一出好大好爽视频| 中文资源天堂在线| 男人的好看免费观看在线视频| 99国产极品粉嫩在线观看| 在线播放国产精品三级| 精品久久久久久久毛片微露脸| 亚洲精品一卡2卡三卡4卡5卡| 操出白浆在线播放| 久久亚洲精品不卡| 亚洲av片天天在线观看| 又大又爽又粗| 亚洲乱码一区二区免费版| 国产又黄又爽又无遮挡在线| 精品国产乱码久久久久久男人| 日韩中文字幕欧美一区二区| 久久久国产成人精品二区| 97超视频在线观看视频| 一个人看视频在线观看www免费 | 又爽又黄无遮挡网站| 日本精品一区二区三区蜜桃| 在线十欧美十亚洲十日本专区| 午夜精品在线福利| 美女免费视频网站| 中文字幕久久专区| 亚洲黑人精品在线| 此物有八面人人有两片| 欧美大码av| 观看美女的网站| 麻豆久久精品国产亚洲av| 91麻豆av在线| 国产午夜精品久久久久久| 精品福利观看| 精品无人区乱码1区二区| 日韩欧美一区二区三区在线观看| 亚洲av片天天在线观看| 亚洲五月婷婷丁香| 啦啦啦免费观看视频1| 99久久精品一区二区三区| 亚洲国产欧美一区二区综合| 欧美又色又爽又黄视频| 久久伊人香网站| 黄色视频,在线免费观看| 成年女人看的毛片在线观看| 精品午夜福利视频在线观看一区| 国产一区二区在线av高清观看| 黄色日韩在线| 国产97色在线日韩免费| 亚洲精品美女久久久久99蜜臀| 国产精品国产高清国产av| 熟女人妻精品中文字幕| 在线永久观看黄色视频| 少妇的逼水好多| 国模一区二区三区四区视频 | 亚洲男人的天堂狠狠| 国产精品爽爽va在线观看网站| 国产精品久久久久久久电影 | 免费看光身美女| 久99久视频精品免费| 一进一出好大好爽视频| 老司机深夜福利视频在线观看| 久久九九热精品免费| 亚洲精品在线观看二区| 99热这里只有是精品50| 51午夜福利影视在线观看| 国内少妇人妻偷人精品xxx网站 | 亚洲国产欧美人成| av福利片在线观看| 精品无人区乱码1区二区| 国产91精品成人一区二区三区| 欧美乱码精品一区二区三区| 好男人在线观看高清免费视频| 国产不卡一卡二| 高潮久久久久久久久久久不卡| 欧美zozozo另类| 午夜免费成人在线视频| 夜夜躁狠狠躁天天躁| 91久久精品国产一区二区成人 | avwww免费| 欧美乱码精品一区二区三区| 免费人成视频x8x8入口观看| 亚洲中文字幕日韩| 91久久精品国产一区二区成人 | 亚洲精品粉嫩美女一区| 国产成人精品久久二区二区免费| 国产伦人伦偷精品视频|