• <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 improvement of spray characteristics of aero-engines using gliding arc plasma

    2023-03-15 00:54:26LeiZHANG張磊DengchengZHANG張登成JinluYU于錦祿BingbingZHAO趙兵兵XinyuQU屈新宇YiCHEN陳一andWeidaCHENG程偉達(dá)
    Plasma Science and Technology 2023年3期
    關(guān)鍵詞:張磊

    Lei ZHANG (張磊), Dengcheng ZHANG (張登成), Jinlu YU (于錦祿),2,3,*,Bingbing ZHAO (趙兵兵), Xinyu QU (屈新宇), Yi CHEN (陳一) and Weida CHENG (程偉達(dá))

    1 Aeronautics Engineering College,Air Force Engineering University,Xi’an 710038,People’s Republic of China

    2 School of Power and Energy,Northwestern Polytechnical University,Xi’an 710003,People’s Republic of China

    3 AECC Sichuan Gas Turbine Establishment, Chengdu 610599, People’s Republic of China

    Abstract A gliding arc plasma fuel atomization actuator suitable for aeroengines was designed, and a gliding arc plasma fuel spray experimental platform was built to address the fuel atomization problem in aeroengine combustion chambers.The spray characteristics for different airflows,fuel flows,and discharge voltages were analyzed using laser particle size analysis.The research shows that the fuel atomization effect is improved from the increased airflow.The decreased fuel flow not only reduces the injection pressure of the fuel but also changes the discharge mode of the gliding arc, which affects reductions in the discharge power and inhibits fuel atomization.Gliding arc discharges accelerate the breaking, atomization, and evaporation of fuel droplets while reducing the particle size,which increases the proportion of small droplets.Compared with the working conditions of plasma-assisted atomization without the gliding arc,the D0.5,D0.9,and average particle size of the fuel droplets are reduced by 4.7%, 6.5%, and 4.1%, respectively,when the modulation voltage of the gliding arc power supply is 200 V.

    Keywords: gliding arc discharge, spray characteristics, droplet size distribution, aeroengine

    1.Introduction

    Increasingly strict emission requirements have put focus on pollutants from aeroengines.Therefore, it is necessary to reduce aeroengine emissions while improving their performance.The fuel atomization quality in the combustion chamber of an aeroengine directly impacts the engine performance and pollutant emissions [1].However, at high altitudes and low temperatures, the viscous force and surface tension of the fuel increase[2,3],which creates larger ejected fuel particles.In addition,with increased flight altitude,the air density is reduced; thus, the airflow of the engine decreases(and correspondingly the fuel supply), the fuel pressure in front of the nozzle decreases, the fuel atomization quality worsens,and the fuel particle size increases.This reduces the effective contact area between the fuel and air, which is not conducive to its atomization and evaporation and inhibits combustion reactions [4, 5].This directly leads to increased pollutant emissions from the aeroengine.

    There has been significant research on the atomization cone angle and enhanced spray diffusion to improve the performance of aeroengine sprays while reducing emissions.Wu et al [6] studied the spray characteristics of aviation kerosene (RP-3) under flash and non-flash boiling conditions and found that flash boiling promotes the effective diffusion of fog droplets to reduce their size.In numerical simulations,Boretti et al[7]studied the spray characteristics of air-assisted injectors in direct injection spark ignition(DISI)engines,and Mohammad et al [8] analyzed the influence of cavitation phenomena on spray characteristics under different injection pressures.These studies showed that an increased injection pressure causes a greater spray penetration length and decreased Sauter diameter.Although the fuel spray effect can be improved through a better combustor aerodynamic structure or by increasing the fuel temperature or fuel injection pressure,these conventional methods cause structural changes to the combustion chamber, which requires significant additional research.Therefore, new atomization methods are needed to improve the atomization quality of fuel in the given combustion chamber.

    Plasma-enhanced combustion technology is a new approach that utilizes the thermal, kinetic, and transport effects [9, 10] generated in the plasma discharge process to influence the combustion chemical reactions [11, 12].The excited state particles and high-energy electrons generated in the plasma discharge process collide with fuel molecules,which shows significant advantages in improving the fuel atomization quality [13].A gliding arc is a typical representative of non-equilibrium plasmas, and it has the advantages of both low- and high-temperature plasmas [14, 15].Because the gliding arc has a large fuel contact area, high active particle density, high active particle excitation efficiency, and uniform and stable discharge [16], it has unique advantages in promoting fuel atomization.In addition, the gliding arc has a simple electrode structure that easily integrates with the combustion chamber structure, which gives it unique advantages in practical applications [17].

    To improve the fuel atomization quality using plasmas,Matveev et al [18] designed a plasma fuel nozzle to enhance combustion using transient plasma discharge.The group also researched the effects of nanosecond pulse discharge plasma on the extinguishing limit of a depleted C3H8/air premix and kerosene spray fuel [19, 20].Gomez et al [21] designed a multipoint lean fuel direct injection(LDI)plasma injector and tested it using a simulator with a fan-shaped combustion chamber.The study showed that plasma discharge could reduce the equivalent ratio of depleted fuel from 1% to 24%.Huang et al [22] proposed a gliding arc plasma airblast fuel injector and studied the spray characteristics of aviation kerosene under different working conditions.Research has shown that gliding arc plasma discharge gives more uniform fuel atomization, smaller droplet sizes, and larger spray cone angles.

    Current research to improve fuel atomization quality via plasma is in the exploratory stage.It is difficult to apply plasma to aeroengines because of its complex structure and poor working conditions.In this study, a gliding arc plasma fuel atomization exciter was designed, and a fuel spray experimental platform was established to address the problem of fuel atomization in the aeroengine combustor.The macroscopic structure of the spray field was analyzed using laser particle size analysis.Finally,the effects of gliding arc plasma discharge, airflow, and fuel flow on fuel droplet size distribution and fuel particle size were studied.

    2.Experimental system and devices

    2.1.Experimental system

    Figure 1 shows a schematic of the gliding arc plasma fuel spray experimental system.The device consists of a power supply, gas supply system, fuel supply system, measurement system, gliding arc plasma fuel atomizing exciter, spray experimental section, and other components.The spray section is composed of an optical glass chamber with light transmittance on all sides to ensure the tightness and light transmittance of the spray field.The entire experimental process was performed under no-light conditions to ensure the accuracy of the particle size measurements.

    Figure 1.Schematic of the gliding arc experimental system.

    The gliding arc plasma is excited using a plasma power supply system consisting of a low-temperature millisecond pulse power supply (CTP-2000K, Nanjing Suman Plasma Technology Co.) and a voltage regulator.The power supply changes the output voltage of the voltage regulator into a high-voltage sinusoidal AC signal.The input of the power supply is 220 V, the modulation voltage is 0-250 V, the output voltage peak of the power supply is 0-30 kV, and the power supply frequency is 1-100 kHz with a power of 500 W.

    The air supply system is composed of a screw air compressor (BK132-8, displacement 24 m3min-1, nominal exhaust pressure 0.8 MPa), pressure storage tank (volume of 4 m3), cold and dry machine (LY-D200ah, 28.5 m3min-1),inner cone flowmeter (DYNZ16-8001E12), intelligent flow integrator, and other components.Thus, the system can accurately measure the volumetric airflow.

    2.2.Measurement system

    The measurement system consists of a Winner318 spray laser particle size analyzer and a CNC 3D coordinate displacement mechanism.The Winner318 spray laser particle size analyzer is installed on the cantilever claw end of the 3D coordinate displacement mechanism, and its movement trajectory is controlled with a software program.This completes the particle size measurements for the entire spray area.

    The Winner318 spray laser particle size analyzer consists of a host with an optical system and a signal acquisition and processing system.The optical system is composed of a Nd:YAG laser with a power of 40 mW and wavelength of 532 nm, along with a collimating lens.The host is connected to a computer as a display to output the experimental data.The laser beam is filtered, expanded, filtered, collimated, and then irradiated onto the test area.The particles to be measured in the test area scatter the incident laser irradiation.The intensity and spatial distribution of the scattered spectrum are related to the size and distribution of the measured particle group.After converging again through the lens, the scattered light is received by a photoelectric array detector located on the rear focal plane of the lens.After being converted into an electrical signal, it is amplified, digitized, and sent to the computer through the communication port.After inversion and data processing, the volume distribution curve, characteristic particle size, and Sauter average particle size of the measured particle group are output.The adjustable measuring range of the equipment is 15-711 μm,the sampling frequency is 3 Hz, and the device has a 3D fine-tuning detection unit.Multiple automatic scattering correction functions are designed based on the spray concentration, which gives accuracy and repeatability errors of less than 3%.The instrument can capture continuous dynamic testing of the samples.

    In the spray experiments, a high-speed CCD camera(Mini UX50), high-voltage probe (Tektronix P6015A),current probe (Tektronix TCP0030), and oscilloscope (Tektronix DPO4104B) were used to collect the arc motion characteristics and voltage and current signals from the gliding arc discharge.The frame rate of the high-speed CCD camera was set to 10 000 fps and the exposure time was 100 μs,which quickly captures the motion of the gliding arc.

    2.3.Gliding arc plasma fuel atomizing excitation

    A gliding arc plasma fuel atomization exciter was designed to promote the application of gliding arc plasma-enhanced combustion technology in aeroengines and improve fuel atomization quality using gliding arc plasma[15].The design combines the gliding arc plasma actuator with the dome of the aeroengine combustion chamber and uses the rotating airflow at the combustion chamber dome to drive the arc movement.A schematic is shown in figure 2.The three-dimensional rotating gliding arc exciter is composed primarily of a cathode fuel nozzle,anode Venturi tube,cyclone,and bell mouth.The minimum breakdown distance between the cathode fuel nozzle and anode Venturi tube is 10.9 mm, and the outlet angle of the anode venturi tube is 37°.The anode Venturi tube is connected to the high-voltage alternating current through a high-voltage cable, and the cathode fuel nozzle and combustion chamber are grounded together.The anode Venturi tube and cathode fuel nozzle generate a plasma arc under the excitation of a strong electric field to form a three-dimensional rotating gliding arc plasma discharge under the action of rotating airflow in the cyclone.The fuel jet from the cathode fuel nozzle passes through the gliding arc discharge area, further promoting fuel atomization.

    Figure 2.Schematic of the gliding arc plasma fuel atomization actuator.

    2.4.Experimental conditions

    The experimental conditions are shown in table 1.Aviation kerosene (RP-3) was selected as the fuel during the experiment.Three-way valves were used to connect the air source system,argon system,and spray laboratory module to ensure that the plasma arc did not ignite the fuel after atomization in the experimental section.The mass flow controller was used to detect the makeup amount of argon, and the vortex flowmeter was used to detect the total airflow into the experimental section.As long as the argon flow is 1/3 of the total gas flow, the proportion of argon replenishment can be guaranteed.The added argon enters the experimental section together with the airflow of the gas supply system to protect and prevent combustion.All experiments were measured under the conditions of stable air and fuel flows.

    Table 1.Experimental conditions.

    2.5.Operating characteristics of gliding arc plasma fuel atomizing exciter

    A high-speed camera recorded the sliding of the plasma arc along the Venturi tube under atmospheric pressure, as shown in figure 3.The gliding arc plasma discharge was driven by a rotating airflow of 20 m3h-1.The movement of the gliding arc plasma was captured under breakdown, tensile sliding,and extinction.Among them,figures 3(a)and(l)show the arc breakdown moment, and figures 3(b)-(k) are images of the arc movement and extension process.The arc first breaks down at the minimum distance between the two electrodes under the excitation of a high-voltage electric field.Driven by the airflow, the arc stretches and twists while its root slides between the two poles.When the arc develops to a certain length, the voltage between the two electrodes cannot maintain the arc development and is extinguished.When the arc is extinguished, the active components in the arc channel are blown downstream by the rotating airflow.At this time, the condition of the forming arc channel downstream is better than the arc breakdown condition at the minimum spacing.Therefore, the arc does not necessarily continue to break down at the minimum distance, and the arc breakdown location becomes random.According to figures 3(a) and (l),when the arc enters the subsequent development cycle, the breakdown position of the arc differs.

    The voltage and current signals in the gliding arc discharge process, as collected synchronously with the arc movement from figure 3,are shown in figure 4.Time t0is the arc breakdown moment at which the voltage has a pulse peak of approximately 5 kV and the instantaneous current reaches about 15 A.The arc shape corresponds to figure 3(a), where the arc has higher brightness and a shorter length.After arc breakdown, the voltage and current develop a sinusoidal pattern, and the voltage increases with arc stretching.Once the arc develops to a certain length (arc length shown in figures 3(h)-(k)),the voltage and current have repeated pulse peaks(as shown in the voltage and current waveforms within the period from t2to t3in figure 4).This is because the arc stretching increases the channel resistance,and the stability of the current is not sufficient to maintain its development,which causes it to dissipate.However,the conditions to form the arc channel are optimal, and a peak voltage and current pulse appear to maintain arc stretching.When the power supply is insufficient to maintain continuous development of the arc, it extinguishes, and a new arc channel is formed.As illustrated by t4in figures 3(l) and 4, a new arc channel is formed, and the arc movement enters the next cycle.

    Figure 3.Arc sliding process over time.

    Figure 4.Gliding arc discharge voltage and current.

    Figure 5.Gliding arc discharge power.

    Figure 6.Distribution structure of the spray field.(a)Cloud diagram of the spray field distribution, (b) particle size distribution on the radial section.

    Figure 5 shows variations in the discharge power for the gliding arc discharge process.The analysis shows that the gliding arc discharge power first increases and then decreases with the inlet airflow.However, under a 120 V modulation voltage, the gliding arc discharge power increases with the inlet airflow, and the trend of an initial increase followed by a decrease appears in some areas.However, within the airflow range involved in the spray experiments (15, 20, and 25 m3h-1), an increased airflow also increases the discharge power of the gliding arc because a greater gas flow alters the discharge mode of the gliding arc [23].This enhances both the local turbulence degree and the heat and mass transfer between the plasma arc and the environment.Therefore, more energy is needed to develop the plasma arc further.With a greater modulation voltage from the power supply, the discharge power increases significantly because the energy input from the power supply to the two electrodes increases with the modulation voltage, which enhances the discharge power.

    3.Results and discussion

    3.1.Analysis of macroscopic characteristics of the spray field

    The distribution structure of the fuel spray field is shown in figure 6 for a 0.45 g s-1fuel supply flow and no swirling airassisted atomization.Figure 6(a) indicates that the fuel spray field is cone-shaped, and the overall particle size distribution ranges from 30 to 130 μm.The fuel particle size is larger within the spray cone angle and has a concentrated distribution.In figure 6(a),the blue area represents fuel particle sizes less than 40 μm and is defined as outside the fuel spray cone angle as marked with the black dotted line,which is only used for the macroscopic description of particle size distribution.Figure 6(b) shows that the particle size distribution curves along different axial sections have symmetric concave distributions.The local minimum of the fuel particle size is at the center point, and the curve becomes steeper toward the fuel nozzle.This result agrees with the distribution law of the fuel particle size in the spray field structure.

    The influence of the fuel supply flow change on the fuel particle size distribution is considered when there is no airassisted atomization,as shown in figure 7.The Wfin the figure represents the fuel supply flow, Δ is the fuel spray non-uniformity coefficient, given by Δ = (SMDmax-SMDave) /SMDave[17], and smaller values indicate a more uniform fuel spray.With a smaller fuel flow, the fuel particle size increases in the mainstream spray area (red area in figure 6(a)), mainly because the fuel supply flow decreases, and the fuel supply pressure in front of the nozzle decreases.Thus, the fuel supply pressure decreases, the kinetic energy of the injected fuel decreases, the friction and shear between the fuel and air after leaving the nozzle decreases,and the fuel particle size increases.In the range of Z=-10 to -30 mm, the fuel spray nonuniformity coefficient increases for smaller fuel supply flows,which is primarily the result of the increased fuel particle size.In the Z=-40 to-60 mm range, the fuel supply flow decreases, but the spray nonuniformity coefficient also decreases.Combined with the fuel particle size distribution curve,this is mainly due to smaller fluctuations in the particle size, even though the particle size is still large.

    The influence of the airflow on the fuel particle size distribution when the fuel supply flow is 0.57 g s-1is shown in figure 8.The Wain the figure represents the airflow.The increased gas flow reduces the particle size of the fuel droplets and improves the uniformity of the spray field, especially in the spray cone angle.This is because the fuel ejected from the nozzle forms a mixture of fuel and gas under the action of air.The air in the mixture expands rapidly under the disturbance of high-speed airflow [6], which accelerates the atomization of the fuel droplets.Figures 7(a)and(b)show that after air is injected, the fuel particle size curve becomes convex instead of concave within the cone angle of the fuel spray.The absolute maximum fuel particle size appears at the center point, X=0 mm.

    Figure 7.Effects of fuel flow on the fuel particle size distribution.(a)Z=-10 mm,(b)Z=-20 mm,(c)Z=-30 mm,(d)Z=-40 mm,(e) Z=-50 mm, (f) Z=-60 mm.

    After air is introduced,the local minimum particle size of fuel is formed at X=-20 mm and 20 mm.This is because an annular vortex is formed downstream after the air passes through the swirler; a schematic of the ‘vortex structure’ is shown in figure 8(g)[17].The inner region of the red ellipsis of the vortex center in figure 8(g) is defined as the vortex core, where the left side of the vortex core is defined as the outer side of the vortex,and the right side of the vortex core is defined as the inner side of the vortex.The inner side of the vortex is a low-pressure area where the airflow velocity is small.The fuel droplets on the outer side of the vortex are brought to the inner side of the vortex; thus, the absolute maximum fuel particle size appears on the inner side of the vortex.However, at the vortex core, the streamlines are relatively sparse, and the fuel molecules are affected by the airflow trajectory,which forms the local minimum of the fuel particle size.

    Figure 8.Effect of airflow on the fuel particle size distribution.(a) Z=-10 mm, (b) Z=-20 mm, (c) Z=-30 mm, (d) Z=-40 mm,(e) Z=-50 mm, (f) Z=-60 mm, (g) schematic of vortex [17].

    Figure 9.Effect of the modulation voltage on the fuel particle size distribution.(a) Z=-10 mm, (b) Z=-20 mm, (c) Z=-30 mm,(d) Z=-40 mm, (e) Z=-50 mm, (f) Z=-60 mm.

    The influence of the fuel supply flow on the fuel particle size distribution is shown in figure 9 after coupling between the air and gliding arc discharge(airflow is 20 m3h-1and the power supply modulation voltage is 200 V).Within the cone angle of the fuel spray,the fuel supply decreases,the particle size of the fuel droplets increases, and the uniformity of the spray field deteriorates.This is because when the temperature and pressure fields at the outlet of the fuel nozzle are relatively constant, the fuel supply reduces, the pressure difference before the fuel nozzle becomes smaller, there is less kinetic energy for the injected fuel, and the aerodynamic effect between the fuel and air is relatively weak.In addition, as the fuel supply flow is reduced,the proportion of the breakdown mode in the discharge mode of the gliding arc decreases, which leads to a lower discharge power for the gliding arc[23].This weakens the ability of the gliding arc to process the fuel, deteriorating the fuel atomization quality.Most droplets outside the cone angle of the fuel spray escape in the mainstream spray area, and its change law differs from the above law.Figures 8(b)-(f) show that the fuel spray inhomogeneity coefficient does not increase with a lower fuel supply flow but achieves a global minimum of 0.57 g s-1.However, the particle size distribution curve indicates that the particle size global minimum is at 0.76 g s-1.This maximizes the gliding arc discharge mode of the breakdown mode when the fuel supply flow rate reaches the absolute maximum, which increases the local turbulence (gliding arc discharge produces local disturbances to the discharge area)[24].The aerodynamic factors increase the local nonuniformity of the spray.

    3.2.Droplet size-frequency distribution

    The particle size of the fuel spray is an important factor that affects the fuel evaporation rate.The droplet size-frequency distribution is used to analyze the particle size of the fuel spray and quantitatively describe the influence of the gliding arc plasma discharge on the particle size distribution.The droplet size-frequency distribution is given byφB/Δd[25],whereφBis the volume fraction in each size interval andΔdis the size interval.

    Applying the above formula provides the droplet sizefrequency distribution shown in figure 10 when the fuel supply flow is 0.57 g s-1,and the inlet airflows are 0,20,and 25 m3h-1.The ‘w/o plasma’ labels in the figure indicate plasma excitation without the gliding arc.The droplet size frequency presents a biased normal distribution.When the inlet airflow is 0 m3h-1, most droplets are distributed between 90 and 130 μm,and the 110-120 μm range accounts for the largest proportion at 4.41%μm-1.After swirling air is introduced (20 m3h-1), the fuel injected from the nozzle forms a fuel-gas mixture, and the air expands rapidly under the disturbance of the high-speed airflow, which accelerates the fuel droplet breakage.The inlet airflow increases,and the spray field of the aerodynamic force is enhanced.The shear force between the fuel droplets and air friction increases,accelerating fuel droplet fragmentation.Therefore, the proportion of the droplet frequency over 120-130 μm decreases,and the absolute maximum droplet frequency proportion tends to smaller particle sizes.When the inlet airflow increases to 25 m3h-1, the absolute maximum droplet sizefrequency ratio is reduced to 3.03% μm-1with a range of 100-110 μm.

    The influence of the gliding arc plasma on the droplet size-frequency distribution is shown in figure 11, where U0represents the modulation voltage.After the gliding arc plasma-assisted atomization is applied, the frequency proportion of the larger fuel particle sizes is reduced, and the peak value of the droplet size-frequency distribution curve moves toward smaller particle sizes.The working condition with an airflow of 20 m3h-1is taken as an example.Under non-gliding arc plasma-assisted atomization, the fuel particle size is concentrated in the 90-130 μm range,and the absolute maximum droplet size-frequency ratio appears in the 110-120 μm range.Droplets in the 120-130 μm range disappear after the gliding arc plasma-assisted atomization,and the absolute maximum droplet size-frequency ratio shifts to smaller particle sizes.The thermal effects generated during gliding arc plasma discharge instantly heat the discharge area[9, 10] and promote the atomization and evaporation of fuel droplets.At the same time,the gliding arc discharge enhances the eddy current intensity and air velocity downstream of the cyclone[24,26].This increases the interaction force between the fuel droplets and air, and promotes fuel atomization.Therefore,when the gliding arc plasma modulation voltage is 160 V,the fuel particle size is concentrated in the 70-120 μm range,the absolute maximum fuel droplet size-frequency ratio shifts to the 100-110 μm range, and the absolute maximum droplet size-frequency ratio decreases to 3.81% μm-1.

    Figure 10.Effect of airflow on the droplet size-frequency distribution.

    3.3.Cumulative droplet size distribution

    The cumulative distribution curve of the droplet sizes is an important parameter to describe the droplet size distribution.Its expression is [25]

    whereQrepresents the total volume fraction contained in droplets with sizes smaller than D,VDrepresents the total volume of droplets with sizes smaller than D, andVrepresents the volume of any droplet.

    Figure 12 shows the cumulative distribution curves of the droplet sizes at different airflows with a fuel supply flow rate of 0.57 g s-1and power supply modulation voltage of 200 V.The airflow affects the slope and position of the cumulative droplet size distribution curve.A larger airflow gives a greater slope for the cumulative distribution curve,which produces a higher position.

    Figure 11.Effect of gliding arc on the droplet size-frequency distribution.(a) Air flow is 20 m3 h-1, (b) air flow is 25 m3 h-1.

    Figure 12.Effects of airflow on the cumulative distribution of the droplet sizes.

    A larger curve slope gives a smaller absolute maximum value particle size.The higher the curve is, the larger the proportion of small fuel droplets.The analysis of section 3.1 shows that the aerodynamic force plays an important role in the fuel atomization process, and a larger airflow provides stronger interaction forces between the fuel droplets and air.The junction analysis of section 2.5 indicates that a larger airflow gives a higher gliding arc discharge power.Thus,the processing capacity of the gliding arc is greater,and the fuel atomization quality improves.

    To explore the gliding arc discharge effect on the cumulative droplet size distribution, we calculate the cumulative droplet size distribution curve under different working conditions in figure 13.The red chain line represents no gliding arc plasma, the dotted blue line represents a modulation voltage of 160 V,and the solid green line represents a power modulation voltage of 200 V.With increased plasma modulation voltage, the cumulative droplet size distribution curve shifts upward,and its slope increases.That is to say,the absolute maximum particle of fuel particle size becomes smaller as the supply modulation voltage increases.Figure 13(c)indicates that when the airflow is 20 m3h-1and the fuel supply flow is 0.57 g s-1,the cumulative droplet size distribution curve reaches its peak at 122.59 μm without gliding arc plasma excitation.When the gliding arc plasma modulation voltage is 200 V, the position where the cumulative droplet size distribution curve reaches a peak shifts to 111.66 μm,and the absolute maximum droplet size is reduced by 10.93 μm, proving that the gliding arc plasma excitation can effectively reduce the fuel particle size.

    D0.5represents the droplet size when Q is 50%, which characterizes the median diameter; a smaller value gives a higher spray uniformity.Its value decreases with the modulation voltage, indicating that the droplet’s overall atomization effect is improved.The D0.5value is 109.01 μm when there is no gliding arc plasma excitation and 100.951 μm when the gliding arc plasma modulation voltage is 200 V.This indicates that the atomization effect of the spray field improves after applying the gliding arc plasma.When the gliding arc discharges, it instantly increases the temperature of the discharge area to thousands of Kelvin, which directly impacts the viscous force between the fuel molecules and surface tension.Thus, the fuel droplets more easily atomize and evaporate.In addition, the gliding arc discharge process results in high local pressure in its vicinity,which causes local turbulence pulsation and accelerates the fuel droplet breakdown.Thus, the fuel atomization quality is improved.

    Compared with figures 13(b)-(e), the slope of the cumulative distribution curve of droplet size decreases, and the value of D0.5increases when the fuel supply flow decreases.This means that a smaller fuel flow gives a larger maximum fuel particle size and worsens the atomization uniformity.Without gliding arc plasma excitation, the value of D0.5is 103.76 μm when the flow rate is 0.76 g s-1, and D0.5is 111.49 μm when the flow rate decreases to 0.3 g s-1.When the outside conditions are the same as the inlet airflow phase,the temperature and pressure fields at the nozzle outlet are constant,the fuel supply flow decreases,the pressure drop in front of the fuel nozzle decreases,the kinetic energy of the injected fuel is small, and the interaction force between fuel and air decreases.The fuel supply quantity is small, the gliding arc discharge power decreases,and the gliding arc on the fuel atomization effect is abated.Thus, the fuel atomization quality deteriorates.

    Figure 13.Cumulative droplet size distribution for various working conditions.(a)Wa=15 m3 h-1,Wf=0.34 g s-1,(b)Wa=20 m3 h-1,Wf=0.76 g s-1, (c) Wa=20 m3 h-1, Wf=0.57 g s-1, (d) Wa=20 m3 h-1, Wf=0.45 g s-1, (e) Wa=20 m3 h-1, Wf=0.30 g s-1,(f) Wa=25 m3 h-1, Wf=0.57 g s-1.

    3.4.Average particle size and change of representative particle size

    The effects of gliding arc plasma discharge on the representative particle size (D0.5, D0.9) and average particle size (D32) of the fuel spray were studied for an airflow of 20 m3h-1.D0.9represents the droplet size when Q is 90%and characterizes the upper 10% of droplet sizes; the larger the value indicates, the more the number of large-sized droplets.D32is the average particle size,which represents the Sauter mean diameter(SMD)of the droplets.Figure 14(a) shows the variations in D0.5, representing the spray uniformity.Figure 14(b)shows the variations in average particle size, and figure 14(c) shows variations in D0.9.Applying the gliding arc plasma excitation decreases the trend of the average fuel spray particle size and representative particle size(D0.5and D0.9).This decreasing trend becomes more evident when the modulation voltage of the gliding arc plasma is increased.When the fuel flow rate is 0.76 g s-1, the average particle size of fuel droplets without the gliding arc plasma excitation is 93.55 μm.The average fuel droplet particle size is 89.69 μm when the gliding arc plasma modulation voltage is 200 V,which is a reduction of 4.1%.D0.5decreases from 103.76 to 98.91 μm, a reduction of 4.7%, and D0.9decreases from 113.21 to 105.81 μm, a reduction of 6.5%.These results show that the gliding arc discharge has significant advantages in decreasing the average particle size of the fuel spray, which decreases the absolute maximum fuel droplet size and improves the fuel spray uniformity.This is because the ‘a(chǎn)erodynamic effect’produced in the gliding arc discharge process increases the flow field turbulence and accelerates fuel droplet fragmentation.At the same time, the ‘temperature rise effect’ produced in the gliding arc discharge process promotes fuel droplet atomization and evaporation.

    Figure 14.Effect of gliding arc plasma discharge on the mean and representative particle sizes.(a) D0.5, (b) D32, (c) D0.9.

    The D0.5, D0.9, and average particle size of the fuel spray tend to increase for smaller fuel flow rates.When applying the gliding arc plasma,the particle size decreases significantly at 0.76 and 0.57 g s-1compared with 0.30 g s-1.For D0.9,a fuel supply of 0.30 g s-1reduces the particle size by 6.3%after applying the gliding arc plasma excitation, and fuel supplies of 0.76 and 0.57 g s-1give 9.1% and 6.5% reductions, respectively.For the average particle size, when the fuel supplies are 0.76, 0.57, and 0.30 g s-1, the reductions are 4.3%, 4.1%, and 3.6%, respectively.There are two reasons for this.On one hand, the fuel supply flow and front nozzle pressure decrease, and the preliminary atomization effect worsens.On the other hand, the decreased fuel supply reduces the proportion of the breakdown mode in the gliding arc discharge process, and the discharge energy in the breakdown mode is greater than that in the stable sliding mode[23].Thus,the fuel supply volume is reduced.The power of the gliding arc discharge process decreases, so the processing capacity decreases,and the fuel particle size increases.

    4.Conclusions

    A gliding arc plasma fuel atomization exciter was proposed to improve the fuel atomization quality in the dome of an aeroengine combustor.The structure of the spray field was analyzed using laser particle size analysis.The effects of the airflow,fuel supply, and plasma excitation on the fuel particle size distribution were studied.The main conclusions are as follows:

    (1) Under the excitation of a strong electric field, a 3D gliding arc discharge was formed between the fuel nozzle and Venturi,and the discharge power of the gliding arc increased with the modulation voltage.When the modulation voltage was 160 and 200 V, the gliding arc discharge power first increased and then decreased with the gas flow.

    (2) The temperature rise and pneumatic effect from the gliding arc discharge accelerated the crushing, atomization, and vaporization of fuel droplets, while the average and representative particle sizes of fuel droplets were reduced.Thus, the cumulative droplet size distribution curve slope increased, the small droplet size proportion increased, and the atomization effect improved significantly.Compared with the working condition of plasma-assisted atomization without a gliding arc, when the modulation voltage of the gliding arc power supply was 200 V, the D0.5, D0.9, and average particle sizes of the fuel droplets were reduced by 4.7%, 6.5%, and 4.1%, respectively.

    (3) The increased airflow not only changed the structure and enhanced the aerodynamic action of the spray field but also increased the power of the gliding arc discharge, which improved the fuel atomization effect.The fuel supply flow decreased and reduced the injection pressure in front of the nozzle, worsening the fuel atomization effect.At the same time,reduced fuel flow also affected the discharge mode of the gliding arc, reduced its discharge power, worsened the fuel atomization effect,and increased the particle size of fuel atomization.When the modulation voltage was 160 V, the D0.9,D0.5,and average particle sizes at a fuel supply flow of 0.30 g s-1increased by 3.6%,4.7%,and 4.9%,respectively,compared with a fuel supply flow of 0.76 g s-1.

    Acknowledgments

    This work was supported by National Natural Science Foundation of China (Nos.91741112 and 52276142).

    猜你喜歡
    張磊
    Spin transport characteristics modulated by the GeBi interlayer in Y3Fe5O12/GeBi/Pt heterostructures
    張磊治療反流性食管炎經(jīng)驗(yàn)
    Temperature and doping dependent flat-band superconductivity on the Lieb-lattice?
    風(fēng)雨中逆行的抗“疫”巾幗戰(zhàn)士——記呼吸科副主任張磊
    北極光(2020年1期)2020-07-24 09:04:06
    THE GLOBAL ATTRACTOR FOR A VISCOUS WEAKLY DISSIPATIVE GENERALIZED TWO-COMPONENT μ-HUNTER-SAXTON SYSTEM?
    “口”“ㄙ”偏旁混用趣談
    “好聲音”冠軍張磊:哦,我的田螺姑娘
    幸福(2016年6期)2016-12-01 03:07:57
    什么是四輪驅(qū)動(dòng)?
    車(chē)迷(2015年6期)2015-03-20 02:43:54
    配型
    張磊老師的大醫(yī)情懷和大家風(fēng)范
    69人妻影院| 国产黄色视频一区二区在线观看| 99久久人妻综合| 久久亚洲国产成人精品v| 一区二区三区高清视频在线| 最近最新中文字幕免费大全7| 一二三四中文在线观看免费高清| 女人被狂操c到高潮| 波多野结衣巨乳人妻| 最近中文字幕2019免费版| av免费观看日本| 我要看日韩黄色一级片| 国产视频内射| kizo精华| 午夜福利高清视频| 国产成人精品婷婷| 午夜爱爱视频在线播放| 看黄色毛片网站| 国产精品一区二区性色av| 亚洲国产精品sss在线观看| 狂野欧美激情性xxxx在线观看| 在线免费十八禁| 国产在线一区二区三区精| 人妻制服诱惑在线中文字幕| 免费观看a级毛片全部| 青春草亚洲视频在线观看| 乱系列少妇在线播放| 亚洲av成人精品一区久久| av福利片在线观看| 97超视频在线观看视频| 老司机影院成人| 观看免费一级毛片| 三级国产精品欧美在线观看| 免费不卡的大黄色大毛片视频在线观看 | 婷婷六月久久综合丁香| 男女视频在线观看网站免费| 天堂影院成人在线观看| 寂寞人妻少妇视频99o| 亚洲综合色惰| 国语对白做爰xxxⅹ性视频网站| 九九爱精品视频在线观看| www.av在线官网国产| 久久国产乱子免费精品| 国产免费视频播放在线视频 | 久久久久久伊人网av| 国模一区二区三区四区视频| 久久精品国产鲁丝片午夜精品| 99久久九九国产精品国产免费| 黄片wwwwww| 国产高清有码在线观看视频| 国产一级毛片在线| 欧美日韩综合久久久久久| 搡女人真爽免费视频火全软件| 一二三四中文在线观看免费高清| 欧美成人午夜免费资源| 成年版毛片免费区| 爱豆传媒免费全集在线观看| 成人性生交大片免费视频hd| 尾随美女入室| 在线观看一区二区三区| 插阴视频在线观看视频| 久久精品熟女亚洲av麻豆精品 | 18禁动态无遮挡网站| 美女被艹到高潮喷水动态| 国产女主播在线喷水免费视频网站 | 久久草成人影院| 久久精品国产亚洲av涩爱| 亚洲乱码一区二区免费版| 青春草亚洲视频在线观看| 韩国高清视频一区二区三区| 看免费成人av毛片| 日日撸夜夜添| 国产精品美女特级片免费视频播放器| 亚洲无线观看免费| 国产午夜精品久久久久久一区二区三区| a级毛片免费高清观看在线播放| 精品欧美国产一区二区三| 成人一区二区视频在线观看| 久久99精品国语久久久| 亚洲国产精品成人综合色| 亚洲精品国产av蜜桃| 国产精品一区二区三区四区免费观看| 七月丁香在线播放| 久久精品夜色国产| 91狼人影院| 亚洲成人精品中文字幕电影| 久久99热这里只频精品6学生| 男人狂女人下面高潮的视频| 国产毛片a区久久久久| 免费播放大片免费观看视频在线观看| 日本欧美国产在线视频| 亚洲三级黄色毛片| 欧美日韩在线观看h| 搡老妇女老女人老熟妇| 在线观看人妻少妇| 精品人妻偷拍中文字幕| 国产黄片视频在线免费观看| 在线观看人妻少妇| 午夜福利网站1000一区二区三区| 一级毛片久久久久久久久女| 亚洲精品456在线播放app| 亚洲av二区三区四区| 欧美日韩精品成人综合77777| 两个人视频免费观看高清| 亚洲精品乱码久久久v下载方式| 看十八女毛片水多多多| 亚州av有码| 成年人午夜在线观看视频 | 国产精品国产三级国产av玫瑰| 高清欧美精品videossex| 欧美一级a爱片免费观看看| av在线蜜桃| 色综合站精品国产| 亚洲久久久久久中文字幕| av天堂中文字幕网| 丰满乱子伦码专区| 精品久久久久久成人av| 大陆偷拍与自拍| 亚洲精品色激情综合| 精品一区二区三区视频在线| 少妇熟女aⅴ在线视频| 久久精品夜夜夜夜夜久久蜜豆| 丰满少妇做爰视频| kizo精华| 亚洲av.av天堂| 久久久色成人| 欧美高清成人免费视频www| 深夜a级毛片| 69av精品久久久久久| 亚洲av在线观看美女高潮| 肉色欧美久久久久久久蜜桃 | 国产精品美女特级片免费视频播放器| av专区在线播放| 久久久久久久久久人人人人人人| 亚洲成人中文字幕在线播放| 日本wwww免费看| 亚洲18禁久久av| 欧美日韩国产mv在线观看视频 | 亚洲国产成人一精品久久久| 欧美xxⅹ黑人| 日本av手机在线免费观看| 99久久人妻综合| 日本黄大片高清| 久久久精品94久久精品| 国产精品久久久久久久久免| 美女脱内裤让男人舔精品视频| 国产男人的电影天堂91| 久久精品国产亚洲网站| 欧美日韩综合久久久久久| 日韩av免费高清视频| 非洲黑人性xxxx精品又粗又长| 日本熟妇午夜| 亚洲国产精品专区欧美| 51国产日韩欧美| 久久久成人免费电影| 婷婷色综合www| 久久这里有精品视频免费| av播播在线观看一区| 一二三四中文在线观看免费高清| 国产高清三级在线| av福利片在线观看| 欧美日韩精品成人综合77777| 男女边吃奶边做爰视频| 国产免费视频播放在线视频 | 国产一级毛片七仙女欲春2| 黄片无遮挡物在线观看| 国产一区亚洲一区在线观看| 卡戴珊不雅视频在线播放| 国产精品一二三区在线看| 国产午夜精品一二区理论片| 1000部很黄的大片| 亚洲熟妇中文字幕五十中出| 国产老妇伦熟女老妇高清| 国产探花在线观看一区二区| 亚洲av日韩在线播放| 在线观看一区二区三区| 少妇裸体淫交视频免费看高清| 久久精品夜色国产| 欧美潮喷喷水| 欧美日韩一区二区视频在线观看视频在线 | 中文在线观看免费www的网站| 九草在线视频观看| 哪个播放器可以免费观看大片| 亚洲av电影在线观看一区二区三区 | 超碰av人人做人人爽久久| 免费黄网站久久成人精品| av在线蜜桃| 综合色丁香网| 天堂av国产一区二区熟女人妻| 精品久久久久久久久亚洲| 亚洲精品国产av成人精品| 18+在线观看网站| 亚洲欧洲国产日韩| 天堂俺去俺来也www色官网 | 精品久久久精品久久久| 国产伦一二天堂av在线观看| 九色成人免费人妻av| 嘟嘟电影网在线观看| 夜夜看夜夜爽夜夜摸| 舔av片在线| 内地一区二区视频在线| 街头女战士在线观看网站| av国产久精品久网站免费入址| 欧美潮喷喷水| 丰满人妻一区二区三区视频av| 日韩一区二区三区影片| 亚洲精品一二三| 高清午夜精品一区二区三区| 色尼玛亚洲综合影院| 成人午夜精彩视频在线观看| 亚洲欧美一区二区三区黑人 | 国产成人精品一,二区| 国产精品一区二区在线观看99 | 亚洲国产精品成人综合色| 亚洲精品456在线播放app| 日韩成人伦理影院| 神马国产精品三级电影在线观看| 男女边吃奶边做爰视频| 国产美女午夜福利| 成人一区二区视频在线观看| 91久久精品电影网| 国产av不卡久久| 国产片特级美女逼逼视频| 成人一区二区视频在线观看| 国产高潮美女av| 国产男人的电影天堂91| 欧美日韩在线观看h| 丰满少妇做爰视频| 精品久久久久久久人妻蜜臀av| 亚洲av中文av极速乱| 老司机影院毛片| 女人久久www免费人成看片| 少妇被粗大猛烈的视频| 国产黄色小视频在线观看| 国产美女午夜福利| 国产成人午夜福利电影在线观看| 人妻夜夜爽99麻豆av| 一级毛片电影观看| 欧美3d第一页| 夫妻性生交免费视频一级片| 欧美日韩精品成人综合77777| 在线免费观看的www视频| av国产久精品久网站免费入址| 干丝袜人妻中文字幕| 国产成人a区在线观看| 国内精品美女久久久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产爱豆传媒在线观看| 国产精品一区二区三区四区久久| 成人性生交大片免费视频hd| 天堂影院成人在线观看| 国产成人精品久久久久久| 色综合亚洲欧美另类图片| 亚洲精品影视一区二区三区av| 女的被弄到高潮叫床怎么办| 国产亚洲精品av在线| 久久这里只有精品中国| 日日摸夜夜添夜夜爱| 伊人久久国产一区二区| 人体艺术视频欧美日本| 久久久久久久久久黄片| 中文乱码字字幕精品一区二区三区 | 久久久久免费精品人妻一区二区| av国产免费在线观看| 日本一二三区视频观看| 可以在线观看毛片的网站| 亚洲精品,欧美精品| 91在线精品国自产拍蜜月| 欧美日韩亚洲高清精品| 国产白丝娇喘喷水9色精品| 五月玫瑰六月丁香| 亚洲国产色片| 国产午夜精品一二区理论片| 天天躁夜夜躁狠狠久久av| 观看美女的网站| 亚洲国产精品国产精品| 国产av码专区亚洲av| 七月丁香在线播放| 亚洲国产色片| 亚洲欧美中文字幕日韩二区| 神马国产精品三级电影在线观看| 国产亚洲午夜精品一区二区久久 | 性插视频无遮挡在线免费观看| 午夜精品国产一区二区电影 | 日日啪夜夜爽| 亚洲国产高清在线一区二区三| 国产精品三级大全| 一级毛片黄色毛片免费观看视频| 免费无遮挡裸体视频| 嫩草影院精品99| 国产免费又黄又爽又色| 日日撸夜夜添| 久久精品国产鲁丝片午夜精品| 日本色播在线视频| 免费看日本二区| 日日啪夜夜撸| 97精品久久久久久久久久精品| 国产成人精品久久久久久| 精品久久久久久久末码| 国产精品一区二区在线观看99 | 一级av片app| 久久久亚洲精品成人影院| 午夜免费激情av| 成人亚洲精品一区在线观看 | 日韩伦理黄色片| 成年免费大片在线观看| 久久99蜜桃精品久久| 国产成人aa在线观看| 欧美xxⅹ黑人| 肉色欧美久久久久久久蜜桃 | 亚洲国产成人一精品久久久| 精华霜和精华液先用哪个| 韩国av在线不卡| 国产永久视频网站| 成人av在线播放网站| 肉色欧美久久久久久久蜜桃 | 最近最新中文字幕大全电影3| 免费播放大片免费观看视频在线观看| 中文字幕免费在线视频6| 国产精品福利在线免费观看| 麻豆精品久久久久久蜜桃| 中文字幕av成人在线电影| 2021少妇久久久久久久久久久| 亚洲成人av在线免费| 一区二区三区四区激情视频| 不卡视频在线观看欧美| 一夜夜www| 中文资源天堂在线| 亚洲av一区综合| 赤兔流量卡办理| 精品久久国产蜜桃| 久久这里有精品视频免费| 亚洲av成人av| 久久这里有精品视频免费| 亚洲熟妇中文字幕五十中出| 精品国内亚洲2022精品成人| 国产av不卡久久| 亚洲在线观看片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 免费少妇av软件| 国产 亚洲一区二区三区 | 国产一级毛片在线| 亚洲熟女精品中文字幕| 国产v大片淫在线免费观看| 国产精品无大码| 亚洲一区高清亚洲精品| 久久国产乱子免费精品| 欧美极品一区二区三区四区| 亚洲欧洲国产日韩| 天堂影院成人在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | av专区在线播放| 国产黄a三级三级三级人| 国产日韩欧美在线精品| 亚洲内射少妇av| 免费人成在线观看视频色| 又粗又硬又长又爽又黄的视频| 天堂网av新在线| xxx大片免费视频| 99热这里只有是精品50| 国产精品人妻久久久影院| 十八禁国产超污无遮挡网站| 久久久国产一区二区| 狂野欧美白嫩少妇大欣赏| 亚洲欧美精品自产自拍| 成人av在线播放网站| 国产精品综合久久久久久久免费| 日韩 亚洲 欧美在线| 天天躁日日操中文字幕| 岛国毛片在线播放| 18禁在线无遮挡免费观看视频| h日本视频在线播放| 一二三四中文在线观看免费高清| 国产男女超爽视频在线观看| 99久久精品国产国产毛片| 欧美xxⅹ黑人| 国产美女午夜福利| 亚洲精品,欧美精品| 日本黄色片子视频| 汤姆久久久久久久影院中文字幕 | 国产午夜福利久久久久久| 午夜福利视频1000在线观看| 永久免费av网站大全| 亚洲国产精品成人久久小说| 日韩人妻高清精品专区| 久久精品国产亚洲网站| 哪个播放器可以免费观看大片| 极品教师在线视频| 日韩制服骚丝袜av| 观看免费一级毛片| 国产黄a三级三级三级人| 天堂av国产一区二区熟女人妻| 久久久精品94久久精品| 国产成人一区二区在线| 免费观看性生交大片5| 亚洲最大成人中文| 欧美成人a在线观看| 久久久久久久久久人人人人人人| 亚洲av免费高清在线观看| 国产免费一级a男人的天堂| 亚洲欧美中文字幕日韩二区| 国产极品天堂在线| 一级av片app| 少妇被粗大猛烈的视频| 日韩欧美国产在线观看| 一边亲一边摸免费视频| 亚洲精品一二三| 少妇裸体淫交视频免费看高清| 日本熟妇午夜| 成人无遮挡网站| 人体艺术视频欧美日本| 亚洲国产色片| 只有这里有精品99| 久久久久久久大尺度免费视频| 热99在线观看视频| 国产精品麻豆人妻色哟哟久久 | 日韩国内少妇激情av| 欧美zozozo另类| 亚洲成色77777| 男插女下体视频免费在线播放| 成人毛片60女人毛片免费| 最近的中文字幕免费完整| 国产精品国产三级专区第一集| 又爽又黄无遮挡网站| 又大又黄又爽视频免费| 国产大屁股一区二区在线视频| 中文精品一卡2卡3卡4更新| 日韩欧美一区视频在线观看 | 亚洲成人av在线免费| 高清欧美精品videossex| 一区二区三区高清视频在线| 色视频www国产| 国产一区二区在线观看日韩| 在线观看人妻少妇| 丰满乱子伦码专区| 天天一区二区日本电影三级| 一级毛片我不卡| 成人毛片60女人毛片免费| 最新中文字幕久久久久| 熟妇人妻不卡中文字幕| 免费电影在线观看免费观看| 又黄又爽又刺激的免费视频.| 国产成人精品婷婷| av在线亚洲专区| av在线播放精品| 国产老妇女一区| 2021少妇久久久久久久久久久| 久久热精品热| 如何舔出高潮| 亚洲精品乱久久久久久| 人人妻人人看人人澡| 超碰av人人做人人爽久久| 国产精品.久久久| 一本久久精品| 国内少妇人妻偷人精品xxx网站| 久久99热6这里只有精品| 国产不卡一卡二| 日韩欧美国产在线观看| 久久久精品94久久精品| 国产毛片a区久久久久| 国产在视频线精品| 亚洲精品日韩在线中文字幕| 午夜福利成人在线免费观看| 国产精品国产三级国产av玫瑰| 国产伦精品一区二区三区四那| 久久综合国产亚洲精品| 黄色日韩在线| 亚洲av免费高清在线观看| 亚洲精品日韩在线中文字幕| 亚洲人成网站在线播| 亚洲av免费在线观看| 亚洲人成网站在线观看播放| 啦啦啦中文免费视频观看日本| 欧美最新免费一区二区三区| 人妻少妇偷人精品九色| 亚洲精品视频女| 国产精品蜜桃在线观看| 国产成年人精品一区二区| 亚洲精华国产精华液的使用体验| 91精品伊人久久大香线蕉| 亚洲电影在线观看av| 在线免费观看不下载黄p国产| 欧美日韩精品成人综合77777| 观看免费一级毛片| 国产高清国产精品国产三级 | 亚洲三级黄色毛片| 视频中文字幕在线观看| 夜夜看夜夜爽夜夜摸| 国产中年淑女户外野战色| 欧美bdsm另类| 亚洲在线观看片| 99九九线精品视频在线观看视频| 三级男女做爰猛烈吃奶摸视频| 久99久视频精品免费| 深夜a级毛片| 老师上课跳d突然被开到最大视频| 亚洲国产欧美人成| 欧美不卡视频在线免费观看| 午夜精品在线福利| 只有这里有精品99| 欧美 日韩 精品 国产| 丰满人妻一区二区三区视频av| 秋霞伦理黄片| 日日干狠狠操夜夜爽| 天堂中文最新版在线下载 | 国产在视频线在精品| 国国产精品蜜臀av免费| 黄色配什么色好看| 在现免费观看毛片| 欧美日韩在线观看h| 老女人水多毛片| 精品一区在线观看国产| 女人久久www免费人成看片| 国产片特级美女逼逼视频| 久久国产乱子免费精品| 亚洲人成网站高清观看| 午夜亚洲福利在线播放| 国产免费福利视频在线观看| 国产精品伦人一区二区| 久久久久久久久久黄片| 99re6热这里在线精品视频| 美女cb高潮喷水在线观看| 2021天堂中文幕一二区在线观| 极品少妇高潮喷水抽搐| 国产精品精品国产色婷婷| 大陆偷拍与自拍| 免费在线观看成人毛片| 七月丁香在线播放| 国产午夜精品论理片| 午夜免费激情av| 国产三级在线视频| 丝袜喷水一区| 99九九线精品视频在线观看视频| 全区人妻精品视频| 久久99热6这里只有精品| www.色视频.com| 18+在线观看网站| 97热精品久久久久久| 久久精品人妻少妇| 国产乱人偷精品视频| 激情 狠狠 欧美| 最近中文字幕2019免费版| 免费在线观看成人毛片| 免费电影在线观看免费观看| 亚洲,欧美,日韩| 99久久精品国产国产毛片| 国产一区二区三区av在线| 国产成人精品一,二区| 日韩成人av中文字幕在线观看| 亚洲欧美精品自产自拍| 欧美成人精品欧美一级黄| 天堂√8在线中文| 国产精品一区二区三区四区免费观看| 熟女人妻精品中文字幕| 亚洲国产av新网站| 看黄色毛片网站| 亚洲国产最新在线播放| 自拍偷自拍亚洲精品老妇| 日本一本二区三区精品| 少妇裸体淫交视频免费看高清| 亚洲欧美精品自产自拍| 18禁动态无遮挡网站| 男女边吃奶边做爰视频| 老司机影院毛片| 韩国高清视频一区二区三区| 国产一区亚洲一区在线观看| 99re6热这里在线精品视频| 中文乱码字字幕精品一区二区三区 | 国产成人aa在线观看| 九九久久精品国产亚洲av麻豆| a级一级毛片免费在线观看| 色播亚洲综合网| 亚洲av国产av综合av卡| 老女人水多毛片| 国产黄片视频在线免费观看| 国内少妇人妻偷人精品xxx网站| 亚洲图色成人| 18禁在线播放成人免费| 麻豆久久精品国产亚洲av| 久久草成人影院| 亚洲熟妇中文字幕五十中出| 欧美不卡视频在线免费观看| 一级黄片播放器| 日韩三级伦理在线观看| 人妻制服诱惑在线中文字幕| 亚洲人与动物交配视频| 七月丁香在线播放| 不卡视频在线观看欧美| 国产美女午夜福利| 国产 亚洲一区二区三区 | 色5月婷婷丁香| 亚洲最大成人手机在线| 日韩精品有码人妻一区| 国产亚洲精品av在线| 一本一本综合久久| 亚洲av成人精品一二三区| 亚洲无线观看免费| 精品一区二区三区视频在线| 亚洲伊人久久精品综合| 熟妇人妻久久中文字幕3abv| 久99久视频精品免费| av免费在线看不卡| 人妻制服诱惑在线中文字幕| 免费黄色在线免费观看| 亚洲av免费在线观看| 国产白丝娇喘喷水9色精品| 久久99热6这里只有精品| 国产国拍精品亚洲av在线观看| 黄色日韩在线| 男人狂女人下面高潮的视频| 国国产精品蜜臀av免费| 日本免费在线观看一区| 99热这里只有精品一区| 美女xxoo啪啪120秒动态图| 欧美三级亚洲精品|