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

    Dynamic evolution of low-viscosity fuel particle distribution driven by constant flow

    2023-09-05 08:48:06ZhongKunYang楊仲坤GaoJunAn安高軍XiMengXu徐曦萌ZheZheng鄭哲YongXuWang王永旭LiFengXie解立峰DanZhang張丹andBinLi李斌
    Chinese Physics B 2023年7期
    關(guān)鍵詞:張丹李斌

    Zhong-Kun Yang(楊仲坤), Gao-Jun An(安高軍), Xi-Meng Xu(徐曦萌), Zhe Zheng(鄭哲),Yong-Xu Wang(王永旭), Li-Feng Xie(解立峰), Dan Zhang(張丹), and Bin Li(李斌),?

    1School of Chemistry and Chemical Engineering,Nanjing University of Science and Technology,Nanjing 210000,China

    2Institute of Systems Engineering,AMS,PLA,Beijing 100000,China

    Keywords: low-viscosity fuel,spray and atomization,particle distribution,laser diffraction

    1.Introduction

    With the rapid development of global modernization and industrialization, the dependence on fuel and demand for it have increased.Low-viscosity fuels have been widely used as alternative energy sources and fuel additives for internal combustion engines owing to their inherent advantages.[1–3]Compared with high-viscosity fuels, they can provide good fuel economy.[4]Moreover, they exhibit significant advantages in improving emission quality, reducing fuel consumption, increasing machine power, and reducing failure and maintenance rates.[5,6]The fuel spray and atomization characteristics are the main factors affecting the combustion and emission performances of internal combustion engines.[7]Good dispersion atomization can affect the fuel distribution in the burner to a certain extent, thereby improving the ignition,[8]combustion,[9]and emission characteristics[10,11]of the fuel.Moreover, the formation of fuel–air mixture, which is an important factor in controlling the combustion rate to reduce emissions, is influenced by the fuel atomization characteristics.Reitz and Duraisamy[12]and Biswaset al.[13]found that a shorter mixing time of fuel and air can improve engine performance to some extent.Effective fuel combustion depends on an effective atomization process.[14]This is determined by the physicochemical qualities of the fuel and the injection settings.[15]

    Several studies have been conducted to determine the spray properties of gasoline and diesel fuels;[16,17]however,some disagreement about the macroscopic behavior still exists.In the process of fuel injection and atomization, owing to the strong aerodynamic force,the fuel is broken and atomized almost instantaneously; therefore,it is necessary to have a clear and accurate understanding of the microscopic characteristics of the physical process of the spray and to conduct a comprehensive study of its atomization process and cloud development.[18]The droplet size after spraying is a key parameter for evaluating the atomization performance from a microscopic perspective.[19,20]Parket al.[21]studied the atomization characteristics of gasoline and diesel fuel blends and showed that the fuel droplet size decreases,the ignition delay time increases,and the NOxand soot emissions decreases.

    Hanet al.[22]also conducted an experimental study of the spray characteristics of diesel and gasoline fuel blends,but they mainly focused on the effects of the proportion of gasoline on volume and mass injection rates and spray cone angle.In their study Kimet al.[23]obtained the conclusion that narrower spray cone angle has better performance in terms of maximum combustion pressure,maximum heat release rate and ignition delay.It is advantageous under the early combustion strategy.However, the spray penetration decreases and the spray width increases, which promotes the mixing of fuel and gas.This was mainly attributed to the blending of low viscosity fuel additives.[24]Linet al.[25]combined the microscopic measurement means of a phase Doppler particle analyzer to study the characteristics of a liquid spray under the action of high-speed airflow.They found that the local droplet size of the distribution of clouds significantly affects the overall cloud spray characteristics.By determining the time-averaged properties of the spray,the researchers converted the droplet diameters in the center of the axial flow into the surface at the measurement points on the cross-section.Kannaiyan and Sadr[26]used a laser diagnostic measurement technique for spray experiments on jet fuel under ambient atmospheric conditions and reported the effect of particle size parameters on the microscopic spray performance.Recently,they compared the microscopic spray characteristics of jet fuel in spray chambers.[27]The droplet states in different regions of the spray process were derived from the measurements at the different spray axial positions.Changet al.[28]investigated droplet behavior under ultra-high injection pressure conditions by using a particle image analysis detection technique.Parameters, such as the Sauter mean diameter (SMD/D32) and droplet number, were used to evaluate the spray atomization characteristics.[29]Most of studies have been performed in constant-volume containers.Some researchers have also studied the volume of the test container,arguing that the constant volume of the closed chamber affects the atomization characteristics of the liquid.[30]Some researchers also used the principle of maximum entropy and the law of conservation of mass to deduce and study the droplet size and distribution,and they used models to predict the droplet formation characteristics of diesel sprays.[31–34]From all the findings they concluded that the spray performance is closely related to the microscopic behavior of the droplets in the spraying process and cloud development,[35]which would serve as the key motivation for this study.

    In contrast,fuels can form mixtures with air in the event of a leakage during production, storage, and transportation.The mixture has a risk of explosion when exposed to an ignition source.Especially in high-pressure situations, such as broken pipes and containers,fuel injection,accompanied with high-speed air flow,may make spray explosions more likely to occur.[36–38]However,there is little research on low-viscosity fuels.Liuet al.[39]emphasized the necessity of investigating the spray particle size.They investigated the effect of the diameter ofn-hexane on the explosion characteristics of its mixture with air.Both the explosion pressure and temperature,which represent the explosion hazard, were inversely associated with the particle size of the cloud.Lyuet al.[40]found that the variation in the SMD has a remarkable effect on the explosion characteristics.Particle size also affects the ignition sensitivity.Recently, Wuet al.[41]investigated the atomization and spray explosion properties of methanol based on a 20L spherical container.The variations in the spray particle size measured by the laser diffraction technique are consistent with the changes of explosion pressure.

    Therefore, from the perspective of improving combustion performance or preventing explosion accidents,it is worthy of attention to evaluate the microscopic atomization performance by the droplet size after spraying.[42]At present,studies mainly focus on measurements under fixed-test conditions,with research and exploration of the combined effects of multiple factors lacking.Existing studies have paid insufficient attention to the spray characteristics of low-viscosity fuels, so it is necessary to investigate their atomization characteristics and cloud development status in depth and detail.To this end, in this work we study and compare the combined effects of the fuel mass concentration (ratio of fuel mass to spray chamber volume,c=200, 400, 600, 800, and 1000 g/m3)and injection pressure(pressure of the driving airflow at the fuel inlet,pi=0.50, 0.75, 1.00, 1.25, 1.50, and 2.00 MPa)on low-viscosity fuels under ambient conditions of 1 atm (1 atm=1.01325×105Pa) and 20?C.The fuel droplet size parameters are measured by using laser diffraction techniques, and the droplet formation characteristics of different types of low-viscosity fuels are compared at several axial and radial measurement positions downstream of the nozzle outlet.

    This study fill the gaps in the atomization characteristics and spray development mode of low-viscosity fuel, and provides references for improving combustion performance and avoiding accidental disasters.

    2.Experimental setup and methods

    2.1.Experimental setup

    The spray tests were conducted in a 20-L spherical container as shown in Fig.1.The inner diameter of the spray chamber was 340 mm, and the thickness of the wall was 10 mm.The device was equipped with two quartz optical windows(each with an observation diameter of 100 mm),a pneumatic atomization system,cooling circulation system,and terminal control system.The centers of the two optical windows were aligned on the same level as the center of the spray chamber,which was used for better optical and laser diagnostic measurements.Owing to the presence of optical windows and material limitations,the maximum pressure bearing limit of the spray chamber was 2 MPa.The nozzle was installed at the same level as the center of the spray chamber,while the two optical windows were orthogonal to the vertical axis of the nozzle (see Fig.1).The nozzle was connected to the sample chamber by an explosion-proof solenoid valve.The sample chamber was pressurized with high-pressure air and injected through the nozzle at the required pressure.The pneumatic atomization system was equipped with a pressure sensor to monitor the injection pressure.

    Fig.1.Schematic diagram of the experimental setup.

    2.2.Optical/laser diagnostic measurement

    A laser diffraction (LD) system from Malvern Panalytical Spraytec was used for point-by-point measurements of the spray particle and droplet size distributions in real time.The LD system used for the experimental measurements is illustrated in Fig.2.It was equipped with a highly stable He–Ne laser, a transmitter, a receiver, an optical bench, and a signal processing software.The system used the laser diffraction principle to measure the particle size by measuring the intensity of the diffracted light as the laser beam passed through the spray.Subsequently, the size of the droplets that created the diffraction pattern was analyzed and calculated based on Mie theory[43–45]and Fraunhofer scattering optical model.For the LD system with a red parallel laser of wavelengthλ=633 nm,a quartz glass with light transmission wavelength in the range 633±10 nm and transmission≥98%was used to eliminate interference from other wavelengths to achieve the optimal measurement accuracy.[37,39]The operating parameters of the LD system are presented in Table 1.

    Fig.2.The LD system.

    Table 1.Operating parameters of the LD system.

    2.3.Materials

    Considering the viscosity characteristics of the fuel, a low-viscosity fuel with viscosity less than 3 mPa·s was used for the experiments in this study.1-butanol and cyclooctane were used as fuel substitutes and additives,which can increase the calorific value of fuel and improve the combustion efficiency, providing better fuel economy than traditional fuel economy.Also,they can effectively reduce the carbon smoke emission in the combustion process under the same quality,which is strategically important in energy saving and emission reduction.The properties of the experimental materials are presented in Table 2.

    Table 2.Properties of experimental materials.

    2.4.Methods

    The experiments were conducted under the following environmental conditions: the initial ambient pressure in the spray chamber was 0.10 MPa, and the initial ambient temperature and fuel temperature were maintained at 293 K by using the cooling circulation system.The particle size measurement experiments were performed according to the following steps.(i)The spray chamber and quartz optical glasses were cleaned.(ii) Fuel was added into the sample chamber.(iii)The control system was initialized.The air inlet solenoid valve was opened,and high-pressure air was injected into the sample chamber.The air inlet solenoid valve was closed when the pressure reached the preset spray pressure.(iv)The spray solenoid valve was opened (continuous opening time was set to 100 ms),and the fuel was sprayed from the nozzle together with the high-speed airflow.(v)The Spraytec LD system was automatically triggered off, which tested and collected data.(vi)The chamber was cleaned and flushed repeatedly with dry air.Each spray measurement was performed after the environmental and operating conditions had reached a stable and uniform state.

    The particle size measurements in this study were performed downstream of the nozzle outlet at different axial (x)and vertical radial (z) positions, as shown in Fig.3.Because the laser measurement range was limited, the optical bench was adjusted to obtain measurements at different positions.The radial orientation on one side was selected for the measurement because the spray was axially symmetric about the spray direction.

    Fig.3.Details of particle size measurement locations.

    2.5.Experimental uncertainty analysis

    The experimental uncertainty in the results of the spray characteristics primarily arises from the test apparatus, fuel addition,initial ambient pressure,and initial ambient temperature.The LD system is based on Mie theory, where large and small particles are diffracted against the laser at small and large angles, respectively.Thus, the particle size was determined by the intensities of the diffracted light at different angles.The measurement of the ambient background before the experiment contributed to the accuracy of the experimental data.However, the uncertainty of the particle size measurement was confirmed to be under 5%by using a monodisperse droplet generator.Samples were added by using a pipette with an accuracy of 1% and volume rang of 0 mL–50 mL in this study.Thus,the fuel content in the sample chamber was estimated to have a maximum uncertainty of±0.5 mL.The initial ambient pressure was measured by using a pressure gauge with an accuracy of 1.6% in the pressure range of 0 MPa–1 MPa.The initial ambient temperature was measured by using a digital electronic thermometer with an accuracy of 1%in the temperature range of 263 K–363 K.Therefore, the pressure and temperature measurements in this study displayed maximum errors of±16 kPa and±1 K,respectively.

    3.Results and discussion

    The droplet size distribution of a spray is a critical indicator in the evaluation of atomization characteristics.To characterize the fineness and uniformity of the fuel spray, the arithmetic mean droplet diameter(D10)and Sauter mean diameter(SMD/D32), which are of significance for the combustion industry,are quantified and compared.[46]

    3.1.Variations of spray particle size under effects of mass concentration and injection pressure

    Experiments on 1-butanol spray parameters, including droplet diameter and size distribution,are conducted with variations in fuel mass concentration and injection pressure.Figure 4 presents a comparison among the particle sizes of the 1-butanol spray as a function of mass concentration at various injection pressures.As shown in Fig.4(a), the diameter of the 1-butanol droplet increases with the concentration increasing, and the variation inD10presents a good linear fit.The statistical information of the linear fit curves is shown in Table 3.This trend is consistent for all the cases investigated under this condition.Atpi=0.50 MPa, the diameter of 1-butanol at a concentration of 1000 g/m3is the largest, 1.85 times larger than that at 200 g/m3.The increase in the concentration weakens the driving force of the spray airflow per unit mass of fuel, causing the spray velocity at the nozzle exit to decrease.Further analysis shows that the decrease in the turbulence intensity leads the droplet diameter to increase.This implies that fuel with lower concentration has better atomization.As the injection pressure increases, the droplet size decreases for the same concentration 1-butanol.The maximum difference inD10betweenpi= 0.50 MPa and 2.00 MPa is observed at the concentration of 1000 g/m3,which is reduced by 77.55%.This phenomenon can be explained based on dimensionless analysis,where the augment of the injection pressure increases the Mach number(Ma)of the airflow under the condition of constant ambient pressure, resulting in a larger driving force.The resulting increase in the Weber number(We) promotes the generation of small droplets of the spray.Furthermore, figure 4(a) shows that with the increase of the injection pressure, the effect of the change in concentration has less significant effect onD10.As the injection pressure rises from 0.50 MPa to 2.00 MPa,the slope of the fitted curve ofD10versusconcentration drops from 0.00864 to 0.00163,which is a reduction of 81.1%.The differences inD10among the sprays at different concentrations decrease as the injection pressure increases.Therefore, a higher fuel concentration requires a larger driving force to obtain smaller droplets, so as to promote the atomization effect.Likewise, a lower injection pressure is not sufficient to promote the fragmentation of the fuel droplets.It can be observed that under a low-pressure difference, the effect of fuel concentration on the droplet diameter is more significant.

    Figure 4(b)shows the comparison ofD32of the 1-butanol spray with the rise of concentration under various injection pressures, whereD32represents the more general average diameter of the fuel spray.Aside from atpi=2.00 MPa, the same variation trends for particle size as that forD10are observed.It is noteworthy that the minimum value of theD32curve forpi=2.00 MPa is observed at the maximum concentration, which is contrary to the trend of the other curves.In particular,at a low concentration,D32atpi=2.00 MPa is higher than at other injection pressures.The reason for this is the fuel breaking into small droplets in the atomization process,which is also known as secondary atomization.However,the motion direction of the droplets in the spraying process is not definite, and the probability of collisions between the droplets is high.Thus, some small droplets reaggregate into big droplets during the dispersion, which is called the aggregation effect.Particularly, the combined effect of low concentration and highMaairflow creates an extremely strong turbulence intensity,which causes a more unstable spray process,resulting in an increased probability of collision between droplets and an obvious aggregation effect.This theory is further confirmed by comparing the droplet size distribution with cumulative number percentage of 1-butanol atpi=2.00 MPa.This is consistent with the conclusion concluded by Zhanget al.[47]that the augment of injection pressure further increases the irregularity of the spray.

    Fig.4.Effects of mass concentration (c) on arithmetic mean droplet diameter(D10,left)and Sauter mean diameter(D32,right)of 1-butanol at different injection pressures(pi)(x=100 mm;z=0 mm).

    Figure 5 shows the droplet size distribution and cumulative number percentage for the casesc= 200 g/m3and 1000 g/m3.At low concentration, the droplet size exhibits a shift toward smaller diameters, which demonstrates the decreasing trend ofD10.However, the collisional aggregation of droplets also results in the appearance of a small fraction of larger diameter droplets, which affects theD32value to some extent.Thus, the low concentration ofD32appears to rise rather than drop atpi=2.00 MPa.The display in Fig.5 further supports the above analysis.

    Table 3.Statistical information for linear fit curves.

    Fig.5.Comparison of droplet size distribution and cumulative number percentage of 1-butanol spray for c=200 g/m3 and 1000 g/m3 at pi=2.00 MPa(x=100 mm;z=0 mm).

    The variation curves of the particle size with injection pressure for different concentrations of butanol and water spray are shown in Fig.6.The arithmetic mean droplet diameter (D10) is described as a function of the mass concentration and injection pressure as shown in Fig.6(a).Both 1-butanol and water exhibit the best atomization effect forpi=2.00 MPa,andD10drops to below 4μm.Thus,aerodynamic forces remain dominant in the entire process of dispersion and atomization.Droplet fragmentation and atomization are improved with the increase ofMaof the airflow.Moreover,the values ofD10of 1-butanol decrease by 77.55%and 46.04%with the injection pressure and concentration increasing, respectively, and by 85.35% with the coupling effect increasing, while for water, the values ofD10of 1-butanol decrease by 60.86%, 28.99%, and 59.27%, correspondingly.The low concentration and high injection pressure promote the fragmentation and atomization of the spray.However,the droplet particle sizes are relatively similar for different concentrations at the same injection pressure, thereby further demonstrating that the injection pressure is the main factor affecting the fuel spray,reducing the particle size to its limit,and weakening the concentration effect.However, the decreasing trend ofD10gradually decrease with injection pressure increasing,indicating that the injection pressure significantly contributes to the spray effect at high concentrations.Water exhibits a smallerD10as shown in Fig.6(a).This is because water has a lower viscosity than 1-butanol,and it is easier to break into small droplets after experiencing the same high-speed airflow,achieving a significant atomization effect.Furthermore,asMaof the jet airflow increases, a smooth change in the particle size is observed earlier for water than for 1-butanol.This phenomenon is more obvious at low concentrations.

    Figure 6(b) shows the comparison among the values of SMD(D32)which is a function of the mass concentration and injection pressure for 1-butanol and water sprays.Majority of the data in the graph maintain the same trend as that forD10.The droplet size distribution and cumulative number percentage explain this phenomenon more clearly as shown in Fig.7(a).It can be observed that the droplet distribution of the fuel spray atpi=2.00 MPa shifts toward smaller sizes.TheD10andD32decrease at higher injection pressures,which is mainly attributed to the contribution of the small droplets.This further confirms that a high injection pressure significantly promotes the droplet breakup.However, for 1-butanol droplets at low concentrations and high injection pressure,D32appears to become larger,and the maximum value is observed at the minimum concentration.This is because the influence of low concentration and highMalead theWeof the high-speed airflow to increase.Thus,the increase in the turbulence intensity results in the small droplets having a higher probability of aggregating into larger droplets.Figure 7(b)displays a visual illustration of the theory.The droplet size distribution and cumulative number percentage at low concentrations are presented.The high injection pressure promotes droplet breakup,resulting in the increase of the number of small droplets,which contributes toD10decreasing.Droplet aggregation is also observed,and the increase ofD32is mainly caused by the contribution of the little fraction of large droplets.Similarly,D32of the water spray demonstrates an increasing trend with injection pressure increasing.It should be noted thatD32values of water forc=600 g/m3and 200 g/m3show increasing trends atpi=1.5 MPa and 1.25 MPa,respectively,compared with the variation observed for 1-butanol atpi=2.00 MPa.The difference in viscosity between butanol and water may be a key factor influencing this phenomenon.The lower viscosity of water enables the aggregation effect in the spray process to be more prominent and to occur easily at smaller mass concentrations.Based on the dimensionless analysis, it is suggested that the high-speed airflow along with highWesignificantly promotes the droplet aggregation of low-viscosity liquids.The effect of fuel viscosity is analyzed and discussed in later sections.TheD32of the water spray presents an increasing trend;however,it shows a drop atpi=2.00 MPa,indicating that the highWeof the airflow improves the probability of collision between the droplets.However,it is further confirmed that the highWealso promotes the breakup of droplets.This can also be observed in Fig.7(b).

    Fig.6.Effects of injection pressure(pi)on(a)arithmetic mean droplet diameter (D10) and (b) Sauter mean diameter (D32) of 1-butanol and water at different mass concentrations(c)(x=100 mm and z=0 mm).

    Fig.7.Comparisons among droplet size distributions and cumulative number percentages of 1-butanol spray for different injection pressures at c=1000 g/m3 (a)and 200 g/m3 (b)(x=100 mm and z=0 mm).

    3.2.Variations of spray particle size in axial direction

    Experiments are conducted to compare the droplet diameters and size distributions of three liquid sprays subjected to different injection pressures at different axial distances.Figure 8 shows comparisons of the particle sizes of the test liquids under different injection pressures at axial locations downstream of the nozzle exit(x=50 mm,100 mm,and 150 mm).TheD10andD32both decrease with injection pressure increasing,and their downward trends gradually decrease.This trend is consistent for all the axial locations considered in this work.The magnitude of the particle size reduction is attributed to the increase of the injection pressure.The comparison of the data with respect toD10shows that the droplet size near the nozzle exit is smaller than the droplet size far from the nozzle exit, and the droplet size increases fromx= 50 mm to 150 mm.This indicates that the increase in axial distance results in enhanced evaporation of tiny droplets,and hence larger droplet sizes are observed.In addition,the low viscosity of the test fuel facilitates the disaggregation of the fuel droplets,thus making it easier to produce smaller droplets.In this process,the probability of collision between the droplets increases and the small droplets tend to collide and re-agglomerate into larger droplets, which was not observed in the experiments with liquids of higher viscosity, however, which is demonstrated more distinctly by the differences in the droplet size distribution and cumulative number percentage in the axial direction, as shown in Fig.9.The tendency of the droplet size distribution to move toward larger diameters is the main reason for the increase of the diameter in the axial direction.Comparing with the location near the nozzle exit,the spray away from the nozzle exit not only reduces the small sized droplets, but also increases the number of the larger sized droplets, which responsible for the trend seen in Fig.8.

    Fig.8.Comparisons of arithmetic mean droplet diameter(D10,top)and Sauter mean diameter(D32,bottom)of fuel sprays for different values of injection pressure(pi)at axial locations,x=50 mm,100 mm,and 150 mm,downstream of the nozzle exit(c=1000 g/m3;z=0 mm).

    At locations close to the nozzle exit (x= 50 mm–100 mm), the difference in droplet diameter decreases as the injection pressure increases.The minimum difference is observed atpi=2.00 MPa.However, it is noteworthy that the difference in droplet size on the spray axis (x= 100 mm–150 mm) is different in opposite directions.As the injection pressure increases, the difference augments.The droplet diameter increases along the spray axis due to the aggregation effect.However,the droplet size difference caused by droplet aggregation gradually decreases as the injection pressure increases.It is worth emphasizing that a higher injection pressure promotes breakage of the droplets.This indicates that the increase in injection pressure weakens the effect of small droplet aggregation and promotes droplet fragmentation.In the“Waterc=200 g/m3”group in Fig.6(b),the droplet size increases due to droplet aggregation and reaches a peak atpi=1.75 MPa.The droplet diameter shows a decrease when the injection pressure reaches 2 MPa, when droplet fragmentation dominates.Moreover, the smaller mass concentration also results in a larger airflow driving force per unit mass of fuel, which also has a beneficial effect on droplet fragmentation.Similarly,near the nozzle exit,the higher-velocity airflow breaks the aggregated large droplets into smaller droplets,weakening the effect of aggregation and reducing the difference.However, a significant difference is observed downstream of the nozzle.Figures 9(a) and 9(b) show comparisons of changes of droplet size distribution and cumulative number percentage in the axial direction at low and high injection pressures,respectively.At low injection pressures,the difference in size distribution in the region ofx= 50 mm–100 mm is evident,indicating a significant aggregation effect.However,a different trend is observed at higher injection pressures.AlthoughD10andD32forx=100 mm are generally higher than those forx=50 mm,the droplet size distribution is more concentrated in the small diameter range.It is observed that the increase in injection pressure inhibits droplet aggregation and promotes breakup.However, the size distribution in the region ofx=100 mm–150 mm tends to have a large diameter range, and aggregation occurs again.It is observed that the droplet breakup and aggregation effects coexist and suppress each other in the whole spraying process,which affects the spraying characteristics of the fuel.Except for the water droplets,the curves quantified and compared forD32in Fig.8 were similar to those forD10.When the injection pressure increases to 1.75 MPa, the droplet particle size decreases near the nozzle exit(x=50 mm–100 mm).Furthermore, the trend is more prominent as the injection pressure increases, owing to the significant contribution of the majority of the small droplets.This indicates that droplet breakup and atomization dominate at this stage.The beneficial effects of the low viscosity of the liquid and highMaairflow on the fuel spray are further demonstrated.However, the droplet diameter atx=150 mm increases slightly, indicating that the droplet fragmentation consumes the airflow driving force,and the aggregation effect continues to dominate as the axial distance increases.Atx=50 mm–100 mm,the increase of small droplets further promotes the aggregation of droplets.

    Fig.9.Comparison of droplet size distribution and cumulative number percentage of 1-butanol spray for different axial locations at pi =0.50 MPa(a)and 2.00 MPa(b)(c=1000 g/m3;z=0 mm).

    Comparisons ofD32between liquids with different viscosities under various injection pressures at axial locations,x= 50 mm, 100 mm, and 150 mm, are shown in Fig.10.The variation inD32is influenced by the viscosity of the test liquid.At lower injection pressures, the experimental results have larger error bars.The main reason is that the low injection pressure produces a lower airflow driving force and a small and uneven force per unit mass of fuel at the moment of airflow driving, so any slight fluctuation in injection pressure will affect the spraying effect.This is also illustrated by the significant effect of lower injection pressure on droplet diameter.As shown in Table 2,the viscosity of the tested liquids is in the order of 1-butanol>cyclooctane>water;therefore,larger droplet sizes are observed in the experiments for the samples with higher viscosities.The same trend is observed in the axial direction downstream of the nozzle exit.However,the breakup and atomization of droplets and the aggregation effect affecte the variations in the droplet size with increasing axial distance in the whole spray development process.Less viscous fuels produce smaller droplets,which is consistent with the previous analysis.Further analysis demonstrates that the less viscous fuel is more easily broken into small droplets after experiencing the high-speed airflow,and the atomization effect is more evident.In the spraying process, viscosity affects the droplet formation characteristics.In addition, the increase of injection pressure has a more significant effect on the reduction in droplet diameter for high-viscosity liquids.

    The effect of droplet aggregation on the fuel spray characteristics is significant, and collision-induced fragmentation and aggregation have a significant effect on the droplet size distribution.The droplet position distribution will have an effect on the subsequent spray combustion, including ignition energy and ignition position.Through the distribution of droplets at different positions,the optimal ignition conditions are set up, and the spray position or ignition state is changed appropriately to achieve the best combustion state and improve the combustion efficiency.The work done in the confined space is enhanced, and the heating effect is obvious.In the open environment,the velocity of the flow field increases,and the thrust force is enhanced.Therefore,the aggregation effect of droplets is an important factor when exploring the spray characteristics and flow field state.

    Fig.10.Effects of viscosity on Sauter mean diameter(D32)for different injection pressures(pi)at axial locations,x=50 mm(a),100 mm(b),and 150 mm(c),downstream of nozzle exit(c=1000 g/m3;z=0 mm).

    3.3.Variations of spray particle size in radial direction

    Experiments are conducted to compare the droplet diameter and size distribution of three liquid sprays subjected to different injection pressures at different radial distances.Figure 11 shows the comparisons of the particle sizes of the test liquids under various injection pressures at radial locationsz=0 mm, 25 mm, and 50 mm for the axial locationx=100 mm.The trend of the decreasing of size of the droplet at different radial locations with injection pressure increasing is consistent with that reported in previous study.TheD10andD32decrease from the spray axis(z=0 mm)toward the spray boundary,and the decrease is gradually reduced.After the entrainment effect is rapidly enhanced,the cloud diffuses in the radial direction in the spraying process, increasing the number of smaller droplets.Furthermore, owing to the effect of air turbulence and entrainment, the large sized droplets move axially with the airflow,whereas the small droplet population deviates outward.The fuel at the spray boundary is broken into a fewer number of droplets.The density of droplets in the peripheral region is less than that in the axis region,resulting in insignificant droplet aggregation.This results in the droplet size decreasing in the radial direction.The droplet size distribution of 1-butanol in the radial direction atc=1000 g/m3andpi=0.5 MPa is compared as shown in Fig.12(a).It can be clearly observed that the size distribution of the droplet population at the spray boundary shifts toward a smaller diameter range.The decreasing trends ofD10andD32presented in Fig.11 are mainly due to the increase in the number of smaller droplets.This is also present in other cases.

    Atz=50 mm, the decrease in the droplet size becomes smaller than that atz=0 mm and 25 mm.The main reason for this phenomenon is that the variation of the spray from the spray axis to the spray boundary is due to radial expansion and the distance from the spray axis, which leads to the decrease of the turbulent intensity of the airflow, and turbulent intensity is an important parameter affecting droplet breakup and atomization.Furthermore, in all the cases investigated,the difference in droplet diameter in the radial direction gradually narrows as the injection pressure increases.This indicates that the injection pressure significantly promotes more intense droplet breakup and atomization near the spray axis compared with that in spray boundary region.It is noteworthy that an increasing trend in the droplet diameter in the radial direction is observed at higher injection pressures.When the injection pressure is high, the droplet diameter near the spray axis is smaller than that at the spray boundary, which is particularly evident in the case of low-viscosity fuels.It is clearly observed in the spray of water droplets in Fig.11.The stronger droplet fragmentation near the spray axis region increases the outward deflection of smaller droplets and the droplet density at the spray boundary.The larger spray pressure also strengthens the radial turbulence intensity,leading to a significant aggregation effect.This indicates that the promotion of breakup via the high-Maairflow is significant near the spray axis.To some extent, the trend of droplet reduction in the radial direction can be reversed to achieve good atomization in the axial flow region.This can be more intuitively explained by comparing the changes in the droplet size distribution and cumulative number percentage in the radial direction under low and high injection pressures as shown in Figs.12(a)and 12(b),respectively.Although the high injection pressure also promotes the size distribution atz=50 mm to shift toward smaller diameters, the effect is subtle.By contrast,the size distribution atz=0 mm exhibits a significant shift.Therefore,the low viscosity of the liquid and highMaairflow have equally beneficial effects in the radial direction.

    Fig.11.Comparisons of arithmetic mean droplet diameter (D10, top) and Sauter mean diameter (D32, bottom) of fuel sprays for different injection pressures(pi)at radial location z=0 mm,25 mm,and 50 mm for c=1000 g/m3,x=100 mm.

    Fig.12.Comparisons of droplet size distribution and cumulative number percentage of 1-butanol spray for different radial locations at pi =0.50 MPa(a)and 2.00 MPa(b)for c=1000 g/m3,x=100 mm.

    Fig.13.Effects of viscosity on Sauter mean diameter(D32)for different injection pressures(pi)at radial locations z=0 mm(a),25 mm(b),and 50 mm(c)for c=1000 g/m3,x=100 mm.

    Comparison ofD32between liquids with different viscosities under different injection pressures at radial locationsz=0 mm,25 mm,and 50 mm is shown in Fig.13.The droplet diameter variation follows the same trend as the fuel viscosity,with a lower viscosity liquid possessing a better breakup and atomization in the radial position,which is consistent with that observed in the spray axis area.Change in the droplet size during spray development does not change the effect of the intrinsic viscosity of the liquid on the atomization effect.Moreover,the effect of the injection pressure on the difference between the droplet diameters is more evident in the sprays of highviscosity liquids.As shown in Figs.10 and 13,low-viscosity fluids appear to change smoothly with injection pressure increasing at an earlier stage compared with the high-viscosity fluids.It can be seen from Fig.13 that the effect of liquid viscosity on droplet particle size varies at different radial positions.The differences in radial droplet size among different viscosity fuels are mainly in the region close to the spray axis,while the difference at the spray boundary is not significant.Owing to air entrainment,the larger droplet sizes move axially with the air stream,and the increase in spray pressure mainly promotes the breakup of droplets close to the spray axis.This is especially true for the larger viscosity liquids, where the larger density droplet groups move along the axial direction.Meanwhile, smaller droplets shifted toward the spray boundary are affected by the weaker turbulence intensity, and the effect of fuel viscosity on particle size difference is much less.Although the increased injection pressure also enhances the radial driving force,the overall low density and small particle size of droplets at the spray boundary make the particle size difference of different viscosity fuels insignificant.

    4.Conclusions

    The effects of mass concentration and injection pressure on the spray characteristics of fuels in a 20-L spray chamber are investigated,which are quantified and compared for three different test fuels.To this end,test experiments are conducted at seven injection pressures(pi=0.50,0.75,1.00,1.25,1.50,1.75,and 2.00 MPa)by using high-speed airflow as the driving force.The characteristics of the spray droplets are measured by using Malvern Panalytical Spraytec at three axial locations(x=50, 100, and 150 mm) and three radial locations (z=0,25, and 50 mm)downstream of the nozzle exit.From the results,the main conclusions are obtained below.

    (I) The increase of mass concentration has an inhibitory effect on atomization,and the variation in the particle size with mass concentration can be fitted linearly.Moreover,the mass concentration significantly affects the particle size of the spray droplet at low injection pressures.

    (II) The increase of injection pressure enhancesMaand driving force of the airflow,which promotes droplet breakup;however,the effect gradually diminishes.Similarly,the injection pressure significantly affects the spraying effect at high mass concentration.

    (III) In comparison with the mass concentration, the injection pressure dominates the spraying process and has a greater influence on the spraying effect, reducing the droplet diameter to its limit and weakening the effect of the mass concentration.

    (IV) Although the low mass concentration and high injection pressure promote droplet breakup, they also increase the probability of collision between the small droplets and increases the likelihood of the aggregation effect.This effect is more prominent in low-viscosity fuels, which was not previously observed in higher-viscosity liquid sprays.

    (V) In the axial direction, the droplets tend to increase owing to the aggregation effect.Moreover,because of the air turbulence and entrainment effect, the large droplets develop in the axial direction, whereas the small droplet population deviates outward.Therefore,the droplet diameter tends to decrease in the radial direction; however, the trend is reversed when the injection pressure increases to a certain level.

    (VI) The high-speed airflow significantly promotes the breakup of droplets in the region close to the nozzle exit and spray axis and suppresses the aggregation effect to a certain extent,where the difference in droplet diameter gradually decreases with injection pressure increasing.The effect is minimal at the spray boundary and far from the nozzle outlet.

    (VII) Droplet fragmentation and aggregation effects coexist and mutually inhibit in the whole spraying process.Nevertheless, the trend of the viscosity of the fuel on the droplet size does not change.The low-viscosity fuels exhibit better atomization, and smoother changes are observed with injection pressure increasing at an earlier stage compared with the high-viscosity fluids.

    Acknowledgement

    Project supported by the Young Scientists Fund of National Natural Science Foundation of China (Grant No.11802136).

    猜你喜歡
    張丹李斌
    World Wetlands Day
    The Wizard of Oz
    The Wizard of Ozby L. Frank Baum
    介詞的時間搭配
    The Wizard of Ozby L. Frank Baum
    Factory Life in the 1800's
    Probing Nonclassicality of Two-Mode SU(2)Generator Based on Quantum Fisher Information?
    李斌:蔚來為未來而來
    金橋(2018年4期)2018-09-26 02:25:08
    Application of Communicative Approach to Junior English Teaching
    嚇人奶奶,新年快樂
    99国产精品一区二区三区| 亚洲欧美日韩卡通动漫| 久久中文看片网| 高清在线国产一区| 免费人成视频x8x8入口观看| 国产黄a三级三级三级人| av中文乱码字幕在线| 亚洲电影在线观看av| 亚洲欧美激情综合另类| 午夜免费观看网址| 国产一区二区在线av高清观看| 国产精品av视频在线免费观看| 国产精品香港三级国产av潘金莲| 久久精品人妻少妇| 午夜福利在线观看免费完整高清在 | 国产亚洲精品久久久久久毛片| netflix在线观看网站| 99国产精品一区二区蜜桃av| 欧美三级亚洲精品| 中文亚洲av片在线观看爽| 午夜福利成人在线免费观看| 国产一区二区激情短视频| 免费在线观看影片大全网站| 久久久久精品国产欧美久久久| 亚洲欧美精品综合久久99| 国产精品女同一区二区软件 | 俄罗斯特黄特色一大片| 精品一区二区三区视频在线观看免费| 麻豆一二三区av精品| 99久久无色码亚洲精品果冻| 老鸭窝网址在线观看| 久久久久九九精品影院| 制服丝袜大香蕉在线| 免费观看人在逋| 日韩欧美 国产精品| 丰满人妻熟妇乱又伦精品不卡| svipshipincom国产片| 老司机福利观看| 中国美女看黄片| 搡老岳熟女国产| 日韩欧美 国产精品| 成人国产一区最新在线观看| 每晚都被弄得嗷嗷叫到高潮| 一级黄色大片毛片| 脱女人内裤的视频| 国产乱人伦免费视频| 成熟少妇高潮喷水视频| 国产麻豆成人av免费视频| 少妇的丰满在线观看| 中文字幕精品亚洲无线码一区| 99久久精品国产亚洲精品| 国产精品美女特级片免费视频播放器 | 99国产精品99久久久久| 久久久久国产精品人妻aⅴ院| 全区人妻精品视频| 国产极品精品免费视频能看的| av中文乱码字幕在线| 1024香蕉在线观看| АⅤ资源中文在线天堂| 69av精品久久久久久| 日本五十路高清| 日韩 欧美 亚洲 中文字幕| 久久99热这里只有精品18| 亚洲av成人不卡在线观看播放网| 18禁观看日本| 久久精品夜夜夜夜夜久久蜜豆| 18禁观看日本| 美女免费视频网站| 男女那种视频在线观看| 别揉我奶头~嗯~啊~动态视频| 国产黄色小视频在线观看| 成年人黄色毛片网站| 狂野欧美白嫩少妇大欣赏| 美女高潮的动态| www国产在线视频色| 久久性视频一级片| 久久九九热精品免费| 夜夜躁狠狠躁天天躁| 精品久久久久久成人av| 午夜福利成人在线免费观看| 日韩欧美在线二视频| 国产精品 欧美亚洲| 可以在线观看毛片的网站| 欧美丝袜亚洲另类 | 精品久久久久久成人av| 后天国语完整版免费观看| 免费看美女性在线毛片视频| 天堂动漫精品| 免费大片18禁| 国产亚洲av嫩草精品影院| 国产视频内射| 怎么达到女性高潮| 麻豆成人av在线观看| 亚洲真实伦在线观看| 国产精品久久久av美女十八| 久久这里只有精品中国| 琪琪午夜伦伦电影理论片6080| 欧美性猛交╳xxx乱大交人| 制服丝袜大香蕉在线| 色播亚洲综合网| 高清毛片免费观看视频网站| 国产亚洲精品久久久com| www国产在线视频色| 久久久久九九精品影院| 长腿黑丝高跟| 欧美日本亚洲视频在线播放| av国产免费在线观看| 亚洲,欧美精品.| 我要搜黄色片| 十八禁人妻一区二区| 国产精品亚洲av一区麻豆| 久久国产精品人妻蜜桃| 成人永久免费在线观看视频| 免费在线观看影片大全网站| 日韩免费av在线播放| 丰满的人妻完整版| 热99在线观看视频| 亚洲电影在线观看av| 欧美乱码精品一区二区三区| 亚洲午夜精品一区,二区,三区| 亚洲av五月六月丁香网| 国产伦在线观看视频一区| 757午夜福利合集在线观看| 在线a可以看的网站| 在线观看免费视频日本深夜| 1000部很黄的大片| 日韩欧美 国产精品| 欧美+亚洲+日韩+国产| 99精品在免费线老司机午夜| 99精品久久久久人妻精品| 色综合亚洲欧美另类图片| 日本黄色视频三级网站网址| 婷婷亚洲欧美| 白带黄色成豆腐渣| 久久国产乱子伦精品免费另类| 亚洲18禁久久av| 99国产综合亚洲精品| 51午夜福利影视在线观看| 成人无遮挡网站| 最新在线观看一区二区三区| 国内精品久久久久精免费| 一个人看的www免费观看视频| 午夜两性在线视频| 五月伊人婷婷丁香| 天天添夜夜摸| 国产成人aa在线观看| 久久热在线av| 亚洲精品一卡2卡三卡4卡5卡| 成人三级黄色视频| 成人特级av手机在线观看| 国产精华一区二区三区| 色尼玛亚洲综合影院| 亚洲色图av天堂| 亚洲国产看品久久| 日本免费a在线| 黄色片一级片一级黄色片| 啦啦啦观看免费观看视频高清| 亚洲国产欧美网| 亚洲专区国产一区二区| 国产午夜福利久久久久久| 亚洲自拍偷在线| 久久久久久久久久黄片| 亚洲七黄色美女视频| 日韩有码中文字幕| 婷婷丁香在线五月| 小蜜桃在线观看免费完整版高清| 黄色日韩在线| 国产av在哪里看| 婷婷六月久久综合丁香| 国产午夜精品久久久久久| 国产成人影院久久av| 成人鲁丝片一二三区免费| 窝窝影院91人妻| 欧美日韩黄片免| 草草在线视频免费看| 国产精品自产拍在线观看55亚洲| 亚洲黑人精品在线| 夜夜爽天天搞| 成人午夜高清在线视频| 两人在一起打扑克的视频| 国产欧美日韩精品一区二区| 三级国产精品欧美在线观看 | 可以在线观看的亚洲视频| 欧美一级a爱片免费观看看| 黑人巨大精品欧美一区二区mp4| 五月伊人婷婷丁香| 日本a在线网址| 免费看日本二区| 成人永久免费在线观看视频| av中文乱码字幕在线| 久久久精品欧美日韩精品| 在线观看免费午夜福利视频| 看免费av毛片| 国产亚洲av高清不卡| 国产精品久久久久久精品电影| 亚洲七黄色美女视频| 老司机福利观看| 动漫黄色视频在线观看| 色综合站精品国产| 国产在线精品亚洲第一网站| 日韩 欧美 亚洲 中文字幕| 1024手机看黄色片| 婷婷精品国产亚洲av| 久久精品综合一区二区三区| 欧美中文综合在线视频| 国产野战对白在线观看| av女优亚洲男人天堂 | 18禁裸乳无遮挡免费网站照片| 一级毛片高清免费大全| 日韩欧美三级三区| 欧洲精品卡2卡3卡4卡5卡区| 此物有八面人人有两片| 色播亚洲综合网| 看片在线看免费视频| 久久精品aⅴ一区二区三区四区| 给我免费播放毛片高清在线观看| 亚洲午夜精品一区,二区,三区| 床上黄色一级片| 亚洲欧美精品综合久久99| 免费在线观看影片大全网站| 日韩高清综合在线| 757午夜福利合集在线观看| 国产黄a三级三级三级人| 欧美极品一区二区三区四区| 午夜福利在线在线| 一本久久中文字幕| 国产黄片美女视频| 一个人免费在线观看电影 | 欧美性猛交╳xxx乱大交人| 最新美女视频免费是黄的| 婷婷精品国产亚洲av在线| www.www免费av| 国产午夜精品论理片| 男女床上黄色一级片免费看| 国产久久久一区二区三区| 久久精品人妻少妇| 亚洲色图 男人天堂 中文字幕| 日本熟妇午夜| 极品教师在线免费播放| 真人做人爱边吃奶动态| 深夜精品福利| 久久伊人香网站| 男人舔奶头视频| 99riav亚洲国产免费| 欧美日韩国产亚洲二区| 久久久久久大精品| 午夜精品一区二区三区免费看| 一a级毛片在线观看| 亚洲成人久久爱视频| 一二三四在线观看免费中文在| 露出奶头的视频| 级片在线观看| 国产精品香港三级国产av潘金莲| 久久久国产成人精品二区| 成人av一区二区三区在线看| 青草久久国产| 国产精品98久久久久久宅男小说| av在线天堂中文字幕| 亚洲欧美日韩高清在线视频| 国产欧美日韩一区二区三| 国产真实乱freesex| 18禁观看日本| 嫩草影视91久久| 久久久国产精品麻豆| 亚洲成av人片免费观看| 国产欧美日韩精品一区二区| 日本免费一区二区三区高清不卡| 村上凉子中文字幕在线| 国产一区二区在线av高清观看| 国产成人aa在线观看| 亚洲无线观看免费| 午夜福利欧美成人| 亚洲精华国产精华精| 一本精品99久久精品77| 免费在线观看影片大全网站| 国产精品日韩av在线免费观看| aaaaa片日本免费| 岛国在线观看网站| 国产一区二区在线观看日韩 | 中文字幕人妻丝袜一区二区| 操出白浆在线播放| 成人永久免费在线观看视频| avwww免费| 亚洲色图 男人天堂 中文字幕| 搡老熟女国产l中国老女人| 男女做爰动态图高潮gif福利片| 91在线精品国自产拍蜜月 | 亚洲中文日韩欧美视频| 午夜视频精品福利| 国产午夜精品久久久久久| 国产成人aa在线观看| 免费一级毛片在线播放高清视频| av国产免费在线观看| a在线观看视频网站| 黑人欧美特级aaaaaa片| 色哟哟哟哟哟哟| 久久久久久久午夜电影| 99久久成人亚洲精品观看| av国产免费在线观看| 亚洲欧美激情综合另类| 国产精品久久久久久亚洲av鲁大| 操出白浆在线播放| 精品电影一区二区在线| 国产精品美女特级片免费视频播放器 | 久久久水蜜桃国产精品网| 色尼玛亚洲综合影院| 18禁黄网站禁片免费观看直播| 精品一区二区三区四区五区乱码| 日韩精品青青久久久久久| 美女cb高潮喷水在线观看 | 国产真实乱freesex| 日韩中文字幕欧美一区二区| 亚洲aⅴ乱码一区二区在线播放| 欧美高清成人免费视频www| 国产欧美日韩精品亚洲av| 国产 一区 欧美 日韩| 在线观看一区二区三区| 亚洲男人的天堂狠狠| 天堂动漫精品| 国产高清三级在线| 午夜免费观看网址| 欧美大码av| 久久婷婷人人爽人人干人人爱| 特大巨黑吊av在线直播| 99re在线观看精品视频| 熟女少妇亚洲综合色aaa.| 久久人人精品亚洲av| 一进一出抽搐动态| 日韩欧美国产一区二区入口| 一进一出抽搐动态| 97碰自拍视频| 国产成人av教育| 日韩欧美 国产精品| 禁无遮挡网站| 精品国产亚洲在线| 国产亚洲av嫩草精品影院| 国产精品亚洲一级av第二区| 性色avwww在线观看| 亚洲欧美精品综合一区二区三区| 欧美激情在线99| 欧美乱色亚洲激情| 亚洲七黄色美女视频| 国产伦一二天堂av在线观看| 欧美一区二区国产精品久久精品| 国产欧美日韩精品一区二区| 成人特级av手机在线观看| 日韩大尺度精品在线看网址| 香蕉av资源在线| 亚洲国产精品久久男人天堂| 久久精品aⅴ一区二区三区四区| 美女扒开内裤让男人捅视频| 亚洲国产看品久久| 国产久久久一区二区三区| 男女下面进入的视频免费午夜| 色视频www国产| 一进一出抽搐gif免费好疼| 欧美乱妇无乱码| 国产午夜福利久久久久久| 桃红色精品国产亚洲av| 九九在线视频观看精品| 一本久久中文字幕| 999精品在线视频| 久久久色成人| av福利片在线观看| 日本 欧美在线| 免费无遮挡裸体视频| 18禁美女被吸乳视频| 日韩欧美在线乱码| 国产精品自产拍在线观看55亚洲| 国产精品99久久久久久久久| 国产成人精品久久二区二区91| 两性午夜刺激爽爽歪歪视频在线观看| 欧美日韩乱码在线| 曰老女人黄片| 性色avwww在线观看| 欧美日韩综合久久久久久 | 亚洲天堂国产精品一区在线| 五月玫瑰六月丁香| 亚洲天堂国产精品一区在线| 精品乱码久久久久久99久播| 国内精品久久久久精免费| 亚洲熟妇中文字幕五十中出| 波多野结衣巨乳人妻| 久久精品综合一区二区三区| 9191精品国产免费久久| 一夜夜www| 国内精品一区二区在线观看| 精品久久久久久久末码| 亚洲精品在线美女| 一级作爱视频免费观看| 男人舔奶头视频| 欧美色欧美亚洲另类二区| 日韩精品青青久久久久久| 久久久久性生活片| 午夜免费成人在线视频| 欧美一级毛片孕妇| 在线永久观看黄色视频| 高潮久久久久久久久久久不卡| 五月玫瑰六月丁香| 亚洲中文日韩欧美视频| 欧美日韩乱码在线| 一个人看的www免费观看视频| 中国美女看黄片| 亚洲中文字幕日韩| 久99久视频精品免费| 亚洲av成人不卡在线观看播放网| 美女高潮的动态| 久久国产乱子伦精品免费另类| 99热6这里只有精品| 亚洲精品一区av在线观看| 美女高潮的动态| 国产真人三级小视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 91麻豆精品激情在线观看国产| 欧美一级毛片孕妇| 国产欧美日韩一区二区三| 国产在线精品亚洲第一网站| 十八禁网站免费在线| 国产伦精品一区二区三区四那| 综合色av麻豆| 成人高潮视频无遮挡免费网站| 母亲3免费完整高清在线观看| 亚洲欧美精品综合一区二区三区| 欧美成狂野欧美在线观看| 久久久久久人人人人人| 日本撒尿小便嘘嘘汇集6| 美女大奶头视频| 欧美黄色淫秽网站| 人人妻,人人澡人人爽秒播| 亚洲欧美日韩高清专用| 国产精品野战在线观看| av欧美777| 两人在一起打扑克的视频| 日本黄色视频三级网站网址| 日韩av在线大香蕉| 免费看光身美女| 97超视频在线观看视频| 色视频www国产| 免费在线观看影片大全网站| 午夜福利在线观看吧| 亚洲专区字幕在线| 九色成人免费人妻av| 国产综合懂色| 国产精品av视频在线免费观看| 久久草成人影院| 午夜免费成人在线视频| 亚洲成人久久爱视频| 久久精品91蜜桃| www.www免费av| 老熟妇乱子伦视频在线观看| 久久久国产成人免费| 一二三四在线观看免费中文在| 性色avwww在线观看| 1000部很黄的大片| 亚洲人成伊人成综合网2020| 97人妻精品一区二区三区麻豆| 亚洲男人的天堂狠狠| 97超视频在线观看视频| 99精品在免费线老司机午夜| 黄色视频,在线免费观看| 欧美乱码精品一区二区三区| 欧美在线一区亚洲| 国内精品美女久久久久久| 精品一区二区三区四区五区乱码| 色噜噜av男人的天堂激情| 久久午夜综合久久蜜桃| 午夜福利免费观看在线| 激情在线观看视频在线高清| 午夜福利视频1000在线观看| 久久久久亚洲av毛片大全| 色综合欧美亚洲国产小说| 色吧在线观看| 亚洲av成人一区二区三| 一个人观看的视频www高清免费观看 | 99re在线观看精品视频| 亚洲在线自拍视频| 99国产综合亚洲精品| 国产精品99久久99久久久不卡| 男人和女人高潮做爰伦理| 成人国产一区最新在线观看| 校园春色视频在线观看| 国产成年人精品一区二区| 亚洲午夜理论影院| 中亚洲国语对白在线视频| 999精品在线视频| 成年免费大片在线观看| 国产欧美日韩一区二区三| 国产av麻豆久久久久久久| 一夜夜www| 日韩中文字幕欧美一区二区| 免费在线观看成人毛片| 日本 av在线| 脱女人内裤的视频| 搡老妇女老女人老熟妇| 成年女人看的毛片在线观看| 两人在一起打扑克的视频| 久久婷婷人人爽人人干人人爱| 脱女人内裤的视频| 欧美成人性av电影在线观看| 999精品在线视频| 成人鲁丝片一二三区免费| 亚洲激情在线av| 中文字幕av在线有码专区| 亚洲国产精品成人综合色| 亚洲精品色激情综合| netflix在线观看网站| 国产精品 国内视频| 老熟妇仑乱视频hdxx| 性欧美人与动物交配| 香蕉av资源在线| 日韩欧美在线乱码| 欧美xxxx黑人xx丫x性爽| 男人舔女人的私密视频| 法律面前人人平等表现在哪些方面| 长腿黑丝高跟| 99精品在免费线老司机午夜| 亚洲国产精品久久男人天堂| 国内少妇人妻偷人精品xxx网站 | 国产高清videossex| 亚洲美女视频黄频| 麻豆av在线久日| 桃色一区二区三区在线观看| 亚洲 国产 在线| x7x7x7水蜜桃| 亚洲 国产 在线| 精品无人区乱码1区二区| 97碰自拍视频| 国产成人一区二区三区免费视频网站| 国产精品影院久久| e午夜精品久久久久久久| 国产三级在线视频| 成人三级做爰电影| 成人三级黄色视频| 国产v大片淫在线免费观看| 亚洲国产精品sss在线观看| 婷婷丁香在线五月| 男女下面进入的视频免费午夜| 国产成人啪精品午夜网站| 热99在线观看视频| 校园春色视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 熟女人妻精品中文字幕| 亚洲专区中文字幕在线| 狂野欧美激情性xxxx| 午夜免费成人在线视频| 在线观看日韩欧美| 每晚都被弄得嗷嗷叫到高潮| 午夜精品在线福利| 成年女人看的毛片在线观看| 免费观看的影片在线观看| 欧美日韩黄片免| 午夜精品久久久久久毛片777| 好男人电影高清在线观看| 99久久国产精品久久久| 欧美又色又爽又黄视频| 久久久国产成人精品二区| 亚洲国产看品久久| 久久久久精品国产欧美久久久| 美女高潮喷水抽搐中文字幕| 99精品在免费线老司机午夜| a级毛片a级免费在线| 搞女人的毛片| 久久亚洲真实| 日本黄色视频三级网站网址| 日韩有码中文字幕| 国产男靠女视频免费网站| 久久精品国产99精品国产亚洲性色| avwww免费| 久久热在线av| 亚洲国产欧美人成| 国产伦在线观看视频一区| 久久久久久久精品吃奶| 成人特级黄色片久久久久久久| 欧美日韩福利视频一区二区| 亚洲真实伦在线观看| 一卡2卡三卡四卡精品乱码亚洲| 丰满的人妻完整版| 怎么达到女性高潮| 成在线人永久免费视频| 精品国产亚洲在线| 欧美极品一区二区三区四区| a级毛片a级免费在线| 亚洲电影在线观看av| 精品日产1卡2卡| 久久香蕉国产精品| 黑人巨大精品欧美一区二区mp4| 欧美黄色淫秽网站| 国产成人福利小说| 三级男女做爰猛烈吃奶摸视频| 欧美3d第一页| 老汉色av国产亚洲站长工具| 亚洲最大成人中文| 午夜激情欧美在线| 亚洲av免费在线观看| 国产精品一区二区三区四区免费观看 | 欧美日韩乱码在线| 午夜福利高清视频| 中文在线观看免费www的网站| 午夜亚洲福利在线播放| 18禁观看日本| 久久国产精品人妻蜜桃| 国产成人精品久久二区二区91| 日韩人妻高清精品专区| 久久久久精品国产欧美久久久| 黄频高清免费视频| 久久精品人妻少妇| 色尼玛亚洲综合影院| 国产精品女同一区二区软件 | 精品日产1卡2卡| 国产一区二区激情短视频| 一级黄色大片毛片| 91麻豆av在线| 国产精品爽爽va在线观看网站| 亚洲中文字幕一区二区三区有码在线看 | 波多野结衣高清无吗|