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

    Effect of carbon dioxide concentration on the combustion characteristics of boron agglomerates in oxygen-containing atmospheres

    2023-11-18 08:41:24LianDUANZhixunXIAYunchaoFENGBinbinCHENJiaruiZHANGLikunMA

    Lian DUAN, Zhixun XIA, Yunchao FENG, Binbin CHEN, Jiarui ZHANG, Likun MA

    Research Article

    Effect of carbon dioxide concentration on the combustion characteristics of boron agglomerates in oxygen-containing atmospheres

    Lian DUAN, Zhixun XIA, Yunchao FENG, Binbin CHEN, Jiarui ZHANG, Likun MA

    College of Aerospace Science and Engineering, National University of Defense Technology, Changsha 410073, China

    In ramjet combustion chambers, carbon dioxide (CO2) produced by the combustion of carbonaceous fuel enters the chamber together with boron agglomerates. In order to investigate the effect of CO2concentration present in an oxygen-containing atmosphere on the combustion characteristics and oxidation mechanisms of boron agglomerates, we used a laser ignition system, an X-ray diffractometer (XRD), and a thermogravimetric-differential scanning calorimetry (TG-DSC) combined thermal analysis system. Single-particle boron was tested in the laser-ignition experiments as the control group. The ignition experiment results showed that with a fixed O2concentration of 20%, when the particle temperature reaches the melting point of boron, increasing CO2content causes the combustion process of boron agglomerates to transition from single-particle molten droplet combustion to porous-particle combustion. Furthermore, XRD analysis results indicated that the condensed-phase combustion products (CCPs) of boron particles in a mixed atmosphere of O2and CO2contained B4C, which is responsible for the porous structure of the particles. At temperatures below 1200 ℃, the addition of CO2has no obvious promotion effect on boron exothermic reaction. However, in the laser-ignition experiment, when the oxygen concentration was fixed at 20% while the CO2concentration increased from 0% to 80%, the maximum temperature of boron agglomerates rose from 2434 to 2573 K, the self-sustaining combustion time of single-particle boron decreased from 396 to 169 ms, and the self-sustaining combustion time of boron agglomerates decreased from 198 to 40 ms. This study conclusively showed that adding CO2to an oxygen-containing atmosphere facilitates boron reaction and consumption pathways, which is beneficial to promoting exothermic reaction of boron agglomerates at relatively high temperatures.

    Boron combustion; Amorphous boron; Boron-containing propellant; Solid fuel ramjet

    1 Introduction

    Boron is regarded as an ideal energetic additive for solid propellants because of its high volumetric and gravimetric heating values (Fry, 2004), so it is widely used in solid-fuel ramjets and scramjets (Lv et al., 2017; Chen et al., 2018; Zhang et al., 2020). In a ramjet gas generator, boron particles are likely to agglomerate into large particles due to low melting point and high viscosity of the oxide layer on the surface of boron particles (Meerov et al., 2015; Liu et al., 2017). The combustion of the carbonaceous formulation in a boron-based fuel-rich propellant produces a large amount of CO2(Liu et al., 2015), which affects the combustion characteristics of boron agglomerates. In addition, direct combustion heating is often used in ramjet ground simulation tests, to increase the air temperature during simulation of the total temperatures of air flows entering the secondary combustion chamber in flight state (Roux et al., 2014). When the working medium of the heater uses an air/alcohol or air/kerosene mixture, the level of CO2in the air entering the combustion chamber exceeds 10% (Li et al., 2007).

    Researchers have performed a considerable amount of work on the effects of fluorine-containing atmosphere (Krier et al., 1996; Zhou et al., 1998; Ulas et al., 2001) and water vapor (Smolanoff et al., 1996; Foelsche et al., 1999; Yoshida and Yuasa, 2000) on the ignition and combustion characteristics of boron particles, and improved the understanding of the energy-release characteristics of such particles. However, few studies have addressed the ignition and combustion characteristics of boron in a CO2-containing atmosphere. Experiments conducted by DiGiuseppe and Davidovits (1981) revealed that the reaction between boron and CO2produced BO and CO, with an exothermic heat of (67±5) kcal/mol (1 kcal=4.186 kJ) and a relatively slow reaction rate. Burkholder et al. (1993) conducted a condensation reaction of pulsed laser-evaporated boron atoms with CO2in excess Ar, and infrared spectrum analysis of the reaction products showed that OBCO was generated in addition to BO and CO. Thus, they concluded that OBCO was the main product of the reaction between boron and CO2. Chin et al. (2003) calculated the reaction between B and CO2by ab initio and density functional methods, and found that the reaction first generated a dynamically unstable intermediate, BOCO, which then decomposed into BO and CO, and finally formed the most stable thermodynamic substance: OBCO. Li et al. (2014) analyzed the effect of CO2concentration on the thermal oxidation characteristics of boron particles by thermogravimetry, and found that appropriate concentrations of CO2(10% or 30%) in a limiting O2(10%) atmosphere facilitated the onset of oxidation of boron particles and improved ignition performance.

    A CO2atmosphere can affect the thermal oxidation characteristics of boron particles, but the specific reasons for this effect are yet to be elucidated. Especially about the effects of boron agglomerates, little is known. In this study, we directly observed the structural and morphological evolution of boron agglomerates during combustion, and analyzed the combustion process and spectral data. The temperature evolution of the particle surface was obtained through a fitting procedure using the Planck law. The microscopic morphology and rough elemental composition of condensed-phase combustion products (CCPs) were analyzed with a scanning electron microscope (SEM) and energy dispersive spectrometer (EDS), and the phase-crystal structures were analyzed by X-ray diffractometer (XRD). Additionally, we employed thermogravimetric-differential scanning calorimetry (TG-DSC) to study the oxidation process. The results of this study provide a reference for the energy-release mechanism of boron agglomerates when O2and CO2coexist.

    2 Experimental procedures

    2.1 Materials

    We carried out our comparative experiments on amorphous boron agglomerates and single-particle crystalline boron. The boron agglomerates were prepared by drying a slurry of the micron-sized amorphous boron in water (Mi et al., 2013). The initial amorphous boron was available commercially from Yingkou Liaobin Fine Chemicals Co., Ltd., China, with a purity of 96% and mean diameter of 2.55 μm. The agglomerate density was measured by the Archimedes drainage method to be 1.27 g/cm3. The single-particle crystalline boron was obtained from Zhongnuo Advanced Material Technology Co., Ltd., China, with a purity of 99.9%. The boron agglomerates and single-particle boron between 355 and 400 μm were screened by a sieve shaker. The samples used in the test are shown in Fig. 1.

    Fig. 1 Amorphous boron agglomerate (a) and single-particle crystalline boron (b)

    2.2 Apparatus and methods

    A laser-ignition test system was used to investigate the effect of CO2concentration on the combustion characteristics of boron agglomerates in oxygen-containing atmospheres. A schematic diagram of the experimental setup is given in Fig. 2. The same system was successfully used in an earlier ignition study on boron agglomerates (Duan et al., 2022). All combustion tests were conducted under constant laser power and laser duration, and the laser loading time was set to 250 ms. We measured the actual power emitted to the sample surface in different atmospheres with a power meter (F150A-BB-26, Ophir, Israel) as (85.0±0.3) W.

    Fig. 2 Schematic diagram of the laser-ignition experimental system

    A monochrome high-speed camera (Dimax HS4, PCO, Germany), a color high-speed camera (FASTCAM SA-X2, Photron, Japan), and a fiber-optic spectrometer (HR2000+ES, Oceanoptics, USA) were employed to monitor the combustion process of the sample through high-purity quartz-glass windows. Flame images were photographed using the color high-speed camera, with an exposure time of 200 μs and a frame rate of 5000 frame/s, while the particle-surface morphology evolution was recorded with the monochrome high-speed camera, with an exposure time of 900 μs and a frame rate of 1000 frame/s. The fiber spectrometer measured the spectrum within the wavelength range of 200?–?1100 nm, while the integration time was fixed at 1 ms. Before each experiment, the spectrometer recorded the background radiation spectrum. During the experiment, the data collected by the spectrometer automatically deducted the background radiation spectrum.

    In all ignition tests, the burner pressure was kept at 1.013×105Pa, the total gas-flow rate was precisely controlled by a mass flowmeter to be constant at 10 mL/min, the flow rate of O2was fixed at 2 mL/min, and the flow ratios of Ar and CO2were varied to set the concentrations of CO2in the burner as 0%, 20%, 40%, 60%, or 80%. In addition, the operating conditions of a pure CO2atmosphere were set as the control. At least 10 repetitive experiments were performed under each condition to ensure the reliability of the experimental results.

    After the ignition test, we analyzed the microscopic morphology and rough elemental composition of the CCPs with an SEM (JEOL, JSM-7900F, Japan) combined with EDS. The phase crystal structures of CCPs were analyzed with an XRD instrument (D8 Advance, Bruker, Germany).

    For thermo-analytical measurements, we used a Netzsch STA 449F3 simultaneous thermal analyzer. A total of 8.5?–?10.0 mg boron agglomerates were packed in alumina crucibles. Three flowmeters controlled the total flow rate into the furnace of the analyzer to be 100 mL/min, and different atmospheres were obtained by changing the flow ratios of O2, CO2, and Ar. The test atmospheres were set to 20% O2+80% Ar, 20% O2+40% CO2+40% Ar, 20% O2+80% CO2, and pure CO2. A well-proportioned gas flow was continuously fed into the furnace for 10 min before each programmed temperature rise to ensure that the gas environment was well controlled. In addition, the test data of an empty crucible in various atmospheres were obtained for baseline correction.

    3 Results and discussion

    3.1 Surface-morphology analysis

    Figs. 3a?–?3c illustrate the temporal evolution of the surface morphology of boron agglomerates when the oxygen concentrations were kept at 20% while the CO2concentration was 0%, 40%, and 80%, respectively. The image parameters such as contrast, brightness, and other settings were the same. We recorded the moment when the laser started to radiate to the particle as 0 ms. It can be seen from Fig. 3a that, in the oxygen-containing atmosphere without CO2, the particles gradually melted into a liquid state. After 250 ms, the laser stopped loading, and the particles underwent self-sustained combustion for a while. This was similar to the process observed in our previous study (Duan et al., 2022). As can be seen from Fig. 3b, in the oxygen-containing atmosphere with 40% CO2concentration, no melting trace of particles could be observed in the first 150 ms of laser loading. At 200 ms, partial reflection appeared on the particle surface, and it changed into a semi-molten state at 250 ms.

    Fig. 3 Surface morphologies of boron agglomerates with different CO2 concentrations: (a) 20% O2+80% Ar; (b) 20% O2+40% CO2+40% Ar; (c) 20% O2+80% CO2

    Fig. 3c shows that in the oxygen-containing atmosphere with 80% CO2concentration, the particle-combustion process did not exhibit a melting state, and there was no obvious liquid film covering the particle surface. The upper surface of the particle collapsed downward and the surface morphology was irregular. At 300 ms, the particles had cooled and the surface had dimmed. Comparing the images shown in Figs. 3a?–?3c, we found that with a fixed oxygen concentration of 20%, the higher the CO2concentration, the lower the melting degree of the particles, the more irregular the particle morphology, and the faster the particle-cooling rate after the laser had stopped loading.

    The light source was placed opposite to the camera as a background light so that we could observe the morphological changes in the particle profile during combustion. Figs. 4a and 4b show the combustion processes of single-particle crystalline boron in atmospheres of 20% O2+80% Ar and 20% O2+80% CO2, respectively. As seen from Fig. 4a, in the mixed atmosphere of oxygen and argon, the contour lines of single-particle boron gradually became smooth, and the particles gradually melted into a hemispherical shape. Fig. 4b shows that in the oxygen-containing atmosphere with 80% CO2concentration, the smoke overflowed from the particle surface to the environment at about 50 ms, and a large amount of smoke was generated around the crystalline boron particles during the subsequent combustion process. After flameout, the thickness of deposition products accumulated around the particles was greater.

    Fig. 4 Particle-profile evolution of single-particle boron: (a) 20% O2+80% Ar; (b) 20% O2+80% CO2

    3.2 Analysis of the combustion process

    The ignition and combustion processes of boron agglomerates in atmospheres with different CO2concentrations are illustrated in Fig. 5. The time of laser initiation (0 ms) was used as the reference time. After laser loading, a faint green flame appeared above the particles in both the 20% O2+80% Ar atmosphere and 20% O2+80% CO2atmosphere, caused by light emission from the gaseous intermediate BO2(Yuasa and Isoda, 1991).

    Fig. 5 Combustion process of boron agglomerates in different atmospheres: (a) 20% O2+80% Ar; (b) 20% O2+80% CO2. References to color refer to the online version of this article

    Comparing Figs. 5a and 5b, it becomes clear that in the 20% O2+80% Ar atmosphere, there was an ejection phenomenon during the combustion process of boron agglomerates, as shown in Fig. 5a at 248.8 ms. However, the boron agglomerates did not melt into a hemispherical shape in an atmosphere of 20% O2+80% CO2, and the upper surface of the particles was highly irregular during combustion. We speculated that this was due to the reaction of boron and CO2generating other substances in the CO2-rich atmosphere, which did not melt easily or suppressed the exotherm of the reaction system. After 250 ms, the laser stopped loading, and the particles were in self-sustained combustion for a while. Then, the reaction rate gradually decreased and the flame was extinguished. Fig. 5 shows that the flameout rate was faster in an atmosphere containing CO2.

    3.3 Combustion-spectrum analysis

    Fig. 6 illustrates the combustion spectra of boron agglomerates at the moment of maximum spectral intensity in a pure CO2atmosphere, pure Ar atmosphere, and 20% O2+80% CO2atmosphere. During the combustion process of boron agglomerates in the oxygen-containing atmosphere, the spectral signals of gas-phase intermediates BO and BO2were emitted, causing the wavelength curves to contain multiple distinct characteristic peaks (Yoshida and Yuasa, 2000; Song et al., 2021). In a pure CO2atmosphere, there were no emission peaks for BO and BO2in the spectral curve, indicating that the reaction between boron and CO2did not produce the gas-phase intermediate products BO and BO2. In a pure Ar atmosphere, the particle temperature was raised by laser radiation, and spectral fluctuation was detected by the spectrometer. However, the spectral intensity was relatively low because boron particles did not undergo exothermic reaction in the Ar atmosphere.

    Fig. 6 Combustion spectra with the maximum intensity for boron agglomerates in different atmospheres

    Both the ignition and combustion processes of boron particles in the oxygen-containing atmosphere emitted the spectral signal of BO2, and the strongest spectrum intensity corresponded to a wavelength of 547.6 nm. Fig. 7 shows that after laser loading, the spectral intensity signals increased rapidly for atmospheres of 20% O2+80% CO2and 20% O2+80% Ar, and then decreased slowly, due to the surface reaction and phase transition during the combustion process changing the emissivity of the particle surface, as well as the solidified B2O3around the particle absorbing some of the radiation (Yuasa and Isoda, 1991; Millot et al., 2002). Curve B has two troughs at about 50 and 190 ms, which may be related to other substances produced by the reaction of boron and CO2. After the laser has been continuously loaded for 250 ms and then deactivated, the spectral intensity of curve B dropped sharply.

    Fig. 7 Temporal variation of spectral intensity at a wavelength of 547.6 nm

    The duration of availability of the 547.6 nm spectral signal was defined as the combustion time (Song et al., 2021). The laser turn-on time was recorded as 0 ms. The combustion time (c) can be further divided into combustion time with laser loading (lc) and self-sustaining combustion time after laser loading (sc). Under these experimental conditions, the longer the particle-combustion time, the slower the decrease in chemical reaction rate after laser deactivation, and thus the self-sustaining combustion lasted longer.

    Fig. 8 shows the self-sustaining combustion time of single-particle crystalline boron and amorphous boron agglomerates at different CO2concentrations when the oxygen concentration was 20%. The self-sustained combustion time of both single-particle boron and boron agglomerates decreased significantly with increasing CO2concentration. When the CO2concentration was increased from 0% to 80%, the self-sustaining combustion time of boron agglomerates decreased from 198 to 41 ms, and that of single-particle boron decreased from 397 to 170 ms. The higher the CO2concentration, the faster the chemical reaction rate decreased, and therefore the shorter the self-sustaining combustion time became. This suggests that when boron is burned in an oxygen-containing atmosphere, the higher the CO2concentration, the more sensitive the chemical reaction rate is to the temperature changes.

    Fig. 8 Self-sustaining combustion time of single-particle boron and boron agglomerates with different CO2 concentrations

    3.4 Temperature analysis

    The temperature evolution of the particle surface was obtained through a fitting procedure using the Planck law (Mi et al., 2013). The calculation method and error analysis are illustrated in the electronic supplementary materials (ESM).

    Fig. 9 shows the temporal variation in particle temperature when the oxygen concentration was fixed at 20% and the CO2concentrations were 0%, 40%, and 80%. The time interval of temperature data was 10 ms. After laser activation, the particle temperature in all three atmospheres rapidly increased to above the boiling point of B2O3(2316 K), so the boron oxide layer on the particle surface could be quickly removed by evaporation. Within 250 ms of laser loading, when the CO2concentration was increased from 0% to 80%, the maximum surface temperature of the boron agglomerates rose from 2434 to 2573 K, which is between the melting points of boron (2350 K) and B4C (2623 K). After 250 ms, the particle temperature in the 20% O2+80% CO2atmosphere decreased rapidly, while in the 20% O2+80% Ar atmosphere, the high temperature was maintained for a period.

    Fig. 9 Temperature variation of the boron agglomerates with different CO2 concentrations

    Fig. 10 FE-SEM images of CCPs: (a) CCP of boron agglomerate in 20% O2+80% CO2 atmosphere; (b) CCP of single-particle boron in 20% O2+80% CO2 atmosphere

    3.5 Analysis of condensed combustion products

    It is helpful to explore the combustion process of boron agglomerates in typical atmospheres by analyzing the residual CCPs after particle flameout. Field emission scanning electron microscope (FE-SEM) images of the CCPs of amorphous boron agglomerates in the 20% O2+80% CO2atmosphere, as well as single-particle crystalline boron in the 20% O2+80% CO2atmosphere, are shown in Figs. 10a and 10b, respectively. Fig. 10a shows that the structure was curly and the surface morphology was very rough; they lost their initial granular structure and contained a large number of pores with obvious wrinkles. The rough surface increased the contact area between the oxidizing gas and the particle surface. The clear crystal structure was discernible by magnifying the particle surface. The surface of CCPs of single-particle boron lost its initial flat morphology and also became rough, with an obvious crystal structure (Fig. 10b). We have previously observed the microscopic morphology of the CCPs of boron agglomerates in pure oxygen and mixed atmospheres of oxygen and argon (Duan et al., 2022), and found that there were no wrinkles or similar crystal structures on the particle surface, which is significantly different from the morphologies observed here.

    Subsequently, we analyzed the micro-area elementary composition of the CCPs of single-particle crystalline boron in an atmosphere of 20% O2+80% CO2. The results of this analysis are shown in Fig. 11. We observed that the distribution of C was dense in the wrinkled regions, while the distribution of O was dense in the smooth regions. Elemental composition analysis of the wrinkled regions in Figs. 10a and 10b revealed that the predominant elements were B and C.

    Fig. 11 Reference image of SEM and EDS elemental maps corresponding to elements C (green points), O (blue points), and B (red points). References to color refer to the online version of this figure

    To explore the crystal structures formed by the reaction of boron and CO2in oxygen-containing atmospheres, we placed the compacted amorphous boron stack in an atmosphere with an oxygen concentration of 20% and a CO2concentration of 80%, for laser ignition. The power level and duration of the laser were consistent with the above ignition experiments. Then, the CCPs were taken for XRD analysis. The obtained XRD pattern is shown in Fig. 12. The intense peaks at around 2=20°, 22°, 23°, 35°, and 38° corresponded to boron carbide (B4C) (Jain and Anthonysamy, 2015), indicating that boron can react with CO2to form B4C in an oxygen-containing atmosphere. Additionally, peaks for B2O3were identified at around 2=14° and 28° (Hashim et al., 2021), while few low-intensity crystal boron peaks were observed due to the incomplete reaction of the tested sample (Sun et al., 2018).

    Fig. 12 XRD pattern of CCPs in 20% O2+80% CO2 atmosphere

    3.6 Low-temperature exothermic reaction

    TG and DSC traces for amorphous boron in 20% O2+80% Ar, 20% O2+40% CO2+40% Ar, 20% O2+80% CO2, and pure CO2atmospheres at a heating rate of 20 ℃/min are shown in Fig. 13.

    Fig. 13 TG and DSC traces of amorphous boron in different atmospheres heated at 20 ℃/min to a maximum temperature of 1200 ℃

    Oxidation of boron began at around 550 ℃ in all atmospheres. In an atmosphere with a fixed oxygen concentration of 20%, the oxidation rate increased with the increase of CO2content after reaching approximately 750 ℃. We speculate that after the oxide layer on the outer surface of the boron particles melts, O2and CO2diffuse simultaneously through the liquid oxide layer and react with the boron core. Therefore, the higher the CO2content, the greater the mass gain of the sample and the more advanced the exothermic peak. At temperatures above 850 ℃, the mass gain of the boron particles tended to be gentle. This is likely due to the accumulation of reaction products on the particle surface, which thickens the liquid oxide layer and inhibits the reaction between the gases and the boron core.

    In an atmosphere of pure CO2, the oxidation process continued and eventually surpassed the mass gain of boron in O2-CO2mixtures. However, the heat release was comparatively low. This is most likely because boron reacts with CO2to form BOCO during this period, with an exothermic heat of about 25 kcal/mol (Chin et al., 2003), which is far lower than the exothermic heat of 146 kcal/mol of the reaction between boron and oxygen (King, 1973).

    3.7 Discussion

    The above experimental results demonstrated that B4C in a ramjet combustion chamber can form not only by the reaction of boron with hydrocarbons generated by the decomposition of carbon-based compounds (such as hydroxyl terminated polybutadiene in propellant (Liu et al., 2015)) but also by the reaction of boron with atmospheric CO2at high temperature. B4C can further react with oxygen to form the thermodynamically stable compound B2O3(Jain and Anthonysamy, 2015). B2O3is the main product of B4C combustion (Li and Qiu, 2007), and the relevant reaction is expressed by Eq. (1). During the continuous loading of the laser, the surface temperature of particles was higher than the boiling point of B2O3(2316 K), and the B2O3evaporated from the surface and condensed around the particles to form obvious smoke, as shown in Fig. 4b. Our conjecture is that at high temperature, CO2in an oxygen-containing atmosphere increases the combustion efficiency of boron particles, due to the continuous generation and consumption of B4C.

    The calorific value of B4C is 52 MJ/kg, which is similar to that of boron (59 MJ/kg). However, the melting point of B4C is about 300 K higher than that of boron; hence, it will be relatively difficult to melt into liquid phase. The temperature in this experiment reached the melting point of boron but not that of B4C. The latter adhered to molten boron and changed the physicochemical properties and combustion state of the particles.

    Based on the above analysis, a schematic diagram of the structural evolution of boron agglomerates during combustion is shown in Fig. 14. When the agglomerated boron is heated, individual boron particles tend to bond together due to the low melting point and high surface tension of their oxide layer, resulting in the volume of the agglomerate shrinking, as shown in Fig. 14b. At a fixed O2concentration of 20%, as the CO2concentration increased from 0% to 80%, the maximum surface temperature of the boron agglomerates rose from 2434 to 2573 K, which was between the melting points of B and B4C. In an atmosphere with low or no CO2content, after the agglomerate's temperature reaches the melting point of boron, it gradually melts into droplets, as shown in Fig. 14c. These droplets lose their original agglomeration morphology, and the oxidizing gases cannot diffuse inside through the pores. Droplets with larger particle sizes tend to have longer combustion times in practical application, resulting in incomplete combustion (Sun et al., 2020). In an atmosphere with a CO2content of more than 40%, the amount of boron carbide produced is greater than the amount consumed by its reaction with oxygen, causing it to dominate the agglomerate and further transforming it into a porous structure, as shown in Fig. 14d. Gases such as O2and CO2can diffuse into the interior of the agglomerate through the pores, leading to an accelerated rate of boron consumption. When the CO2concentration is between 20% and 40%, the agglomerate structure is in a semi-molten state during combustion, and the local area of the particle surface may be rough due to the reaction with CO2.

    Fig. 14 Schematic of the internal structural evolution of boron agglomerates during heating and combustion in atmospheres with different CO2 concentrations: (a) original boron agglomerate; (b) shrunken agglomerate; (c) molten droplet; (d) porous particle

    After the laser stopped radiation heating, we noted that the higher the concentration of CO2became, the shorter the self-sustaining combustion time of single-particle boron and boron agglomerates was. With a fixed oxygen concentration of 20%, when CO2concentration was increased from 0% to 80%, the self-sustaining combustion time of single-particle boron decreased from 396 to 169 ms, and that of boron agglomerates decreased from 198 to 40 ms. Therefore, it seems that the higher the CO2concentration, the greater the content of generated B4C, and the higher the content of B2O3generated by the reaction between B4C and oxygen. After the laser stopped loading, in the absence of continuous heating from any external heat source, the temperature gradually decreased. When the particle surface temperature was lower than the boiling point of B2O3(2316 K), B2O3generated by surface reaction accumulated on the surface of B4C, forming a glassy liquid film that acted as a barrier to the further diffusion of oxygen into B4C and to the release of CO2from B4C (Li and Qiu, 2007), such that the reaction rate decreased rapidly and the flame was quickly extinguished. In the atmosphere without CO2, the self-sustaining combustion time of particles was longer. This is because the B2O3produced by the reaction of boron and oxygen was formed by further oxidation of the intermediate products in the gas phase (Burkholder and Andrews, 1991; Zhou, 1998). The relevant reaction pathway is presented in Eqs. (2)?–?(4) (Burkholder and Andrews, 1991; Liang et al., 2017). Therefore, when the particle temperature drops below the boiling point of boron oxide, the B2O3in the environment will not completely diffuse and cover the particle surface, and the particle surface can still be in direct contact with oxygen in the environment for an exothermic reaction, thereby self-sustaining combustion for a relatively long time.

    To the best of our knowledge, the mechanism model of boron ignition and combustion established by Yetter et al. (1991), Zhou (1998), and Zhou et al. (1999), based on the detailed reaction mechanism and molecular dynamics theory, did not consider the reaction process after boron reacts with CO2to generate CO and BO. Therefore, the path and thermodynamic parameters of the reaction on the particle surface that generates B4C need to be further explored.

    4 Conclusions

    In this work, we investigated the effect of CO2concentration on the combustion characteristics of amorphous boron agglomerates by changing the flow ratios of CO2to Ar while maintaining a constant oxygen concentration of 20%. The main findings are summarized as follows:

    1. Boron agglomerates can react with CO2in an oxygen-containing atmosphere at high temperature to form B4C, which increases the reaction and consumption pathways of boron and changes the physicochemical properties and combustion state of particles. With a fixed oxygen concentration of 20%, when particle temperatures reach the melting point of boron, an increase in CO2content causes the combustion process of boron agglomerates to transition from single-particle molten droplet combustion to porous-particle combustion.

    2. At temperatures below 1200 ℃, the addition of CO2has no obvious promotion effect on a boron exothermic reaction. However, in the laser ignition experiment, the maximum surface temperature of boron agglomerates increased from 2434 to 2573 K when the oxygen concentration was fixed at 20% and the CO2concentration increased from 0% to 80%. Therefore, adding CO2to an oxygen-containing atmosphere is beneficial to promoting the exothermic reaction of boron agglomerates at relatively high temperatures.

    3. With a fixed oxygen concentration of 20%, when the CO2concentration was increased from 0% to 80%, the self-sustaining combustion time of single-particle boron decreased from 396 to 169 ms, and that of boron agglomerates decreased from 198 to 40 ms. The higher the CO2concentration, the more quickly the chemical reaction rate decreases, and therefore the shorter the self-sustaining combustion time will be. This is presumably due to the fact that when the particle temperature is lower than the boiling point of B2O3, B2O3produced by B4C oxidation forms a glassy liquid film that covers the particle surface and acts as a barrier to hinder the further combustion of particles.

    Acknowledgments

    This work is supported by the National Natural Science Foundation of China (No. 52006240) and the Hunan Provincial Natural Science Foundation of China (Nos. 2020JJ4665 and 2021JJ30775).

    Author contributions

    Lian DUAN designed the research and wrote the first draft of the manuscript. Yunchao FENG and Binbin CHEN processed the corresponding data. Likun MA and Jiarui ZHANG helped to organize the manuscript. Zhixun XIA revised and edited the final version.

    Conflict of interest

    Lian DUAN, Zhixun XIA, Yunchao FENG, Binbin CHEN, Jiarui ZHANG, and Likun MA declare that they have no conflict of interest.

    Burkholder TR, Andrews L, Bartlett RJ, 1993. Reaction of boron atoms with carbon dioxide: matrix and ab initio calculated infrared spectra of OBCO., 97(14):3500-3503. https://doi.org/10.1021/j100116a010

    Chen BB, Xia ZX, Huang LY, et al., 2018. Characteristics of the combustion chamber of a boron-based solid propellant ducted rocket with a chin-type inlet., 82-83:210-219. https://doi.org/10.1016/j.ast.2018.08.035

    Chin CH, Mebel AM, Hwang DY, 2003. Theoretical study of the reaction mechanism of boron atom with carbon dioxide., 375(5-6):670-675. https://doi.org/10.1016/S0009-2614(03)00964-3

    DiGiuseppe TG, Davidovits P, 1981. Boron atom reactions. II. Rate constants withO2, SO2, CO2, and N2O., 74(6):3287-3291. https://doi.org/10.1063/1.441534

    Duan L, Xia ZX, Chen BB, et al., 2022. Ignition and combustion characteristics of boron agglomerates under different oxygen concentrations., 197:81-90. https://doi.org/10.1016/j.actaastro.2022.05.020

    Foelsche RO, Burton RL, Krier H, 1999. Boron particle ignition and combustion at 30?–?150 atm., 117(1-2):32-58. https://doi.org/10.1016/S0010-2180(98)00080-7

    Fry RS, 2004. A century of ramjet propulsion technology evolution., 20(1):27-58. https://doi.org/10.2514/1.9178

    Hashim SA, Islam M, Kangle SM, et al., 2021. Performance evaluation of boron/hydroxyl-terminated polybutadiene-based solid fuels containing activated charcoal., 58(2):363-374. https://doi.org/10.2514/1.A34820

    Jain A, Anthonysamy S, 2015. Oxidation of boron carbide powder., 122(2):645-652. https://doi.org/10.1007/s10973-015-4818-3

    King MK, 1973. Boron particle ignition in hot gas streams., 8(5-6):255-273. https://doi.org/10.1080/00102207308946648

    Krier H, Burton RL, Pirman SR, et al., 1996. Shock initiation of crystalline boron in oxygen and fluorine compounds., 12(4):672-679. https://doi.org/10.2514/3.24088

    Li HP, Ao W, Wang Y, et al., 2014. Effect of carbon dioxide on the reactivity of the oxidation of boron particles., 39(4):617-623. https://doi.org/10.1002/prep.201300180

    Li XP, Ge LH, Luan XT, 2007. Applications of gas generator in ramjet direct-connect test facility., 33(3):14-19 (in Chinese). https://doi.org/10.3969/j.issn.1672-9374.2007.03.003

    Li YQ, Qiu T, 2007. Oxidation behaviour of boron carbide powder., 444(1-2):184-191. https://doi.org/10.1016/j.msea.2006.08.068

    Liang DL, Liu JZ, Zhou YN, et al., 2017. Ignition and combustion characteristics of molded amorphous boron under different oxygen pressures., 138:?118-128. https://doi.org/10.1016/j.actaastro.2017.05.019

    Liu LL, He GQ, Wang YH, et al., 2015. Chemical analysis of primary combustion products of boron-based fuel-rich propellant., 5(123):101416-101426. https://doi.org/10.1039/C5RA13693H

    Liu LL, He GQ, Wang YH, et al., 2017. Factors affecting the primary combustion products of boron-based fuel-rich propellants., 33(2):333-337. https://doi.org/10.2514/1.B36134

    Lv Z, Xia ZX, Liu B, et al., 2017. Preliminary experimental study on solid-fuel rocket scramjet combustor., 18(2):106-112. https://doi.org/10.1631/jzus.A1600489

    Meerov D, Monogarov K, Bragin A, et al., 2015. Boron particles agglomeration and slag formation during combustion of energetic condensed systems., 72:85-88. https://doi.org/10.1016/j.phpro.2015.09.024

    Mi XC, Goroshin S, Higgins AJ, et al., 2013. Dual-stage ignition of boron particle agglomerates., 160(11):2608-2618. https://doi.org/10.1016/j.combustflame.2013.06.004

    Millot F, Rifflet JC, Sarou-Kanian V, et al., 2002. High-temperature properties of liquid boron from contactless techniques., 23(5):1185-1195. https://doi.org/10.1023/A:1019836102776

    Roux JA, Choi J, Shakya N, 2014. Parametric scramjet cycle analysis for nonideal mass flow rate., 28(1):166-171. https://doi.org/10.2514/1.T4217

    Smolanoff J, Sowa-Resat M, ?apicki A, et al., 1996. Kinetic parameters for heterogenous boron combustion reactions via the Cluster Beam approach., 105(1-2):68-79. https://doi.org/10.1016/0010-2180(95)00155-7

    Song QG, Cao W, Wei X, et al., 2021. Laser ignition and combustion characteristics of micro- and nano-sized boron under different atmospheres and pressures., 230:111420. https://doi.org/10.1016/j.combustflame.2021.111420

    Sun YL, Ren H, Du FZ, et al., 2018. Preparation and characterization of sintered B/MgB2as heat release material., 759:100-107. https://doi.org/10.1016/j.jallcom.2018.05.038

    Sun YL, Ren H, Jiao QJ, et al., 2020. Oxidation, ignition and combustion behaviors of differently prepared boron-magnesium composites., 221:11-19. https://doi.org/10.1016/j.combustflame.2020.07.022

    Ulas A, Kuo KK, Gotzmer C, 2001. Ignition and combustion of boron particles in fluorine-containing environments., 127(1-2):1935-1957. https://doi.org/10.1016/S0010-2180(01)00299-1

    Yetter RA, Rabitz H, Dryer FL, et al., 1991. Kinetics of high-temperature B/O/H/C chemistry., 83(1-2):43-62. https://doi.org/10.1016/0010-2180(91)90202-M

    Yoshida T, Yuasa S, 2000. Effect of water vapor on ignition and combustion of boron lumps in an oxygen stream., 28(2):2735-2741. https://doi.org/10.1016/S0082-0784(00)80694-3

    Yuasa S, Isoda H, 1991. Ignition and combustion of small boron lumps in an oxygen stream., 86(3):216-222. https://doi.org/10.1016/0010-2180(91)90101-G

    Zhang H, Wang NF, Wu ZW, 2020. Effect of fuel grain configuration on the thrust of a solid-fuel scramjet., 106:106145. https://doi.org/10.1016/j.ast.2020.106145

    Zhou W, 1998. Numerical Study of Multi-Phase Combustion: Ignition and Combustion of an Isolated Boron Particle in Fluorinated Environments. PhD Thesis, Princeton University, Princeton, USA.

    Zhou W, Yetter RA, Dryer FL, et al., 1998. Effect of fluorine on the combustion of “clean” surface boron particles., 112(4):507-521. https://doi.org/10.1016/S0010-2180(97)00129-6

    Zhou W, Yetter RA, Dryer FL, et al., 1999. Multi-phase model for ignition and combustion of boron particles., 117(1-2):227-243. https://doi.org/10.1016/S0010-2180(98)00079-0

    Electronic supplementary materials

    Section S1

    題目:含氧氣氛中的二氧化碳濃度對硼團聚體燃燒特性的影響

    作者:段煉,夏智勛,馮運超,陳斌斌,張家瑞,馬立坤

    機構(gòu):國防科技大學,空天科學學院,中國長沙,410073

    目的:在沖壓發(fā)動機燃燒室中,含碳燃料燃燒產(chǎn)生的CO2與硼團聚體一起進入燃燒室。本文旨在探究在O2濃度固定的情況下,CO2濃度對團聚硼燃燒特性的影響規(guī)律,以深入認識在O2和CO2共存的情況下,團聚硼的燃燒模式以及硼的反應和能量釋放路徑。

    創(chuàng)新點:1. 直接觀測燃燒過程中,團聚硼的表面形貌的變化過程; 2. 通過凝相燃燒產(chǎn)物分析測試獲得團聚硼在O2和CO2共存的氣氛中的反應消耗路徑。

    方法:1. 通過激光點火實驗,直接觀察在不同CO2濃度的氣氛下,團聚硼燃燒過程中的表面形貌、顆粒結(jié)構(gòu)和火焰形貌的演變過程(圖3~5);2. 分析得到凝相燃燒產(chǎn)物的表面形貌、元素成分以及晶體結(jié)構(gòu)(圖10~12); 3. 通過熱重-差示掃描量熱法獲得在不同CO2濃度的氣氛中,團聚硼的低溫氧化過程。

    結(jié)論:1. 硼團聚體在高溫含氧氣氛中可與CO2反應生成B4C,增加硼的反應和消耗途徑,改變了顆粒的物理化學性質(zhì)和燃燒狀態(tài);2. 在O2濃度固定為20%的氣氛中,當顆粒溫度達到硼的熔點時,CO2含量的增加使得硼團聚體的燃燒模式由單顆粒液滴燃燒轉(zhuǎn)變?yōu)槎嗫最w粒燃燒;3. 當顆粒溫度低于B2O3的沸點時,B4C氧化生成的B2O3形成玻璃態(tài)液膜覆蓋在顆粒表面,阻礙顆粒進一步燃燒。

    關(guān)鍵詞:硼燃燒;無定形硼;含硼推進劑;固體燃料沖壓發(fā)動機

    https://doi.org/10.1631/jzus.A2200468

    https://doi.org/10.1631/jzus.A2200468

    ? Zhejiang University Press 2023

    Oct. 4, 2022;

    Feb. 22, 2023;

    Aug. 7, 2023

    黄片wwwwww| 中文字幕熟女人妻在线| 欧美色视频一区免费| 18禁动态无遮挡网站| 免费观看精品视频网站| av福利片在线观看| av免费在线看不卡| 午夜免费激情av| 99久国产av精品| 91狼人影院| 精品一区二区三区人妻视频| 亚洲av成人av| 人妻少妇偷人精品九色| 国产午夜精品一二区理论片| 最近中文字幕2019免费版| 91av网一区二区| 欧美性感艳星| 91精品国产九色| 尾随美女入室| 国产 一区精品| 久久久精品欧美日韩精品| 成人av在线播放网站| 十八禁国产超污无遮挡网站| 日本黄大片高清| 日韩欧美三级三区| 精品欧美国产一区二区三| 中文字幕人妻熟人妻熟丝袜美| 伦理电影大哥的女人| 91av网一区二区| 观看美女的网站| 人妻少妇偷人精品九色| 99久久精品一区二区三区| 亚洲欧美日韩东京热| 成年版毛片免费区| 国产av码专区亚洲av| 亚洲激情五月婷婷啪啪| 91久久精品国产一区二区三区| 精品一区二区三区视频在线| 波多野结衣巨乳人妻| 男的添女的下面高潮视频| 国产成人精品婷婷| 国产淫片久久久久久久久| 男女那种视频在线观看| 欧美日韩一区二区视频在线观看视频在线 | 一个人免费在线观看电影| a级毛色黄片| 日韩大片免费观看网站 | 国产一区有黄有色的免费视频 | 又爽又黄a免费视频| 国产精品野战在线观看| 国产高潮美女av| 中文精品一卡2卡3卡4更新| 最近手机中文字幕大全| 联通29元200g的流量卡| 婷婷色麻豆天堂久久 | 国产精品.久久久| www日本黄色视频网| 免费观看性生交大片5| 国产精品一二三区在线看| videos熟女内射| 成人一区二区视频在线观看| 久久久久久久久久久免费av| 99久久中文字幕三级久久日本| 六月丁香七月| 国产一区有黄有色的免费视频 | 国产三级在线视频| 久久6这里有精品| 亚洲av不卡在线观看| 天天一区二区日本电影三级| 精品久久久久久电影网 | 国产乱人视频| 欧美性感艳星| 蜜桃亚洲精品一区二区三区| 国产白丝娇喘喷水9色精品| 久久亚洲精品不卡| 99久久精品国产国产毛片| 国产精品麻豆人妻色哟哟久久 | 国产av不卡久久| 国产淫片久久久久久久久| 国产伦精品一区二区三区四那| 22中文网久久字幕| 亚洲av.av天堂| 久久精品夜夜夜夜夜久久蜜豆| 在线免费十八禁| 国内精品宾馆在线| 久久久久九九精品影院| 18禁动态无遮挡网站| 国产伦在线观看视频一区| 亚洲精品aⅴ在线观看| 日本wwww免费看| 亚洲熟妇中文字幕五十中出| 日本av手机在线免费观看| 丰满少妇做爰视频| 亚洲欧美一区二区三区国产| 亚洲av不卡在线观看| 一级毛片aaaaaa免费看小| www.av在线官网国产| 身体一侧抽搐| av免费观看日本| 国产国拍精品亚洲av在线观看| 两个人视频免费观看高清| 热99re8久久精品国产| 中文字幕人妻熟人妻熟丝袜美| 尾随美女入室| 我的老师免费观看完整版| 九九久久精品国产亚洲av麻豆| 色哟哟·www| 男女国产视频网站| 97在线视频观看| 男女国产视频网站| 日日干狠狠操夜夜爽| 亚洲国产精品成人久久小说| 嘟嘟电影网在线观看| 欧美潮喷喷水| 亚洲av二区三区四区| 欧美又色又爽又黄视频| 色吧在线观看| 69人妻影院| 麻豆精品久久久久久蜜桃| 亚洲美女搞黄在线观看| 在线天堂最新版资源| 亚洲av福利一区| 欧美xxxx性猛交bbbb| 欧美日韩国产亚洲二区| 国产精品久久久久久精品电影| 亚洲精品国产av成人精品| 国产探花在线观看一区二区| 国产三级中文精品| 99国产精品一区二区蜜桃av| 老师上课跳d突然被开到最大视频| 日本黄大片高清| 淫秽高清视频在线观看| 亚洲va在线va天堂va国产| 国产亚洲精品久久久com| 精品99又大又爽又粗少妇毛片| 久久这里有精品视频免费| 波多野结衣高清无吗| 综合色丁香网| 男女那种视频在线观看| 亚洲自拍偷在线| 美女内射精品一级片tv| 欧美潮喷喷水| 久久亚洲精品不卡| 最近的中文字幕免费完整| 欧美人与善性xxx| 搡女人真爽免费视频火全软件| 免费观看的影片在线观看| 午夜视频国产福利| 日韩精品有码人妻一区| 水蜜桃什么品种好| 看十八女毛片水多多多| 男女国产视频网站| 国产久久久一区二区三区| 国产极品精品免费视频能看的| 97在线视频观看| 极品教师在线视频| 纵有疾风起免费观看全集完整版 | 国产国拍精品亚洲av在线观看| 国产视频内射| 网址你懂的国产日韩在线| 黑人高潮一二区| 热99re8久久精品国产| 久久鲁丝午夜福利片| 精品一区二区免费观看| 一级黄片播放器| 女人久久www免费人成看片 | 波多野结衣高清无吗| 日本av手机在线免费观看| 精品久久久久久久久久久久久| .国产精品久久| 天堂√8在线中文| 日本黄色片子视频| 三级国产精品欧美在线观看| 两个人视频免费观看高清| 亚洲av成人精品一区久久| 久久欧美精品欧美久久欧美| 亚洲欧美清纯卡通| 亚洲欧美日韩东京热| 免费不卡的大黄色大毛片视频在线观看 | 热99re8久久精品国产| 岛国毛片在线播放| 少妇的逼好多水| 亚洲精品456在线播放app| 精品不卡国产一区二区三区| 日本av手机在线免费观看| 国产 一区 欧美 日韩| 亚洲精品自拍成人| 亚洲婷婷狠狠爱综合网| 亚洲国产日韩欧美精品在线观看| 日韩精品有码人妻一区| 三级毛片av免费| 免费无遮挡裸体视频| 国产精品一区二区三区四区久久| 免费黄网站久久成人精品| 大话2 男鬼变身卡| 91久久精品国产一区二区成人| 99热这里只有是精品50| eeuss影院久久| 中文乱码字字幕精品一区二区三区 | 大香蕉久久网| 国产成人精品婷婷| 高清在线视频一区二区三区 | 老司机影院成人| 亚洲性久久影院| 国产成人免费观看mmmm| 久久草成人影院| 99久久成人亚洲精品观看| 日本熟妇午夜| 久久人人爽人人爽人人片va| 极品教师在线视频| 久久欧美精品欧美久久欧美| 97热精品久久久久久| 国产又色又爽无遮挡免| 亚洲美女视频黄频| 99热全是精品| 国产午夜福利久久久久久| 18禁在线播放成人免费| 久久人妻av系列| 村上凉子中文字幕在线| 性色avwww在线观看| 尾随美女入室| 亚洲精品色激情综合| 麻豆av噜噜一区二区三区| 国产av码专区亚洲av| 成人特级av手机在线观看| 亚洲人与动物交配视频| 国产精品野战在线观看| 人妻夜夜爽99麻豆av| 亚洲自拍偷在线| 舔av片在线| 国产69精品久久久久777片| 成人午夜精彩视频在线观看| 精品国产一区二区三区久久久樱花 | 国产精品野战在线观看| ponron亚洲| 我要搜黄色片| 色尼玛亚洲综合影院| 在线天堂最新版资源| 国产精品国产三级国产av玫瑰| 亚洲国产最新在线播放| 久久久久久久国产电影| 成年女人看的毛片在线观看| 国产成人免费观看mmmm| 国产亚洲精品久久久com| 国产色婷婷99| 2022亚洲国产成人精品| 亚洲五月天丁香| 亚洲人成网站在线播| www.av在线官网国产| 久久精品影院6| 麻豆乱淫一区二区| 久久鲁丝午夜福利片| 婷婷六月久久综合丁香| 岛国在线免费视频观看| 午夜精品国产一区二区电影 | 国产黄色视频一区二区在线观看 | 欧美性感艳星| av线在线观看网站| av.在线天堂| 少妇裸体淫交视频免费看高清| 国产av一区在线观看免费| 22中文网久久字幕| 成人午夜高清在线视频| 久久人妻av系列| 国产精品福利在线免费观看| 91在线精品国自产拍蜜月| 免费黄网站久久成人精品| 日本免费一区二区三区高清不卡| 一区二区三区四区激情视频| 国产免费又黄又爽又色| 亚洲精品亚洲一区二区| 99久久无色码亚洲精品果冻| 免费观看精品视频网站| 亚洲中文字幕一区二区三区有码在线看| 亚洲熟妇中文字幕五十中出| 免费黄色在线免费观看| 国产黄片美女视频| 精品久久久久久久久av| 午夜激情福利司机影院| 天堂av国产一区二区熟女人妻| 国产三级在线视频| 国产免费视频播放在线视频 | 波野结衣二区三区在线| 国产精品国产三级国产av玫瑰| 欧美xxxx黑人xx丫x性爽| 看片在线看免费视频| 精品久久久久久电影网 | 欧美最新免费一区二区三区| 欧美日韩在线观看h| 91精品伊人久久大香线蕉| 亚洲av成人精品一区久久| 国产三级中文精品| 女人十人毛片免费观看3o分钟| 亚洲天堂国产精品一区在线| 久久久久久九九精品二区国产| 99九九线精品视频在线观看视频| 晚上一个人看的免费电影| 欧美区成人在线视频| 亚洲国产欧美人成| 老司机福利观看| 人妻制服诱惑在线中文字幕| 深夜a级毛片| 美女内射精品一级片tv| 国产伦在线观看视频一区| 麻豆国产97在线/欧美| av在线老鸭窝| 国产欧美另类精品又又久久亚洲欧美| 久久精品国产亚洲av天美| 国产 一区精品| 久久精品国产亚洲网站| 蜜臀久久99精品久久宅男| 人妻少妇偷人精品九色| 欧美一区二区精品小视频在线| 午夜福利视频1000在线观看| 亚洲中文字幕日韩| 插阴视频在线观看视频| 小说图片视频综合网站| 午夜福利在线观看吧| 亚洲最大成人中文| 精品免费久久久久久久清纯| 中国美白少妇内射xxxbb| 精品久久国产蜜桃| 亚洲丝袜综合中文字幕| 国产一区二区在线av高清观看| 观看免费一级毛片| 在线免费观看的www视频| 一卡2卡三卡四卡精品乱码亚洲| 国产真实伦视频高清在线观看| 国产美女午夜福利| 亚洲中文字幕一区二区三区有码在线看| 成人无遮挡网站| 三级经典国产精品| 97热精品久久久久久| 网址你懂的国产日韩在线| 一区二区三区四区激情视频| 欧美不卡视频在线免费观看| 国产成人a∨麻豆精品| 国内精品一区二区在线观看| 色播亚洲综合网| 亚洲,欧美,日韩| 一级二级三级毛片免费看| 午夜免费男女啪啪视频观看| 综合色av麻豆| 免费av观看视频| 九九热线精品视视频播放| 国产伦精品一区二区三区四那| АⅤ资源中文在线天堂| 日韩一区二区视频免费看| 变态另类丝袜制服| 可以在线观看毛片的网站| 亚洲av福利一区| 国产成人精品一,二区| 97超视频在线观看视频| 成人综合一区亚洲| 中文精品一卡2卡3卡4更新| 亚洲中文字幕一区二区三区有码在线看| 日韩成人av中文字幕在线观看| 中国美白少妇内射xxxbb| 亚洲最大成人av| 欧美变态另类bdsm刘玥| 欧美性猛交╳xxx乱大交人| 在线播放无遮挡| 国产亚洲最大av| 精品久久久久久久人妻蜜臀av| 国产一区二区在线av高清观看| 亚洲丝袜综合中文字幕| 国产视频首页在线观看| 日韩欧美在线乱码| 少妇裸体淫交视频免费看高清| 伦精品一区二区三区| 免费观看人在逋| 91狼人影院| 免费大片18禁| 伊人久久精品亚洲午夜| 女人十人毛片免费观看3o分钟| 久久鲁丝午夜福利片| 亚洲欧美日韩卡通动漫| 国产高清不卡午夜福利| 亚洲最大成人手机在线| 日本免费一区二区三区高清不卡| 韩国av在线不卡| 日本五十路高清| 韩国av在线不卡| 只有这里有精品99| 蜜桃亚洲精品一区二区三区| 亚洲欧美成人综合另类久久久 | 欧美一区二区精品小视频在线| www.色视频.com| 日韩欧美 国产精品| 久久精品综合一区二区三区| 干丝袜人妻中文字幕| 波野结衣二区三区在线| 我的女老师完整版在线观看| 亚洲国产色片| 国产精品乱码一区二三区的特点| 亚洲熟妇中文字幕五十中出| 成年av动漫网址| 亚洲精品亚洲一区二区| 午夜免费激情av| 伦理电影大哥的女人| 最近视频中文字幕2019在线8| 99热这里只有是精品50| 午夜a级毛片| 青春草国产在线视频| 99视频精品全部免费 在线| 午夜亚洲福利在线播放| 啦啦啦韩国在线观看视频| 久久婷婷人人爽人人干人人爱| 免费在线观看成人毛片| 我的老师免费观看完整版| 国产精品久久视频播放| 一个人看视频在线观看www免费| 爱豆传媒免费全集在线观看| 亚洲天堂国产精品一区在线| АⅤ资源中文在线天堂| 欧美日本视频| 久久6这里有精品| 三级男女做爰猛烈吃奶摸视频| 午夜免费激情av| 两个人的视频大全免费| 国产精品国产三级国产专区5o | 最近的中文字幕免费完整| 人人妻人人澡欧美一区二区| 久久久久久久国产电影| 亚洲精品乱久久久久久| 成人午夜精彩视频在线观看| 精品免费久久久久久久清纯| 亚洲成人av在线免费| 在现免费观看毛片| 日韩av在线免费看完整版不卡| 成年免费大片在线观看| 国产 一区精品| 日本三级黄在线观看| 国产精品国产三级国产av玫瑰| 成人一区二区视频在线观看| 国产精品女同一区二区软件| 伦理电影大哥的女人| 一本久久精品| 亚洲精品自拍成人| 亚洲欧美精品综合久久99| 国产精品久久久久久精品电影| 成人美女网站在线观看视频| 我要看日韩黄色一级片| 黄色配什么色好看| 久久国内精品自在自线图片| 久久久亚洲精品成人影院| 干丝袜人妻中文字幕| 亚洲美女视频黄频| 精品国内亚洲2022精品成人| 偷拍熟女少妇极品色| 亚洲国产精品合色在线| 久久国内精品自在自线图片| 汤姆久久久久久久影院中文字幕 | 亚洲人成网站高清观看| 国产高清国产精品国产三级 | 午夜福利视频1000在线观看| 午夜日本视频在线| 精品99又大又爽又粗少妇毛片| 美女黄网站色视频| 麻豆一二三区av精品| 亚洲av一区综合| 伦理电影大哥的女人| 2021少妇久久久久久久久久久| 淫秽高清视频在线观看| 日韩国内少妇激情av| 成人无遮挡网站| 成人毛片a级毛片在线播放| 麻豆av噜噜一区二区三区| 免费播放大片免费观看视频在线观看 | .国产精品久久| 1024手机看黄色片| 国产午夜精品久久久久久一区二区三区| 亚洲美女搞黄在线观看| 变态另类丝袜制服| 久久久久久久国产电影| 亚洲av成人av| 亚洲在久久综合| 天天躁日日操中文字幕| 国产精品一区二区性色av| 麻豆精品久久久久久蜜桃| 亚洲欧美成人综合另类久久久 | 国产黄色小视频在线观看| 欧美日韩综合久久久久久| 久久久久网色| 亚洲精品日韩av片在线观看| 男女那种视频在线观看| 一级二级三级毛片免费看| 亚洲经典国产精华液单| 精品一区二区三区视频在线| av免费观看日本| 日本一本二区三区精品| 激情 狠狠 欧美| 国语自产精品视频在线第100页| 国产老妇伦熟女老妇高清| 国产亚洲5aaaaa淫片| av线在线观看网站| 97在线视频观看| 国产精品一区二区三区四区免费观看| 好男人在线观看高清免费视频| 国产三级在线视频| 中文在线观看免费www的网站| 九九热线精品视视频播放| 色综合色国产| 在线免费十八禁| 日本午夜av视频| 最近最新中文字幕大全电影3| 91久久精品国产一区二区三区| 白带黄色成豆腐渣| 乱系列少妇在线播放| 少妇被粗大猛烈的视频| 国产免费又黄又爽又色| 最近中文字幕高清免费大全6| 成年女人看的毛片在线观看| 女人被狂操c到高潮| 99久久中文字幕三级久久日本| 少妇的逼水好多| 99国产精品一区二区蜜桃av| 欧美精品一区二区大全| 亚洲人成网站在线播| 精品久久久久久久久av| 白带黄色成豆腐渣| 日韩 亚洲 欧美在线| 黄片wwwwww| 99热精品在线国产| 丝袜美腿在线中文| 精品人妻偷拍中文字幕| 91精品国产九色| 欧美精品国产亚洲| 成人漫画全彩无遮挡| 亚洲成人av在线免费| 99热全是精品| 七月丁香在线播放| 91精品一卡2卡3卡4卡| 亚洲精品影视一区二区三区av| 日韩av不卡免费在线播放| 亚洲不卡免费看| 97在线视频观看| 黄色日韩在线| 日韩一区二区三区影片| 丰满人妻一区二区三区视频av| 在线播放无遮挡| 插阴视频在线观看视频| 村上凉子中文字幕在线| 国产免费又黄又爽又色| 亚洲aⅴ乱码一区二区在线播放| 国产精品综合久久久久久久免费| 欧美高清性xxxxhd video| 精品久久国产蜜桃| 一个人看的www免费观看视频| 少妇人妻精品综合一区二区| 久久久久久久久大av| 国产又色又爽无遮挡免| 亚洲综合精品二区| 日本与韩国留学比较| 国产精品麻豆人妻色哟哟久久 | 中文字幕免费在线视频6| 成年女人永久免费观看视频| 黄色日韩在线| 国产精品一二三区在线看| 日韩三级伦理在线观看| 黄色配什么色好看| 最近的中文字幕免费完整| 99热网站在线观看| 欧美+日韩+精品| 亚洲四区av| 久久精品国产自在天天线| av专区在线播放| 欧美成人a在线观看| 日产精品乱码卡一卡2卡三| 男女那种视频在线观看| av视频在线观看入口| 91午夜精品亚洲一区二区三区| 又爽又黄a免费视频| 青青草视频在线视频观看| www.av在线官网国产| 日韩高清综合在线| 日韩欧美 国产精品| 国模一区二区三区四区视频| 麻豆成人午夜福利视频| 国产三级中文精品| 久久精品国产99精品国产亚洲性色| 少妇被粗大猛烈的视频| 久久久a久久爽久久v久久| 国产一区二区在线av高清观看| 国产成人精品一,二区| 欧美性感艳星| 波多野结衣巨乳人妻| 两个人视频免费观看高清| 免费观看在线日韩| 一级黄片播放器| 在线观看一区二区三区| 天堂网av新在线| 长腿黑丝高跟| 秋霞伦理黄片| 国产亚洲91精品色在线| 免费人成在线观看视频色| 男女国产视频网站| 国产伦一二天堂av在线观看| 最近中文字幕高清免费大全6| 亚洲国产精品合色在线| 男女视频在线观看网站免费| 一级黄色大片毛片| 午夜福利视频1000在线观看| 亚洲av日韩在线播放| 能在线免费看毛片的网站| 全区人妻精品视频| 日韩亚洲欧美综合| 国产精品综合久久久久久久免费| 最近2019中文字幕mv第一页| 波多野结衣高清无吗| 国产69精品久久久久777片| 精品久久久久久久久亚洲| 久久精品国产亚洲av天美| 国产精品,欧美在线| 联通29元200g的流量卡|