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

    Collapse of cavitation bubbles near air bubbles *

    2020-12-16 02:19:38JingLuoWeilinXuRuiLi

    Jing Luo, Wei-lin Xu, Rui Li

    State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu 610065,China

    Abstract: The cavitation erosion is one of the common damage modes in water engineering. The study of the interaction among cavitation bubbles, air bubbles and the wall is of great significance for understanding the mechanism of the air entrainment to alleviate the cavitation and to enhance the cavitation erosion mitigation effect of aeration. By using the high-speed camera, the regular patterns of the collapse of cavitation bubbles in the vicinity of the wall and the air bubbles are studied in this paper. It is shown that in the vicinity of air bubbles, the cavitation bubbles may only collapse towards or from air bubbles, while under the dual impacts of air bubbles and the wall, the direction of collapse of the cavitation bubbles depends on the combined vector of the impact forces of the air bubbles and the wall. When the air bubbles are very close to the cavitation bubbles, the air bubbles will be penetrated and stretched by the cavitation bubbles, when the distance between them is short enough, the cavitation bubbles and the air bubbles will connect to form gas-type cavitation bubbles, and the collapse strength will be significantly decreased, when there are two air bubbles near the cavitation bubbles, they may even be penetrated successively by the cavitation bubbles. Thus it can be concluded that during the air entrainment to alleviate the cavitation, the number of air bubbles is more important than the total air concentration.

    Key words: Cavitation bubble, air bubble, high-speed photography, air entrainment to alleviate cavitation

    Introduction

    The cavitation is common in many engineering applications, such as large cavitation bubbles generated by underwater explosions, and millimeterscale cavitation bubbles that produce the cavitation erosion in ships, hydraulic machinery and hydraulic energy dissipation structures, as well as ultrasonicinduced millimeter-scale cavitation bubbles used in the biomedicine, the material preparation process and the chemical production to accelerate chemical reaction rates. Wang et al.[1]gave a very comprehensive review of the bubble dynamics and its applications, and pointed out the research difficulties and research directions in this field.

    Philipp and Lauterborn[2]studied the collapse of a cavitation bubble near the wall and it was suggested that the jet plays a less important role in the cavitation damage, while the bubble-wall dimensionless distanceγ(γ=h/Rmax, whereRmaxis the maximal radius of the cavitation bubble,his the distance between the center of the cavitation bubble and the wall)determines whether the cavitation damage takes place or not. It would cause the cavitation erosion whenγis equal to or less than 2, but the jet can be only observed when the cavitation bubble is extremely close to the wall (withγequal to or less than 0.7).Furthermore, there were many other intensive researches of the mechanism of the cavitation erosion[3-5]. Shaw et al.[6]studied the interaction of a laser-generated bubble with a nearby solid boundary,and recorded the pressure increases for the cavities created with sizes between 0.56 and 1.5 times the maximum bubble radius from the solid boundary by a series of pressure transducer traces. However, in order to generate the cavitation bubble for different dimensionless distances, the maximum bubble radius should be varied, while the distance between the cavitation bubble and the wall should be fixed. With the high-speed photographic equipment with a camera speed of up to 108frames per second used to study the evolution properties of the cavitation bubbles around a wall, Lindau and Lauterborn[7]discovered the annular form during the collapse of a cavitation bubble, the formation of a contrajet and the shock wave.

    In practical engineering applications, one sees often not isolated cavitation bubbles, but ones surrounded by other cavitation bubbles, particles and air bubbles. With the coupled effects between the cavitation bubbles,the multiple cavitation bubbles at different nascent moments and different sized cavitation bubbles, the study becomes extremely difficult[8-9]. Rungsiyaphornrat et al.[10]studied the interactions of explosion bubbles formed by underwater explosive of multiple underwater weapons at the same time or with some delay. Fong et al.[11]studied the coupling effect of two cavitation bubbles with different nascent moments by the method of a low voltage discharge. Luo et al.[12]found that the two cavitation bubbles collapse towards each other and then gradually merge, and the merged collapse body quickly moves to the wall. Pain et al.[13]studied jets in quiescent bubbles caused by a nearby oscillating bubble, and found that the jet in the initially stationary bubble can reach a velocity up to 250 m/s. Luo et al.[14]found that the low pressure formed in the opposite direction to the cavitation bubble jet stretches the air bubble in the low-pressure area into a step-like shape.Goh et al.[15]studied the interaction between a cavitation bubble and a non-oscillating air bubble attached to a horizontal polyvinyl chloride plate submerged in the de-ionized water, and found the ratio of the oscillation time and the wall-attached bubble oscillation time is an important parameter for determining the jet direction. The size of the air bubbles plays an outstanding role in the cavitation damage reduction[16]. In addition, the interaction between a cavitation bubble and a particle was researched by adopting a high-voltage discharge technology to induce the cavitation bubble[17]and a low-voltage underwater discharge technology to induce a cavitation bubble[18]. Wu and Gou[19]determined the critical size of the sand particles in the liquid by means of the special vibratory device, and it was found to be related to various effects on the cavitation damage. The bubble-particle interactions were found to be strongly dependent on two dimensionless parameters, the particle-bubble size ratio and the particle-liquid density ratio[20]. Four different interactions of a laser-induced bubble and a free-settling particle were studied by Wu et al.[21], a force balance model was built and it was shown that several non-dimensional factors involving the particle radius, the maximum bubble radius, the initial separation distance , the fluid viscosity and the particle and fluid densities influence the particlebubble dynamics. The interaction between the ice and the cavitation bubble was studied by Cui et al.[22], and the direction of the micro-jet and the propagation rate of the shock waves were captured.

    With the development of the computer technology, the boundary element method (BEM) is widely used in numerical simulations of cavitation bubble dynamics for its inherent nature of reducing the dimension of the problem, and saving computing time, and also for its own unique advantages in dealing with three-dimensional and complex boundary surfaces (such as the cavitation bubble) or the free surface. Klaseboer et al.[23]simulated the ring stage of the cavitation bubble and the vortex surfaces, to study the dynamic characteristics of the cavitation bubble and the interaction with walls for different distances between the bubble and the wall, and the results are in agreement with the experimental ones. Wang and Khoo[24], and many others carried out numerical simulations for the cavitation bubbles.

    The air bubbles are ubiquitous in the water flow in actual engineering. However, the behavior of a cavitation bubble near one or more air bubbles with a wall is not well studied. The focus of this paper is to predict the direction of the collapse of a cavitation bubble under the dual impacts of an air bubble and a wall. In the micro-level, this paper studies the cavitation bubble collapse near the air bubbles, the interaction among the cavitation bubble, the air bubbles and the wall, and the interacted process of the cavitation bubble with double air bubbles near a wall,as a study of the aeration mitigation cavitation.

    1. Experimental setup and method

    As shown in Fig. 1, several complex systems are used to investigate the interactions among cavitation bubbles, air bubbles and the walls, including a single-pulsed high voltage discharge cavitation bubble induced system, an air bubble release system and a high-speed dynamic data acquisition and analysis system.

    Fig. 1 (Color online) The schematic diagram of experimental devices

    When a single-pulsed high voltage discharge cavitation bubble system works, the voltage of the discharge device could reaches 0 kv-50 kv when the electric current passes through a step-up transformer,with the protected water resistance and the rectifier regulator of silicon stack. Once the voltage reaches or exceeds the air breakdown voltage threshold, the circuit is connected and the tungsten electrodes discharge, to immediately generate cavitation bubbles.By adjusting the distance between the discharge devices and the resistance value in the circuit, the desirable cavitation bubbles can be obtained to meet the experiment requirements. The plexiglass is used as the rigid wall material in the interaction experiments of the cavitation bubbles with the wall. Since the air bubbles are required to appear around the cavitation bubbles in a very short time for the study of the interaction of the cavitation bubble with the air bubbles, an air bubble release system is designed to be compatible with the cavitation bubbles. The air bubbles in the water are automatically pushed by the air injection pump to the pinholes (of 0.5 mm in inner diameter) pre-fixed in the water tank (0.4 m×0.3 m×0.4 m), which is filled with deionized water. In the experiment, the air bubble size is controlled by the diameter of the pinholes. In the high-speed data acquisition and analysis system, a certain area of the air bubble rising path is selected to monitor the change of the pixel RGB. If the pixel value changes in the selected area, the high-speed data acquisition and analysis system orders the pre-charged capacitor to discharge, and the high-speed camera to take pictures.The position of the electrode determines the position of the cavitation bubble, and the area in which the pre-selected air bubble passes controls the position of the air bubble and the angle between the air bubble,the cavitation bubble and the wall. Through coordinated work of above-mentioned several systems,we can have air bubbles around the cavitation bubbles.The water temperature keeps constant 22°C during the experiment for avoiding the evident influence of the outside temperature to the evolution of the cavitation bubbles.

    Since the time range of the whole evolution is within a millisecond level, a high-speed dynamic data acquisition and analysis system is needed to study the very fast evolution process. The high-speed dynamic data acquisition and analysis system consists of a high speed camera, the lens assembly, a synchronization device, a lighting equipment, and computer components. The high-speed camera is the main equipment in the whole system, which is the MotionPro Y3-classic high-speed camera (Integrated Design Tools Inc., USA, Maximum data acquisition rate: 105fps). The cavitation bubbles and the air bubbles are video taped with a telephoto microscope lens Navitar (Zoom 6000, Navitar, USA) combined with the high-speed cameras. The high-speed camera’s fast filming speed, rigorous illumination condition for a short exposure time, combined with the illumination requirements of the microscope make more severe requirements to the lighting equipment when these two facilities are used together. Therefore,these experiments need a cold light source (Halogen,150 W) with better lighting and smaller heat.

    2. Collapse of cavitation bubble near an air bubble

    The detailed interaction between the air bubble and the cavitation bubble at different distances can be realized by keeping the constant distance of the discharge needle-plate in the circuit while changing the distance between the port that releases the air bubble and the electrodes. The voltage and the resistance of the four sets of experiments remain unchanged as shown in Fig. 2, Theoretically, the maximum radius of the cavitation bubbles in the four sets of experiments are the same, however, the maximum radius of the cavitation bubbles in groups c and d vary a lot compared with groups a and b. It is caused by the distance between the air and the cavitation bubbles, and the fact that the cavitation bubble has merged with the air bubbles in the expansion phase (see Figs. 2(c3) and 2(d3)). During the fusion of the two bubbles, the cavitation bubble contains a large amount of non-condensable gas, to form a gas-type cavitation bubble. In Fig. 2, the maximum expansion radiusRmaxof the cavitation bubble is selected as the characteristic parameter of the cavitation bubble,ε, the ratio of the initial air bubble radius to the characteristic parameter of the cavitation bubble, is the dimensionless air bubble radius,dis the minimum distance of the center of the cavitation bubble to the surface of the air bubble,δis the ratio of the nearest distance between the cavitation bubble center and the air bubble surface to the characteristic parameter of the cavitation bubble.

    The first row of pictures in Fig. 2 show the air bubble positions from the electrode center before the cavitation bubble is incipient. In the second row of pictures, the bright white light is the water discharge electrode when the air bubbles in group a is farther away from the electrode center as compared with those in group d. The light intensity is strongest in the electrode center, and weaker outward. The bubbles in groups a and b have expanded to the maximum volume before the line 4, while those in groups c and d reach their maximum volume at the line 5 and line 6,respectively. At this stage, the bubbles in groups a and b farther away from the electrode centers are penetrated by the main shock wave generated by the expansion of the cavitation bubbles. Compared with that in group b, the bubble in group a have later breakdown time point (respecitvely Figs. 2(a3) and 2(b2)). For the bubbles in groups c and d, due to the close distance between the cavitation bubble and the air bubble, the cavitation bubble in group c breaks in the expansion phase much earlier, while the bubble in group d is fully integrated with the air bubble, to become a gas-type cavitation bubble. In group c, some air bubbles are split in the wake of the main shock wave and for the same reason, some of them, as precursors, form some breakdown split bubbles, as shown in c4, c5and c6.

    Fig. 2 (Color online) The interaction between the air bubbles and the cavitation bubbles at different distances (Filming speed:19 700 fps, Exposure time: 48 μs, Frame size: 17.33 mm×7.36 mm)

    After the contraction phase, the space is quickly released due to the cavitation bubble contraction and rapidly filled with the surrounding water, meanwhile,the air bubbles move slightly under the flow effect of the cavitation bubbles, as shown in group b.Compared with group a, the air bubbles of group b move more apparently to the center of the cavitation bubbles. The air bubbles are in dumbbell shape under the double actions of the collapse shock and the backfill flow, and the air bubbles in group a is farther away from the center of the cavitation bubbles as compared with the bubbles in group b, with a more obvious dumbbell shape. During the collapse process of the bubbles in groups a and b, the cavitation bubbles do not move to the air bubble, and the bubbles of both groups generate resilient bubbles after the complete collapse of the cavitation bubbles.Compared with the bubbles in groups a and b, the cavitation bubbles in groups c and d have longer expansion-collapse period. Under the same experimental conditions, this phenomenon is due to the fact that the air bubbles are closer to the cavitation bubble center, part or even all of the non-condensable gases of the air bubbles are mixed into the cavitation bubbles to form the gas-type cavitation bubbles in the expansion or collapse phase. As can be seen from the photos, the more non-condensable gas is mixed in the expansion phase, the longer the expansion-collapse period of the cavitation bubble will be.

    Through the study of the interaction between the air bubbles and the cavitation bubbles at different distances, it can be found that with the gradually decreased distance between the cavitation bubbles and the center of the air bubbles, the air bubbles are more likely to break down in the expansion phase of the cavitation bubble. When the distance between the cavitation bubble and the air bubble is short enough,the cavitation bubble absorbs the air bubble to form a gas-type cavitation bubble. By generating the cavitation bubbles by means of underwater lowvoltage discharge, and changing the position of the cavitation bubbles to interact with the hemispherical bubbles adhered to the underwater plates, Goh et al.[15]found that when the bubble-wall dimensionless distance is in the range 1.5<γ<4.5, the direction of the jet flow has the tendency of facing to or deviating from the collapse. In this paper, the cavitation bubbles in Fig. 2(a) are collapsed away from the air bubble,while the cavitation bubbles in Fig. 2(b) are collapsed in situ. The energy of the cavitation is mainly concentrated in the first two expansion and contraction processes. Therefore, the critical bubblewall distance of the fusion in these two stages needs to be considered in combination with the shape and the size of the cavitation bubbles, In the experimental method of this paper, since the air bubble is free-moving, the distance between the air bubble and the cavitation bubble cannot be precisely controlled.As a result, we fail to obtain the critical condition for the fusion of free-moving air bubbles and cavitation bubbles.

    2.1 Effect of air bubble on the collapse direction of cavitation bubble

    A single cavitation bubble collapsing near a solid boundary is a well-studied phenomenon[1-6]. Under certain conditions, a liquid jet directed towards the boundary is formed during the collapse of the cavitation bubble. A large number of studies indicate that the wall strongly attracts the cavitation bubble in the collapse phase.

    Figure 3 shows the interaction of the cavitation bubble with an air bubble: The black area in each picture on the left side is the air bubble, and the cavitation bubble is on the right side. The maximum radiiRmaxof the cavitation bubbles, respectively, are 2.76 mm and 3.35 mm, the dimensionless sizesεof the air bubbles are 0.39 and 0.28, and the minimum distancesδbetween the centers of the cavitation bubble and the surface of the air bubble are 2.29 and 2.35.

    The pictures a12-a25and b13-b25in Fig. 3 show two cases of the collapse directions of the cavitation bubble, collapsing away from and to the air bubble,and it can be said that the air bubbles have the repulsion and attraction tendency for the collapse of the cavitation bubbles.

    Fig. 3 (Color online) Repulsion and attraction of air bubbles on the cavitation bubbles (Filming speed: 197 00 fps, Frame size:17.33 mm×7.36 mm)

    Under certain conditions, two different situations can be considered regarding to the bubble-bubble interaction: the cavitation bubble collapses towards and away from the air bubble. Furthermore, the rigid wall has an attraction to the collapse of the cavitation bubbles, while the air bubbles show certain repulsion and attraction to the collapse of the cavitation bubbles.By generating cavitation bubbles by means of underwater low-voltage discharge, and changing the position of the cavitation bubbles to interact with the hemispherical bubbles adhering to the underwater plates, Goh et al.[15]found that the dimensionless cavitation oscillation timeT′ is an important parameter that determines the collapse direction of the cavitation bubble (T′ is the ratio of the time in the first cycle of the cavitation bubble to the inverse of the resonance frequency of the air bubble). WhenT′> 1 .2, the cavitation bubble collapses toward the air bubble. WhenT′<1.0, the cavitation bubble collapses away from the air bubble and the wall.WhenT′≈1.0, the cavitation bubble collapses and splits into two parts, to produce a jet toward the air bubble and a jet away from the air bubble. The interaction between the cavitation bubble and the air bubble studied in this paper has two main features different from those in the references. First, with the experimental method used in this paper, the motion of the air bubbles in the water is not affected by the outside world. Second, the shape of the air bubble interacting with the cavitation bubble is not a regular sphere. The air bubble generated by the experimental method used in the Ref. [15] is a hemispherical and still air bubble, similar to that in this paper whenθ=0 (θis the angle between the center line of two bubbles and the normal to the boundary). Therefore,these sections will study the interactions, respectively,when the air bubble locates between the cavitation bubble and the wall, and when the cavitation bubble locates between the air bubble and the wall.

    In Fig. 4, a1and a2show the collapse directions of the cavitation bubble when the cavitation bubble locates between the wall and the air bubbles, a3, a4show the collapse directions of the cavitation bubbles when the air bubble locates between the wall and the cavitation bubble. As can be seen in experimental results, the cavitation bubbles move away from the air bubbles pushed by the repulsion forces of the air bubbles, meanwhile, the wall has attraction to the cavitation bubbles, which eventually collapse in the combined vector direction under the combined action of the two forces. The picture b1, b2in Fig. 4 present the collapse direction of the cavitation bubbles when the cavitation bubbles locate between the wall and the air bubbles, while picture b3, b4present the collapse direction of the cavitation bubbles when the air bubbles locate between the wall and the cavitation bubble. Similar to the experimental results in group a,the ultimate collapse direction of the cavitation bubbles depends on the combined vector direction when the air bubbles have an intensified attraction to the cavitation bubbles and the rigid wall. At present,our experiments are limited by the cavitation bubble induction methods and measurement techniques. The attraction of the wall towards the cavitation bubbles and the attraction or rejection of the air bubbles to the cavitation bubbles need further exploration, but the mechanism of the collapse direction of the cavitation bubbles under the action of the air bubbles and the walls is well studied.

    2.2 Effect of air bubble on the collapse strength of cavitation bubble

    Fig. 4 (Color online) Vector synthesis of the cavitation bubbles,air bubbles and the wall

    Figure 5 shows the high-speed photographs of the cavitation bubbles which absorb and merge with the air bubbles in the interaction process under the same experimental conditions. In Fig. 5, the dimensionless radius of the air bubbles relative to the cavitation bubblesεranges from 0.21 to 0.22, the maximum radiusRmaxof the cavitation bubble varies from 3.57 mm to 4.06 mm, and the nearest dimensionless distance between the cavitation bubble center and the air bubble surfaceδranges from 0.72 to 1.01.The first column in Fig. 5 is the relative position of the electrode center and the air bubbles before the nascent cavitation bubbles. The cavitation bubbles are formed at the fourth column in each group experiment,while the cavitation bubbles in this process are in the expansion stage until the sixth column. At this moment, the air bubbles are broken down by the main shock generated by the cavitation bubble expansion.For instance, the air bubble in group b, due to the smallest dimensionless radius of the air bubble in the four groups of experiments, is not only broken down,but also generates a split part, (see Fig. 5(b6)). With a further expansion of the cavitation bubbles, the intensity of the main shock wave is reduced, and the outline shape of the air bubbles begins to grow due to the strong compressibility of the non-condensable gas inside the air bubbles, as shown in pictures from column 7 to column 9. When the cavitation bubbles are in the contraction phase as shown in pictures from column 10 to column 17, and the surrounding water fills the space released by the cavitation bubbles at once due to the rapid contraction, and the air bubbles are stretched by the action of the backfill. The main part of the air bubbles begins to grow from the tabular shape in the buffer stage into the slender shape, and with a further contraction of the cavitation bubbles,the water in the jets from the bubbles originally broken down moves faster to the center of the cavitation bubbles than the periphery water. At this moment, the air bubbles take tabular shapes at the original bubble breakdown position, mainly due to the uneasily stretched nature of the peripheral water with no air bubbles inside on one hand. And on the another hand, the strong compressibility of the noncondensable gas of the air bubbles makes themselves easily to be pulled to the center of the cavitation bubbles to take the tabular form. As can be seen from Fig. 4, sinceδin the group b is greater than that of the other three groups in the experiments, under the action of the backfill flow, the tabular shapes of the air bubbles develops significantly slower than the other three groups, which have merged with the cavitation bubbles in groups b, c and d in the later collapse.

    From the expansion and the collapse of the cavitation bubbles and the morphology evolution of the air bubbles, it can be found that the surface of the air bubbles near the cavitation bubbles becomes concave under the action of the main shock wave,then the air bubbles break down to form a jet (see Fig.6(a)). Similar observations were made by Pain et al.[13]and Goh et al.[15]. The air bubble is stretched by the backfill water in the collapse process of the cavitation bubbles, and its surface near the center of the cavitation bubble becomes stepped (see Fig. 6(b)).From the shape analysis of the air bubbles at different stages, it can be found that the air bubbles are stretched into a stepped shape, which largely buffers the contraction and the collapse of the cavitation bubbles. On the other hand, the cavitation bubbles merge with the air bubbles, to form gas-type cavitation bubbles.

    In the high-speed flow, the cavitation bubbles near the wall merge with the air bubbles during the development process when they are close enough, to have an incomplete collapse. The merged cavitation bubbles, containing non-condensable gas, have different collapse characteristics, as compared with other cavitation bubbles.

    Fig. 5 (Color online) The breakdown jet and the backfill elongational flow (Filming speed: 19 700 fps, Exposure time: 48 μs,Frame size: 17.33 mm×7.36 mm)

    Fig. 6 (Color online) Air bubble is penetrated and elongated

    Figure 7 shows the merging process when the air bubble is between the cavitation bubble and the wall.The distance between the center of the cavitation bubble and the wall is 5.17 mm and 5.48 mm, and the dimensionless distanceγ(the ratio of the distance from the center of the cavitation bubble to the wall to the characteristic parameter of the cavitation bubble)between the centers and the wall is 1.68 and 1.78, both within the attraction range of the wall to the cavitation bubbles. Moreover, the bubble-boundary angleθis selected as the measure of the geometrical arrangement of a cavitation bubble and an air bubble near a rigid boundary. In Fig. 7, with the gradual increased volume of the cavitation bubble, the air bubble breaks down, to form a breakdown jet, as shown in a3and b3.The breakdown jet moves to the wall, and eventually rolls out along the flat wall by the blocking effect of the wall. In the contraction process of the cavitation bubble, the peripheral water quickly moves to the center of the cavitation bubbles, due to the breakdown air bubbles in contact with the wall, while the wall shows a strong adsorption to the air bubbles, then the air bubbles are stretched by the peripheral backfill flow, as shown in a8and b8. The cavitation bubble merges with the deformed air bubble in the collapse phase, to form a gas-type cavitation bubble, (see a12and b12in Fig. 7).

    Figure 8 shows the merging process when the cavitation bubble locates between the air bubble and the wall. In the experiments, a constant distance is kept between the center of the cavitation bubble and the wall, and the size of the air bubbles is kept as 1.03 mm and 0.92 mm, and the dimensionless distances between the centers of the cavitation bubbles and the wall are 2.28 and 2.59.

    In the two groups of experiments, the air bubble breaks down under the action of the main shock wave radiated by the cavitation bubble in the expansion phase, to form a breakdown jet flow. The length of the jet flow gradually decreases in the contraction phase of the cavitation bubble, and the air bubble is absorbed by the cavitation bubble in the collapse phase, with each cavitation bubble forming a gas-type cavitation bubble. Since the gas-type cavitation bubble contains a large amount of non-condensable gas, the cavitation bubble will not completely collapse, which would greatly relieve for the wall the destructive effects of the cavitation bubble collapse.

    Figure 9 is the scatter plot of the maximum radiusRmaxof the cavitation bubbles with regard to dBmax(dBmaxis the peak noise level of the collapse of the cavitation bubble) in cases with and without the air bubbles. By adjusting the position of the air bubble releasing port over the holder in the experiment, the air bubbles can move freely to the water surface by the buoyant effect. Therefore, the cavitation bubble would interact with the air bubble in any position on the floating trajectory. After fixing the air bubble releasing port, the cavitation bubbles are generated by multiple discharging, in the experiments of the cavitation bubbles and the air bubbles of different distances.

    In Fig. 9, it is seen that the noise level peak dBmaxradiated by the collapsed cavitation bubbles continues to rise with the gradual increase of the maximum radiusRmaxof the cavitation bubbles, but its growth trend gradually slows down. On the other hand, compared with the cases of no air bubbles, the noise peak level radiated by the collapsed cavitation bubbles gradually gets weaker when there is an air bubble around the cavitation bubble. Especially for the cavitation bubbles with the maximum radiusRmaxfrom 1.00 mm to 3.00 mm, the presence of the air bubbles will attenuate the noise peak level dBmaxmost significantly, and the impact of the air bubbles to the noise peak level gets softened with the gradually increased radiusRmaxof the cavitation bubble.

    Fig. 7 (Color online) The merger when the air bubble is between the cavitation bubble and the wall (Filming speed: 19 700 fps,Exposure time: 48 μs, Frame size: 17.33 mm×7.36 mm)

    Fig. 8 (Color online) The merger when the cavitation bubble is between the air bubble and the wall (Filming speed: 19 700 fps,Exposure time: 48 μs, Frame size: 17.33 mm×7.36 mm)

    Figure 9 can be explained as follows: Firstly, the stretched air bubble, when in its contraction phase,could attenuate the powerful contraction of the cavitation bubble to some degree due to its compressibility.Secondly, when the air bubble and the cavitation bubble are close enough, the cavitation bubble merges with the air bubble into a gas-type cavitation bubble containing a large amount of non-condensable gas.The gas-type cavitation bubbles will have a non-complete collapse and largely decrease the collapse intensity. What is more, the air bubble size affects the buffering effect when the air bubble does not merge with the cavitation bubble. While forming gas-type bubbles, different sizes of the air bubble mean different amounts of non-condensable gas,which leads to varying collapse strength.

    3. Collapse of cavitation bubble near two air bubbles

    Figure 10 is the high-speed photographs of the interaction process of the cavitation bubbles with double air bubbles. The size of the air bubbles and the cavitation bubbles, and the relative positional relationship of the two are shown in Fig. 10, with the subscript 1 representing the air bubbles on the left side,and the subscript 2 representing the air bubbles on the right side.

    In the expansion phase of the cavitation bubbles,whether the cavitation bubble locates on the same side of the air bubble or not, the double air bubbles near the cavitation bubbles will have an inner concave shape to varying degrees under the effect of the main shock, and the closer the air bubble to the cavitation bubble, the greater the inner concave level will be. In the contraction phase, the air bubble near the cavita-tion bubble forms a breakdown jet, and with further development, it will be detached from the air bubbles to form a single tiny bubble, as shown in a8and b6in Fig. 10.

    Fig. 9 (Color online) The maximum radius Rmax of cavitation bubbles and the corresponding peak values of collapse noises

    By generating the cavitation bubbles by means of underwater low-voltage discharge, Pain et al.[16]studied the jet phenomenon of the air bubbles wrapped inside the silicone oil under the action of the cavitation bubbles, and found that, the dimensionless distance from the center of the cavitation bubble to the center of the air bubble must be less than 2 to generate jets under the action of the cavitation bubbles.However, in Fig. 10(a), the distance from the air bubble to the closest cavitation bubble is greater than 2.33 (the bubble-bubble center distance is used as the reference, and the minimum distance from the cavitation bubble center to the air bubble surface is used in this paper), the jet can still be observed in the experiments in this paper. The main reason may be the difference in the method of obtaining the air bubbles in the experiment. The method used in the reference is to wrap the air bubble in silicone oil droplets, and in this paper a free-floating air bubble is generated in the water. Figure 11 is the high-speed photographs of the interaction process of the cavitation bubbles with double air bubbles near the wall. The size of the bubbles and the relative positional relationship of the two are shown in Fig. 11, with the subscript 1 representing the air bubbles on the left side and the subscript 2 representing the air bubbles on the right side.

    Fig. 10 (Color online) The interaction of cavitation bubbles with double air bubbles (Filming speed: 19 700fps, Exposure time:48 μs, Frame size: 26.1 4mm×5.84 mm)

    Fig. 11 (Color online) The interaction of cavitation bubbles near the wall with double air bubbles at the same side. (Filming speed:19 700 fps, Exposure time: 48μs, Frame size: 13.86 mm×7.36 mm)

    As can be seen from Fig. 11, in the expansion stage of the cavitation bubbles, the two air bubbles experience a contraction deformation. As the radius of the cavitation bubble is further increased, the air bubbles near the cavitation bubbles first generate the jet, followed by the air bubbles near the wall with the jet flow being eventually developed towards the wall,as shown in a5, b5and c4. After the cavitation bubbles entering into the collapse phase, the space occupied by the cavitation bubbles is released, and the backfill flow is rapidly formed by the surrounding water. At this moment, the air bubbles are stretched into irregular shapes, while the air bubbles close to the wall is still attached to the wall.

    4. Conclusions

    With the experimental system of using high voltage pulse sparks to generate cavitation bubbles and a high-speed dynamic data acquisition and analysis system, a large number of experiments are conducted to study the interactions among the cavitation bubbles, the air bubbles and the wall. From the experimental results the following conclusions can be made:

    (1) In the vicinity of the air bubbles, the cavitation bubbles may only collapse towards or from the air bubbles, while under the dual impacts of the air bubbles and the wall, the cavitation bubbles are subjected also to the attraction of the wall , along with the attraction or repulsion of the air bubbles, the final direction of the collapse of the cavitation bubbles depends on the combined vector of such two forces.

    (2) If the cavitation bubbles are very close to the air bubbles, the cavitation bubbles in the expansion stage will penetrate the air bubbles, and the cavitation bubbles in the collapse stage will stretch the air bubbles. When the distance between the two is short enough, the cavitation bubbles and the air bubbles will connect to form gas-type cavitation bubbles, and then the collapse noise peak level will be significantly decreased, suggesting a significantly weakened collapse strength.

    (3) By observing the interaction of two air bubbles and cavitation bubbles near the wall, it is found that the cavitation bubbles in the development process may penetrate the two air bubbles successively from the near to the distant ones.

    The above results show that during the aeration and cavitation erosion mitigation, the number of air bubbles is more important than the total air concentration. These results have potential theoretical values for the anti-erosion design in hydraulic engineering, the erosion resistance of blades in hydraulic machinery,and the better use of the cavitation effect in ultrasonic cleaning technology.

    成年av动漫网址| 精品福利永久在线观看| 多毛熟女@视频| 久久99蜜桃精品久久| 午夜影院在线不卡| 亚洲成人av在线免费| 日本色播在线视频| av国产精品久久久久影院| 美女国产高潮福利片在线看| 日韩精品免费视频一区二区三区 | a级毛片在线看网站| 一本色道久久久久久精品综合| 亚洲色图综合在线观看| 啦啦啦啦在线视频资源| 精品人妻在线不人妻| 欧美xxⅹ黑人| 国产亚洲欧美精品永久| 美国免费a级毛片| 老熟女久久久| 亚洲经典国产精华液单| 少妇人妻久久综合中文| 日韩欧美精品免费久久| 一个人免费看片子| 99国产综合亚洲精品| 久久久久精品久久久久真实原创| 国产免费视频播放在线视频| 内地一区二区视频在线| 久久久久国产网址| 草草在线视频免费看| 国产免费视频播放在线视频| 成年av动漫网址| 久久鲁丝午夜福利片| 国产精品偷伦视频观看了| 又粗又硬又长又爽又黄的视频| 一本久久精品| 日本欧美视频一区| 成人毛片a级毛片在线播放| 亚洲精品色激情综合| 一级毛片我不卡| 纯流量卡能插随身wifi吗| 久久久久久久久久成人| 国产黄色免费在线视频| 久久热在线av| 国产乱来视频区| 国产午夜精品一二区理论片| 欧美精品人与动牲交sv欧美| 成人影院久久| 婷婷成人精品国产| 国产在线视频一区二区| 69精品国产乱码久久久| 国产1区2区3区精品| 国产精品久久久久久久久免| www.av在线官网国产| 亚洲色图 男人天堂 中文字幕 | 日本欧美视频一区| 纵有疾风起免费观看全集完整版| 欧美另类一区| 欧美精品亚洲一区二区| 欧美变态另类bdsm刘玥| 欧美精品一区二区免费开放| 少妇的逼好多水| 久久精品国产鲁丝片午夜精品| 大香蕉久久网| 亚洲精品国产av蜜桃| 蜜桃国产av成人99| 成人二区视频| 亚洲国产毛片av蜜桃av| 国产亚洲精品第一综合不卡 | 欧美精品人与动牲交sv欧美| 极品少妇高潮喷水抽搐| 女的被弄到高潮叫床怎么办| 丝袜喷水一区| 久久精品国产亚洲av涩爱| 国产在线视频一区二区| 插逼视频在线观看| 亚洲天堂av无毛| 亚洲av福利一区| 久久久久久久精品精品| 国产精品嫩草影院av在线观看| 三上悠亚av全集在线观看| 欧美激情 高清一区二区三区| 国产av码专区亚洲av| 免费在线观看黄色视频的| 中文精品一卡2卡3卡4更新| 卡戴珊不雅视频在线播放| 少妇 在线观看| 亚洲第一av免费看| 少妇精品久久久久久久| 人妻系列 视频| 亚洲欧洲精品一区二区精品久久久 | 国产精品三级大全| 如日韩欧美国产精品一区二区三区| 久久国产精品大桥未久av| 久久毛片免费看一区二区三区| 在线亚洲精品国产二区图片欧美| 婷婷色综合www| 中文字幕最新亚洲高清| 亚洲欧洲精品一区二区精品久久久 | 丁香六月天网| 久久久久久久亚洲中文字幕| 少妇的逼水好多| 久久精品人人爽人人爽视色| 超色免费av| 精品国产一区二区久久| 久久午夜综合久久蜜桃| 日本欧美国产在线视频| 国产无遮挡羞羞视频在线观看| 国产日韩欧美在线精品| 赤兔流量卡办理| 亚洲成人av在线免费| 欧美日韩精品成人综合77777| 色哟哟·www| 精品国产国语对白av| a级毛色黄片| 十八禁高潮呻吟视频| 成年女人在线观看亚洲视频| 国产探花极品一区二区| 精品福利永久在线观看| 少妇人妻久久综合中文| 国产成人av激情在线播放| 精品一区二区三区视频在线| 女的被弄到高潮叫床怎么办| 18+在线观看网站| 99视频精品全部免费 在线| 午夜福利网站1000一区二区三区| 一级,二级,三级黄色视频| 中文欧美无线码| 久久久久视频综合| 啦啦啦中文免费视频观看日本| 在现免费观看毛片| 国产精品99久久99久久久不卡 | 中文天堂在线官网| 搡女人真爽免费视频火全软件| 久久热在线av| 国产精品久久久久久精品电影小说| 亚洲少妇的诱惑av| 国产毛片在线视频| 国产福利在线免费观看视频| 国产精品无大码| 日韩三级伦理在线观看| 成人影院久久| 午夜免费鲁丝| 在线观看国产h片| 亚洲,一卡二卡三卡| 亚洲精品久久午夜乱码| 久久久久网色| 国产欧美日韩一区二区三区在线| 国产一区有黄有色的免费视频| 欧美激情国产日韩精品一区| 老司机影院毛片| 99精国产麻豆久久婷婷| 欧美精品一区二区免费开放| 亚洲成国产人片在线观看| 免费观看性生交大片5| 51国产日韩欧美| 久久精品国产鲁丝片午夜精品| 美女中出高潮动态图| 国产黄频视频在线观看| 最近最新中文字幕大全免费视频 | 成人国语在线视频| 国产高清国产精品国产三级| 亚洲国产最新在线播放| 欧美 日韩 精品 国产| 宅男免费午夜| 日本黄大片高清| 人人妻人人澡人人看| 国产精品女同一区二区软件| 久久久久精品性色| 妹子高潮喷水视频| 美女中出高潮动态图| 色5月婷婷丁香| 免费观看在线日韩| 69精品国产乱码久久久| 乱码一卡2卡4卡精品| 香蕉精品网在线| 另类亚洲欧美激情| 永久免费av网站大全| 中国三级夫妇交换| 一二三四中文在线观看免费高清| 亚洲国产色片| tube8黄色片| 国产成人精品在线电影| 超色免费av| 777米奇影视久久| 中文字幕av电影在线播放| 日韩 亚洲 欧美在线| 精品国产一区二区三区四区第35| 国产成人精品福利久久| 亚洲精品久久久久久婷婷小说| 欧美国产精品一级二级三级| 中文字幕精品免费在线观看视频 | 男女边摸边吃奶| 亚洲精品色激情综合| 国产视频首页在线观看| 午夜激情av网站| 午夜免费观看性视频| 99视频精品全部免费 在线| 日韩制服骚丝袜av| 成人漫画全彩无遮挡| 97人妻天天添夜夜摸| 久久国产精品男人的天堂亚洲 | 天美传媒精品一区二区| 亚洲精品第二区| 考比视频在线观看| 1024视频免费在线观看| 激情五月婷婷亚洲| 成年人午夜在线观看视频| 亚洲第一av免费看| 精品国产露脸久久av麻豆| 少妇被粗大猛烈的视频| 日韩成人av中文字幕在线观看| 超碰97精品在线观看| 免费看不卡的av| 交换朋友夫妻互换小说| 亚洲一级一片aⅴ在线观看| 欧美人与性动交α欧美软件 | 26uuu在线亚洲综合色| 天美传媒精品一区二区| 成人二区视频| 大香蕉久久网| av片东京热男人的天堂| 免费观看无遮挡的男女| 亚洲精品中文字幕在线视频| 人人妻人人爽人人添夜夜欢视频| 在线观看免费日韩欧美大片| 日韩成人av中文字幕在线观看| 欧美成人精品欧美一级黄| 五月开心婷婷网| 亚洲美女黄色视频免费看| 国产日韩一区二区三区精品不卡| 欧美人与性动交α欧美软件 | 看非洲黑人一级黄片| 中国国产av一级| 99久久精品国产国产毛片| 一个人免费看片子| 欧美日韩综合久久久久久| 99久久中文字幕三级久久日本| 午夜老司机福利剧场| 丝瓜视频免费看黄片| 岛国毛片在线播放| 久久精品国产自在天天线| 日本爱情动作片www.在线观看| 国产精品久久久久久久久免| 在线天堂中文资源库| 免费在线观看黄色视频的| 午夜福利视频精品| a级毛色黄片| 人妻一区二区av| 成人黄色视频免费在线看| 久久人人爽av亚洲精品天堂| 日韩一区二区三区影片| a级片在线免费高清观看视频| 欧美日韩成人在线一区二区| 精品第一国产精品| 蜜桃国产av成人99| 国产免费现黄频在线看| 亚洲综合色惰| videossex国产| 精品亚洲成a人片在线观看| 国产av码专区亚洲av| 美女国产视频在线观看| 一二三四中文在线观看免费高清| 熟妇人妻不卡中文字幕| 国产成人精品在线电影| av福利片在线| 啦啦啦中文免费视频观看日本| 纯流量卡能插随身wifi吗| 亚洲中文av在线| 看非洲黑人一级黄片| 9191精品国产免费久久| 国产乱人偷精品视频| 国产一级毛片在线| 久久久久精品人妻al黑| 日本与韩国留学比较| 少妇猛男粗大的猛烈进出视频| videos熟女内射| 亚洲一码二码三码区别大吗| 中文精品一卡2卡3卡4更新| av女优亚洲男人天堂| 亚洲国产精品999| 少妇被粗大猛烈的视频| 咕卡用的链子| 欧美97在线视频| freevideosex欧美| 日日爽夜夜爽网站| 成人综合一区亚洲| 亚洲av中文av极速乱| 国产精品三级大全| 国产精品一区www在线观看| 亚洲人成网站在线观看播放| 国产黄频视频在线观看| 精品少妇久久久久久888优播| 亚洲图色成人| 人人妻人人澡人人爽人人夜夜| xxx大片免费视频| 国产精品99久久99久久久不卡 | av国产久精品久网站免费入址| 日韩视频在线欧美| 一级片免费观看大全| 久久人人爽av亚洲精品天堂| 内地一区二区视频在线| 日本-黄色视频高清免费观看| 精品久久国产蜜桃| 黄网站色视频无遮挡免费观看| 久久鲁丝午夜福利片| 18禁裸乳无遮挡动漫免费视频| 亚洲精品日本国产第一区| a级片在线免费高清观看视频| 26uuu在线亚洲综合色| 99九九在线精品视频| 天天影视国产精品| 欧美日韩成人在线一区二区| av在线播放精品| 亚洲伊人色综图| 欧美xxⅹ黑人| 精品第一国产精品| 高清在线视频一区二区三区| 日本免费在线观看一区| 精品久久国产蜜桃| 999精品在线视频| 香蕉丝袜av| √禁漫天堂资源中文www| 亚洲国产欧美日韩在线播放| 久久久国产精品麻豆| 777米奇影视久久| 热re99久久精品国产66热6| 波野结衣二区三区在线| av又黄又爽大尺度在线免费看| 狂野欧美激情性bbbbbb| 欧美人与善性xxx| 亚洲精品色激情综合| 亚洲欧美日韩另类电影网站| 亚洲成av片中文字幕在线观看 | 亚洲欧美一区二区三区黑人 | 国精品久久久久久国模美| 老女人水多毛片| 日韩人妻精品一区2区三区| 青青草视频在线视频观看| 精品人妻一区二区三区麻豆| 日韩av在线免费看完整版不卡| 久久韩国三级中文字幕| 一个人免费看片子| 蜜臀久久99精品久久宅男| 大片免费播放器 马上看| 国产精品久久久av美女十八| 久久免费观看电影| 久久人人爽av亚洲精品天堂| 黑人巨大精品欧美一区二区蜜桃 | 国产精品嫩草影院av在线观看| 欧美激情 高清一区二区三区| 视频区图区小说| 亚洲精品日本国产第一区| 赤兔流量卡办理| 欧美bdsm另类| 精品少妇内射三级| 精品一品国产午夜福利视频| 日韩中字成人| 卡戴珊不雅视频在线播放| 久久狼人影院| 免费在线观看黄色视频的| 一本久久精品| a级片在线免费高清观看视频| 精品一区二区免费观看| 欧美人与性动交α欧美软件 | 男女高潮啪啪啪动态图| 国产精品久久久久久久久免| 精品熟女少妇av免费看| 校园人妻丝袜中文字幕| 国产成人精品一,二区| 国产极品天堂在线| 啦啦啦视频在线资源免费观看| 天天影视国产精品| 亚洲美女黄色视频免费看| 美国免费a级毛片| av视频免费观看在线观看| 亚洲美女黄色视频免费看| 春色校园在线视频观看| 国产又色又爽无遮挡免| 少妇人妻 视频| 老司机影院成人| 欧美老熟妇乱子伦牲交| 草草在线视频免费看| 51国产日韩欧美| 成人亚洲精品一区在线观看| 一级a做视频免费观看| 女人被躁到高潮嗷嗷叫费观| 波野结衣二区三区在线| 一级毛片电影观看| 自线自在国产av| 伊人亚洲综合成人网| 国产综合精华液| 新久久久久国产一级毛片| 18在线观看网站| 午夜老司机福利剧场| 亚洲国产精品999| av又黄又爽大尺度在线免费看| 欧美 日韩 精品 国产| 纵有疾风起免费观看全集完整版| 亚洲,一卡二卡三卡| 成年动漫av网址| 菩萨蛮人人尽说江南好唐韦庄| 曰老女人黄片| 午夜福利,免费看| 国产精品人妻久久久影院| 国产 精品1| 国产成人精品福利久久| 国产亚洲午夜精品一区二区久久| 美女视频免费永久观看网站| 国产av码专区亚洲av| www.熟女人妻精品国产 | 国产精品一国产av| 在现免费观看毛片| 国产成人aa在线观看| 91精品三级在线观看| 国产激情久久老熟女| 国产日韩欧美亚洲二区| 这个男人来自地球电影免费观看 | 99精国产麻豆久久婷婷| 中文字幕最新亚洲高清| 一级爰片在线观看| 欧美bdsm另类| 欧美成人午夜精品| 免费av不卡在线播放| 伦理电影免费视频| 免费黄频网站在线观看国产| 中文字幕人妻熟女乱码| 亚洲国产色片| 国产熟女欧美一区二区| 色视频在线一区二区三区| 丁香六月天网| 中文字幕人妻熟女乱码| 97人妻天天添夜夜摸| 国产高清不卡午夜福利| 国产综合精华液| 80岁老熟妇乱子伦牲交| 99久久中文字幕三级久久日本| 女的被弄到高潮叫床怎么办| 一级片免费观看大全| av电影中文网址| 成年人免费黄色播放视频| 亚洲精品国产色婷婷电影| 国产色婷婷99| 草草在线视频免费看| 另类精品久久| 22中文网久久字幕| 一级毛片我不卡| 不卡视频在线观看欧美| 一级片'在线观看视频| 午夜福利视频在线观看免费| 国产熟女午夜一区二区三区| 欧美精品亚洲一区二区| 欧美最新免费一区二区三区| 国产av国产精品国产| 午夜免费观看性视频| 最近最新中文字幕免费大全7| 国产爽快片一区二区三区| 最新的欧美精品一区二区| 人妻少妇偷人精品九色| 熟女电影av网| 丰满迷人的少妇在线观看| 国产又爽黄色视频| 国产成人91sexporn| 另类精品久久| 亚洲精品自拍成人| 五月天丁香电影| 精品国产露脸久久av麻豆| 黄色毛片三级朝国网站| 亚洲欧美成人精品一区二区| 有码 亚洲区| 国产午夜精品一二区理论片| 999精品在线视频| 免费观看在线日韩| 国产成人a∨麻豆精品| 免费看光身美女| 久久久久人妻精品一区果冻| 精品视频人人做人人爽| 人人妻人人澡人人看| 飞空精品影院首页| 免费人妻精品一区二区三区视频| 水蜜桃什么品种好| 免费女性裸体啪啪无遮挡网站| 亚洲高清免费不卡视频| 黄色配什么色好看| 综合色丁香网| 亚洲性久久影院| 久久热在线av| 国产成人午夜福利电影在线观看| 久久人人爽人人爽人人片va| 午夜福利视频在线观看免费| 国产成人免费观看mmmm| 久久精品国产鲁丝片午夜精品| 免费在线观看黄色视频的| 亚洲国产欧美日韩在线播放| 久久99蜜桃精品久久| 欧美xxxx性猛交bbbb| 精品人妻偷拍中文字幕| 高清欧美精品videossex| 免费高清在线观看日韩| 久久 成人 亚洲| 九草在线视频观看| 精品第一国产精品| 国产精品不卡视频一区二区| 九色成人免费人妻av| 美女中出高潮动态图| 大码成人一级视频| 亚洲精品一区蜜桃| 视频区图区小说| 18禁观看日本| 哪个播放器可以免费观看大片| 麻豆精品久久久久久蜜桃| 天天躁夜夜躁狠狠久久av| av不卡在线播放| 免费女性裸体啪啪无遮挡网站| 99久久精品国产国产毛片| 亚洲av电影在线进入| 日日爽夜夜爽网站| 久久精品国产亚洲av涩爱| 黄片播放在线免费| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 中文字幕另类日韩欧美亚洲嫩草| 在线观看免费高清a一片| 在线天堂中文资源库| 久久青草综合色| 捣出白浆h1v1| av在线播放精品| 国产在线视频一区二区| 精品少妇久久久久久888优播| 欧美激情国产日韩精品一区| 综合色丁香网| kizo精华| 欧美性感艳星| 熟女av电影| 亚洲国产精品一区三区| 国产av一区二区精品久久| 国产精品欧美亚洲77777| 亚洲国产看品久久| 国产亚洲午夜精品一区二区久久| 午夜免费观看性视频| 午夜久久久在线观看| 午夜影院在线不卡| 午夜激情久久久久久久| 久久精品熟女亚洲av麻豆精品| 一级片'在线观看视频| 成年动漫av网址| 欧美另类一区| 亚洲欧美成人精品一区二区| 久久综合国产亚洲精品| 亚洲精品美女久久av网站| 日日爽夜夜爽网站| 夫妻性生交免费视频一级片| 日韩伦理黄色片| 精品亚洲乱码少妇综合久久| 久久久久久伊人网av| 精品亚洲成国产av| 日韩中字成人| 九九爱精品视频在线观看| 久久精品国产自在天天线| 狠狠精品人妻久久久久久综合| 亚洲国产av影院在线观看| av黄色大香蕉| 五月玫瑰六月丁香| 美女中出高潮动态图| 97在线视频观看| 欧美激情 高清一区二区三区| 丰满迷人的少妇在线观看| 国产探花极品一区二区| 欧美日韩视频高清一区二区三区二| 尾随美女入室| 黄色怎么调成土黄色| 成人免费观看视频高清| 精品国产国语对白av| 国产精品久久久久久久电影| 在现免费观看毛片| 国产永久视频网站| 国产精品国产三级国产专区5o| 国产极品天堂在线| 丝袜脚勾引网站| 国产精品 国内视频| 国产日韩欧美视频二区| 亚洲色图 男人天堂 中文字幕 | 两个人免费观看高清视频| 夫妻午夜视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲美女黄色视频免费看| 狠狠精品人妻久久久久久综合| 亚洲国产精品专区欧美| 一个人免费看片子| 日本猛色少妇xxxxx猛交久久| 午夜免费男女啪啪视频观看| 人人妻人人添人人爽欧美一区卜| 免费大片黄手机在线观看| 国产成人免费无遮挡视频| 99视频精品全部免费 在线| 久久久久网色| 国产极品天堂在线| 中文字幕免费在线视频6| 一本—道久久a久久精品蜜桃钙片| 国产精品麻豆人妻色哟哟久久| 午夜久久久在线观看| 中国美白少妇内射xxxbb| 亚洲精品久久午夜乱码| 欧美最新免费一区二区三区| 大香蕉97超碰在线| 久久免费观看电影| 成人二区视频| 日韩制服骚丝袜av| 成人综合一区亚洲| 中文字幕最新亚洲高清| 国产探花极品一区二区| 欧美日韩国产mv在线观看视频| 高清黄色对白视频在线免费看| 久久久久国产精品人妻一区二区| 国产精品偷伦视频观看了| 18禁观看日本| 精品久久蜜臀av无|