• <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

    x7x7x7水蜜桃| 午夜福利一区二区在线看| 老司机靠b影院| 天天躁日日躁夜夜躁夜夜| 久久影院123| 高潮久久久久久久久久久不卡| 欧美精品啪啪一区二区三区| 涩涩av久久男人的天堂| 老司机影院毛片| 亚洲一卡2卡3卡4卡5卡精品中文| 久久影院123| 国产精品国产高清国产av | 女警被强在线播放| 精品第一国产精品| 少妇 在线观看| 亚洲精品在线观看二区| 亚洲国产看品久久| 成人亚洲精品一区在线观看| 国产成人精品无人区| 成人黄色视频免费在线看| 成人黄色视频免费在线看| 最近最新中文字幕大全电影3 | 岛国毛片在线播放| 黄色片一级片一级黄色片| 看免费av毛片| 久久精品国产a三级三级三级| 亚洲av成人av| 黄片大片在线免费观看| 淫妇啪啪啪对白视频| 一边摸一边做爽爽视频免费| 午夜福利,免费看| 中国美女看黄片| 国产男女内射视频| 99国产精品一区二区三区| 久久人人97超碰香蕉20202| 天堂√8在线中文| 久久精品亚洲av国产电影网| 午夜成年电影在线免费观看| 亚洲欧美精品综合一区二区三区| 国产一卡二卡三卡精品| tocl精华| 一区二区三区国产精品乱码| 新久久久久国产一级毛片| 露出奶头的视频| 国产午夜精品久久久久久| 18禁国产床啪视频网站| 男人操女人黄网站| 男男h啪啪无遮挡| 国产精品综合久久久久久久免费 | 欧美日韩黄片免| 日韩欧美一区视频在线观看| 免费人成视频x8x8入口观看| svipshipincom国产片| 下体分泌物呈黄色| 国产不卡av网站在线观看| 中亚洲国语对白在线视频| 美女扒开内裤让男人捅视频| 在线观看免费视频日本深夜| 久久国产精品男人的天堂亚洲| 免费观看人在逋| 十分钟在线观看高清视频www| 久久精品国产综合久久久| 男人舔女人的私密视频| 久久人妻福利社区极品人妻图片| 亚洲综合色网址| 欧美 亚洲 国产 日韩一| 精品久久久久久,| 国产有黄有色有爽视频| 国产精品免费大片| 亚洲专区国产一区二区| 亚洲成av片中文字幕在线观看| 亚洲人成电影观看| 人妻丰满熟妇av一区二区三区 | 大陆偷拍与自拍| 老司机午夜十八禁免费视频| av不卡在线播放| 丁香欧美五月| 亚洲性夜色夜夜综合| 两人在一起打扑克的视频| 欧美乱码精品一区二区三区| 国产不卡一卡二| 欧美在线一区亚洲| 黄片播放在线免费| 亚洲少妇的诱惑av| av网站免费在线观看视频| 亚洲精品国产一区二区精华液| 美女国产高潮福利片在线看| 18禁国产床啪视频网站| 日本a在线网址| 欧美成人免费av一区二区三区 | 久久中文字幕人妻熟女| 久久久国产成人精品二区 | 新久久久久国产一级毛片| 中文字幕另类日韩欧美亚洲嫩草| 欧洲精品卡2卡3卡4卡5卡区| 在线观看66精品国产| 热99久久久久精品小说推荐| 欧美人与性动交α欧美精品济南到| 脱女人内裤的视频| 国产亚洲精品第一综合不卡| 麻豆av在线久日| 国产高清国产精品国产三级| 男人舔女人的私密视频| 亚洲一码二码三码区别大吗| 久久久久精品国产欧美久久久| 国产视频一区二区在线看| 99热只有精品国产| 国产欧美日韩综合在线一区二区| 搡老熟女国产l中国老女人| 亚洲国产毛片av蜜桃av| 久久精品成人免费网站| 久久国产亚洲av麻豆专区| 婷婷精品国产亚洲av在线 | 亚洲中文字幕日韩| 亚洲精品国产精品久久久不卡| 人人妻人人澡人人看| 久久精品国产综合久久久| 久久国产精品男人的天堂亚洲| 亚洲国产欧美一区二区综合| av天堂久久9| 免费在线观看影片大全网站| 国产成人精品在线电影| 色尼玛亚洲综合影院| 建设人人有责人人尽责人人享有的| 99re6热这里在线精品视频| 亚洲九九香蕉| 精品人妻熟女毛片av久久网站| 国产亚洲欧美在线一区二区| 露出奶头的视频| 久久久国产成人免费| 少妇的丰满在线观看| 国产精品永久免费网站| 男男h啪啪无遮挡| 一本综合久久免费| 久久精品亚洲精品国产色婷小说| 老司机靠b影院| 一个人免费在线观看的高清视频| 美女扒开内裤让男人捅视频| 黄色 视频免费看| 欧美日韩中文字幕国产精品一区二区三区 | 咕卡用的链子| 国产真人三级小视频在线观看| 中国美女看黄片| 欧美最黄视频在线播放免费 | 亚洲精品久久成人aⅴ小说| 十八禁人妻一区二区| av国产精品久久久久影院| 国产亚洲精品第一综合不卡| 国产精品久久久久久人妻精品电影| 在线播放国产精品三级| 亚洲伊人色综图| 亚洲一码二码三码区别大吗| 成人18禁高潮啪啪吃奶动态图| 12—13女人毛片做爰片一| 精品久久久久久久毛片微露脸| 精品午夜福利视频在线观看一区| 身体一侧抽搐| 久久99一区二区三区| www.熟女人妻精品国产| 女人被躁到高潮嗷嗷叫费观| 高清黄色对白视频在线免费看| 伦理电影免费视频| 成人av一区二区三区在线看| 天天操日日干夜夜撸| 丰满人妻熟妇乱又伦精品不卡| 每晚都被弄得嗷嗷叫到高潮| 电影成人av| 夜夜躁狠狠躁天天躁| www.自偷自拍.com| 1024视频免费在线观看| 欧美最黄视频在线播放免费 | 在线免费观看的www视频| 一进一出抽搐动态| 国产精品九九99| 欧美黄色淫秽网站| 丝袜在线中文字幕| av网站免费在线观看视频| 正在播放国产对白刺激| tube8黄色片| 一边摸一边抽搐一进一小说 | 亚洲国产精品一区二区三区在线| 丁香六月欧美| 国产野战对白在线观看| 久久午夜亚洲精品久久| 成年人午夜在线观看视频| 亚洲视频免费观看视频| 首页视频小说图片口味搜索| 久久久精品区二区三区| 香蕉丝袜av| 男女高潮啪啪啪动态图| 欧美av亚洲av综合av国产av| 高清毛片免费观看视频网站 | 色老头精品视频在线观看| 80岁老熟妇乱子伦牲交| 久久这里只有精品19| 一级毛片高清免费大全| 久久精品熟女亚洲av麻豆精品| 国产欧美日韩一区二区精品| 黄色成人免费大全| 叶爱在线成人免费视频播放| 久久久久精品国产欧美久久久| 天天躁夜夜躁狠狠躁躁| 在线观看免费日韩欧美大片| 日韩免费高清中文字幕av| 国产精品免费一区二区三区在线 | 国产一区二区三区视频了| 亚洲九九香蕉| 最近最新中文字幕大全电影3 | 久久这里只有精品19| 一级a爱视频在线免费观看| 国产欧美亚洲国产| 亚洲avbb在线观看| av中文乱码字幕在线| 亚洲精品自拍成人| 99国产精品99久久久久| 亚洲男人天堂网一区| 亚洲成国产人片在线观看| 中文亚洲av片在线观看爽 | 亚洲专区国产一区二区| 国产成人精品在线电影| 高潮久久久久久久久久久不卡| 在线视频色国产色| 国产麻豆69| 一级a爱视频在线免费观看| 免费高清在线观看日韩| 亚洲精品国产色婷婷电影| 丝瓜视频免费看黄片| 国产精品 国内视频| 色播在线永久视频| 免费不卡黄色视频| 亚洲欧美精品综合一区二区三区| 国产一区二区三区在线臀色熟女 | 久久精品91无色码中文字幕| 老司机靠b影院| 亚洲第一av免费看| av有码第一页| 国产精品秋霞免费鲁丝片| 亚洲色图综合在线观看| 黄色怎么调成土黄色| 免费久久久久久久精品成人欧美视频| 精品久久久精品久久久| 国产欧美亚洲国产| 黄片播放在线免费| 他把我摸到了高潮在线观看| 日本撒尿小便嘘嘘汇集6| 欧美成人午夜精品| 国产成人精品无人区| 久久国产亚洲av麻豆专区| 亚洲精品久久成人aⅴ小说| 国产成人一区二区三区免费视频网站| 精品少妇一区二区三区视频日本电影| 男女下面插进去视频免费观看| av国产精品久久久久影院| 两个人免费观看高清视频| 国产成人一区二区三区免费视频网站| 悠悠久久av| 亚洲黑人精品在线| 一区二区三区国产精品乱码| 久久久精品免费免费高清| 亚洲人成伊人成综合网2020| 身体一侧抽搐| 久久久久久久久久久久大奶| 国产精品免费大片| 亚洲精品成人av观看孕妇| 国产淫语在线视频| 视频在线观看一区二区三区| 91大片在线观看| 精品国产超薄肉色丝袜足j| 村上凉子中文字幕在线| 啦啦啦 在线观看视频| av超薄肉色丝袜交足视频| 亚洲国产中文字幕在线视频| 777久久人妻少妇嫩草av网站| 欧美午夜高清在线| 国产成人免费观看mmmm| 黄色 视频免费看| 国产精品久久视频播放| 国产亚洲精品一区二区www | 久久久国产成人精品二区 | 五月开心婷婷网| 黄色丝袜av网址大全| 免费黄频网站在线观看国产| 久久久国产成人免费| 亚洲色图av天堂| 天天添夜夜摸| 中文字幕av电影在线播放| 亚洲成人国产一区在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲第一av免费看| 夜夜夜夜夜久久久久| 国产成人欧美| 午夜久久久在线观看| 亚洲精品在线观看二区| 精品国产一区二区久久| 国产精品综合久久久久久久免费 | 欧美黄色片欧美黄色片| 伦理电影免费视频| 在线永久观看黄色视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲成人免费av在线播放| 波多野结衣一区麻豆| 亚洲欧美精品综合一区二区三区| 久久亚洲真实| 777久久人妻少妇嫩草av网站| 美女高潮到喷水免费观看| 亚洲五月天丁香| 国产主播在线观看一区二区| 母亲3免费完整高清在线观看| 国产三级黄色录像| 日韩三级视频一区二区三区| e午夜精品久久久久久久| 亚洲中文字幕日韩| 一夜夜www| 精品国产乱子伦一区二区三区| av天堂在线播放| 亚洲av欧美aⅴ国产| 国产精品电影一区二区三区 | 最新的欧美精品一区二区| 91大片在线观看| 午夜福利视频在线观看免费| 欧美人与性动交α欧美精品济南到| 在线观看免费高清a一片| 啦啦啦免费观看视频1| 亚洲一区二区三区欧美精品| 脱女人内裤的视频| 久久国产精品男人的天堂亚洲| 国产精品成人在线| 国产在视频线精品| 国产蜜桃级精品一区二区三区 | 成人av一区二区三区在线看| 免费女性裸体啪啪无遮挡网站| 91大片在线观看| 午夜免费成人在线视频| 午夜福利在线观看吧| 日日夜夜操网爽| 国产淫语在线视频| 一区在线观看完整版| 热99国产精品久久久久久7| 日韩欧美三级三区| 十分钟在线观看高清视频www| 女同久久另类99精品国产91| 一个人免费在线观看的高清视频| 亚洲欧洲精品一区二区精品久久久| 亚洲av成人不卡在线观看播放网| av网站免费在线观看视频| 脱女人内裤的视频| 啦啦啦免费观看视频1| 免费观看a级毛片全部| www日本在线高清视频| 深夜精品福利| 最新在线观看一区二区三区| 在线观看免费视频日本深夜| 黄色视频,在线免费观看| 搡老乐熟女国产| 十八禁高潮呻吟视频| 在线观看午夜福利视频| 9191精品国产免费久久| 日韩欧美一区二区三区在线观看 | 日韩免费高清中文字幕av| 一本一本久久a久久精品综合妖精| 美女 人体艺术 gogo| 精品福利观看| 久久久久久久午夜电影 | 久久久久精品国产欧美久久久| 在线观看免费视频网站a站| 美女扒开内裤让男人捅视频| 在线永久观看黄色视频| 久久热在线av| 日本精品一区二区三区蜜桃| 一级毛片精品| 国产av精品麻豆| 精品卡一卡二卡四卡免费| 19禁男女啪啪无遮挡网站| 精品少妇久久久久久888优播| 国产精品二区激情视频| 免费在线观看日本一区| 美女高潮到喷水免费观看| 在线观看免费高清a一片| 90打野战视频偷拍视频| 亚洲专区字幕在线| 在线看a的网站| 亚洲va日本ⅴa欧美va伊人久久| 狠狠婷婷综合久久久久久88av| 女人爽到高潮嗷嗷叫在线视频| 国产av精品麻豆| 日韩免费高清中文字幕av| 中文字幕人妻丝袜制服| 免费日韩欧美在线观看| 最近最新中文字幕大全免费视频| 一二三四在线观看免费中文在| 午夜福利在线观看吧| 国产亚洲精品久久久久久毛片 | av线在线观看网站| 国产亚洲欧美98| videos熟女内射| 国产区一区二久久| avwww免费| 国内久久婷婷六月综合欲色啪| 欧美黑人精品巨大| 搡老乐熟女国产| 一本一本久久a久久精品综合妖精| 夜夜躁狠狠躁天天躁| 一本一本久久a久久精品综合妖精| 成人三级做爰电影| 亚洲全国av大片| 中国美女看黄片| 91精品国产国语对白视频| 欧美午夜高清在线| 可以免费在线观看a视频的电影网站| 久热爱精品视频在线9| videosex国产| 国精品久久久久久国模美| 一二三四社区在线视频社区8| 一级毛片精品| 91av网站免费观看| 一边摸一边抽搐一进一出视频| 一级,二级,三级黄色视频| 亚洲中文字幕日韩| 99香蕉大伊视频| 老汉色∧v一级毛片| 女警被强在线播放| 亚洲国产毛片av蜜桃av| 国产精品永久免费网站| 人妻一区二区av| 一级毛片精品| 免费女性裸体啪啪无遮挡网站| 久久精品亚洲精品国产色婷小说| 母亲3免费完整高清在线观看| 一二三四在线观看免费中文在| 亚洲人成电影观看| 飞空精品影院首页| 国产成人免费观看mmmm| 久久精品熟女亚洲av麻豆精品| 99热网站在线观看| 18在线观看网站| 王馨瑶露胸无遮挡在线观看| 亚洲七黄色美女视频| 下体分泌物呈黄色| 中文欧美无线码| 999久久久国产精品视频| av片东京热男人的天堂| 国产真人三级小视频在线观看| 免费久久久久久久精品成人欧美视频| 精品乱码久久久久久99久播| 一级毛片高清免费大全| 午夜精品在线福利| 婷婷精品国产亚洲av在线 | 国产男靠女视频免费网站| 99久久综合精品五月天人人| 国产不卡av网站在线观看| 久久人人爽av亚洲精品天堂| 久久香蕉国产精品| 午夜福利影视在线免费观看| 一级,二级,三级黄色视频| 欧美日韩黄片免| 无人区码免费观看不卡| 美女午夜性视频免费| 久久久久精品国产欧美久久久| 另类亚洲欧美激情| 国产精品一区二区免费欧美| 亚洲av日韩精品久久久久久密| 又黄又粗又硬又大视频| 日韩人妻精品一区2区三区| 亚洲中文字幕日韩| 精品无人区乱码1区二区| 巨乳人妻的诱惑在线观看| 婷婷丁香在线五月| 在线永久观看黄色视频| 搡老乐熟女国产| 国产午夜精品久久久久久| 一本综合久久免费| 18在线观看网站| 999久久久国产精品视频| 国产成人欧美在线观看 | 成人18禁高潮啪啪吃奶动态图| 日韩三级视频一区二区三区| 欧美 亚洲 国产 日韩一| 在线观看午夜福利视频| www.自偷自拍.com| av网站在线播放免费| 久久精品国产清高在天天线| 中文字幕高清在线视频| 亚洲七黄色美女视频| 亚洲三区欧美一区| 三级毛片av免费| 久久狼人影院| 免费在线观看日本一区| 久久人妻熟女aⅴ| 国产亚洲av高清不卡| 午夜精品久久久久久毛片777| 日本五十路高清| 黑人猛操日本美女一级片| 一本一本久久a久久精品综合妖精| 亚洲av成人一区二区三| 久久天躁狠狠躁夜夜2o2o| 久久婷婷成人综合色麻豆| 国产97色在线日韩免费| 国产91精品成人一区二区三区| 一级毛片高清免费大全| 久久久久久久久久久久大奶| 日韩免费av在线播放| 国产成人欧美在线观看 | 亚洲欧美一区二区三区久久| 精品福利永久在线观看| 一进一出抽搐gif免费好疼 | 欧美性长视频在线观看| 国产精品成人在线| 日韩成人在线观看一区二区三区| 国产高清videossex| 色94色欧美一区二区| 老司机在亚洲福利影院| 操出白浆在线播放| 日韩 欧美 亚洲 中文字幕| 亚洲精品国产色婷婷电影| 国产精品.久久久| 下体分泌物呈黄色| 国产欧美日韩综合在线一区二区| avwww免费| 亚洲精品av麻豆狂野| 天堂中文最新版在线下载| 亚洲三区欧美一区| 久久久久久亚洲精品国产蜜桃av| 91麻豆精品激情在线观看国产 | 色在线成人网| 国产在线精品亚洲第一网站| 欧美日韩乱码在线| 亚洲熟女毛片儿| 男人操女人黄网站| 国产真人三级小视频在线观看| 50天的宝宝边吃奶边哭怎么回事| 国内毛片毛片毛片毛片毛片| 久久久精品区二区三区| 少妇的丰满在线观看| 一进一出好大好爽视频| 丁香六月欧美| 国产精品 国内视频| www.自偷自拍.com| 正在播放国产对白刺激| 亚洲精品国产精品久久久不卡| 热re99久久国产66热| 99国产综合亚洲精品| 一边摸一边做爽爽视频免费| 色婷婷av一区二区三区视频| 婷婷精品国产亚洲av在线 | 精品人妻在线不人妻| 国精品久久久久久国模美| 日韩欧美免费精品| 激情视频va一区二区三区| 国产单亲对白刺激| 亚洲精品美女久久av网站| 一边摸一边抽搐一进一出视频| 一区二区三区精品91| 12—13女人毛片做爰片一| 很黄的视频免费| 亚洲国产精品sss在线观看 | 亚洲色图 男人天堂 中文字幕| 一二三四在线观看免费中文在| 91国产中文字幕| 老熟妇仑乱视频hdxx| 亚洲成a人片在线一区二区| 男女午夜视频在线观看| 老司机午夜十八禁免费视频| 丝瓜视频免费看黄片| 首页视频小说图片口味搜索| 久久精品熟女亚洲av麻豆精品| 少妇猛男粗大的猛烈进出视频| 久久亚洲真实| 亚洲少妇的诱惑av| 国产成人av激情在线播放| 男女之事视频高清在线观看| 亚洲精品中文字幕在线视频| 熟女少妇亚洲综合色aaa.| 午夜福利免费观看在线| 亚洲三区欧美一区| 99久久综合精品五月天人人| 亚洲中文av在线| 国产精品1区2区在线观看. | 老熟妇仑乱视频hdxx| 多毛熟女@视频| 俄罗斯特黄特色一大片| 999精品在线视频| 欧美中文综合在线视频| 亚洲性夜色夜夜综合| 超碰97精品在线观看| 国产精品免费大片| 亚洲性夜色夜夜综合| 不卡一级毛片| 法律面前人人平等表现在哪些方面| 老司机影院毛片| 久久久久久免费高清国产稀缺| 婷婷丁香在线五月| 美女高潮喷水抽搐中文字幕| 王馨瑶露胸无遮挡在线观看| 免费人成视频x8x8入口观看| 亚洲人成电影免费在线| 建设人人有责人人尽责人人享有的| 久久性视频一级片| 亚洲专区中文字幕在线| 可以免费在线观看a视频的电影网站| 亚洲片人在线观看| 国产aⅴ精品一区二区三区波| 高清视频免费观看一区二区| xxxhd国产人妻xxx| 国产精品九九99| 丝袜美足系列| 亚洲精品中文字幕一二三四区| 国产精品一区二区精品视频观看| 黄片大片在线免费观看| 老汉色∧v一级毛片| 免费在线观看亚洲国产| 亚洲全国av大片| 女人被狂操c到高潮| 国产xxxxx性猛交| 国产日韩一区二区三区精品不卡| 90打野战视频偷拍视频|