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

    Experimental study on propagation characteristics of rotating detonation wave with kerosene fuel-rich gas

    2022-08-30 14:33:32JiaxiangHanQiaodongBaiShijianZhangFangWangChunshengWeng
    Defence Technology 2022年8期

    Jia-xiang Han,Qiao-dong Bai,Shi-jian Zhang,Fang Wang,Chun-sheng Weng

    National Key Laboratory of Transient Physics,Nanjing University of Science and Technology,Nanjing,210094,China

    Keywords:Rotating detonation wave Kerosene fuel-rich gas Initiation process Propagation mode

    ABSTRACT In this study,kerosene fuel-rich gas produced by the combustion in the gas generator was used as the fuel and oxygen-rich air was used as the oxidant to investigate the propagation characteristics of the rotating detonation wave(RDW).The initiation of the kerosene fuel-rich gas and propagation process of the RDW were analyzed.The influences of the oxygen content in the oxidizer,kerosene mass flow rate of the gas generator,and temperature of the kerosene fuel-rich gas on the propagation process of the RDW were studied.The experimental results revealed that the propagation velocity of the RDW could be improved by increasing the three parameters mentioned above with the kerosene mass flow rate as the strongest factor.The minimum oxygen content that could successfully initiate and maintain the stable propagation of the RDW was 32%,achieving the RDW velocity of 1141.9 m/s.The RDW mainly propagated as two-counter rotating waves and a single wave when the equivalent ratios were 0.62-0.79 and 0.85-0.87,respectively.The highest RDW velocity of 1637.2 m/s was obtained when the kerosene mass flow rate,oxygen content,and equivalent ratio were 74.6 g/s,44%,and 0.87,respectively.

    1.Introduction

    Continuous rotating detonation engine(CRDE) is a new type of engine that utilizes a rotating detonation wave(RDW) to generate thrust in a rotating detonation combustor (RDC).After an RDW passes through combustible mixtures,the high-temperature and high-pressure products generated expand,and are discharged from the end of the RDC at a high speed,thereby generating thrust.CRDEs have broad application prospects and offer advantages,such as high thermal cycle efficiency,fast heat release rate,compact structure,and wide working range.Therefore,CRDEs have attracted extensive attention in the field of hypersonic propulsion technology.

    The formation and propagation characteristics of the RDW with hydrogen as the fuel have been widely studied[1-3].The effects of the oxygen volume fraction,mass flow rate,and equivalent ratio on the propagation of continuous RDW have been studied in a circular RDC [4].Moreover,experimental studies on the combustion characteristics and working principle under poor fuel conditions have been conducted [5].The effects of the combustion chamber configurations on the RDW propagation mode [6],relationship between high-frequency tangential instability and continuous rotating detonation[7],and flow field structure in a stable working mode [8] have been investigated.In addition,there have been numerical studies on the H/air RDW flow field and propagation characteristics of the shock wave at the head of the combustor[9],and mechanism of chemical dynamics of a CRDE [10].

    Kerosene,which is the optimum fuel used in hypersonic propulsion system,has the advantages of high calorific value,satisfactory combustion performance,convenience of use and maintenance,non-toxicity,and strong controllability,among others.However,the ignition of liquid kerosene remains a challenge because of the difficulty of atomization and uneven mixing process.Rotating detonation using kerosene as the fuel has been experimentally investigated by adding hydrogen,oxygen,or heated air to the reactants to improve the activity of combustible mixtures,thereby minimizing the challenges associated with liquid hydrocarbon fuel initiation.Stable rotating detonation was successfully achieved using liquid kerosene as the fuel with hydrogen [11],syngas additives[12],and excess oxygen[13,14].Further,preheated liquid hydrocarbon fuel with heated air[15,16]was also applied to facilitate rotating detonation.The effects of the total pressure,equivalence ratio,and total inlet temperature on the working performance of kerosene/air rotating detonation have also been studied [17-19].Furthermore,the reactivity of liquid hydrocarbon fuel has been studied,and improved through vaporization,cracking,or precombustion to mitigate the challenges associated with its initiation.

    A mixture of ethylene-acetylene-hydrogen gas has been used as a substitute of precombustion-cracked kerosene to investigate the rotating detonation characteristics [20-22].The detonation delay time and propagation characteristics of an RDW have also been studied using methane-hydrogen-carbon dioxide as the fuel and air as the oxidant[23].Furthermore,the effects of the fuel ratio[24],refueling method,and refueling distance [25] on the distribution of pyrolysis products have been studied to realize kerosene/air rotating detonation without adding hydrogen and oxygen.The influence of the channel and oxidizer injection slot width on the rotating detonation characteristics of the precombustion-cracked kerosene has also been investigated [26-28].An experimental study has been carried out on the process of rotating detonation combustion using the scheme of the secondary detonation combustion of kerosene fuel-rich gas (KFRG) and oxygen-rich air after the primary combustion of liquid kerosene [29].Additionally,the influence of the nozzle configuration on the rotating detonation characteristics of pyrolysis gas in kerosene precombustion was studied [30],and a two-dimensional numerical simulation was conducted to study the initiation process of kerosene,which was partially cracked by preheating and subsequent ignition [31].

    Owing to the complex composition of KFRG,its rotating detonation propagation characteristics are different from those of existing conventional fuels;thus,a fundamental understanding of the rotating detonation propagation characteristics of KFRG is required.To date,only a few studies have investigated the influence of the width of the inlet ring,width of the combustion chamber,equivalent ratio,and oxygen content of the oxidizer on the rotating detonation using KFRG as the fuel.Further,research on the characteristics of RDW of KFRG is currently in its initial stage.

    In this study,a precombustion method of kerosene was adopted to form KFRG,thereby reducing the difficulty of initiating liquid kerosene.The initiation of KFRG and propagation characteristics of the RDW were experimentally studied.The influences of the oxygen mass fraction,kerosene mass flow rate,and gas temperature on the propagation of the RDW of KFRG were further explored.The results were compared with those of the rotating detonation of liquid kerosene.This paper reports the results of the investigation of multi-component mixtures in the rotating detonation phenomenon.The experimental results can be used in the research of combined air turbine rocket detonation engine for future studies.

    2.Experimental apparatus and measurement system

    The experimental system and experimental equipment photo for the rotating detonation of KFRG are shown in Fig.1.It includes a fuel/oxidant supply system,control system,data acquisition system,hot jet ignition system,KFRG generator,and RDC.The data acquisition system comprised a static pressure sensor,highfrequency dynamic pressure sensor (PCB),turbine flowmeter,and temperature sensor.

    The KFRG generator comprised an injection section,combustion chamber,and nozzle.Kerosene and precombustion oxygen were injected into the combustion chamber,burned to form the KFRG,and subsequently passed through the nozzle into the KFRG collection chamber.The KFRG pressure in the combustion chamber was controlled by changing the nozzle throat diameter.Jet fuel no.3 was used in the experiment.

    The inner and outer diameters of the RDC were 94 mm and 110 mm,respectively.The width and length of the combustion chamber were 8 and 180 mm,respectively.The KFRG was injected into the rotating combustion chamber through 80 circular holes uniformly distributed on the inner wall with a total injection area of 40.19 mm.The combustible mixture was ignited by a hot jet.

    High-frequency dynamic pressure sensors were installed on the wall surface of the hot jet tube (PCB0) and RDC (PCB1,PCB2,and PCB3).The pressure sensors installed on the circumferential plane were aligned with the ignition position and evenly distributed at an angle of 120,as shown in Fig.2.The sampling frequency was set to 500 kHz.

    Static pressure sensors were installed on the oxygen gas collecting chamber,combustion chamber,KFRG collecting chamber,air collecting chamber,and outer wall surface of the RDC,labelled as P1,P2,P3,P4,and P5,respectively.Temperature sensors were installed on the gas generator combustion chamber,and KFRG collecting chamber,labelled as T1 and T2,respectively,as shown in Fig.1.

    Fig.1.(a) Schematic of the rotating detonation experimental system for KFRG;(b).Photo of experimental equipment.

    Fig.2.Installation position of the high-frequency dynamic pressure sensors.

    The time sequence of the experiments is shown in Fig.3.The data acquisition system was operational throughout the experimental process.After the stabilization of the process in the gas generator,the hydrogen/oxygen mixture in the hot jet tube was ignited to produce high-temperature and high-pressure jet flow into the RDC.The KFRG/oxygen-rich air mixture was subsequently ignited to form an RDW.

    Fig.3.Time sequence of the experiments.

    In this study,the mass fraction of oxygen in the oxidizer,kerosene mass flow rate,and gas temperature were used as the variables.The control variable method was used to design the experimental conditions.A total of 14 groups of experimental conditions were set to analyze the propagation characteristics of the RDW of the KFRG under various conditions.The ambient temperature and environmental pressure were 288 K and 1 atm,respectively.

    3.Results and discussion

    Table 1,Table 2,and Table 3 present the influence of the oxygen mass fraction in the oxidizer on the RDW velocity,kerosene mass flow rate on the RDW wave velocity and mode,and KFRG temperature on the RDW formation and propagation velocity,respectively.The detailed discussion of the results is presented in the subsequent sections.

    Table 1.Experimental parameters and results with varying oxygen mass fraction.

    Table 2.Experimental parameters and results with varying kerosene mass flow rate at 44% oxygen mass fraction.

    Table 3.Experimental parameters and results with varying gas temperature at 32% oxygen mass fraction.

    3.1.Initiation of the RDW of KFRG

    The RDW initiation of the KFRG was analyzed for case 2.Fig.4 shows the static pressure history of the gas generator and RDC.Fig.4(a) shows the pressure history of the gas generator oxygen chamber(P1)and combustor(P2).When oxygen was injected into the gas generator at 1.1 s,P1 gradually increased and reached an equilibrium point of 0.74 MPa at 1.8 s.At 2.6 s,the high-energy spark plug ignited,resulting in an ignition pressure peak owing to the fuel residue.At 3.2 s,P2 increased owing to the large amount of KFRG generated in the gas generator.At 4.0 s,P2 and P1 exhibited steady state pressures of 1.09 and 1.25 MPa,respectively.At 4.388 s,the hot jet tube was ignited,and the RDC operation began.The pressure of the rotating detonation in the RDC increased the P2 to 1.12 MPa.As the kerosene supply was terminated at 5.1 s,P2 slightly increased owing to the combustion of the residual kerosene,and subsequently sharply decreased.

    Fig.4(b) shows the pressure history of the KFRG collecting chamber(P3),air collecting chamber(P4),and RDC(P5).The trend of P3 was similar to that of P2.Air was injected into the RDC at 2.6 s,and oxygen was added to the air chamber at 3.1 s.Subsequently,oxygen-rich air was formed,thereby P4 reached 0.39 MPa.At 4.388 s,hot jet was injected into the RDC.Owing to the back pressure of the rotating detonation,P4 and P5 reached 0.49 and 0.22 MPa,respectively.At 5.1 s,kerosene supply was shut down,thereby rapidly decreasing P4 and P5.

    Fig.4.Static pressure history under case 2: (a) P1 and P2;and (b) P3-P5.

    Fig.5 shows the high-frequency dynamic pressure history obtained at PCB0,PCB1,and PCB2 during the ignition process of the rotating detonation.In this study,the high-frequency dynamic pressure data acquisition was triggered 4 s after the signal acquisition,that is,the time corresponding to 0 s of the high-frequency dynamic pressure data corresponds to 4 s of the overall experimental signal acquisition data.The hot jet tube was ignited at 387.4 ms.After 0.9 ms,the hydrogen and oxygen mixture in the hot jet tube was ignited to produce a thermal jet.PCB0 captured the ignition pressure peak of 1.02 MPa.After Δt,PCB1 recorded a peak pressure of 0.99 MPa at 388.3 ms owing to the hot jet generated in the hot jet tube and injected into the RDC.After the injection of the hot jet into the RDC,the pressure peak oscillated irregularly.Here,deflagration combustion occurred in the combustor,and the RDW was not immediately formed.After t,under the combined action of the shock wave reflection,friction,and turbulence in the combustion chamber wall,and continuous injection of fresh KFRG/air combustible mixtures,the deflagration wave developed into a detonation wave.PCB1 captured a pressure peak of 0.72 MPa at 390.2 ms.Subsequently,an RDW was formed in the RDC and gradually stabilized,from which the deflagration to detonation transition(DDT)process was completed,indicating that tlasted approximately 2 ms.

    Fig.5.High-frequency dynamic pressure of PCB0,PCB1,and PCB2 during the ignition process.

    Fig.6 shows the RDW pressure at the ignition process in case 6,case 8,and case 9.The time for the formation of the RDW was identified as the delay time from the formation of the hot jet to the formation of a stable RDW,which was equal to 1.1 ms,1.6 ms,and 2.4 ms for case 6,case 8,and case 9,respectively.As the kerosene mass flow rate was increased,the tneeded for RDW formation from the DDT gradually increased because of the increased product volumes of the KFRG.As the mass ratio of the fuel injected into the RDC increased,the mass proportion of the oxidizer decreased,and the overall reactivity of the combustible premixed gas decreased.Thus,the DDT process required a longer time for completion.

    Fig.6.Rotating detonation delay time in case 6,case 8,and case 9.

    3.2.RDW propagation process of KFRG

    The analysis of the RDW propagation process of the KFRG revealed the significant effect of the equivalence ratio on its propagation mode.With the equivalence ratio of 0.62-0.79,a mixed propagation mode dominated by two-counter rotating waves and supplemented by a single wave was noted.With the equivalence ratio of 0.85-0.87,a single-wave RDW propagation mode was observed.The characteristics of the rotating detonation under different propagation modes are analyzed as follows.

    3.2.1.Mixed RDW propagation mode

    A mixed RDW propagation mode indicates that single-wave and two-counter RDW propagation modes alternately appear in one experimental cycle.Based on case 2,the propagation of the RDW of the KFRG was analyzed with an equivalence ratio of 0.76.In the single-wave propagation mode,the RDW rotated and propagated in one direction within the RDC,and the pressures at PCB1 and PCB2 alternately and regularly increased.In the two-counter RDW propagation mode,two RDWs propagated towards each other in the RDC,and the two waves collided twice per cycle.The PCB pressure near the collision point became similar to that in the single-wave mode,and the PCB at other positions captured two pressure peaks within one cycle.Fig.7 shows the high-frequency dynamic pressure with respect to time.The dynamic RDW pressure peak captured by PCB1 constantly fluctuated owing to two reasons: the uneven composition of KFRG in RDC,and changes in the RDW propagation between a single-wave and two-counter rotating wave modes.Under a single-wave propagation mode,the pressure peak was relatively stable and fluctuated at approximately 1.2 MPa.Under the two-counter rotating wave propagation mode,and the collision point at PCB1,a high-pressure peak was generated.However,for a collision points different from PCB1,the pressure peak was low.

    Fig.8 illustrates the pressure at PCB1 and PCB2 in the singlewave propagation mode.At 486.6-493.95 ms,the RDW propagated for 25 cycles in the RDC with the pressure peak in the order of a→b→c→d,which is a typical codirectional rotating mode.As viewed from the RDC exit end,i.e.,counterclockwise,the rotating direction was from PCB1 to PCB2.According to the time and distance of the RDW propagating for 25 cycles,its average propagation velocity was 1191 m/s,as calculated using Eq.(1) and Eq.(2).

    where L is the total propagation distance of the RDW,K is the number of detonation wave propagation periods,and D is the outer diameter of the RDC,D=0.11 m.

    The transient propagation and average frequency of the RDW was calculated using Eq.(3) and Eq.(4),respectively.

    where f is the dominant frequency and n is the number of cycles.In Fig.8,the transient frequency of a Δt cycle was 3846 Hz,and the average dominant frequency of the RDW propagation for 25 cycles was 3486 Hz.

    Based on the analysis of the RDW pressure peak sequence,calculated wave velocity,and dominant frequency,RDW followed a single-wave propagation mode during this period.The number of waves can be verified by Eq.(5) [32]:

    where n is the number of RDW fronts;θ is the radian angle between PCB1 and PCB2 with a value of 2π/3;and tand tare the interval between two pressure peaks in one cycle.Here,n=1,that is,the number of wavefronts was 1,and the RDW propagated following a single-wave mode.The average pressure peak of the RDW propagating in the single-wave mode was 0.94 MPa.

    Fig.9 shows the diagrams obtained by the fast Fourier transform(FFT)and short-time Fourier transform(STFT)of the pressure data at PCB1.The sampling time corresponds to the entire RDC operation.The dominant frequency of the RDW propagating in RDC was 3229.9 Hz,as seen in Fig.9(a).There was another peak at 6513 Hz,which was close to the secondary harmonic frequency of 6459.8 Hz.The average propagation velocity of the RDW was estimated as 1115.6 m/s,as follows:

    Fig.7.High-frequency dynamic pressure curve of PCB1.

    Fig.8.High-frequency dynamic pressure at PCB1 and PCB2 in single-wave propagation mode.

    Fig.9.Frequency domain analysis diagram of RDW at PCB1: (a) FFT results;(b) STFT results.

    where f is the dominant frequency.A relatively low RDW velocity was noted owing to the uneven mixing of kerosene and oxygen in the KFRG generator,which resulted in the uneven composition of the KFRG generated by the combustion.Moreover,this can be ascribed to the uneven mixing of the KFRG and oxidizer,resulting in the RDW velocity and pressure loss.

    By magnifying the curve in the red box of the FFT result in Fig.9(a),a slightly lower peak can be observed at 3337.5 Hz,whose frequency was higher than the dominant frequency of 3229.9 Hz.Using Eq.(6),the average propagation velocity was 1152.8 m/s,which was equivalent to the RDW velocity in the single-wave mode with a difference of 3%.This indicates the presence of multiple propagation modes in the RDC.As the single-wave frequency was higher than the double-wave frequency,the dominant frequency corresponded to that of the double wave mode,and the subfrequency is that of the single-wave mode.Moreover,according to the power spectrum density,the detonation wave in this working condition mainly rotates following a double-wave mode with a small number of single-wave modes.After detonation,the dominant and secondary harmonic frequencies were approximately 3200 and 6400 Hz,respectively,as seen in Fig.9(b)consistent with the above analysis.

    Under all conditions in this experiment,there are mainly two propagation modes of the RDW formed by KFRG,stable singlewave,and two-counter rotating waves modes.The stable single propagation wave mode is described above.The pressure peaks captured by PCB1 and PCB2 appear alternately,and the rotating direction is determined according to the occurrence order.The number of wave heads in the RDC was verified according to Eq.(5).In the two-counter rotating propagation wave mode,the pressure peak near the collision point(PCB1)was similar to that in a singlewave mode;however,the peak pressure was significantly higher.Two pressure peaks were measured at other positions(PCB2,PCB3)in each cycle,which were lower than that in the single-wave mode.

    As can be seen from Fig.10,the propagation mode of the RDW is constantly transformed because of initiation,and oscillating peak pressure captured by PCB1 and PCB2 in the working period.When the peak value at PCB1 was higher,the two-counter rotating wave propagation mode and collision point were located near PCB1.Similarly,if the pressure peak at PCB2 was higher,the two-counter rotating wave propagation mode and collision point were located near PCB2.If the two pressure peak values were similar,a singlewave propagation mode was achieved.According to Eq.(5) and the PCB pressure,there was no two-codirectional RDW propagation mode,that is,two RDWs rotate and propagate in the same direction within the RDC,in this case.The two-counter detonation rotating wave propagation mode appeared frequently during the working period of the RDC,and the collision point was located near PCB1.Meanwhile,the single-wave mode propagation appeared less frequently.As shown in Fig.11,the RDW comprised a two-counter rotating wave propagation between 433.5 ms and 439 ms,which produced a high-pressure peak at the collision point of two-counter rotating waves.In addition,lower and similar pressure peaks were obtained by the RDW at the other positions.PCB1 captured a highpressure peak of 2.12 MPa,and PCB2 captured two low pressure peaks of approximately 0.3 MPa.Therefore,the collision point of the two-counter rotating waves was near PCB1.

    Fig.10.High-frequency dynamic pressure curves at PCB1 and PCB2.

    Fig.11.High-frequency dynamic pressure curves at PCB1 and PCB2: two-counter rotating wave propagation mode.

    Fig.12.Propagation of two counter rotating waves when the collision point is located in PCB1: (a) PCB1 and PCB2 pressure curve;(b) Two-counter RDW propagation process.

    The propagation process of two counter rotating waves is shown in Fig.12.In Fig.12(a),the counter rotating waves Aand Bwhich were formed after the previous cycle,collided at PCB1 at 444.7 ms when PCB1 captured a pressure peak of 2.3 MPa,which was significantly higher than that of the normal single wave.As shown in Fig.12(b),the leading shock waves of the two RDW fronts transmitted into the detonation product region of the other waves after their collision.Here,the wave degenerated into shock waves Aand B,and the wave velocity and pressure peak values correspondingly decreased.After the collision,the two transmitted shock waves continually propagated in the original direction.When fresh combustible mixtures were injected into the RDC,transmitted shock wave Breformed into RDW B.When RDW Breached PCB2 after traversing the circumference for 120,the pressure peak was 0.3 MPa.Subsequently,RDWs Aand Bcollided again at the position opposite to PCB1.Transmitted shock wave Agradually increased to Aand reached PCB2 after traversing the circumference for 60.Here,the pressure peak was slightly lower than Bbecause the transmitted shock wave was not instantly strengthened into a stable RDW after the collision process of the two waves.After the double-wave collision,the transmitted shock wave entered the post-wave reaction zone of the opposite RDW;thus,the pressure peak value of RDW in the two-counter rotating wave propagation mode was lower than that of single wave.With the exhaustion of the detonation products,the height of the fresh combustible mixture layer gradually increased.The RDW propagated to this region prior to the filling of the fresh fuel layer in the double-wave propagation mode.This consistently decreased the fuel layer height in the double-wave mode than that in the singlewave mode.Additionally,this resulted in a poor mixing process.

    In a double-wave collision process,the collision point consumes the energy of the detonation wave;thus,the two-counter rotating wave propagation mode had a weak detonation state.In engineering applications,the collision propagation of RDWs should be avoided.After RDWs Aand Bcollided at 0,RDW Bpassed through PCB2 after rotating for 120,and RDW Apassed through PCB2 after rotating for 240.The two waves collided again at 0,thereby ending the double-wave propagation in the RDC for one cycle.In Fig.12(a),within one cycle,PCB1 captured the pressure peak during collision,and PCB2 captured the pressure peak when two counter rotating waves passed successively.The time intervals t,t,and tof the three pressure peaks were similar at 92,104,and 92 μs,respectively.Every interval corresponds to the time required for the RDW to rotate for 120of the RDC circumference.The blue box in the figure shows that the collision point was offset away from PCB1;hence,PCB1 captured two pressure peaks in a short time interval.In the two-counter rotating wave propagation process,the collision points of the two counter rotating waves constantly shifted.Most of the collision points in this working condition were stable at PCB1 owing to the differences in the intensity and velocity of the two waves after collision.

    The propagation mode of the RDW changed from two-counter rotating wave mode to single-wave mode,as shown in Fig.13.From 478 to 479.5 ms,the collision point of RDW continuously changed at PCB1 and PCB2,and the propagation subsequently became unstable and attenuated into a deflagration wave.After 2 ms,the RDW was reformed as fresh fuel was injected into the RDC.At 481.8 ms,the two waves collided at PCB2,and PCB1 captured two lower pressure peaks,which were the two counter RDWs after the collision.In the next cycle,PCB1 captured the pressure peak of 1.2 MPa at 482.3 ms for 2/3rd of the cycle,indicating that the detonation wave has rotated for 240after collision at PCB2 and has passed through PCB1.Thereafter,a pressure peak was captured for each PCB in each cycle,indicating a single counterclockwise-rotating wave in RDC observed from the exit of the combustion chamber.Consequently,the propagation of the RDW changed from the two-counter rotating wave mode to an unstable single wave in a counterclockwise rotating mode.The propagation gradually stabilized after several cycles at 486.3 ms.During the shift of the propagation mode and after the collision of the two-counter rotating waves,a portion of the RDW entered the hot jet tube,where there was no fresh combustible mixture after ignition,eventually causing detachment and attenuation of the RDW.The attenuated RDW was completely annihilated when it collided with a strong wave that rotated counterclockwise.

    The RDW propagated in a single-wave mode,as shown in Fig.14.The average pressure peak of the RDW was 0.94 MPa,which was lower than the peak pressure at the collision point in the twocounter rotating wave mode and higher than the peak pressure of each RDW.According to the time interval between the two peak points,the velocity range of the RDW propagation was calculated to be 987-1439 m/s.The velocity of the RDW remained unstable and fluctuated considerably in each cycle,which may be attributed to the uneven composition of the KFRG,and uneven mixing of the fuel and oxidant injection at the RDC head.

    As shown in Fig.15,the RDW propagation mode switched from the single-wave mode to two-counter rotating wave mode at 493-506 ms.At 493 ms,the RDW propagated following singlewave mode that rotated counterclockwise as observed from the end of combustion chamber.The RDW propagated in the RDC at 494 ms,as magnified in the black box.When the RDW passed through PCB2,a pressure peak of 0.95 MPa was collected,followed by three low-pressure peaks of approximately 0.2 MPa.The product gas with high temperature behind the detonation front ignited the fresh combustible mixtures and formed three deflagration waves,which induced the RDW propagation modal conversion.In the subsequent two cycles,the three deflagration waves were gradually superimposed and strengthened to form a new RDW.The propagation mode switched into the two-counter rotating wave mode,and the pressure peak decreased.The collision point constantly varied during the modal conversion process,and the position of the collision was determined according to the high-pressure peak generated during the collision of the two waves.At 495.9 ms,the collision point was at PCB2.At 499.2 ms,the collision point changed to PCB1 until at 504 ms,when the collision point became stable and the mode conversion was terminated;the RDW mode changed from the single-wave to a stable two-counter rotating wave mode.

    Fig.13.High-frequency dynamic pressure curves at PCB1 and PCB2: two-counter rotating wave mode to single-wave mode.

    Fig.14.High-frequency dynamic pressure/wave velocity curves at PCB1 and PCB2: single-wave mode.

    Fig.15.High-frequency dynamic pressure curves at PCB1 and PCB2: single-wave mode to two-counter rotating wave mode.

    When the RDW passed through the location of the hot jet tube,multiple deflagration waves were generated owing to the RDC wall reflection and fresh combustible mixtures.During the propagation,these deflagration waves were constantly superimposed to form a reverse reflected shock wave,which gradually developed into a reverse RDW.After a period of unstable collision propagation,the single-wave mode switched to a relatively stable two-counter rotating wave mode.At 494.5 ms,the peak pressure of 1.25 MPa was captured by PCB1,and two low-pressure peak signals were captured at PCB2 in the two cycles before and after it.Therefore,the RDW propagating in the opposite direction to the original single wave was generated near PCB1,where the hot jet tube was located during injection.This confirmed the correctness of the above viewpoint.

    As shown in Fig.16,the RDW propagation in the RDC at 564-581 ms underwent five stages: two-counter rotating wave,unstable mode,single-wave,unstable mode,and two-counter rotating wave modes.Throughout the RDC cycle,the RDW propagation mode was constantly switching.The instability of the RDW propagation occurred during the entire RDC operation was partly related to the decoupling and reinitiation phenomena shown in Fig.17.After the decoupling phenomenon,reinitiation was realized in a relatively short time.The detonation and reinitiation can be attributed to the detonation wave attenuated to the transmitted shock wave and high pressure at the collision area of the two waves,which blocked the injection of fresh combustible mixtures,after the collision of two waves.Furthermore,the unstable operation of the gas generator led to the uneven composition of the KFRG and unsteady mass flow rate of gas.Therefore,the transmitted shock wave in the RDC cannot be immediately enhanced into an RDW.Instead of propagating as a deflagration wave for several periods with the normal injection of the fresh combustible mixtures into the RDC,the wave underwent a DDT process and gradually reinforced to form RDW,thereby maintaining self-sustenance and stable propagation.As the generation of the KFRG and injection of the fuel oxidant in the RDC are dynamic processes,the RDW propagation process appeared as modal switching or a detonation discontinuity in the unstable mode.

    Fig.16.High-frequency dynamic pressure curves at PCB1 and PCB2: mode switching.

    Fig.17.High-frequency dynamic pressure curves at PCB1 and PCB2: detonation discontinuity.

    To eliminate errors,the experiment was repeated for case 2.The experimental results were stable under the same working conditions.The main frequencies of the RDW in the three experiments were 3249 Hz,3229.9 Hz,and 3221 Hz,respectively,and the calculated wave velocities were 1122.2 m/s,1115.6 m/s,and 1112.5 m/s,respectively.The RDW propagation modes were mixed mode with two-counter rotating waves as the main mode and single-wave as the auxiliary mode.The median value of the wave velocity was taken as the RDW propagation velocity for the followup analysis,and the fluctuations of the RDW velocity were within 0.6%.The results obtained by the experiments were stable and reproducible.

    Table 4 shows the product composition and mass fraction of the KFRG after the enriched combustion.The data were calculated by the CEA software.The theoretical C-J velocity of the KFRG was calculated to be 1967.2 m/s.According to Eq.(2) and Eq.(6),the calculated RDW propagation velocity was 60.5% and 56.7% of the theoretical value,respectively.The RDW had the highest wave velocity in the single-wave mode,which greatly deviated from the calculated and theoretical wave velocity.The reasons for this deviation are as follows.First,the theoretical calculated value had anerror owing to the complex composition of KFRG.Second,KFRG had an uneven distribution owing to the combustion of the KFRG generator.Third,the uneven mixing of KFRG and oxidizer resulted in a loss in the propagation velocity and pressure of the RDW.

    Table 4.Theoretical results of the KFRG composition.

    3.2.2.Single-wave propagation mode

    Based on case 10,the RDW propagated stably in a single-wave propagation mode during the experiment,and the pressure peak appeared alternately and regularly.Fig.18(a) shows the frequency diagram obtained by the FFT analysis at PCB3 with the sampling time covering the entire working process of the RDC.The dominant frequency of the RDW propagating in the RDC was 4740 Hz.According to Eq.(6),the calculated average propagation velocity of the RDW was 1637.2 m/s.Meanwhile,the theoretical RDW propagation velocity was 2048 m/s;the experimental value was 79.9% of the theoretical value,which was significantly higher than the RDW velocity in case 2.Fig.18(b) shows the STFT analysis at PCB3.The dominant frequency of the RDW propagation was approximately 4800 Hz.According to dynamic pressure and Eq.(5),RDW propagated in a single-wave mode under this condition,and the wave velocity was significantly higher than that of the two-counter rotating wave mode.

    Fig.18.Frequency domain analysis diagram of the RDW in case 10: (a) FFT results;(b) STFT results.

    Fig.19.High-frequency dynamic pressure/wave velocity curve of PCB1 and PCB2.

    Fig.19 shows the RDW dynamic pressure and velocity at 735-748 ms.The wave velocity distribution was in the range of 1644.8-1918.9 m/s,with an average value of 1768.7 m/s.This value was higher than that obtained in the FFT analysis,indicating the highly stable propagation of the RDW during this period.The average pressure peak pressure at PCB2 was 0.84 MPa.

    3.3.Influence of different factors on the RDW propagation

    3.3.1.Influence of the oxygen mass fraction in the oxidizer on the RDW propagation

    The experiments were conducted under different conditions to investigate the influence of the oxygen mass fraction in the oxidizer on the RDW propagation.Table 1 shows the propagation velocity of the RDW with different oxygen mass fractions in the oxidizer.The RDW propagated mainly through the two-counter rotating wave mode,as indicated by the dominant frequency obtained by the FFT analysis of the dynamic pressure shows.The RDW wave velocity in this mode was calculated by Eq.(6).As shown in Table 1,the precombustion oxygen,kerosene,and air mass flow rates were fixed,whereas the oxygen mass fraction in the oxidizer was varied from 36% to 44%,and the equivalence ratio was decreased from 0.91 to 0.65.Stable continuous rotating detonation was achieved in the RDC.As the oxygen mass fraction was increased,the average propagation velocity of the RDW increased from 903.2 to 1166.8 m/s.

    Fig.20 shows the dynamic pressure at PCB1 in case 1-case 3(mass flow of 36%,40%,and 44%,respectively).When the oxygen mass fraction was 36%,the delay time from the ignition of the hot jet tube to the RDW in the RDC was 20.3 ms,which was significantly higher than that when the oxygen mass fraction was 40%(2 ms) and 44% (1.6 ms).Therefore,under a low oxygen mass fraction,the initiation delay increased,suggesting the increased initiation difficulty.Under an oxygen mass fraction of 36%,the average pressure peak was only 0.3 MPa,indicating an unstable detonation mode.Thus,an increase in the oxygen mass fraction improved the activity of the combustible mixtures,which accelerated the detonation combustion reaction rate and reduced the loss of RDW velocity.The experimental results confirmed the results of existing studies [27].

    3.3.2.Influence of the jet fuel mass flow rate on the RDW propagation

    Table 2 shows the experimental parameters and results of the rotating detonation with the change in the kerosene mass flow rate for cases 4-10 with a fixed precombustion oxygen mass flow rate of 12 g/s,and oxygen mass fraction of 44%.As the mass flow rate of kerosene was increased from 53 to 74.6 g/s,the equivalence ratio also increased from 0.62 to 0.87,resulting in a stable RDW in the RDC.Based on the analysis of the PCB pressure and FFT results,when the equivalence ratio was 0.62-0.79,the RDW propagation mode was dominated by the two-counter rotating wave mode supplemented by single wave.When the equivalence ratio was 0.85-0.87,the RDW followed a single-wave propagation mode.The dominant frequency of the mixed and single-wave modes was employed in the FFT analysis for the wave velocity calculation of the RDW using Eq.(6).The RDW velocity distribution ranged from 1148.1 to 1637.2 m/s.As shown in Fig.21,the RDW velocity increased with the fuel mass flow rate.Under the single-wave mode,the RDW wave velocity increased substantially.When the velocity of the RDW increased with the equivalence ratio,the increase in the kerosene mass flow rate improved the velocity of RDW,which is consistent with a previous report [30].

    Fig.20.High-frequency dynamic pressure curve at PCB1 in case 1-case 3.

    Fig.21.RDW velocity and temperature of the KFRG as a function of kerosene mass flow rate.

    Fig.22.Temperature change of the KFRG collecting chamber.

    With the increase in the kerosene mass flow rate,the temperature of the KFRG decreased slowly.This is attributed to fixed amount of heat released in the oxygen and kerosene reaction.As the kerosene mass flow rate was increased,the temperature of the product decreased.Simultaneously,the pyrolysis reaction of kerosene,which was known to be endothermic,decreased the temperature of the reaction product.Owing to the minimal temperature variation range of the KFRG,the influence of temperature can be ignored for case 4-case 10.

    3.3.3.Influence of the KFRG temperature on the RDW propagation

    Table 3 shows the experimental parameters and results of the RDW under different precombustion oxygen mass flow rates for cases 11-14.The temperature of the KFRG increased with the precombustion oxygen mass flow rate.In these cases,the mass flow rates of the kerosene and oxidizer were fixed,and the oxygen mass fraction was 32%.The precombustion oxygen mass flow rate increased from 12 to 21.6 g/s,the equivalence ratio decreased from 0.93 to 0.89,and the KFRG temperature increased from 480 to 725 K.

    Fig.22 shows the change in temperature of the KFRG collecting chamber over time.At 2.6 s,the gas generator began to ignite.The oxygen was injected into the gas generator and reacted with the residual fuel in the combustor,causing the gas temperature in the collecting chamber to increase at 3.1 s.The kerosene was injected,and fuel-rich combustion was initiated in the gas generator.When the kerosene initially entered the gas generator combustion chamber,the temperature of the gas in the collecting chamber slightly decreased owing to the heat absorbed by the evaporation and pyrolysis of kerosene.The temperature of the gas gradually stabilized with the gas generator in a steady state.At 5.1 s,the kerosene injection ceased,and the oxygen supply was turned off sequentially.The oxidant volume increased,and the temperature briefly rose in the combustion chamber.After the oxygen was turned off,the temperature gradually decreased.

    Fig.23.Pressure change curves of combustor and gas collecting chamber as a function of the precombustion oxygen mass flow rate.

    Fig.23 shows the pressure curve of the combustor and KFRG collecting chamber with the change of the precombustion oxygen mass flow rate when the KFRG generator was operating.The fuel and oxidizer used by the gas generator were kerosene and precombustion oxygen,respectively.During the KFRG generator operation,the change in the kerosene mass flow rate had a weak influence on the pressure of the combustor and KFRG collecting chamber.In particular,the pressure change was only in the range of 0.8-0.9 MPa,and its influence on the gas temperature was weak.Unlike the influence of the kerosene mass flow rate on the gas generator,the change in the precombustion oxygen mass flow rate greatly affected the internal pressure.As the precombustion oxygen mass flow rate was increased from 12 to 21.6 g/s,the pressure in the combustion chamber increased from 1.18 to 2.16 MPa,and the pressure in the gas collecting chamber increased from 0.88 to 1.56 MPa.

    As the temperature of KFRG was increased by increasing the precombustion oxygen mass flow rate,the composition of the KFRG varied.With the increase in temperature,the mass fraction of C decreased,which was difficult to initiate.The mass fraction of Hincreased,which was easy to initiate,thereby promoting the increase of the RDW velocity.Precombustion oxygen mass flow rate also has a great influence on the composition of KFRG from the analysis by the CEA software,as shown in Fig.24.The components with highest mass fraction in KFRG were C,whose mass fraction decreased with the increase of the precombustion oxygen mass flow rate,and CH,which also exhibited a decreasing law.CO,CO,H,and HO exhibited an increasing trend.The mass fraction of CO exhibited the most pronounced increase.The theoretical temperature of the kerosene combustion gas increased from 894 to 1000 K.

    Fig.24.Theoretical composition and temperature of the product as a function of the precombustion oxygen mass flow rate.

    Fig.25.RDW velocity as function of the KFRG temperature.

    Fig.25 shows the change in the RDW velocity with the gas temperature.In case 11,the KFRG failed to initiate,whereas in case 12-case 14,the KFRG achieved a stable self-sustaining RDW propagation.The minimum oxygen mass fraction of the RDC initiation was 32%,which was lower than the minimum oxygen mass fraction (39.2%) needed for liquid kerosene initiation [33].According to the PCB and FFT analysis,RDW propagation followed a single-wave mode.Based on Eq.(6),the wave velocity increased from 967.1 to 1141.9 m/s,indicating that the KFRG temperature had a significant influence on the RDW initiation and propagation.With the increase in the KFRG temperature,the activity of the combustible medium increased,the rate of chemical reaction accelerated,and the RDW wave velocity increased.Thus,temperature increase has a gain effect on the propagation of RDW.

    4.Conclusions

    In this study,the effects of the basic parameters on the rotating detonation KFRG were investigated.The following conclusions were drawn from the study.

    (1) The use of KFRG as the fuel resulted in complex RDW propagation modes.The continuous switching of the single-wave and two-counter rotating wave modes was observed in most experimental conditions,and stable single-wave mode was noted in a few experimental conditions.When the oxygen mass fraction was increased from 36 to 44%,the propagation velocity of the RDW increased from 903.2 to 1166.8 m/s.When the mass flow rate of kerosene was increased from 53 to 74.6 g/s,the propagation velocity of the RDW increased from 1148.1 to 1637.2 m/s.When the equivalence ratio was within 0.62-0.79,the RDW propagation mode predominantly followed the two-counter rotating wave mode.Meanwhile,with the equivalence ratio of 0.85-0.87,the RDW followed a single-wave propagation mode.When the KFRG temperature was 480 K,the KFRG did not initiate the RDW formation.When the temperature increased from 480 to 725 K,the propagation velocity of the RDW increased from 967.1 to 1141.9 m/s.

    (2) When the kerosene mass flow rate was 57.1 g/s,the gas generator oxygen mass flow rate was 12 g/s,the oxidant mass flow rate was 594.6 g/s,the oxygen mass fraction was 40%,and the equivalence ratio was 0.76.Thus,a stable RDW with a dominant two-counter rotating wave mode and auxiliary single-wave mode was obtained.The main frequency and velocity of the RDW were 3229.9 Hz and 1115.6 m/s,respectively.Correspondingly,when the kerosene mass flow rate was 74.6 g/s,the gas generator oxygen mass flow rate,oxidant mass flow rate,oxygen mass fraction,and equivalence ratio were 12 g/s,646.6 g/s,44%,and 0.87,respectively,thereby a stable single-wave propagation mode was obtained.The main frequency of this RDW was 4740 Hz,and its propagation velocity was 1637.2 m/s.

    (3) The initiation of KFRG cannot be achieved with low oxygen mass fraction and KFRG temperature of KFRG.The propagation velocity of the RDW can be improved by increasing the oxygen fraction in oxygen-rich air,kerosene mass flow rate,and temperature of KFRG.Compared with liquid kerosene,KFRG can achieve rotating detonation at a lower oxygen mass fraction of 32%.

    This study investigated the RDW propagation mode of KFRG,and its factors.Further research need to be carried out to study the influence of the geometric structure of the CRDE on RDW,which will lay a foundation for the engineering research of the rotating detonation engine.Further,the KFRG composition should be determined,and its effects on the RDW propagation should be analyzed in detail.The RDW propagation modes should be investigated with high-speed photography.Moreover,it is necessary to optimize the structure of the CRDE and explore the limit working range of it,the detonation methods at lower oxygen levels and even air conditions should be investigated.

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    国产伦理片在线播放av一区| 汤姆久久久久久久影院中文字幕 | 免费观看的影片在线观看| 超碰97精品在线观看| 久久久精品免费免费高清| 国产老妇女一区| 日本与韩国留学比较| 亚洲欧洲日产国产| 中文乱码字字幕精品一区二区三区 | 小蜜桃在线观看免费完整版高清| 国产乱人视频| 亚洲国产欧美人成| 亚洲人与动物交配视频| 亚洲av成人精品一区久久| 三级男女做爰猛烈吃奶摸视频| 亚洲成人中文字幕在线播放| 久久6这里有精品| 国产精品一二三区在线看| 99久久精品一区二区三区| 亚洲第一区二区三区不卡| www.色视频.com| 丝袜美腿在线中文| 日本三级黄在线观看| 免费大片黄手机在线观看| 中文欧美无线码| 激情 狠狠 欧美| 久99久视频精品免费| 天堂网av新在线| 精品国产露脸久久av麻豆 | 韩国高清视频一区二区三区| 久久亚洲国产成人精品v| 肉色欧美久久久久久久蜜桃 | 大片免费播放器 马上看| 成人亚洲精品av一区二区| 高清视频免费观看一区二区 | 99热这里只有精品一区| 欧美区成人在线视频| 日韩精品有码人妻一区| 国产精品久久视频播放| 熟妇人妻不卡中文字幕| 一级毛片黄色毛片免费观看视频| 亚洲av成人av| 亚洲国产精品成人久久小说| 80岁老熟妇乱子伦牲交| 2021天堂中文幕一二区在线观| 97在线视频观看| 简卡轻食公司| 国产伦理片在线播放av一区| 久久国产乱子免费精品| 又大又黄又爽视频免费| 国产精品三级大全| 国产黄色视频一区二区在线观看| 91精品一卡2卡3卡4卡| 国产高清不卡午夜福利| 免费av毛片视频| 午夜福利在线观看吧| 如何舔出高潮| 免费无遮挡裸体视频| 少妇被粗大猛烈的视频| 日韩欧美精品免费久久| 极品少妇高潮喷水抽搐| 三级国产精品欧美在线观看| 精品一区在线观看国产| 亚洲一区高清亚洲精品| 国产黄片美女视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国语对白做爰xxxⅹ性视频网站| 亚洲av在线观看美女高潮| 看十八女毛片水多多多| 免费电影在线观看免费观看| 亚洲国产精品成人久久小说| 人妻一区二区av| 久久99精品国语久久久| 69人妻影院| 久久草成人影院| 特级一级黄色大片| 亚洲欧美成人综合另类久久久| 日韩av在线大香蕉| 国产色爽女视频免费观看| 波野结衣二区三区在线| 777米奇影视久久| 国产精品嫩草影院av在线观看| 国产成人a区在线观看| 国内精品美女久久久久久| 亚洲国产欧美在线一区| 亚洲精品成人av观看孕妇| 看黄色毛片网站| 亚洲精品久久久久久婷婷小说| 男女边吃奶边做爰视频| 啦啦啦啦在线视频资源| 十八禁网站网址无遮挡 | 久久精品国产鲁丝片午夜精品| 三级男女做爰猛烈吃奶摸视频| av免费在线看不卡| 最近的中文字幕免费完整| 国产在线男女| 国产 一区精品| 免费观看性生交大片5| 国产精品蜜桃在线观看| 欧美区成人在线视频| 国产欧美另类精品又又久久亚洲欧美| 少妇人妻精品综合一区二区| 麻豆成人av视频| 3wmmmm亚洲av在线观看| 国产亚洲精品久久久com| 精品一区二区三卡| 黄片无遮挡物在线观看| 亚洲高清免费不卡视频| 欧美潮喷喷水| 亚洲精品自拍成人| 亚洲av一区综合| 18禁动态无遮挡网站| 白带黄色成豆腐渣| 亚洲一级一片aⅴ在线观看| 一本一本综合久久| 亚洲精品一二三| 国产一区二区三区av在线| 国产高清有码在线观看视频| 观看免费一级毛片| 久久久久久九九精品二区国产| 综合色av麻豆| 中文字幕免费在线视频6| 三级男女做爰猛烈吃奶摸视频| 五月玫瑰六月丁香| 国产乱来视频区| 免费看光身美女| 最新中文字幕久久久久| 好男人在线观看高清免费视频| 波多野结衣巨乳人妻| 2022亚洲国产成人精品| 午夜福利视频精品| 日韩欧美一区视频在线观看 | 国产精品蜜桃在线观看| 三级国产精品欧美在线观看| 国产精品女同一区二区软件| 亚洲av不卡在线观看| 国产精品一区二区在线观看99 | 欧美日韩亚洲高清精品| 久久久久久久久久人人人人人人| 久久精品综合一区二区三区| 亚洲av不卡在线观看| 99视频精品全部免费 在线| 久久人人爽人人爽人人片va| 欧美日韩亚洲高清精品| 一区二区三区高清视频在线| 欧美激情在线99| av黄色大香蕉| 久久久午夜欧美精品| 搡老妇女老女人老熟妇| 精品人妻一区二区三区麻豆| 国产精品国产三级专区第一集| 综合色丁香网| 国产大屁股一区二区在线视频| 91久久精品电影网| 丝袜美腿在线中文| 亚洲人成网站在线观看播放| 草草在线视频免费看| 久久国内精品自在自线图片| 特级一级黄色大片| 又粗又硬又长又爽又黄的视频| 午夜激情欧美在线| 色网站视频免费| 一边亲一边摸免费视频| 中文字幕免费在线视频6| 天堂中文最新版在线下载 | 欧美zozozo另类| 老师上课跳d突然被开到最大视频| av网站免费在线观看视频 | 亚洲丝袜综合中文字幕| 日韩伦理黄色片| 美女cb高潮喷水在线观看| 国产伦精品一区二区三区四那| 亚洲精品日韩av片在线观看| 国产精品久久久久久精品电影小说 | 午夜爱爱视频在线播放| h日本视频在线播放| 国产成人91sexporn| av女优亚洲男人天堂| 国产伦一二天堂av在线观看| 亚洲精品色激情综合| 97精品久久久久久久久久精品| 十八禁网站网址无遮挡 | 99久久精品热视频| 熟妇人妻不卡中文字幕| 久久亚洲国产成人精品v| 欧美最新免费一区二区三区| 大话2 男鬼变身卡| 国产乱人视频| 国产精品av视频在线免费观看| 亚洲av电影在线观看一区二区三区 | 亚洲av国产av综合av卡| 男女边吃奶边做爰视频| 亚洲人与动物交配视频| 一级爰片在线观看| 最近中文字幕2019免费版| 久久精品国产自在天天线| 欧美日韩视频高清一区二区三区二| 欧美xxxx性猛交bbbb| 精品一区二区三卡| 久久精品人妻少妇| 国产黄a三级三级三级人| 2021少妇久久久久久久久久久| 乱码一卡2卡4卡精品| 精品国产一区二区三区久久久樱花 | 爱豆传媒免费全集在线观看| 国产成人a∨麻豆精品| 国产精品蜜桃在线观看| 国产一区二区三区综合在线观看 | 18禁在线无遮挡免费观看视频| 成人鲁丝片一二三区免费| 80岁老熟妇乱子伦牲交| 国产永久视频网站| 九色成人免费人妻av| 五月天丁香电影| 99久久精品一区二区三区| 十八禁网站网址无遮挡 | 亚洲真实伦在线观看| 国产国拍精品亚洲av在线观看| 五月天丁香电影| 又爽又黄无遮挡网站| 午夜免费男女啪啪视频观看| 91午夜精品亚洲一区二区三区| 久久久久久久久大av| 色播亚洲综合网| 亚洲av不卡在线观看| 欧美激情久久久久久爽电影| 欧美日韩在线观看h| 日韩制服骚丝袜av| 亚洲国产欧美人成| 最近2019中文字幕mv第一页| 可以在线观看毛片的网站| 免费观看av网站的网址| xxx大片免费视频| 久久久久精品性色| 国产欧美日韩精品一区二区| 男的添女的下面高潮视频| 精品午夜福利在线看| 老司机影院成人| 免费人成在线观看视频色| 一区二区三区高清视频在线| 亚洲精品国产成人久久av| 久久久久网色| 精品国内亚洲2022精品成人| 1000部很黄的大片| 国产黄a三级三级三级人| 国产黄色视频一区二区在线观看| 欧美变态另类bdsm刘玥| 日韩av在线免费看完整版不卡| 久久久精品欧美日韩精品| 免费在线观看成人毛片| av在线亚洲专区| 午夜福利在线在线| 国产高清国产精品国产三级 | 国产成人91sexporn| 97热精品久久久久久| 高清在线视频一区二区三区| 欧美xxⅹ黑人| 色综合亚洲欧美另类图片| 嫩草影院入口| 伊人久久国产一区二区| 日本黄大片高清| 亚洲18禁久久av| 国产老妇伦熟女老妇高清| 午夜免费男女啪啪视频观看| 婷婷色麻豆天堂久久| 偷拍熟女少妇极品色| 免费看日本二区| 欧美不卡视频在线免费观看| 久久久成人免费电影| 午夜福利视频1000在线观看| 成人美女网站在线观看视频| 亚洲精品aⅴ在线观看| 久久久久久久久久黄片| 精品欧美国产一区二区三| 午夜免费观看性视频| 婷婷色麻豆天堂久久| 一个人免费在线观看电影| 亚洲av电影不卡..在线观看| 久久久久久久久久久免费av| 日韩大片免费观看网站| 国产精品熟女久久久久浪| 秋霞伦理黄片| 免费看a级黄色片| 国产黄色小视频在线观看| 亚洲av电影不卡..在线观看| 欧美日韩综合久久久久久| 国产成人精品婷婷| 韩国av在线不卡| 一个人看视频在线观看www免费| 成人国产麻豆网| 禁无遮挡网站| 精品久久久久久久人妻蜜臀av| 久久99精品国语久久久| 卡戴珊不雅视频在线播放| 国产色爽女视频免费观看| 国产女主播在线喷水免费视频网站 | 亚洲欧美一区二区三区国产| 欧美变态另类bdsm刘玥| 婷婷色综合大香蕉| 亚洲国产最新在线播放| 最近中文字幕2019免费版| 国产永久视频网站| 国产国拍精品亚洲av在线观看| 精品一区二区三区视频在线| 国产大屁股一区二区在线视频| 女人久久www免费人成看片| 亚洲怡红院男人天堂| 丝袜美腿在线中文| 亚洲国产成人一精品久久久| 免费av不卡在线播放| 国产精品福利在线免费观看| 69人妻影院| 成年av动漫网址| 精品酒店卫生间| 丰满人妻一区二区三区视频av| 人妻少妇偷人精品九色| 伦精品一区二区三区| 国产不卡一卡二| 国产亚洲av嫩草精品影院| 午夜视频国产福利| 26uuu在线亚洲综合色| 午夜精品在线福利| 十八禁国产超污无遮挡网站| 夜夜看夜夜爽夜夜摸| 色综合亚洲欧美另类图片| 日本一二三区视频观看| 亚洲在线观看片| 精品一区二区三区人妻视频| 日韩av免费高清视频| 国产精品熟女久久久久浪| 中国美白少妇内射xxxbb| 免费播放大片免费观看视频在线观看| 美女cb高潮喷水在线观看| 黄色欧美视频在线观看| 日产精品乱码卡一卡2卡三| 一区二区三区四区激情视频| av线在线观看网站| 免费观看av网站的网址| 久久久精品欧美日韩精品| 女人被狂操c到高潮| 身体一侧抽搐| 老司机影院成人| 精品熟女少妇av免费看| 亚洲四区av| 国产精品久久久久久久电影| 国产亚洲精品久久久com| 熟女电影av网| 日韩欧美精品免费久久| 国产有黄有色有爽视频| 午夜福利视频1000在线观看| 身体一侧抽搐| 精品人妻熟女av久视频| 亚洲精品视频女| 水蜜桃什么品种好| 久久精品人妻少妇| 日韩一区二区三区影片| av播播在线观看一区| 99热这里只有是精品在线观看| 国产黄色免费在线视频| 成年免费大片在线观看| 汤姆久久久久久久影院中文字幕 | 亚洲成色77777| 久久97久久精品| 国产探花极品一区二区| 一个人免费在线观看电影| 国产伦一二天堂av在线观看| 最近中文字幕高清免费大全6| 国产成人一区二区在线| 久久精品久久精品一区二区三区| 色综合色国产| av国产久精品久网站免费入址| 亚洲成色77777| 欧美成人a在线观看| 一个人看的www免费观看视频| 男女视频在线观看网站免费| 国产精品一区二区在线观看99 | 伦精品一区二区三区| 一级毛片久久久久久久久女| 久久鲁丝午夜福利片| 国产av国产精品国产| 国产精品一区二区三区四区免费观看| 91精品伊人久久大香线蕉| 伦理电影大哥的女人| 91久久精品国产一区二区成人| 少妇高潮的动态图| 日韩欧美精品免费久久| 三级国产精品欧美在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲精品视频女| 麻豆乱淫一区二区| 久久国产乱子免费精品| 欧美高清性xxxxhd video| 又大又黄又爽视频免费| videos熟女内射| 99热网站在线观看| 欧美一级a爱片免费观看看| 中文精品一卡2卡3卡4更新| 日韩电影二区| 精品一区在线观看国产| 欧美人与善性xxx| 亚洲av一区综合| 国产成人freesex在线| 国产高清有码在线观看视频| 老司机影院毛片| 男插女下体视频免费在线播放| 你懂的网址亚洲精品在线观看| 久久国内精品自在自线图片| 国产精品综合久久久久久久免费| 午夜激情福利司机影院| 美女大奶头视频| 久久99热6这里只有精品| 高清午夜精品一区二区三区| 一区二区三区乱码不卡18| av.在线天堂| 啦啦啦中文免费视频观看日本| 久久精品熟女亚洲av麻豆精品 | 别揉我奶头 嗯啊视频| 亚洲av男天堂| 免费人成在线观看视频色| 欧美另类一区| 国产精品99久久久久久久久| 国产成人a∨麻豆精品| 国产老妇女一区| 国产精品熟女久久久久浪| 国产爱豆传媒在线观看| 干丝袜人妻中文字幕| av国产免费在线观看| 丝瓜视频免费看黄片| 七月丁香在线播放| 秋霞伦理黄片| 综合色丁香网| 欧美最新免费一区二区三区| 男人和女人高潮做爰伦理| 欧美精品一区二区大全| 建设人人有责人人尽责人人享有的 | 国精品久久久久久国模美| 我的老师免费观看完整版| 黄色一级大片看看| 亚洲自偷自拍三级| 内地一区二区视频在线| 日韩欧美精品v在线| 国产人妻一区二区三区在| 黄色日韩在线| 国产黄色小视频在线观看| 美女被艹到高潮喷水动态| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品中文字幕在线视频 | 成人美女网站在线观看视频| 亚洲国产av新网站| 中文字幕制服av| av播播在线观看一区| 51国产日韩欧美| 日本-黄色视频高清免费观看| 69人妻影院| 色5月婷婷丁香| 精品久久国产蜜桃| 青春草视频在线免费观看| 男人舔女人下体高潮全视频| 久久久久性生活片| 身体一侧抽搐| 日韩 亚洲 欧美在线| 亚洲精品,欧美精品| 国产v大片淫在线免费观看| www.色视频.com| 中文字幕久久专区| 美女xxoo啪啪120秒动态图| 国产精品一区二区三区四区免费观看| 毛片女人毛片| 午夜福利在线在线| 国产久久久一区二区三区| 亚洲av免费在线观看| 99久久精品一区二区三区| 午夜亚洲福利在线播放| 在线天堂最新版资源| 成人av在线播放网站| 91aial.com中文字幕在线观看| 久久精品国产亚洲av天美| av又黄又爽大尺度在线免费看| 国产又色又爽无遮挡免| 一个人观看的视频www高清免费观看| 亚洲成人av在线免费| 亚洲国产精品成人久久小说| 91av网一区二区| 日韩欧美国产在线观看| 国产视频首页在线观看| 亚洲成色77777| 麻豆成人午夜福利视频| 美女大奶头视频| 在线观看人妻少妇| 淫秽高清视频在线观看| 国产精品一区二区三区四区免费观看| 国产精品嫩草影院av在线观看| 免费看a级黄色片| 别揉我奶头 嗯啊视频| 九色成人免费人妻av| 内射极品少妇av片p| 久久久久九九精品影院| 久久久亚洲精品成人影院| 18禁动态无遮挡网站| 亚洲18禁久久av| 亚洲性久久影院| 精品久久久久久久末码| 中文天堂在线官网| 亚洲精品国产av蜜桃| 国产精品人妻久久久影院| 国产伦理片在线播放av一区| 性色avwww在线观看| 91午夜精品亚洲一区二区三区| 午夜福利成人在线免费观看| 国产国拍精品亚洲av在线观看| 99re6热这里在线精品视频| 最近手机中文字幕大全| 国产极品天堂在线| 99久久人妻综合| 精品久久国产蜜桃| 女人久久www免费人成看片| 亚洲国产色片| 看非洲黑人一级黄片| 亚洲三级黄色毛片| 日本爱情动作片www.在线观看| 日本-黄色视频高清免费观看| 韩国高清视频一区二区三区| 精品久久久久久久久亚洲| 亚洲人成网站在线播| 日韩制服骚丝袜av| 亚洲av在线观看美女高潮| 国产午夜精品久久久久久一区二区三区| 青春草国产在线视频| 日本与韩国留学比较| 国产一级毛片在线| 26uuu在线亚洲综合色| 别揉我奶头 嗯啊视频| 免费电影在线观看免费观看| 久久久a久久爽久久v久久| 日韩三级伦理在线观看| 精品久久久噜噜| 欧美精品国产亚洲| 国产欧美日韩精品一区二区| 少妇的逼好多水| 成年女人在线观看亚洲视频 | 国产亚洲最大av| 国产伦理片在线播放av一区| 一级毛片久久久久久久久女| 亚洲最大成人中文| 国产成人a∨麻豆精品| 欧美区成人在线视频| 亚洲四区av| 国产爱豆传媒在线观看| 在线免费观看不下载黄p国产| 久久久午夜欧美精品| 啦啦啦啦在线视频资源| 国产精品无大码| 久久草成人影院| 26uuu在线亚洲综合色| 高清在线视频一区二区三区| 久久久久久久久久黄片| 免费人成在线观看视频色| 激情五月婷婷亚洲| av在线亚洲专区| 69av精品久久久久久| av福利片在线观看| 国语对白做爰xxxⅹ性视频网站| 一个人看的www免费观看视频| 嫩草影院入口| 最近中文字幕2019免费版| 国内精品宾馆在线| 午夜福利网站1000一区二区三区| 国产人妻一区二区三区在| 日韩精品青青久久久久久| 亚洲一区高清亚洲精品| av国产免费在线观看| 直男gayav资源| 高清日韩中文字幕在线| 国产大屁股一区二区在线视频| 又粗又硬又长又爽又黄的视频| 18+在线观看网站| 国产黄色小视频在线观看| h日本视频在线播放| 男插女下体视频免费在线播放| 中文天堂在线官网| 天堂俺去俺来也www色官网 | 一区二区三区四区激情视频| 精品国产一区二区三区久久久樱花 | 国产一区二区三区av在线| 成人综合一区亚洲| 男女视频在线观看网站免费| 久久久午夜欧美精品| 国产精品麻豆人妻色哟哟久久 | 日本-黄色视频高清免费观看| 欧美成人精品欧美一级黄| 97在线视频观看| 插逼视频在线观看| 99久国产av精品| 欧美另类一区| 日本免费在线观看一区| 久久久久久久大尺度免费视频| 男女国产视频网站| 精品熟女少妇av免费看| 国产91av在线免费观看| 舔av片在线| 亚洲av日韩在线播放| 久久久久久久国产电影| 久久久精品欧美日韩精品| 国产免费福利视频在线观看| 亚洲美女搞黄在线观看| 欧美一区二区亚洲| 纵有疾风起免费观看全集完整版 | 国产一区二区三区av在线| 一级毛片aaaaaa免费看小| 能在线免费看毛片的网站| 免费观看a级毛片全部| 中文字幕av在线有码专区| 毛片女人毛片| 久久这里只有精品中国| 亚洲,欧美,日韩|