• <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è)方言詞匯研究
    精品国产超薄肉色丝袜足j| 国产精品麻豆人妻色哟哟久久| 精品久久蜜臀av无| 18在线观看网站| 欧美激情高清一区二区三区 | 美女视频免费永久观看网站| 天堂俺去俺来也www色官网| www日本在线高清视频| 久久精品国产综合久久久| av国产精品久久久久影院| 黑人欧美特级aaaaaa片| 久久人人爽av亚洲精品天堂| 国产一区二区 视频在线| 欧美国产精品一级二级三级| 国产一区亚洲一区在线观看| 一本色道久久久久久精品综合| 日本wwww免费看| 只有这里有精品99| 黄片无遮挡物在线观看| 久久免费观看电影| 色94色欧美一区二区| 女性生殖器流出的白浆| av有码第一页| 欧美日本中文国产一区发布| 热99久久久久精品小说推荐| 一级爰片在线观看| 国产日韩欧美亚洲二区| 男女午夜视频在线观看| 国产精品久久久久久精品古装| 黄网站色视频无遮挡免费观看| 国产亚洲一区二区精品| 一区二区三区激情视频| 激情视频va一区二区三区| 天堂俺去俺来也www色官网| 免费在线观看视频国产中文字幕亚洲 | 咕卡用的链子| 热re99久久国产66热| 久久久国产一区二区| 男女高潮啪啪啪动态图| videossex国产| 在线 av 中文字幕| 18禁动态无遮挡网站| 如何舔出高潮| 超碰97精品在线观看| 九九爱精品视频在线观看| 亚洲第一av免费看| 99re6热这里在线精品视频| 夫妻性生交免费视频一级片| 久久女婷五月综合色啪小说| 日本免费在线观看一区| 欧美亚洲 丝袜 人妻 在线| av片东京热男人的天堂| √禁漫天堂资源中文www| 欧美少妇被猛烈插入视频| 日日摸夜夜添夜夜爱| av女优亚洲男人天堂| www日本在线高清视频| av国产精品久久久久影院| 欧美国产精品一级二级三级| √禁漫天堂资源中文www| 国产人伦9x9x在线观看 | 自线自在国产av| 国产亚洲av片在线观看秒播厂| 在线观看人妻少妇| 成年动漫av网址| 免费看不卡的av| 捣出白浆h1v1| 久久久久久久久久人人人人人人| 国产成人精品久久久久久| 人成视频在线观看免费观看| 青草久久国产| 极品人妻少妇av视频| av片东京热男人的天堂| 国产精品欧美亚洲77777| 久久久久人妻精品一区果冻| 久久韩国三级中文字幕| 春色校园在线视频观看| 欧美日韩亚洲国产一区二区在线观看 | 国产在线免费精品| 大陆偷拍与自拍| 老鸭窝网址在线观看| 男人操女人黄网站| 精品一区二区免费观看| 亚洲精品国产一区二区精华液| 黄片播放在线免费| 欧美日韩国产mv在线观看视频| 亚洲四区av| 狂野欧美激情性bbbbbb| 如何舔出高潮| 国产免费又黄又爽又色| 熟妇人妻不卡中文字幕| 亚洲四区av| 国产乱来视频区| 伊人久久国产一区二区| 人妻 亚洲 视频| 久久久久精品人妻al黑| 人体艺术视频欧美日本| 久久久久久久久久久久大奶| 欧美日韩视频精品一区| 亚洲成国产人片在线观看| 亚洲欧洲精品一区二区精品久久久 | 深夜精品福利| 亚洲成色77777| 欧美激情高清一区二区三区 | 国产精品国产三级专区第一集| 日本wwww免费看| 又粗又硬又长又爽又黄的视频| 国产人伦9x9x在线观看 | 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产一区二区三区av在线| 日日撸夜夜添| 亚洲精品国产色婷婷电影| 午夜福利视频精品| 下体分泌物呈黄色| 在线亚洲精品国产二区图片欧美| 美女中出高潮动态图| 五月开心婷婷网| 久久99一区二区三区| 男人操女人黄网站| 中文字幕色久视频| 久久青草综合色| 国产精品免费视频内射| 亚洲国产精品国产精品| 少妇被粗大的猛进出69影院| 国产白丝娇喘喷水9色精品| 日韩人妻精品一区2区三区| 又大又黄又爽视频免费| 激情五月婷婷亚洲| 午夜日韩欧美国产| 秋霞在线观看毛片| 亚洲国产毛片av蜜桃av| 熟女av电影| 在线看a的网站| 99热国产这里只有精品6| 久久久精品94久久精品| 最近手机中文字幕大全| 国产成人欧美| 国产成人午夜福利电影在线观看| 中文字幕人妻丝袜一区二区 | 18禁观看日本| 午夜福利一区二区在线看| 国产精品.久久久| 天天躁夜夜躁狠狠久久av| 日日撸夜夜添| 久久97久久精品| 少妇 在线观看| 久久久a久久爽久久v久久| 午夜福利,免费看| 亚洲精品aⅴ在线观看| av免费观看日本| av一本久久久久| 老司机亚洲免费影院| 在线亚洲精品国产二区图片欧美| 免费看不卡的av| 国产免费福利视频在线观看| 国产精品一国产av| 久久久国产一区二区| 日韩免费高清中文字幕av| 国产精品蜜桃在线观看| 18禁国产床啪视频网站| 久久国产精品大桥未久av| 超色免费av| 一边亲一边摸免费视频| 午夜老司机福利剧场| 三上悠亚av全集在线观看| a级片在线免费高清观看视频| 老熟女久久久| 日韩 亚洲 欧美在线| www.自偷自拍.com| 精品人妻在线不人妻| 亚洲人成77777在线视频| 亚洲精品自拍成人| 日韩不卡一区二区三区视频在线| 中文字幕最新亚洲高清| 男人爽女人下面视频在线观看| 伊人久久国产一区二区| 欧美激情极品国产一区二区三区| 午夜日本视频在线| 日韩av免费高清视频| 夜夜骑夜夜射夜夜干| 如何舔出高潮| 亚洲人成网站在线观看播放| 大话2 男鬼变身卡| 日韩一卡2卡3卡4卡2021年| 人人澡人人妻人| 亚洲五月色婷婷综合| 男人添女人高潮全过程视频| 成人免费观看视频高清| 黄片无遮挡物在线观看| 久久久精品94久久精品| 99久国产av精品国产电影| 熟妇人妻不卡中文字幕| 捣出白浆h1v1| 亚洲av中文av极速乱| 国产精品免费大片| 一区二区三区乱码不卡18| 男人舔女人的私密视频| 亚洲精品美女久久久久99蜜臀 | 久久久精品区二区三区| 久久国产精品男人的天堂亚洲| 国产精品.久久久| av在线观看视频网站免费| 国产亚洲午夜精品一区二区久久| 国产成人a∨麻豆精品| a 毛片基地| 欧美日韩国产mv在线观看视频| 日本免费在线观看一区| 国产乱人偷精品视频| 一级毛片 在线播放| 亚洲国产色片| 少妇的丰满在线观看| 国产极品天堂在线| 少妇人妻久久综合中文| 国产男女内射视频| 91精品伊人久久大香线蕉| 777久久人妻少妇嫩草av网站| 亚洲av在线观看美女高潮| 精品少妇黑人巨大在线播放| 免费人妻精品一区二区三区视频| 色哟哟·www| 美女国产高潮福利片在线看| 国产xxxxx性猛交| 国产精品av久久久久免费| 亚洲国产av新网站| 国产成人a∨麻豆精品| 飞空精品影院首页| 宅男免费午夜| a级毛片在线看网站| 免费久久久久久久精品成人欧美视频| 五月天丁香电影| 一级爰片在线观看| 日本免费在线观看一区| 亚洲国产av新网站| 久久国产亚洲av麻豆专区| 如何舔出高潮| 久久精品国产亚洲av涩爱| 日本av免费视频播放| 青春草亚洲视频在线观看| 啦啦啦啦在线视频资源| 18禁国产床啪视频网站| 国产精品偷伦视频观看了| 国产一区二区在线观看av| 精品一区在线观看国产| 日韩欧美精品免费久久| 中文字幕精品免费在线观看视频| 久久青草综合色| 最近手机中文字幕大全| 国产一级毛片在线| 国产成人精品福利久久| 一区在线观看完整版| 久久精品久久精品一区二区三区| 久久人人97超碰香蕉20202| av.在线天堂| 男女国产视频网站| 男人舔女人的私密视频| 国产成人精品福利久久| 最近最新中文字幕免费大全7| 欧美国产精品va在线观看不卡| 99久久人妻综合| 看十八女毛片水多多多| 观看av在线不卡| av在线播放精品| 国产精品一二三区在线看| 亚洲视频免费观看视频| 午夜日本视频在线| 久久久精品94久久精品| 亚洲精品国产av成人精品| 女人高潮潮喷娇喘18禁视频| 人体艺术视频欧美日本| 国产成人91sexporn| 亚洲国产最新在线播放| 国产在线一区二区三区精| 美女高潮到喷水免费观看| 欧美日韩精品成人综合77777| 亚洲精品久久午夜乱码| 999精品在线视频| 亚洲第一青青草原| 一区二区三区精品91| 伦理电影免费视频| 日韩中文字幕欧美一区二区 | 久久99精品国语久久久| 老司机影院毛片| 久久av网站| 亚洲精品日韩在线中文字幕| 我的亚洲天堂| 如何舔出高潮| 免费大片黄手机在线观看| 最新中文字幕久久久久| 9191精品国产免费久久| 91成人精品电影| 亚洲欧美一区二区三区黑人 | 91成人精品电影| 欧美日韩一区二区视频在线观看视频在线| 日韩一卡2卡3卡4卡2021年| 人妻人人澡人人爽人人| 亚洲第一av免费看| 伊人亚洲综合成人网| 久久影院123| 欧美日韩亚洲高清精品| 黑丝袜美女国产一区| 在线观看三级黄色| 人妻少妇偷人精品九色| 两个人免费观看高清视频| 久久99精品国语久久久| 久久人人爽av亚洲精品天堂| 亚洲精品在线美女| 五月开心婷婷网| 亚洲av.av天堂| 一本久久精品| 男人添女人高潮全过程视频| 久久久久久久久久久久大奶| 精品一区在线观看国产| 韩国高清视频一区二区三区| 国产在线视频一区二区| 亚洲av成人精品一二三区| 午夜老司机福利剧场| 久久久精品免费免费高清| 最新的欧美精品一区二区| 欧美精品av麻豆av| 啦啦啦啦在线视频资源| 成人免费观看视频高清| 久久久久视频综合| 最近的中文字幕免费完整| 精品人妻熟女毛片av久久网站| 观看av在线不卡| 精品少妇内射三级| 日本欧美视频一区| 日韩av在线免费看完整版不卡| 天美传媒精品一区二区| 国产爽快片一区二区三区| 人体艺术视频欧美日本| 午夜激情久久久久久久| 亚洲精品成人av观看孕妇| 在线天堂中文资源库| 七月丁香在线播放| 国产黄色视频一区二区在线观看| 精品人妻熟女毛片av久久网站| 美女国产视频在线观看| 亚洲国产精品成人久久小说| kizo精华| 美女主播在线视频| 久久国内精品自在自线图片| 免费看不卡的av| 久久精品国产自在天天线| 亚洲精品国产av成人精品| 18禁国产床啪视频网站| 亚洲av成人精品一二三区| 老司机亚洲免费影院| 十八禁网站网址无遮挡| 精品一区二区三区四区五区乱码 | 熟女av电影| 少妇人妻 视频| 国产黄色视频一区二区在线观看| 黄色一级大片看看| 99久久中文字幕三级久久日本| 秋霞在线观看毛片| 欧美最新免费一区二区三区| 亚洲欧洲日产国产| 婷婷成人精品国产| 夫妻性生交免费视频一级片| 日本猛色少妇xxxxx猛交久久| 999精品在线视频| 国产一区二区激情短视频 | 久久综合国产亚洲精品| 99久久人妻综合| 天天操日日干夜夜撸| 啦啦啦视频在线资源免费观看| 欧美av亚洲av综合av国产av | 国产精品成人在线| 午夜日韩欧美国产| 一区二区日韩欧美中文字幕| 青草久久国产| 免费在线观看视频国产中文字幕亚洲 | 人妻系列 视频| 波野结衣二区三区在线| 又粗又硬又长又爽又黄的视频| 欧美精品亚洲一区二区| 免费看不卡的av| 99久久中文字幕三级久久日本| 国产精品一二三区在线看| 久久久久国产网址| 成人国语在线视频| 亚洲欧美精品综合一区二区三区 | 26uuu在线亚洲综合色| 日韩免费高清中文字幕av| 久久毛片免费看一区二区三区| 少妇被粗大猛烈的视频| 亚洲欧美一区二区三区国产| 亚洲中文av在线| 欧美日韩亚洲国产一区二区在线观看 | 成人国产av品久久久| 亚洲一区二区三区欧美精品| 精品少妇内射三级| 日韩欧美一区视频在线观看| 国产精品一国产av| videossex国产| 国产又色又爽无遮挡免| 国产成人免费观看mmmm| 成年人午夜在线观看视频| 在线精品无人区一区二区三| 欧美日韩亚洲国产一区二区在线观看 | 欧美日韩一级在线毛片| 丰满少妇做爰视频| 蜜桃在线观看..| 十分钟在线观看高清视频www| 久久久精品区二区三区| 亚洲一码二码三码区别大吗| 啦啦啦中文免费视频观看日本| 久久精品熟女亚洲av麻豆精品| 亚洲美女视频黄频| 精品少妇一区二区三区视频日本电影 | 久久精品久久精品一区二区三区| 天天影视国产精品| 亚洲熟女精品中文字幕| 国产精品熟女久久久久浪| www.av在线官网国产| 欧美日韩亚洲国产一区二区在线观看 | 久久青草综合色| 十八禁高潮呻吟视频| 青草久久国产| 不卡av一区二区三区| 在线天堂中文资源库| 欧美日韩亚洲国产一区二区在线观看 | 丝袜美足系列| 精品视频人人做人人爽| 久久这里只有精品19| av女优亚洲男人天堂| 99国产精品免费福利视频| 91国产中文字幕| 午夜91福利影院| 久久av网站| 亚洲精品国产一区二区精华液| 999久久久国产精品视频| 亚洲精品日韩在线中文字幕| 高清欧美精品videossex| 在线天堂中文资源库| 欧美另类一区| 波多野结衣av一区二区av| 国产精品一国产av| 国产精品嫩草影院av在线观看| 黄色一级大片看看| 久久久久人妻精品一区果冻| 国产亚洲精品第一综合不卡| 国产熟女午夜一区二区三区| 欧美日韩精品网址| 亚洲成人av在线免费| 久久综合国产亚洲精品| 一区二区日韩欧美中文字幕| 啦啦啦中文免费视频观看日本| 又粗又硬又长又爽又黄的视频| 国产欧美日韩一区二区三区在线| 精品福利永久在线观看| 少妇猛男粗大的猛烈进出视频| 久久av网站| 亚洲欧洲日产国产| 少妇精品久久久久久久| 最黄视频免费看| 国产一级毛片在线| 最新的欧美精品一区二区| 色婷婷久久久亚洲欧美| 看十八女毛片水多多多| 看免费av毛片| 欧美老熟妇乱子伦牲交| 多毛熟女@视频| tube8黄色片| 精品一区二区免费观看| 九色亚洲精品在线播放| 免费观看av网站的网址| 国精品久久久久久国模美| 久久久精品区二区三区| 深夜精品福利| 欧美少妇被猛烈插入视频| 亚洲国产精品国产精品| av网站免费在线观看视频| 亚洲av电影在线观看一区二区三区| 成人午夜精彩视频在线观看| 日韩一卡2卡3卡4卡2021年| 下体分泌物呈黄色| 色94色欧美一区二区| 国精品久久久久久国模美| 伦精品一区二区三区| 久久ye,这里只有精品| 精品一品国产午夜福利视频| 99久久中文字幕三级久久日本| 亚洲精品一区蜜桃| 精品一区二区三区四区五区乱码 | 老女人水多毛片| 欧美日韩精品成人综合77777| 亚洲av综合色区一区| 这个男人来自地球电影免费观看 | 美女视频免费永久观看网站| 汤姆久久久久久久影院中文字幕| 少妇熟女欧美另类| 亚洲精品自拍成人| 久久人人97超碰香蕉20202| freevideosex欧美| 亚洲欧美中文字幕日韩二区| 99国产精品免费福利视频| 精品一区在线观看国产| 亚洲精品久久成人aⅴ小说| 中国三级夫妇交换| av在线老鸭窝| 亚洲成国产人片在线观看| 97在线人人人人妻| 亚洲欧洲国产日韩| 久久久久久久国产电影| 国产 一区精品| 国产精品久久久久久久久免| 国产极品天堂在线| 免费黄色在线免费观看| 精品福利永久在线观看| 亚洲av免费高清在线观看| 中文字幕人妻丝袜制服| 咕卡用的链子| 大陆偷拍与自拍| 亚洲欧美一区二区三区久久| 大香蕉久久网| 午夜福利视频精品| 久久精品人人爽人人爽视色| 9色porny在线观看| 高清视频免费观看一区二区| 成人漫画全彩无遮挡| 久热这里只有精品99| 国产精品欧美亚洲77777| 日韩欧美一区视频在线观看| 亚洲av福利一区| 欧美97在线视频| 日本猛色少妇xxxxx猛交久久| 日韩精品免费视频一区二区三区| 国产成人精品婷婷| h视频一区二区三区| 男女高潮啪啪啪动态图| 亚洲视频免费观看视频| 国产高清国产精品国产三级| 久久亚洲国产成人精品v| 久久人人爽av亚洲精品天堂| 亚洲,欧美,日韩| av国产精品久久久久影院| 老司机影院毛片| 日本欧美国产在线视频| 日韩人妻精品一区2区三区| 少妇熟女欧美另类| 欧美精品一区二区免费开放| 中文字幕制服av| 青草久久国产| 一区二区三区激情视频| 成人黄色视频免费在线看| 亚洲av综合色区一区| 街头女战士在线观看网站| 国产精品久久久久久久久免| 青草久久国产| 如何舔出高潮| 男人操女人黄网站| 黑人巨大精品欧美一区二区蜜桃| 少妇人妻 视频| 欧美av亚洲av综合av国产av | 国产国语露脸激情在线看| 97在线人人人人妻| 日韩中文字幕视频在线看片| 日韩视频在线欧美| 多毛熟女@视频| 亚洲精品一二三| 丰满迷人的少妇在线观看| 国产精品二区激情视频| 校园人妻丝袜中文字幕| 日韩视频在线欧美| av一本久久久久| 国产在线视频一区二区| 男男h啪啪无遮挡| videosex国产| 免费在线观看视频国产中文字幕亚洲 | 高清视频免费观看一区二区| 看免费av毛片| 女人久久www免费人成看片| 性色avwww在线观看| 大话2 男鬼变身卡| 一二三四在线观看免费中文在| 欧美日韩成人在线一区二区| 黄色 视频免费看| 精品国产一区二区久久| 国产日韩欧美亚洲二区| 99热网站在线观看| 久久精品国产鲁丝片午夜精品| 日韩三级伦理在线观看| 成年动漫av网址| 黑人欧美特级aaaaaa片| 亚洲国产日韩一区二区| 国产免费一区二区三区四区乱码| 欧美日韩一区二区视频在线观看视频在线| 国产片特级美女逼逼视频| av在线观看视频网站免费| 女人被躁到高潮嗷嗷叫费观| 日本vs欧美在线观看视频| 日韩电影二区| 美女xxoo啪啪120秒动态图| 在线天堂中文资源库| 久久人人爽av亚洲精品天堂| 波野结衣二区三区在线| www日本在线高清视频| 日产精品乱码卡一卡2卡三| 三上悠亚av全集在线观看| 中国三级夫妇交换| 国产精品香港三级国产av潘金莲 | 亚洲av国产av综合av卡| 97人妻天天添夜夜摸| 国产福利在线免费观看视频| 最近手机中文字幕大全| 91午夜精品亚洲一区二区三区| 免费看av在线观看网站| 少妇人妻久久综合中文| 日韩精品免费视频一区二区三区| 亚洲国产精品一区三区| 黄色一级大片看看| 2021少妇久久久久久久久久久| 亚洲精品日韩在线中文字幕|