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

    Experimental study on the particle flow patterns in a cyclone dipleg with a trickle valve

    2020-06-22 06:05:16ChaoYuYanGuoGangSunJianYiChenYaoDongWeiChunXiLu
    石油科學(xué)通報(bào) 2020年3期

    Chao?Yu Yan · Guo?Gang Sun · Jian?Yi Chen · Yao?Dong Wei · Chun?Xi Lu

    Abstract An experimental apparatus including a dipleg and a trickle valve was established to simulate the operation of a suspended dipleg-trickle valve system of cyclone used in fluid catalytic cracking (FCC) unit. The flow regimes in the dipleg and the discharge modes in the trickle valve were studied by combining the observation of experimental phenomena with the analysis of transient pressure fluctuation. The results show that the flow regimes in the dipleg have two types–the dilute–dense phase coexisting falling flow and the dilute falling flow. Correspondingly, the trickle valve also has two discharge modes–the intermittent periodic dumping discharge and the continuous trickling discharge. The power spectrum density of pressure fluctuation displays that the gas–solids flow in the dipleg-trickle valve system is characterized by a low-frequency pulsation. The coherence coefficient explains the origin and propagation of pressure fluctuation in the system. Eventually, a map describing the flow regimes and discharge modes related to the operation parameters was proposed, which can provide a helpful guidance for the operation of cyclone dipleg-trickle valve system in FCC unit.

    Keywords Cyclone dipleg · Trickle valve · Flow regime · Discharge mode · Pressure fluctuation

    1 Introduction

    Fluid catalytic cracking (FCC) is an important re finery process for converting heavy oil into more valuable gasoline and lighter products. A typical property of FCC unit is that the catalyst particles are required to keep a continuous and steady circulation flow along a designed loop between reactor and regenerator system. For the purpose, the catalyst particles should be efficiently separated from hydrocarbon vapors in the reactor and from flue gas in the regenerator by cyclone separators and be conveyed back to fluidized beds through diplegs. Usually, two stages of cyclone separators(primary and secondary cyclone) in series are installed in the freeboard above the dense bed of disengager and regenerator, respectively, as shown in Fig. 1. A dipleg is a pipe connected with the bottom hopper of a cyclone separator, which is used to convey the catalysts back to the dense bed from a low-pressure top section of the dipleg to a high-pressure bottom region by particle gravity, as well as maintaining a seal to prevent the gas from flowing into the dipleg (gas leakage). If large quantity of gas flows up into the cyclone due to the gas leakage, the separation efficiency of the cyclone will decrease, and the loss of catalyst particles will happen. Therefore, the operation status of the dipleg has a great effect on the cyclone performance.

    Fig. 1 Schematic diagram of FCC unit

    Generally, due to different inlet solids loading in gas–solids mixture flowing through the two stages of cyclones,the operation status of primary cyclone dipleg is usually different from that of the secondary dipleg. For the primary cyclone dipleg, over 99% of the separated particles flow down the primary cyclone dipleg, which has enough momentum to overcome negative pressure and prevent from gas leakage as well (Reddy Karri and Knowlton 2004; Geldart et al. 1993). So, the bottom end of primary cyclone dipleg is often opened without any valve and submerged in the dense bed (submerged dipleg) (Smolders et al. 2001) as shown in Fig. 1.

    However, for the secondary cyclone dipleg, less than 1% of the separated particles flow through the secondary cyclone dipleg, which does not have enough seal ability to prevent gas leakage (Reddy Karri and Knowlton 2004; Geldart et al. 1993). Therefore, a trickle valve, which had been proposed by DuconCompay and widely used in many FCC units at present, is usually needed to attach to the bottom end of the secondary dipleg to minimize the gas leakage during the startup and other running periods of cyclone separator.And the trickle valve of the secondary dipleg can operate at a status of either submerged in the dense bed or suspended at the freeboard region above dense bed surface (Smolders et al. 2001). For more safety, the trickle valve is often submerged in the dense bed of regenerator (here called as submerged dipleg-trickle valve as shown in Fig. 1). However,due to the space limitation caused by the existence of bottom stripper in the disengager, the secondary dipleg-trickle valve cannot be inserted into the dense bed and often only can be suspended in the freeboard (here called as suspended diplegtrickle valve as shown in Fig. 1).

    The production practice of FCC unit shows that some trouble shooting of catalysts loss could ascribe to the malfunction of the suspended dipleg-trickle valve operation(Niccum 2010). The probable malfunction of suspended secondary dipleg-trickle valve in the disengager is more easily happened than that of submerged dipleg-trickle valve in the regenerator by some causes, such as the gas leaking and the erosion of sealing flapper surface (Shaw and Walter 2007).Therefore, for the sake of the reliable design and operation of the suspended secondary cyclone dipleg-trickle valve in FCC disengager, the research on the flow regimes in the dipleg and the discharge modes in the trickle valve should be carried out systematically.

    The FCC cyclone dipleg is really an over flow standpipe(Smolders and Baeyens 1995; Bouma 1998). Although the FCC cyclone dipleg is a very simple vertical pipe, the gas–solids two-phase flow in it is very complicated. Many factors, including the effect of swirling flow at the bottom outlet of cyclone, the negative pressure gradient flow of particles, and the capacity of solids conveying, make the FCC cyclone dipleg either different from the conventional vertical standpipe or the downcomer of circulating fluidized beds(Li et al. 1997). The studies on conventional standpipe and downcomer have many literature reports about the explanation of flow regimes and flow instability (Leung and Wilson 1973; Leung and Jones 1978; Li 2005; Jing and Li 1999;Zhang and Rudolph 1998; Mountziaris and Jackson 1991).However, in spite of the wide application of the cyclone dipleg in FCC and other industry units, little research work has been done about it so far.

    The basic flow parameters (such as particle velocity, particle concentration, and pressure drop) in dipleg have the characteristics of dynamic nonlinear fluctuation, which leads to the complexity of flow regimes and flow instability in dipleg; for example, there are more than one flow regime coexisting in dipleg (Hoffmann and Stein 2002; Cortés and Gil Cortés and Gil 2007). Dries and Bouma (1997) divided the flow regimes into four regions along the axial direction of dipleg, inlet swirling flow, dense flow, fluidized flow, and packed flow. Wang et al. (2000a, b) thought that the dipleg could generally operate in one of three basic flow regimes,packed bed, fluidized bed flow, and a dilute flow called streaming flow. The experimental results by Smolders et al.(2001), Geldart and Radtke (1986) and Wei et al. (2003,2004) showed that the primary cyclone dipleg in FCC was operated under a dense conveying flow regime due to a high solids flux; however, the secondary cyclone dipleg was operated mainly under a dilute–dense coexisting flow regime due to a low solids flux and a high negative pressure. In recent literatures, the study by Wei et al. (2016) showed that the flow regimes in the primary cyclone dipleg could be divided into two types based on the solids flux and negative pressure,the dilute–dense coexisting falling flow, and the dense conveying flow along the dipleg. However, these research results are mainly focused on the submerged dipleg. The study on the suspended dipleg-trickle valve is rarely involved.

    In fact, for the suspended dipleg-trickle valve, the discharge modes in trickle valve are closely related to the flow regimes in dipleg. So, the relationship between the flow regimes in the dipleg and the discharge modes in the trickle valve should be studied simultaneously with a same set of experimental apparatus. However, the research work is very little despite their importance to FCC processes so far. In the literatures that could be retrieved, only several researchers reported their findings. Bristow and Shingles(1989) proposed four discharge modes of trickle valve,the trickling, the dumping, the trickling/dumping, and the flooding. Geldart et al. (1993), Sun et al. (2001) and Reddy Karri and Knowlton (2004) analyzed the operation of trickle valve. Especially, a blowing mode of trickle valve will occur when operating with positive pressure. In these trickle valve operation studies, Bristow and Shingles (1989) and Geldart et al. (1993) used batch solids transfer system that would operate continuously for a certain length of time. However,the investigation by Sun et al. (2001) and Reddy Karri and Knowlton (2004) used continuous solids transfer system,which is more capable of showing the evolution of the discharge modes.

    The gas–solids flow patterns in dipleg are usually de fined by the analysis of flow parameters (axial distributions of the pressure, particle velocity, or concentration) (Geldart et al.1993; Gil et al. 2002), which are mainly affected by operating conditions (solids flux and negative pressure), physical properties of particles (particle density, particle size distribution, etc.), and geometric structure of dipleg (Smolders and Baeyens 1995; Wang et al. 2000a, b; Geldart and Radtke 1986; Knowlton et al. 1986). These flow parameters have the characteristics of dynamic change, which may be attributed to the gas pulsation, the particle agglomeration and dispersion, the vortex generation and shedding, the bubble formation and eruption, the expansion and compression of gas flow impacted by particles, and the friction between particles and bed wall (Bi 2007; Pun?ochá? and Draho? 2005; van der Schaaf et al. 1998; Fan et al. 1981).

    Comparatively speaking, the nonintrusive measurement of pressure fluctuation is much easier than the intrusive measurement of particle velocity and concentration. And the pressure fluctuation signal data contain sufficient fluidization information implied in gas–solids two-phase flow.So, the measurement of pressure fluctuation has been widely used in many studies of gas–solids fluidization. And several types of data processing methods have been applied to analyze the pressure fluctuation in gas–solids two-phase flow,such as the statistical methods (standard deviation, probability density function, etc.) based on time domain (Yang and Leu 2008; van Ommen et al. 2011; Lin et al. 2001) and the Fourier transform methods (power spectrum density, coherence function, etc.) based on frequency domain (Fan et al.1981; van der Schaaf et al. 1999; Kage et al. 2000; Shou and Leu 2005). The characteristic values of the pressure fluctuation can be extracted by using these signal data processing methods, explaining the fluctuation mechanism of gas–solids two-phase flow in the diplegs. The operation of diplegtrickle valve has the characteristics of dynamic change of flow parameters because the slugs of solids flow out of the dipleg and slugs of gas back up into the trickle valve (Sinclair et al. 1983). But the previous literature reports mainly focused on the static analysis of statistical mean values of pressure data and less on the dynamic change analysis.

    Therefore, in this work, an experimental apparatus of the dipleg-trickle valve system was established to simulate the operation of the suspended secondary cyclone dipleg-trickle valve in FCC disengager, which usually operates at a relatively lower solids flux and a higher negative pressure. The transient pressure fluctuation at different axial positions of the whole dipleg-trickle valve system was measured by pressure transducer under the wide operating conditions (solids flux and negative pressure). The flow regimes in the dipleg and the discharge modes in the trickle valve were identi fied by combining the experimental phenomena with the analysis of transient pressure fluctuation by using the data processing methods. The mechanism about the origin and propagation of pressure fluctuation in the system was also further discussed based on a dynamic analysis.

    Fig. 2 Schematic diagram of experimental apparatus

    Fig. 3 Photographs of experimental apparatus

    2 Experiments

    2.1 Experimental apparatus

    Figures 2 and 3 show the schematic diagram and the photographs of the experimental apparatus, respectively. In order to clearly observe the experimental phenomena of trickle valve, a batch solids transfer system was used in the experiment, which could operate continuously for a certain length of time. The experimental apparatus was mainly composed of a top hopper, a dipleg, a trickle valve, and a bottom rectangular collector for particles. The dipleg had an inner diameter of 150 mm and a height of 5000 mm. The trickle valve with a nominal diameter of DN 150 mm was connected with the bottom end of the dipleg and suspended in the rectangular collector, simulating the case of FCC secondary cyclone dipleg-trickle valve system. The trickle valve was mainly composed of a short inclined pipe (or called as the miter bend) and a flapper plate. The top section of flapper plate was connected with the inclined pipe by two rotatable hand rings and two davit bolts as shown in Figs. 2 and 3. The flapper plate can rotate flexibly back and forth around the two-hand rings. The installation angle between the flapper plate and the vertical line at the trickle valve exit was 5° for generating the natural closure force.In order to simulate the operation of industrial trickle valve,the weight and the size of the flapper plate employed in the experiment were same as the industrial flapper plate. For clearly observation of experimental phenomena, the dipleg,the rectangular collector, and the inclined pipe of trickle valve were made of transparent plexiglas. And the other parts of the experimental apparatus including the flapper plate were made of carbon steel.

    The solid particles employed in the experiment was the equilibrium FCC catalyst with a Sauter mean diameter of 67 μm, a particle density of 1500 kg/m3, and a bulk density of 940 kg/m3. The particle size distribution of FCC catalysts is displayed in Fig. 4.

    In order to simulate the work environment with negative pressure gradient for particle conveying and discharging operation along the dipleg-trickle valve system in FCC disengager, the pressurized air from a compressor was introduced into the rectangular collector to keep a certain operation pressure and form the negative pressure gradient (or pressure difference) between inside and outside of the dipleg-trickle valve system during the experiments. The pressure in the space of rectangular collector was measured by a U-tube manometer. First, the catalyst particles were fed into the top hopper and then flowed down the dipleg and fell into the rectangular collector through the trickle valve. The solids flux in the dipleg was controlled by a butter fly valve under the hopper. The solids flux was determined by simultaneously measuring the accumulating height in the dipleg and the falling time of catalyst particles when closing the butter fly valve installed at the bottom section of the dipleg. In these experiments, the solids flux was varied within 50.0 kg/(m2s) based on the cross-sectional area of the dipleg, and the operation range of negative pressure (or pressure difference)between inside and outside of the dipleg-trickle valve system was within 11.0 kPa.

    In the experiments, six pressure taps along the axial direction of dipleg were arrayed based on the zero height at the exit bottom of trickle valve, namely 0.3 m, 0.9 m, 1.4 m, 2.1 m,2.7 m, 3.3 m, as shown in Fig. 2. The transient pressure fluctuation at each tap was measured by pressure transducer. The sampling frequency of pressure transducer was 333 Hz, and the maximum sampling time was 180 s. The transient pressure fluctuation signals in the six measuring levels were collected simultaneously in the experiments, and the experimental phenomena were also recorded by a digital camera to monitor the flow regimes and discharge modes in the dipleg-trickle valve system.

    Fig. 4 Particle size distribution of FCC catalysts

    2.2 Data processing method

    2.2.1 Power spectrum density analysis

    The power spectrum density (PSD) is a method to analyze the energy content of a signal associated with to the frequencies of a physical phenomenon. The power spectrum density function,Pxx, of a time series signal,x(t), is the Fourier transform of its autocorrelation function,Rxx, as written by Eqs. (1) and (2).The function expresses the behavior of a signal in frequency domain rather than in time domain (Fan et al. 1981; van der Schaaf et al. 1999; Kage et al. 2000; Shou and Leu 2005). In the PSD analysis of pressure fluctuation, the frequency corresponding to the maximum energy of transient pressure fluctuation signal is called as the dominant frequency, which re flects the overall pulsation frequency of gas–solids two-phase flow.So, the PSD analysis of signal can extract the peculiar information re flecting the gas–solids flow behavior in frequency domain.

    2.2.2 Coherence analysis

    The coherence analysis has been used to measure the similarity of two signals at a frequency f (Cai et al. 1990). The coherence coefficientcan be de fined as follows.

    where Pxx(f)and Pyy(f)are the auto-power spectrum of signals Px(t)andPy(t), respectively.Pxy(f)is the cross-power spectrum of signals Px(t)andPy(t).Pxx(f),Pyy(f), and Pxy(f)are de fined as follows.

    In Eq. (3), the coherence coefficient,ranges from 0 to 1, corresponding to fully uncorrelated and fully correlated time series at frequency f, respectively. In this work,the coherence coefficienthas been used to discern the different sources and propagation of pressure fluctuations in the dipleg-trickle valve system.

    3 Results and discussion

    3.1 Experimental phenomena

    3.1.1 Flow regimes in the dipleg

    During the experiments, the flow phenomena in the dipleg were recorded by a digital camera. The flow regimes based on the experimental phenomena in the dipleg can be divided into two patterns, namely the dilute–dense phase coexisting falling flow and the dilute falling flow, as shown in Fig. 5,which depends on the solids flux and the negative pressure.

    Fig. 5 Diagram of flow regimes in the dipleg

    At a given negative pressure, when the solids flux is relatively lower, the dilute–dense phase coexisting falling flow appears in the dipleg, as shown in Fig. 5a. There exists a certain height of dense phase in the bottom of the dipleg,and a dilute phase falling flow locates at the section between the inlet of dipleg and the surface of dense phase bed. As the top inlet of the dipleg employed in this work is not connected with a cyclone, but a hopper, there is no swirling flow in the upper dilute phase falling flow region of the dipleg. In the dilute phase falling flow region, the particles flow downward like rain. However, in the bottom dense phase flow region,the particles move downward at a low velocity. There are some bubbles rising in the dense phase flow, which are produced by gas from the rectangular collector through the exit of trickle valve. So the particle phase moves downward and gas phase moves upward in the dense phase flow region,which causes this region resembling a bubbling fluidized bed. The interface between the dilute falling flow and the dense flow moves up and down unsteadily. The dense phase bed in the bottom section of dipleg has the sealing effect preventing the bypass of large quantity of gas outside of the dipleg. Therefore, the height of the dense phase bed is called as the sealing height of dipleg.

    The experimental phenomena show that when the negative pressure is relatively lower, the sealing height will conversely decrease with increasing the solids flux, and the flow regime becomes the dilute falling flow as shown in Fig. 5b.Due to the relatively higher downward particle velocity at the low solids flux, the dilute flow will occur unlike the dense conveying flow reported by Wei et al. (2016) for primary cyclone dipleg with high solids flux. However, at a given solids flux, the increases in negative pressure will increase the sealing height, inducing the regime transition from the dilute falling flow to the dilute–dense phase coexisting falling flow in the dipleg. Therefore, depending on the operating conditions, the two flow regimes in the dipleg can transform from each other.

    3.1.2 Discharge modes in the trickle valve

    Corresponding to the flow regimes in the dipleg, the trickle valve has two discharge modes in the operation process,namely an intermittent periodic dumping discharge and a continuous trickling discharge, as shown in Fig. 6.

    At a certain solids flux and negative pressure, if the flow regime in the dipleg is in the status of dilute–dense phase coexisting falling flow, the trickle valve will work as the mode of intermittent periodic dumping discharge. When the dense phase bed in the dipleg reaches to a certain height with an enough static pressure head to overcome the force produced by the pressure impacting on the outside surface of the flapper plate, the flapper plate will open abruptly from a close status to a large opening one, slapping the exit of trickle valve intensively. The particles accumulated in the dipleg are discharged as a dumping way in a short time, as shown in Fig. 6a. Meanwhile, the gas leakage occurs along the upper semi-circumference of the exit, which flows from the space of rectangular collector into the trickle valve. The particles pass through the inclined pipe of trickle valve as a strati fied flow, as shown as the collection of photographs in Fig. 7.

    Fig. 6 Diagrams of discharge modes of trickle valve

    The gas flows into the dipleg to form as the bubbles,which makes the dense bed have properties of a bubbling fluidized bed. The surface of the dense bed fluctuates up and down intensively. The static pressure head of the dense bed is also gradually decreased at the same time. At that moment, the static pressure head of the dense bed cannot counteract the force impacting on the outside surface of the flapper plate, and the flapper plate will close. The discharge process is completed in one time. Then, the particles start to accumulate in the dipleg again to form the dense height as before. The trickle valve will discharge particles in the next cycle time. Therefore, the trickle valve runs as a mode of intermittent periodic dumping discharge. Generally, this type of discharge is corresponding to the operating conditions with relatively higher negative pressure and lower solids flux, as well as to the dilute–dense phase coexisting falling flow regime in the dipleg.

    When the negative pressure is relatively lower and the solids flux is relatively higher, the trickle valve will operate as a mode of continuous trickling discharge as shown in Fig. 6b, which is closely related to the dilute falling flow regime in the dipleg. For the flow regime as mentioned above, the conveying particles with a relatively higher velocity impact on the inside surface of the flapper plate, generating enough driving force to balance the resistant force impacting on the outside surface of the flapper plate. Then,the flapper plate keeps open to discharge particles as a continuous and trickling flow mode. In the process of discharge,the flapper plate always keeps open with a slight swing back and forth if the operating conditions are not changed. The maximum opening angle of the flapper plate is about 7° in this case as observed in the experiments, which is slightly larger than the installation angle of 5° of the flapper plate. As the particle flow rate in the continuous trickling discharge is very low, the gas from the space of the rectangular collector can leak into the trickle valve around the almost entire exit circumference of trickle valve, which also results in a stratified flow of particles in the inclined pipe of trickle valve.

    3.2 Pressure fluctuation in the dipleg?trickle valve system

    3.2.1 Time mean pressure pro files

    Figure 8 displays the time mean pressure profiles in the dipleg for different flow regimes. At a given negative pressure (7 kPa as shown in Fig. 8), when the solids flux is gradually increased, the distributions of the time mean pressure clearly vary from the non-uniformity to the uniformity along the axial direction of dipleg, which demonstrates that the flow regime in the dipleg transforms from the dilute–dense phase coexisting falling flow to the dilute falling flow.

    When the solids flux is relatively lower, the axial distribution of the time mean pressure appears an obvious in flection point, as shown in Fig. 8a, b. Above the critical height corresponding to the in flection point (or the dilute–dense surface),the distribution of time mean pressure is very uniform in the dilute phase flow. However, below the critical height (or in the dense phase bed), the time mean pressure is increased gradually with the axial position in close to the bottom of the dipleg. This also means that the concentration above the critical height is much lower than the one below the critical height, indicating that the dilute–dense phase coexisting falling flow occurs in the dipleg.

    Figure 8c, d shows that, when the solids flux increases,the axial distribution of time mean pressure becomes relatively more uniform along the dipleg except a slightly higher in the bottom region in close to the exit of trickle valve, due to the resistance of trickle valve. This type of axial distribution of time mean pressure or the axial distribution of particle concentration corresponds to the dilute falling flow.

    Fig. 7 Photographs of intermittent periodic dumping discharge of trickle valve at different time in a discharge cycle

    3.2.2 Transient pressure fluctuation pro files

    The flow behavior of solids in the dipleg-trickle valve system is a dynamic process, which can be described by pressure fluctuation parameter. Figure 9 displays the typical transient pressure fluctuation pro files in the dipleg under different operating conditions. It can be found that the pressure fluctuation pro files, especially in close to the exit of trickle valve,have signi ficant difference under different operating conditions, which can be used to identify the two flow regimes and discharge modes in the dipleg-trickle valve system as mentioned above.

    In Fig. 9a, the pressure fluctuation pro files in close to the exit of trickle valve are composed of zigzag peaks with low frequency and high amplitude, such as for the axial level of 0.3 m and 0.9 m. In the other axial positions far from the exit of the trickle valve (from 1.4 to 3.3 m), the pressure fluctuation pro files change gently with high frequency and low amplitude. The higher the axial position, the lower the pressure.

    The pressure fluctuation pro files in Fig. 9a are corresponding to the dilute–dense phase coexisting falling flow regime for the dipleg and the intermittent periodic dumping discharge mode for the trickle valve. For the pressure fluctuation pro files at the axial level in closest to the exit of the trickle valve (h= 0.3 m), the gradual ascending part of the zigzag peak with time means the accumulating process of particles in the dipleg when the flapper plate is in a closed state. The static pressure head keeps increasing up to a wave peak due to the formation of a dense bed with a certain height. And then, at the wave peak, the force impacting on the inside surface of flapper plate can balance the resultant force acting on the outside surface of flapper plate. The flapper plate starts to open and discharge particles, decreasing the dense bed level and the static pressure head, which is re flected by the descending part of the zigzag peak pro file.It can be found that one-time discharge cycle is about 50 s for the operating condition in Fig. 9a. Therefore, the zigzag peak of pressure fluctuation located at the position in closest to the exit of trickle valve can be employed to identify the dilute–dense phase coexisting falling flow regime and the intermittent periodic dumping discharge in the dipleg-trickle valve system.

    Figure 9b gives the pressure fluctuations corresponding to the dilute falling flow regime for the dipleg and the continuous trickling discharge mode for the trickle valve. It shows that the transient pressure pro file at each axial level has the similarity basically with high frequency and low amplitude due to the stable discharge of particles in this case. However, due to the slight swing back and forth for the flapper plate, the fluctuation extent of pressure pro files at the axial levels (h= 0.3 m, 0.9 m, 1.4 m) in close to the exit of trickle valve is slightly higher than those at other levels (h= 2.1 m, 2.7 m, 3.3 m), but obviously distinguish from those in the dilute–dense phase coexisting falling flow regime and the intermittent periodic dumping discharge mode.

    3.2.3 Effect of operating conditions on pressure fluctuationIn the experiments, it can be found that the operating condition has the obvious effects on the flow regime and the discharge mode in the dipleg-trickle valve system. Figure 10displays the effects of the solids flux and the negative pressure on the pressure fluctuation profiles at the location(h= 0.3 m) in closest to the exit of trickle valve.

    Fig. 8 Time mean pressure pro files for different flow regimes in the dipleg

    As shown in Fig. 10a, at a given solids flux (Gs= 17.4 kg/(m2s)), when the negative pressure is decreased from 5.1 to 2.1 kPa, the pressure pro files gradually change from the zigzag peaks with low frequency and high amplitude into the smooth fluctuations with high frequency and low amplitude.This clearly indicates that the flow regime and the discharge mode are also changed, respectively, transforming from the dilute–dense phase coexisting falling flow into the dilute falling flow regime for the dipleg and from the intermittent periodic dumping discharge into the continuous trickling discharge mode for the trickle valve.

    On the contrary, as shown in Fig. 10b, at a given negative pressure (ΔP= 7.05 kPa), decreasing the solids flux from 42.0 to 10.5 kg/(m2s) will also change the pressure pro files from the smooth fluctuations into the zigzag peaks. In this case, the dilute falling flow will evolve into the dilute–dense phase coexisting falling flow regime for the dipleg. Meanwhile, the continuous trickling discharge will transform into the intermittent periodic dumping discharge mode for the trickle valve.

    Fig. 9 Typical transient pressure fluctuation pro files for different flow regime and discharge mode in the dipleg-trickle valve system

    3.3 Power spectrum density of pressure fluctuation

    The time series of pressure signals can be processed by using the Fourier transform method in frequency domain to obtain the power spectrum density (PSD), which is used to extract the dominant frequency of pressure fluctuation signals, re flecting the dynamic feature of pressure fluctuation.

    Fig. 10 Effects of operating conditions on pressure fluctuation (h = 0.3 m)

    Fig. 11 Power spectrum density pro files of pressure fluctuations in the dipleg

    Figure 11 shows the PSD pro file of the transient pressure fluctuation at each axial level of the dipleg under different operating conditions. Figure 11a displays the PSD pro files for the dilute–dense phase coexisting falling flow regime and the intermittent periodic dumping discharge mode. Figure 11b shows the PSD pro files corresponding to the dilute falling flow regime and the continuous trickling discharge mode.

    The dominant frequency of pressure signals is closely related to the flow regime and the discharge mode. Figure 11a shows that the PSD pro files have a similarity and an obvious wave crest with a dominant frequency of about 0.01 Hz. The power energy corresponding to the dominant frequency at the axial position (h= 0.3 m) in closest to the exit of trickle valve is relatively higher than that at other axial positions. However, in Fig. 11b, each of the PSD profiles has a dominant frequency of about 0.03 Hz and several sub-frequencies. All the PSD pro files overlap well from each other, indicating the consistency of pressure fluctuation propagation in the dipleg.

    Comparing Fig. 11a with b, it can be found that the power energy corresponding to the dominant frequency in the intermittent periodic dumping discharge mode (and the dilute–dense phase coexisting falling flow regime) is much higher than that in the continuous trickling discharge (and the dilute falling flow regime). Therefore, the flow regimes and the discharge modes in the dipleg-trickle valve system can also be identi fied by using the PSD analysis of the transient pressure fluctuation.

    Fig. 12 Distribution of coherence coefficient of pressure fluctuation between two adjacent positions along the axial height of dipleg

    3.4 Coherence analysis of pressure fluctuation

    As mentioned above, the PSD analysis of the transient pressure fluctuation re flects the fluctuation power under a dominant frequency, discriminating the flow regimes and the discharge modes in the dipleg-trickle valve system. However, this method cannot determine the propagation direction of pressure fluctuation which is very important to explain the dynamic features of gas–solids flow in the system. In this work, the coherence function analysis method is used to determine the similarity or difference of pressure fluctuation signals between two measurement positions (as described in Sect. 2.2.2), so as to determine the propagation direction of pressure fluctuation. If the coefficient value of the coherence function is close to zero, this means that the pressure signal is damped markedly by some other sources from one measurement position to another one. On the other hand, the situation is opposite if the value is close to 1, which proves that the propagation is well-transported and attenuation of the pressure wave is not evident. Generally, when the value of the coherence coefficient is above 0.6, it is considered to be similar for two objective signals.

    In this section, the axial distribution of the coherence coefficient is used to explain the propagation of pressure wave under the different flow regimes and discharge modes in the dipleg-trickle valve system. The maximum values of coherence coefficients for each of two group adjacent pressure fluctuation signals along axial positions are displayed in Fig. 12. It is shown that, at a given solids flux [Fig. 12a,Gs= 17.4 kg/(m2s)], the increases in negative pressure will change the axial distribution of coherence coefficient from a uniformity distribution to a sudden change at an axial position. Similarly, as shown in Fig. 12b, at a given negative pressure (ΔP= 4.1 kPa), the decreases in solids flux will induce the same change tendency as shown in Fig. 12a. As mentioned above, the change of the operating conditions will induce the transition of flow regimes and discharge modes in the dipleg-trickle valve system, which contributes to the different axial distribution of the coherence coefficient.

    In Fig. 12a, b, the case for axial uniform distribution of coherence coefficient is corresponding to the dilute falling flow regime in the dipleg and the continuous trickling discharge mode for the trickle valve. It demonstrates that the pressure fluctuation signal has well consistency from the top to the bottom in the whole dipleg-trickle valve system, which has the unattenuated feature of downward propagation of pressure fluctuation similarly originated from the top in fluidized beds as reported by Bi(2007), Pun?ochá? and Draho?(2005), van der Schaaf et al. (1998) and Fan et al. (1981).That is to say, the source of pressure fluctuation originates from the feed at the top inlet of the dipleg and propagates downward to the bottom exit for this case. The down- flowing particles look like a rain falling inside the dipleg. The gas–solids flow in the dipleg is flowing against negative pressure, and the slip velocity is upward for balancing the negative pressure gradient, which results in compressing the gas flowing up the dipleg. The gas is compressed and expanded alternately to induce the non-uniform distribution of the particle concentration and the pressure. Therefore, the pressure fluctuation is attributed to the unstable feed and the compressed gas for this case.

    However, the sudden change of coherence coefficient at an axial position as shown in Fig. 12a, b demonstrates that the pressure fluctuation signal has well consistency only in the upper section or in the bottom section, respectively. This type of axial distribution of coherence coefficient is closely related to the dilute–dense phase coexisting falling flow regime in the dipleg and the intermittent periodic dumping discharge mode for the trickle valve. The axial uniform distribution of coherence coefficient in the upper section is corresponding to the dilute phase falling flow. The propagation of pressure fluctuation originates from the feed at the top inlet and propagates downward along the dipleg. The falling particle clusters interact with the up flowing gas, causing the fluctuation of pressure in this section. Meanwhile, the consistency between two adjacent pressure signals in the lower dense bed reveals that another pressure fluctuation propagates in this section. The pressure fluctuation mainly originates from the discharge oscillation at the bottom exit, and the bubble coalescence and splitting in the bubble passage,which propagates upward in the dense bed. The intermittent periodic dumping discharge mode causes the dynamic variation in pressure in this region. Meanwhile, the bubble eruption at the surface of the dense phase bed can cause a pressure fluctuation to propagate downward along the dense bed with no attenuation similar to that in fluidized beds reported by Bi (2007), Pun?ochá? and Draho? (2005), van der Schaaf et al. (1998) and Fan et al. (1981). So, four types of pressure fluctuation originating, respectively, from the top inlet feed,the bottom discharge oscillation, the rising bubble passage,and the bubble eruption at the bed surface, encounter at an axial position and interact from each other to attenuate individually, leading to the appearance of the transition point of coherence coefficient along the axial distribution.

    Therefore, these distribution features of the coherence coefficient can also be used to demarcate the flow regimes and discharge modes and further to predict the origin and propagation of pressure fluctuation in the dipleg-trickle valve system as shown in Fig. 13. For the case of the intermittent periodic dumping discharge mode and the dilute–dense phase coexisting falling flow regime (Fig. 13a), the energy of global pressure fluctuation mainly originates from the bottom oscillation of particles dumping discharge through the trickle valve, the coalescence and splitting of rising bubble (or void) in passage in the lower dense bed section, the bubble eruption at the surface of the dense bed, the instable feed at the inlet of dipleg, and the compressed gas impacted by clusters in the upper dilute section. The pressure fluctuation from the discharge oscillation and bubble passage propagates upward in the lower dense bed section. And the pressure fluctuation from the instable feed at the dipleg inlet and the compressed gas in the upper dilute section propagates downward.

    Fig. 13 Diagram showing the origin and propagation direction of pressure fluctuation (or pressure wave) in the dipleg-trickle valve system

    However, for the continuous trickling discharge mode and the dilute falling flow regime (Fig. 13b), due to the continuous and steady discharge of trickle valve, the oscillation produced by this type of discharge is relatively weak. The energy of global pressure fluctuation for this case mainly originates from the upper instable feed at the inlet of dipleg and from the compressed gas impacted by the dilute falling particle swarms in the dipleg. The pressure fluctuation propagates downward along the dipleg.

    3.5 Flow regime and discharge mode map

    Based on the experimental phenomena and the data analysis as mentioned above, the operation of the experimental dipleg-trickle valve system can be divided into two types of discharge modes and flow regimes, which are closely related to the operating parameters, the solids flux (Gs), and the negative pressure (ΔP). So, in this work, a map describing the relationship between the operation modes of diplegtrickle valve system and the operating parameters is plotted in Fig. 14. The region above the demarcation line expresses the dilute–dense phase coexisting falling flow regime in the dipleg and the intermittent periodic dumping discharge mode in the trickle valve, respectively. Correspondingly,the region below the demarcation line represents the dilute falling flow regime and the continuous trickling discharge mode, respectively.

    This map is suitable for predicting the operation status of a dipleg-trickle valve system at low solids flux and variable negative pressure, especially for the secondary suspended cyclone dipleg-trickle valve system in the FCC disengager which usually has a very low particle concentration at the secondary cyclone inlet. Therefore, as long as the operating parameters (Gsand ΔP) are known, the flow regimes and the discharge modes in the dipleg-trickle valve system could be determined from the map, which can be applied to predict and monitor its operation status in the FCC disengager.

    Fig. 14 Flow regime and discharge mode map in the dipleg-trickle valve system

    4 Conclusions

    In this work, the flow regimes and discharge modes in the dipleg-trickle valve system were studied experimentally based on the measurement of transient pressure fluctuation.The experimental phenomena show that the flow regimes of particles down the dipleg could be divided into two types,the dilute–dense phase coexisting falling flow, and the dilute falling flow, which closely correlate with the operating parameters (solids flux and negative pressure). The change of operating parameters can cause the transformation of these two flow regimes from each other. Correspondingly,the trickle valve also has two discharge modes, the intermittent periodic dumping discharge, and the continuous trickling discharge. The former corresponds to the dilute–dense phase coexisting falling flow in the dipleg, and the latter corresponds to the dilute falling flow in the dipleg.

    The pro files of transient pressure fluctuation in close to the exit of trickle valve can be used to demarcate the two flow regimes and discharge modes. The zigzag shape of transient pressure fluctuation represents the operation status with the dilute–dense phase coexisting falling flow regime in the dipleg and the intermittent periodic dumping discharge mode for the trickle valve. Meanwhile, the smooth transient pressure fluctuation describes the characteristics of the dilute falling flow in the dipleg and the continuous trickling discharge mode for the trickle valve.

    The coherence coefficient based on frequency domain well explains the origin and propagation of pressure fluctuation in the dipleg-trickle valve system. The transient pressure fluctuation shows that there exists dynamic property of gas–solids flow in the dipleg-trickle valve system. The operating conditions have obvious effects on the pressure fluctuation. The power spectrum density (PSD) based on frequency domain shows that the dominant frequency of pressure fluctuation time series concentrates on the range of less than 0.03 Hz. It means that the gas–solids flow in the dipleg-trickle valve system is characterized by the low frequency pulsation.

    The bottom discharge oscillation, the coalescence and splitting of rising bubbles (or voids) in passage, the bubble eruption at the surface of dense phase bed, the instable feed at the inlet of dipleg, and the compressed gas impacted by clusters jointly contribute to the fluctuation of pressure in the dipleg-trickle valve system. The pressure fluctuation from the discharge oscillation and bubble passage propagates upward, and the one from the instable feed, the bubble eruption, and the compressed gas propagates downward in the dipleg-trickle valve system.

    AcknowledgementsThe authors acknowledge the supports by the National Natural Science Foundation of China (Grant Nos. 21776305,21978322), and the China Petrochemical Joint Foundation (Grant No.U1862202).

    Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons.org/licen ses/by/4.0/.

    在线播放无遮挡| 精品一区二区三区视频在线| 噜噜噜噜噜久久久久久91| 777米奇影视久久| 国产精品国产三级国产av玫瑰| 男人舔奶头视频| .国产精品久久| 中文天堂在线官网| 精品人妻偷拍中文字幕| 日本免费在线观看一区| 日本熟妇午夜| 国产综合懂色| 亚洲在线观看片| av在线老鸭窝| videos熟女内射| 夜夜爽夜夜爽视频| 精品国产露脸久久av麻豆 | 日韩强制内射视频| 精品人妻熟女av久视频| 中文字幕久久专区| 在线天堂最新版资源| 3wmmmm亚洲av在线观看| 丝袜美腿在线中文| 午夜福利成人在线免费观看| 伊人久久精品亚洲午夜| 黄片无遮挡物在线观看| 日韩成人av中文字幕在线观看| 国产黄色视频一区二区在线观看| 精品久久久久久久人妻蜜臀av| 国产成人精品久久久久久| 亚洲欧美成人综合另类久久久| 欧美日韩视频高清一区二区三区二| 最近的中文字幕免费完整| 欧美一级a爱片免费观看看| 在线a可以看的网站| 日韩中字成人| 91在线精品国自产拍蜜月| 狂野欧美白嫩少妇大欣赏| 在线观看人妻少妇| 亚洲国产高清在线一区二区三| 亚洲国产av新网站| 亚洲欧美清纯卡通| 精品久久久噜噜| 精品国产一区二区三区久久久樱花 | 国产亚洲av嫩草精品影院| 日本三级黄在线观看| www.av在线官网国产| 国产精品福利在线免费观看| 一级毛片久久久久久久久女| 又爽又黄无遮挡网站| 中文资源天堂在线| 国产乱人视频| 成人漫画全彩无遮挡| 免费黄色在线免费观看| a级一级毛片免费在线观看| 永久网站在线| av免费在线看不卡| 99久国产av精品国产电影| 最近最新中文字幕大全电影3| av天堂中文字幕网| 亚洲精品色激情综合| 亚洲高清免费不卡视频| 观看免费一级毛片| 成人毛片a级毛片在线播放| 中文精品一卡2卡3卡4更新| 亚洲综合精品二区| 99久久精品一区二区三区| 国产免费一级a男人的天堂| 国产成人免费观看mmmm| 99久久精品国产国产毛片| 直男gayav资源| 九草在线视频观看| 精品久久久久久电影网| 成人无遮挡网站| 99热这里只有精品一区| 黄色日韩在线| 亚洲精品色激情综合| 亚洲va在线va天堂va国产| av线在线观看网站| 插逼视频在线观看| 在线天堂最新版资源| 久久国产乱子免费精品| 少妇裸体淫交视频免费看高清| 人妻系列 视频| 人人妻人人澡人人爽人人夜夜 | 一本久久精品| 国产在视频线在精品| 黄色一级大片看看| 免费av毛片视频| 亚洲高清免费不卡视频| 亚洲在线自拍视频| 国产精品福利在线免费观看| 卡戴珊不雅视频在线播放| 亚洲av电影不卡..在线观看| 人体艺术视频欧美日本| www.色视频.com| 精品久久久久久久末码| 尤物成人国产欧美一区二区三区| 精品不卡国产一区二区三区| 国产精品人妻久久久影院| 五月天丁香电影| 成人一区二区视频在线观看| 国产精品福利在线免费观看| 大片免费播放器 马上看| 丝瓜视频免费看黄片| 人人妻人人澡人人爽人人夜夜 | 一级av片app| 天天一区二区日本电影三级| 精品国产三级普通话版| 免费观看的影片在线观看| 国产精品麻豆人妻色哟哟久久 | 一区二区三区四区激情视频| 精品99又大又爽又粗少妇毛片| 一个人看的www免费观看视频| av国产免费在线观看| 一夜夜www| 男人狂女人下面高潮的视频| 99视频精品全部免费 在线| 久久久久久伊人网av| 久久午夜福利片| 特大巨黑吊av在线直播| 欧美最新免费一区二区三区| 少妇被粗大猛烈的视频| 搡老妇女老女人老熟妇| 51国产日韩欧美| 亚洲综合色惰| 身体一侧抽搐| 黄片wwwwww| www.色视频.com| 七月丁香在线播放| 激情五月婷婷亚洲| 日韩av在线大香蕉| 一区二区三区乱码不卡18| 国产精品国产三级国产专区5o| 高清在线视频一区二区三区| 亚洲性久久影院| 老司机影院毛片| 街头女战士在线观看网站| 欧美精品国产亚洲| 91久久精品国产一区二区成人| 最后的刺客免费高清国语| 韩国av在线不卡| 亚洲精品一区蜜桃| 久久国内精品自在自线图片| 插阴视频在线观看视频| 99热6这里只有精品| 亚洲久久久久久中文字幕| 免费观看精品视频网站| 日日摸夜夜添夜夜添av毛片| 欧美精品一区二区大全| 亚洲av一区综合| 亚洲精品乱码久久久v下载方式| 三级男女做爰猛烈吃奶摸视频| 国产精品国产三级国产av玫瑰| 成人性生交大片免费视频hd| 91精品一卡2卡3卡4卡| 国产精品无大码| 一个人观看的视频www高清免费观看| videos熟女内射| 一区二区三区高清视频在线| 一级毛片 在线播放| 日本爱情动作片www.在线观看| 成人二区视频| 天堂中文最新版在线下载 | 成人欧美大片| 可以在线观看毛片的网站| 成年版毛片免费区| 美女黄网站色视频| 国产精品一二三区在线看| 国产成人精品久久久久久| 最新中文字幕久久久久| 亚洲国产精品成人久久小说| 久久久欧美国产精品| 精华霜和精华液先用哪个| 欧美成人一区二区免费高清观看| 欧美高清成人免费视频www| 国产一区二区在线观看日韩| 男人和女人高潮做爰伦理| 中文字幕av在线有码专区| 少妇人妻一区二区三区视频| 亚洲精品影视一区二区三区av| 国产免费又黄又爽又色| 亚洲国产精品sss在线观看| 天天躁日日操中文字幕| 亚洲在线自拍视频| 精品久久久久久久久亚洲| 91久久精品国产一区二区三区| 日韩三级伦理在线观看| 国产不卡一卡二| 2018国产大陆天天弄谢| 久久久久九九精品影院| 观看免费一级毛片| 亚洲欧美精品自产自拍| 久久久国产一区二区| 99久国产av精品国产电影| 国产成人精品一,二区| 日本-黄色视频高清免费观看| 色吧在线观看| 男人和女人高潮做爰伦理| 非洲黑人性xxxx精品又粗又长| 国产在线男女| 久久久午夜欧美精品| 国产精品爽爽va在线观看网站| 亚洲欧美一区二区三区国产| 国产精品人妻久久久久久| 日韩一本色道免费dvd| 淫秽高清视频在线观看| 免费在线观看成人毛片| 免费观看精品视频网站| 三级国产精品欧美在线观看| 日韩制服骚丝袜av| 91精品伊人久久大香线蕉| 特大巨黑吊av在线直播| 亚洲国产日韩欧美精品在线观看| 干丝袜人妻中文字幕| videos熟女内射| 大陆偷拍与自拍| 欧美日韩在线观看h| 国产精品人妻久久久影院| av国产久精品久网站免费入址| 大又大粗又爽又黄少妇毛片口| 汤姆久久久久久久影院中文字幕 | 欧美3d第一页| 日本-黄色视频高清免费观看| 激情 狠狠 欧美| 亚洲精品乱码久久久v下载方式| 麻豆久久精品国产亚洲av| 一级毛片电影观看| 尤物成人国产欧美一区二区三区| 亚洲欧美精品专区久久| 97热精品久久久久久| 大香蕉97超碰在线| 国产成人精品一,二区| 国产一区二区三区av在线| 毛片女人毛片| h日本视频在线播放| 一级a做视频免费观看| 亚洲av电影不卡..在线观看| 只有这里有精品99| 亚洲国产精品成人久久小说| 在线a可以看的网站| 久久久精品欧美日韩精品| av线在线观看网站| 亚洲熟女精品中文字幕| 狂野欧美白嫩少妇大欣赏| 黄色日韩在线| 日韩成人av中文字幕在线观看| 日韩视频在线欧美| 26uuu在线亚洲综合色| 国产人妻一区二区三区在| 精品少妇黑人巨大在线播放| 91久久精品国产一区二区三区| 精品人妻熟女av久视频| 久久鲁丝午夜福利片| 国产亚洲91精品色在线| 最近手机中文字幕大全| 国产v大片淫在线免费观看| 深爱激情五月婷婷| 国产精品久久久久久精品电影小说 | 国产成人91sexporn| 91在线精品国自产拍蜜月| 一级黄片播放器| 免费少妇av软件| 在线天堂最新版资源| 色综合亚洲欧美另类图片| 午夜激情欧美在线| 1000部很黄的大片| 亚洲精品影视一区二区三区av| 免费少妇av软件| 国产男女超爽视频在线观看| 搡老妇女老女人老熟妇| 久久久久久久亚洲中文字幕| 日韩av在线免费看完整版不卡| 成人午夜高清在线视频| 尾随美女入室| 亚洲av男天堂| 午夜激情欧美在线| 亚洲第一区二区三区不卡| 淫秽高清视频在线观看| 青春草视频在线免费观看| 热99在线观看视频| 成人鲁丝片一二三区免费| 寂寞人妻少妇视频99o| 亚洲18禁久久av| 久99久视频精品免费| 欧美区成人在线视频| 国产高清有码在线观看视频| 久久久欧美国产精品| 免费不卡的大黄色大毛片视频在线观看 | 国产精品1区2区在线观看.| 成人漫画全彩无遮挡| 国产色婷婷99| 久久久久久伊人网av| 欧美3d第一页| 国产亚洲5aaaaa淫片| 国产精品嫩草影院av在线观看| 国产高清不卡午夜福利| 性插视频无遮挡在线免费观看| 人妻一区二区av| 色哟哟·www| 国产高清三级在线| 国产精品一及| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品成人久久久久久| 一级黄片播放器| 亚洲精品aⅴ在线观看| 一区二区三区乱码不卡18| 白带黄色成豆腐渣| 亚洲aⅴ乱码一区二区在线播放| 人妻一区二区av| 国产精品爽爽va在线观看网站| 国产探花在线观看一区二区| 高清av免费在线| 国产一区有黄有色的免费视频 | 91精品国产九色| 久久99精品国语久久久| 久久久久性生活片| 亚洲精品,欧美精品| 中文乱码字字幕精品一区二区三区 | 听说在线观看完整版免费高清| 日韩欧美一区视频在线观看 | 日韩成人伦理影院| 亚洲怡红院男人天堂| 免费看av在线观看网站| 国产黄频视频在线观看| 久久99热6这里只有精品| 在线观看av片永久免费下载| 精品人妻熟女av久视频| 久久精品人妻少妇| 一级黄片播放器| 久久久久久九九精品二区国产| 一本一本综合久久| 美女脱内裤让男人舔精品视频| 欧美成人a在线观看| 精品酒店卫生间| 欧美日韩视频高清一区二区三区二| 亚洲精品一二三| 亚洲国产欧美人成| 18+在线观看网站| 夫妻午夜视频| 少妇的逼水好多| 床上黄色一级片| 欧美一级a爱片免费观看看| 99热网站在线观看| 免费黄频网站在线观看国产| 九九久久精品国产亚洲av麻豆| 狂野欧美白嫩少妇大欣赏| 欧美最新免费一区二区三区| 99热网站在线观看| 晚上一个人看的免费电影| 亚洲精品日韩av片在线观看| 久久韩国三级中文字幕| 久久久久国产网址| av卡一久久| 99热这里只有是精品在线观看| 精品国产露脸久久av麻豆 | 一区二区三区乱码不卡18| 三级经典国产精品| 最后的刺客免费高清国语| 大陆偷拍与自拍| 久久久久久国产a免费观看| 国产探花在线观看一区二区| 色综合站精品国产| 伊人久久国产一区二区| 国产 亚洲一区二区三区 | videossex国产| 国产精品.久久久| 一级片'在线观看视频| 国产日韩欧美在线精品| 亚洲成人精品中文字幕电影| 午夜福利成人在线免费观看| 在线观看av片永久免费下载| 99热这里只有是精品在线观看| 久久久久久久久久人人人人人人| 精品久久久精品久久久| 欧美激情国产日韩精品一区| 青春草亚洲视频在线观看| 国产在线一区二区三区精| 搡女人真爽免费视频火全软件| 国产淫片久久久久久久久| 成人亚洲精品一区在线观看 | 亚洲国产精品国产精品| 国产伦精品一区二区三区视频9| 久久久久精品性色| 国产乱人视频| 大香蕉97超碰在线| 亚洲精品aⅴ在线观看| 欧美成人a在线观看| 三级男女做爰猛烈吃奶摸视频| 美女xxoo啪啪120秒动态图| 日韩成人伦理影院| 一级a做视频免费观看| 国产午夜精品一二区理论片| 久久久久久久久久黄片| 最近最新中文字幕大全电影3| 亚洲国产精品sss在线观看| 亚洲国产成人一精品久久久| 汤姆久久久久久久影院中文字幕 | 亚洲综合精品二区| 欧美精品一区二区大全| 一区二区三区乱码不卡18| 国产免费福利视频在线观看| 亚洲精品日本国产第一区| 国产白丝娇喘喷水9色精品| 国产精品一区二区三区四区久久| 亚洲丝袜综合中文字幕| 免费不卡的大黄色大毛片视频在线观看 | 老女人水多毛片| 91aial.com中文字幕在线观看| 乱人视频在线观看| 国产成人免费观看mmmm| 最近手机中文字幕大全| 精品久久久噜噜| 日本黄色片子视频| 午夜免费观看性视频| 国产成人精品久久久久久| 成人漫画全彩无遮挡| 久久久久久久亚洲中文字幕| 国产高潮美女av| 色综合站精品国产| 精品人妻一区二区三区麻豆| 欧美精品一区二区大全| videossex国产| 日韩视频在线欧美| 久久精品人妻少妇| 午夜免费激情av| 一级毛片aaaaaa免费看小| 日韩不卡一区二区三区视频在线| 久久国内精品自在自线图片| 久久这里有精品视频免费| 久久久久久久久大av| 国产美女午夜福利| 日本-黄色视频高清免费观看| 欧美成人a在线观看| av.在线天堂| 天堂√8在线中文| 80岁老熟妇乱子伦牲交| 日韩制服骚丝袜av| 色尼玛亚洲综合影院| 最近最新中文字幕免费大全7| 精品久久久久久成人av| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日日摸夜夜添夜夜爱| 国产白丝娇喘喷水9色精品| 乱人视频在线观看| 免费无遮挡裸体视频| av在线天堂中文字幕| 中文资源天堂在线| 日韩av在线大香蕉| 人妻制服诱惑在线中文字幕| 亚洲国产最新在线播放| 女人久久www免费人成看片| 春色校园在线视频观看| 精品酒店卫生间| 男插女下体视频免费在线播放| 1000部很黄的大片| 久久久色成人| 午夜福利在线在线| 一级黄片播放器| 少妇被粗大猛烈的视频| 99久国产av精品| 菩萨蛮人人尽说江南好唐韦庄| 一级av片app| 美女国产视频在线观看| 亚洲久久久久久中文字幕| 国产极品天堂在线| 在现免费观看毛片| 欧美最新免费一区二区三区| 男人爽女人下面视频在线观看| 国产成人精品福利久久| 国产男人的电影天堂91| 中文字幕免费在线视频6| av国产免费在线观看| 国产精品日韩av在线免费观看| av.在线天堂| 久久久久久久久久人人人人人人| 免费看a级黄色片| 看非洲黑人一级黄片| 大话2 男鬼变身卡| 国产伦理片在线播放av一区| 国产成人一区二区在线| 中文资源天堂在线| 夫妻性生交免费视频一级片| 成人高潮视频无遮挡免费网站| 亚洲乱码一区二区免费版| 如何舔出高潮| 国产一区二区三区av在线| 女人被狂操c到高潮| 日韩精品青青久久久久久| 欧美xxxx性猛交bbbb| 啦啦啦韩国在线观看视频| 一级毛片黄色毛片免费观看视频| 内射极品少妇av片p| 精品人妻一区二区三区麻豆| 麻豆乱淫一区二区| 美女高潮的动态| 男女那种视频在线观看| 精品一区二区三区人妻视频| 人妻系列 视频| 中文资源天堂在线| 日韩不卡一区二区三区视频在线| 偷拍熟女少妇极品色| 26uuu在线亚洲综合色| 美女主播在线视频| 一级a做视频免费观看| 国语对白做爰xxxⅹ性视频网站| 日韩欧美精品免费久久| 国产成人aa在线观看| 久久久a久久爽久久v久久| 内射极品少妇av片p| av在线观看视频网站免费| 中文乱码字字幕精品一区二区三区 | 免费黄频网站在线观看国产| 亚洲av二区三区四区| 真实男女啪啪啪动态图| 国产久久久一区二区三区| 丝瓜视频免费看黄片| 国产淫语在线视频| 看黄色毛片网站| 日日啪夜夜撸| 秋霞在线观看毛片| 久久久久久九九精品二区国产| 欧美3d第一页| 国内精品美女久久久久久| 青春草国产在线视频| 汤姆久久久久久久影院中文字幕 | 免费不卡的大黄色大毛片视频在线观看 | 18禁动态无遮挡网站| 最近2019中文字幕mv第一页| 免费少妇av软件| 伦理电影大哥的女人| 美女大奶头视频| 18禁裸乳无遮挡免费网站照片| 婷婷色麻豆天堂久久| 国产精品国产三级国产专区5o| 热99在线观看视频| 日本三级黄在线观看| 精品久久久久久久久av| 色尼玛亚洲综合影院| 综合色av麻豆| 免费av不卡在线播放| 一个人观看的视频www高清免费观看| 男的添女的下面高潮视频| 亚洲国产色片| 精品久久久久久久久亚洲| 日日干狠狠操夜夜爽| 成年av动漫网址| 亚洲精品aⅴ在线观看| 免费看a级黄色片| 51国产日韩欧美| 能在线免费观看的黄片| 国产欧美日韩精品一区二区| 日韩不卡一区二区三区视频在线| 中文字幕亚洲精品专区| 天堂√8在线中文| 97精品久久久久久久久久精品| 人人妻人人澡人人爽人人夜夜 | videos熟女内射| 国产午夜精品论理片| 免费无遮挡裸体视频| 国产精品久久久久久精品电影| 丝袜美腿在线中文| 麻豆精品久久久久久蜜桃| 国产精品不卡视频一区二区| 国国产精品蜜臀av免费| 一个人看的www免费观看视频| 一二三四中文在线观看免费高清| 免费大片黄手机在线观看| 成年人午夜在线观看视频 | 日本黄大片高清| 最新中文字幕久久久久| 99久久精品热视频| 亚洲国产欧美在线一区| 日日摸夜夜添夜夜爱| 一级av片app| 中文字幕av在线有码专区| 2022亚洲国产成人精品| 大片免费播放器 马上看| 秋霞伦理黄片| 久久久午夜欧美精品| 亚洲内射少妇av| 毛片一级片免费看久久久久| 免费观看的影片在线观看| 日本一二三区视频观看| 秋霞伦理黄片| 国产极品天堂在线| 国产黄a三级三级三级人| 成人无遮挡网站| 亚洲性久久影院| 亚洲综合精品二区| 日本免费在线观看一区| av国产久精品久网站免费入址| 寂寞人妻少妇视频99o| 一区二区三区四区激情视频| 日韩欧美国产在线观看| 夫妻性生交免费视频一级片| 国产乱人视频| 亚洲成人精品中文字幕电影| 亚洲最大成人av| 水蜜桃什么品种好| 91久久精品国产一区二区三区| 干丝袜人妻中文字幕| 麻豆成人av视频| 日韩中字成人| 亚洲av免费高清在线观看| 国内揄拍国产精品人妻在线| 中文字幕免费在线视频6| 最近视频中文字幕2019在线8| 欧美丝袜亚洲另类| 一个人看视频在线观看www免费| 草草在线视频免费看| 久久99热6这里只有精品| 久久久久免费精品人妻一区二区| 婷婷六月久久综合丁香| 久99久视频精品免费|