• <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è)方言詞匯研究
    日韩中字成人| 人妻系列 视频| 日日啪夜夜爽| 精品酒店卫生间| 婷婷色麻豆天堂久久| 天堂√8在线中文| 免费观看无遮挡的男女| 免费在线观看成人毛片| 日韩成人伦理影院| 91aial.com中文字幕在线观看| 亚洲内射少妇av| 六月丁香七月| 乱码一卡2卡4卡精品| 国产毛片a区久久久久| 亚洲精品自拍成人| 亚洲最大成人av| 亚洲欧美清纯卡通| 国产一区有黄有色的免费视频 | 亚洲精品成人久久久久久| 一级毛片电影观看| 国产美女午夜福利| 国产在线男女| 亚洲av国产av综合av卡| 日韩强制内射视频| 建设人人有责人人尽责人人享有的 | 日韩中字成人| 亚洲婷婷狠狠爱综合网| 国产91av在线免费观看| 亚洲经典国产精华液单| 亚洲伊人久久精品综合| 中国国产av一级| 少妇猛男粗大的猛烈进出视频 | 国国产精品蜜臀av免费| av又黄又爽大尺度在线免费看| 免费电影在线观看免费观看| 亚洲国产欧美在线一区| 国产 一区精品| 免费不卡的大黄色大毛片视频在线观看 | 天堂影院成人在线观看| 国产精品1区2区在线观看.| 哪个播放器可以免费观看大片| 国产一区亚洲一区在线观看| 亚洲av男天堂| 国产在视频线在精品| 国产伦精品一区二区三区四那| 欧美激情国产日韩精品一区| 国产一级毛片七仙女欲春2| 亚洲精品视频女| 国产黄频视频在线观看| 中文欧美无线码| 国产伦精品一区二区三区视频9| 亚洲av福利一区| 激情 狠狠 欧美| 亚洲国产精品成人综合色| 中文乱码字字幕精品一区二区三区 | 久久精品国产亚洲av涩爱| 国产黄频视频在线观看| 91久久精品国产一区二区成人| 免费看日本二区| 18+在线观看网站| 亚洲经典国产精华液单| 国产午夜精品久久久久久一区二区三区| 久久精品久久久久久噜噜老黄| 看非洲黑人一级黄片| 久久久成人免费电影| 亚洲国产欧美在线一区| 夫妻午夜视频| 国产欧美另类精品又又久久亚洲欧美| 久久97久久精品| 午夜精品在线福利| 午夜激情福利司机影院| 欧美xxxx性猛交bbbb| 亚洲精品456在线播放app| 精品亚洲乱码少妇综合久久| 青春草视频在线免费观看| 中文资源天堂在线| 大话2 男鬼变身卡| 午夜免费男女啪啪视频观看| 国产视频内射| 国产 一区 欧美 日韩| 国产成人午夜福利电影在线观看| 久久久久国产网址| 国内少妇人妻偷人精品xxx网站| 国产精品国产三级国产av玫瑰| 久久99热这里只频精品6学生| 午夜激情福利司机影院| 精品人妻一区二区三区麻豆| 可以在线观看毛片的网站| 成人高潮视频无遮挡免费网站| 在线天堂最新版资源| 婷婷色麻豆天堂久久| 自拍偷自拍亚洲精品老妇| 麻豆成人av视频| 男人狂女人下面高潮的视频| 少妇人妻精品综合一区二区| 国产美女午夜福利| 一级av片app| 国产免费又黄又爽又色| 九草在线视频观看| 国内少妇人妻偷人精品xxx网站| 欧美日韩一区二区视频在线观看视频在线 | 老女人水多毛片| 麻豆乱淫一区二区| 国产视频首页在线观看| 亚洲欧美成人综合另类久久久| 最近最新中文字幕免费大全7| 色吧在线观看| 七月丁香在线播放| 久久久久久国产a免费观看| 丰满少妇做爰视频| 国产精品一区www在线观看| 成人av在线播放网站| 欧美3d第一页| 亚洲精品国产av蜜桃| 国产欧美日韩精品一区二区| 国产真实伦视频高清在线观看| 两个人视频免费观看高清| 亚洲图色成人| 国内少妇人妻偷人精品xxx网站| 激情 狠狠 欧美| 一级av片app| 亚洲经典国产精华液单| 亚洲欧洲日产国产| 亚洲精品中文字幕在线视频 | 欧美日韩精品成人综合77777| 18禁在线无遮挡免费观看视频| 特大巨黑吊av在线直播| 中文资源天堂在线| 中文字幕av成人在线电影| 亚洲av.av天堂| 可以在线观看毛片的网站| 免费观看性生交大片5| 免费黄频网站在线观看国产| 狂野欧美白嫩少妇大欣赏| 精品久久久精品久久久| 日韩成人av中文字幕在线观看| 亚洲国产欧美人成| 少妇的逼好多水| 久久久成人免费电影| 免费电影在线观看免费观看| 七月丁香在线播放| 日韩大片免费观看网站| 亚洲精品乱久久久久久| 国产精品av视频在线免费观看| 日日撸夜夜添| 久久久久久久午夜电影| 国产 亚洲一区二区三区 | 3wmmmm亚洲av在线观看| 最新中文字幕久久久久| 国产真实伦视频高清在线观看| 久久精品国产自在天天线| 我的女老师完整版在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产av国产精品国产| 天美传媒精品一区二区| 欧美日韩视频高清一区二区三区二| 日日干狠狠操夜夜爽| 最后的刺客免费高清国语| 国产片特级美女逼逼视频| 精品人妻熟女av久视频| 久久99热6这里只有精品| 午夜精品国产一区二区电影 | 少妇猛男粗大的猛烈进出视频 | 亚洲av电影在线观看一区二区三区 | 成人亚洲精品av一区二区| 国产高清三级在线| 波多野结衣巨乳人妻| 99久久九九国产精品国产免费| 高清午夜精品一区二区三区| 九九在线视频观看精品| 嘟嘟电影网在线观看| 久久精品久久久久久久性| 欧美成人精品欧美一级黄| 蜜桃亚洲精品一区二区三区| 久久人人爽人人片av| 亚洲精品中文字幕在线视频 | 国产精品一区二区三区四区久久| 只有这里有精品99| 亚洲一级一片aⅴ在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲熟妇中文字幕五十中出| 综合色av麻豆| 少妇裸体淫交视频免费看高清| 能在线免费看毛片的网站| 一区二区三区乱码不卡18| 51国产日韩欧美| 国产单亲对白刺激| videos熟女内射| 又黄又爽又刺激的免费视频.| 国产色婷婷99| 免费av观看视频| 大香蕉久久网| 国产av码专区亚洲av| 成人毛片60女人毛片免费| 久久99蜜桃精品久久| 美女主播在线视频| 亚洲av男天堂| 成人漫画全彩无遮挡| 男女视频在线观看网站免费| 日韩亚洲欧美综合| 久久精品久久久久久久性| 国产高清三级在线| 内射极品少妇av片p| 国产精品.久久久| 你懂的网址亚洲精品在线观看| 亚洲精品久久久久久婷婷小说| 少妇被粗大猛烈的视频| 高清日韩中文字幕在线| 一边亲一边摸免费视频| 我要看日韩黄色一级片| 性色avwww在线观看| 超碰97精品在线观看| 亚洲天堂国产精品一区在线| 国产成人精品一,二区| 日本三级黄在线观看| 国产成人a区在线观看| 国产成人a区在线观看| 国产黄片美女视频| 91午夜精品亚洲一区二区三区| 亚洲av免费在线观看| 成人特级av手机在线观看| 嫩草影院精品99| 亚洲精品久久午夜乱码| 亚洲国产欧美人成| 偷拍熟女少妇极品色| 欧美一区二区亚洲| 噜噜噜噜噜久久久久久91| 久久久久久国产a免费观看| 国产精品精品国产色婷婷| 亚洲国产欧美在线一区| 国产精品国产三级国产av玫瑰| 美女黄网站色视频| 亚洲国产欧美在线一区| 一二三四中文在线观看免费高清| 丰满少妇做爰视频| 亚洲综合精品二区| 晚上一个人看的免费电影| 午夜福利成人在线免费观看| 一级av片app| 黄色欧美视频在线观看| 大香蕉久久网| 伦精品一区二区三区| .国产精品久久| 最近中文字幕2019免费版| 三级经典国产精品| 一区二区三区高清视频在线| 久久99精品国语久久久| 最新中文字幕久久久久| 国产一区二区在线观看日韩| 少妇人妻一区二区三区视频| av免费在线看不卡| 有码 亚洲区| 国产免费福利视频在线观看| 岛国毛片在线播放| 精品不卡国产一区二区三区| 身体一侧抽搐| 两个人的视频大全免费| 久久热精品热| 国产黄色免费在线视频| 日韩精品有码人妻一区| 日韩精品有码人妻一区| 日韩欧美三级三区| 成人特级av手机在线观看| 深爱激情五月婷婷| 最近视频中文字幕2019在线8| 久久精品久久精品一区二区三区| 好男人视频免费观看在线| 国模一区二区三区四区视频| 久久韩国三级中文字幕| 国产黄色视频一区二区在线观看| 国产成人aa在线观看| 女的被弄到高潮叫床怎么办| av线在线观看网站| 秋霞伦理黄片| 久久久久免费精品人妻一区二区| 亚洲美女搞黄在线观看| 亚洲精品乱码久久久v下载方式| 直男gayav资源| 久久久久精品性色| 国产亚洲精品久久久com| 日韩av在线大香蕉| 亚洲高清免费不卡视频| 免费看美女性在线毛片视频| 亚洲av一区综合| 亚洲欧美清纯卡通| 日本黄色片子视频| 亚洲av不卡在线观看| 麻豆乱淫一区二区| 建设人人有责人人尽责人人享有的 | 美女黄网站色视频| 国产成人福利小说| 国产麻豆成人av免费视频| 国产又色又爽无遮挡免| 欧美三级亚洲精品| 亚洲精品色激情综合| 日韩欧美精品免费久久| 日韩欧美国产在线观看| 永久网站在线| 狂野欧美激情性xxxx在线观看| 国产人妻一区二区三区在| 亚洲天堂国产精品一区在线| 国产白丝娇喘喷水9色精品| 久久久久精品久久久久真实原创| 日日摸夜夜添夜夜爱| 欧美高清性xxxxhd video| 亚洲最大成人中文| 街头女战士在线观看网站| 日本免费a在线| 成人二区视频| 十八禁网站网址无遮挡 | www.色视频.com| 超碰av人人做人人爽久久| 中文字幕亚洲精品专区| 久久久久久久久久人人人人人人| 丰满人妻一区二区三区视频av| 亚洲精品国产av成人精品| 日本一本二区三区精品| 欧美精品一区二区大全| 2018国产大陆天天弄谢| 狂野欧美白嫩少妇大欣赏| 日日摸夜夜添夜夜添av毛片| 街头女战士在线观看网站| 日本免费a在线| 亚洲精品乱码久久久v下载方式| 国产老妇女一区| 99久久中文字幕三级久久日本| 老司机影院成人| 久久精品熟女亚洲av麻豆精品 | 国产成人福利小说| 午夜激情福利司机影院| 内射极品少妇av片p| 精品久久久久久电影网| 午夜福利视频精品| 精品亚洲乱码少妇综合久久| 欧美xxⅹ黑人| 激情五月婷婷亚洲| 秋霞在线观看毛片| 成人高潮视频无遮挡免费网站| 亚洲伊人久久精品综合| 国内精品宾馆在线| 国产午夜精品一二区理论片| 久久99精品国语久久久| 一区二区三区四区激情视频| 国产久久久一区二区三区| 黄色一级大片看看| 日韩三级伦理在线观看| 在线免费观看不下载黄p国产| 人妻一区二区av| 成人国产麻豆网| 久久国内精品自在自线图片| 精品久久久噜噜| 天天躁夜夜躁狠狠久久av| 国产精品女同一区二区软件| 九九爱精品视频在线观看| 日韩欧美 国产精品| 校园人妻丝袜中文字幕| 老女人水多毛片| 日韩精品青青久久久久久| 午夜福利视频精品| 日韩视频在线欧美| 亚洲激情五月婷婷啪啪| 一夜夜www| 伊人久久国产一区二区| 日韩av在线免费看完整版不卡| 校园人妻丝袜中文字幕| 99热6这里只有精品| 成年版毛片免费区| 久久久久久国产a免费观看| 久热久热在线精品观看| 精品国产露脸久久av麻豆 | 久久久久网色| 欧美激情国产日韩精品一区| 自拍偷自拍亚洲精品老妇| 观看美女的网站| 亚洲av在线观看美女高潮| 午夜福利网站1000一区二区三区| 国产亚洲最大av| 国产成人aa在线观看| 成人高潮视频无遮挡免费网站| av免费在线看不卡| 国产黄色视频一区二区在线观看| 波野结衣二区三区在线| 日韩强制内射视频| 天堂影院成人在线观看| 精品人妻偷拍中文字幕| 国精品久久久久久国模美| 午夜精品一区二区三区免费看| 亚洲精品亚洲一区二区| 一级二级三级毛片免费看| 精品久久久久久成人av| 国产成人freesex在线| 国产精品福利在线免费观看| 欧美日韩国产mv在线观看视频 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 爱豆传媒免费全集在线观看| 午夜爱爱视频在线播放| 亚洲欧美日韩东京热| 丰满人妻一区二区三区视频av| 亚洲成人久久爱视频| 99久久九九国产精品国产免费| 男女边吃奶边做爰视频| 亚洲精品国产av蜜桃| 日本-黄色视频高清免费观看| 狂野欧美激情性xxxx在线观看| 高清毛片免费看| 在线观看av片永久免费下载| 内地一区二区视频在线| 99视频精品全部免费 在线| 国产精品一区二区三区四区久久| 99热全是精品| 美女脱内裤让男人舔精品视频| a级毛片免费高清观看在线播放| 国产精品一二三区在线看| 在线a可以看的网站| 国产精品女同一区二区软件| 国产一区有黄有色的免费视频 | 18禁在线播放成人免费| av在线老鸭窝| 久久国产乱子免费精品| 免费观看的影片在线观看| 日日摸夜夜添夜夜添av毛片| 日日撸夜夜添| 成年免费大片在线观看| 成人亚洲欧美一区二区av| 激情 狠狠 欧美| 高清毛片免费看| 国产三级在线视频| 免费观看的影片在线观看| 亚洲国产日韩欧美精品在线观看| 美女国产视频在线观看| 麻豆av噜噜一区二区三区| 九九在线视频观看精品| 国产高清国产精品国产三级 | 亚洲自偷自拍三级| 建设人人有责人人尽责人人享有的 | 麻豆成人av视频| 亚洲成人精品中文字幕电影| 婷婷六月久久综合丁香| 久久99蜜桃精品久久| 亚洲精品日韩在线中文字幕| 久久久久久久久久成人| 日日摸夜夜添夜夜爱| 777米奇影视久久| 色综合站精品国产| av线在线观看网站| 亚洲欧美精品自产自拍| 六月丁香七月| 国产精品不卡视频一区二区| 国产成年人精品一区二区| 国产视频内射| 亚洲18禁久久av| 日本一二三区视频观看| 亚洲国产日韩欧美精品在线观看| 国产精品一及| 久久草成人影院| 亚洲精品亚洲一区二区| 日日摸夜夜添夜夜爱| av在线亚洲专区| av一本久久久久| 成人美女网站在线观看视频| 少妇丰满av| 嘟嘟电影网在线观看| 亚洲乱码一区二区免费版| 久久久久网色| 欧美三级亚洲精品| 最后的刺客免费高清国语| 人妻制服诱惑在线中文字幕| 乱人视频在线观看| 一级毛片aaaaaa免费看小| 国产精品国产三级国产av玫瑰| 国产v大片淫在线免费观看| 欧美三级亚洲精品| 高清欧美精品videossex| 色综合色国产| 观看美女的网站| 午夜视频国产福利| 爱豆传媒免费全集在线观看| 高清在线视频一区二区三区| 热99在线观看视频| 九九爱精品视频在线观看| 亚洲欧美精品专区久久| 国产精品美女特级片免费视频播放器| 九九在线视频观看精品| 一二三四中文在线观看免费高清| 欧美+日韩+精品| 国产精品爽爽va在线观看网站| 中文字幕久久专区| 亚洲av成人av| 天堂影院成人在线观看| 久久这里只有精品中国| 午夜精品国产一区二区电影 | 国产综合懂色| 国产伦精品一区二区三区视频9| 中文字幕制服av| 蜜臀久久99精品久久宅男| 亚洲精品日本国产第一区| 国产伦理片在线播放av一区| 久久久国产一区二区| 亚洲av成人精品一区久久| 日韩制服骚丝袜av| 麻豆av噜噜一区二区三区| 亚洲成人精品中文字幕电影| 日本黄色片子视频| 国产精品国产三级专区第一集| 亚洲av中文av极速乱| 99久久精品国产国产毛片| 免费观看在线日韩| 日本一二三区视频观看| 在线a可以看的网站| 久久6这里有精品| 亚洲最大成人中文| 一级二级三级毛片免费看| 久久久欧美国产精品| 亚洲综合色惰| av线在线观看网站| 亚洲av成人精品一二三区| 成人亚洲精品av一区二区| 亚洲欧美日韩无卡精品| 日韩一区二区三区影片| 久久精品人妻少妇| 国产精品一区www在线观看| 亚洲精品aⅴ在线观看| 精品一区二区三区视频在线| 亚洲电影在线观看av| 一区二区三区乱码不卡18| 亚洲成人av在线免费| av国产久精品久网站免费入址| 久久久久久久久久久免费av| 18禁在线播放成人免费| 波多野结衣巨乳人妻| 亚洲自偷自拍三级| 真实男女啪啪啪动态图| 菩萨蛮人人尽说江南好唐韦庄| 欧美 日韩 精品 国产| 亚洲熟女精品中文字幕| 一级av片app| 亚洲熟女精品中文字幕| 国产精品伦人一区二区| 久久久午夜欧美精品| 免费在线观看成人毛片| 午夜日本视频在线| 最近的中文字幕免费完整| 两个人的视频大全免费| 亚洲欧美清纯卡通| 欧美日韩综合久久久久久| av在线播放精品| 伦理电影大哥的女人| 成人午夜高清在线视频| 亚洲欧美日韩东京热| 一夜夜www| 国产精品一区二区三区四区免费观看| av一本久久久久| 欧美高清成人免费视频www| 青春草亚洲视频在线观看| 在线播放无遮挡| 免费大片18禁| 天天一区二区日本电影三级| 观看美女的网站| 国内精品一区二区在线观看| 国产免费视频播放在线视频 | 国产精品爽爽va在线观看网站| 亚洲欧美成人精品一区二区| 久久综合国产亚洲精品| 高清视频免费观看一区二区 | 99re6热这里在线精品视频| 国产黄色小视频在线观看| 亚洲在久久综合| 天美传媒精品一区二区| 国产成人精品福利久久| 国模一区二区三区四区视频| av.在线天堂| 99热6这里只有精品| 国产午夜精品久久久久久一区二区三区| 成人鲁丝片一二三区免费| 亚洲精品日韩av片在线观看| 欧美日韩亚洲高清精品| 又大又黄又爽视频免费| 久久精品国产鲁丝片午夜精品| 亚洲欧美成人综合另类久久久| 亚洲av在线观看美女高潮| 美女黄网站色视频| 一级爰片在线观看| 亚洲熟妇中文字幕五十中出| 可以在线观看毛片的网站| or卡值多少钱| 亚洲国产日韩欧美精品在线观看| 免费大片18禁| 亚洲色图av天堂| 能在线免费看毛片的网站| 国内精品美女久久久久久| 亚洲精品成人av观看孕妇| 国产真实伦视频高清在线观看| 久久草成人影院| 久久99精品国语久久久| 少妇人妻精品综合一区二区| 2021天堂中文幕一二区在线观| 女人被狂操c到高潮| 免费高清在线观看视频在线观看| 久久久久精品性色| 国产午夜精品久久久久久一区二区三区| 午夜福利成人在线免费观看| 亚洲怡红院男人天堂| 国产探花在线观看一区二区| 亚洲熟妇中文字幕五十中出| 秋霞伦理黄片| 日日干狠狠操夜夜爽| 亚洲国产高清在线一区二区三| 成人av在线播放网站| av在线观看视频网站免费| 色尼玛亚洲综合影院| 日本猛色少妇xxxxx猛交久久| 日韩国内少妇激情av| 淫秽高清视频在线观看| 久久久久久久久中文| 在线观看免费高清a一片|