首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 453 毫秒
1.
Jürgen Troe 《Molecular physics》2014,112(18):2374-2383
The relationship between rate constants for dissociation and the reverse association reactions and their potential energy surfaces is illustrated. The reaction systems e? + SF6 ? SF6 ? →SF5 ? + F, H + CH3 ?CH4, 2 CF2 ? C2F4, H + O2 →HO2, HO + O ?HO2 ? H + O2, and C + HO →CHO are chosen as representative examples. The necessity to know precise thermochemical data is emphasised. The interplay between attractive and anisotropic components of the potentials influences the rate constants. Spin–orbit and electronic–rotational coupling in reactions between electronic open-shell radicals so far generally has been neglected, but is shown to have a marked influence on low temperature rate constants.  相似文献   

2.
Roaming     
Joel M. Bowman 《Molecular physics》2014,112(19):2516-2528
In this review the roaming pathway in chemical reaction dynamics is described and illustrated. After a brief review and illustration of this pathway, an update on the chemical reactions where this pathway has been reported is given. Then the focus is on two recent systems where new aspects of roaming have been reported, namely roaming-isomerisation in CH3NO2 unimolecular dissociation and multi-electronic-state roaming in NO3 photodissociation. Full-dimensional potentials for these systems have been developed and used in quasiclassical trajectory calculations. Selected results from these studies are reviewed. So-called roaming saddle points and associated intrinsic reaction coordinate paths have been determined for these systems. These will be reviewed and investigated by making relevant comparisons with quasiclassical trajectories. Manifestations of roaming in quantum scattering and resonance calculations of the MgH + H → Mg + H2 reaction are reviewed. That is followed by a review of several statistical theories of roaming. The review concludes with some ‘editorial’ comments about the challenges that roaming presents to approximate theories and comments about possible future research directions.  相似文献   

3.
刘新国  孙海竹  刘会荣  张庆刚 《物理学报》2010,59(11):7796-7802
运用准经典轨线方法,基于RODRIGO势能面,对碰撞能为20kcal/mol时,O++H2及其同位素取代反应的立体动力学性质进行了理论研究,对k-j′两矢量相关和k-k′-j′三矢量相关的分布函数、极化微分反应截面,以及产物转动取向参数进行了详细的讨论.结果表明,O++H2→OH++H,O++DH→OD++H和O++TH→OT++H反应的立体动力学性质对体系的质量因数非常敏感.  相似文献   

4.
Direct dynamics calculations have been performed for three reactions: C3H8 + H → i-C3H7 + H2, C3H8 + H → n-C3H7 + H2, and C2H3 + O2 → HCO + CH2O. The fraction of the population for the radical products that promptly dissociates is computed. The results for C3H8 + H are qualitatively similar to previous results for C3H8 + OH, but the new results exhibit a slightly higher branching fraction for prompt dissociation products, owing to the fact that a greater fraction of the internal energy in the transition state ends up in the radical. For C2H3 + O2 → HCO + CH2O, the fraction of HCO that promptly dissociates is in excess of 99%. Consequently, the main product for C2H3 + O2 at lower temperatures should be written as H + CO + CH2O and not HCO + CH2O. These results are then compared with four previous systems: CH2O + H → HCO + H2, CH2O + OH → HCO + H2O, C3H8 + OH → i-C3H7 + H2O, and C3H8 + OH → n-C3H7 + H2O. Based upon these seven system, several statistical models are presented. The goal of these statistical models is to predict the fraction of the transition state energy that ends up in the rovibrationally excited radical. On average, these statistical models provide an excellent prediction of product energy distribution. Consequently, these models can be used instead of costly trajectory simulations for predicting prompt radical dissociation for larger species.  相似文献   

5.
The characteristics of possible chain explosive hydrogen burning reactions in an oxidizing medium are calculated on the potential energy surface. Specifically, reactions H2 + O2 → H2O + O, H2 + O2 → HO2 + H, and H2 + O2 → OH + OH are considered. Special attention is devoted to the production of a pair of fast highly reactive OH radicals. Because of the high activation threshold, this reaction is often excluded from the known kinetic scheme of hydrogen burning. However, a spread in estimates of kinetic characteristics and a disagreement between theoretical predictions with experimental results suggest that the kinetic scheme should be refined.  相似文献   

6.
The reaction between CH3 and O2(a1?g) is crucial to understand the effects of electronically excited oxygen in plasma-assisted combustion of methane and other hydrocarbons. In the present work, multireference quantum chemical methods were used to investigate the potential energy surface of CH3 + O2(a1?g). The RRKM/master equation simulation was employed to compute the rate coefficients of various pathways to this reaction over the temperature range of 300–2000 K and a pressure range of 0.1–100 atm. Special attention has been paid to the nonadiabatic transition between the excited state and ground state, which directly leads to a quenching channel from CH3 + O2(a1?g) to CH3 + O2(X3g?). This quenching reaction has been overlooked by previous theoretical and kinetic modeling studies. We also conducted kinetic modeling to examine the effect of this reaction on the ignition enhancement of methane oxidation. Although the channel of CH3 + O2(a1?g) quenching to CH3 + O2(X3g?) has nonnegligible rate constants comparing with other reaction channels, modeling result with the inclusion of 5% O2(a1?g) in molecular oxygen shows that the titled reactions shorten the ignition delay time of methane by more than twenty times at 900 K, 1 atm. The ignition enhancement is mainly from the chain branching channels to CH2O + OH and CH3O + O which has been greatly promoted by excess energy from O2(a1?g). The present study uncovers the kinetic mechanism of this nonadiabatic reaction and provides reasonable rate coefficients for further kinetic modeling of plasma-assisted combustion of methane and other hydrocarbons.  相似文献   

7.
The temporal variation of chemiluminescence emission from OH?(A2 Σ +) and CH?(A2 Δ) in reacting Ar-diluted H2/O2/CH4, C2H2/O2 and C2H2/N2O mixtures was studied in a shock tube for a wide temperature range at atmospheric pressures and various equivalence ratios. Time-resolved emission measurements were used to evaluate the relative importance of different reaction pathways. The main formation channel for OH? in hydrocarbon combustion was studied with CH4 as benchmark fuel. Three reaction pathways leading to CH? were studied with C2H2 as fuel. Based on well-validated ground-state chemistry models from literature, sub-mechanisms for OH? and CH? were developed. For the main OH?-forming reaction CH+O2=OH?+CO, a rate coefficient of k 2=(8.0±2.6)×1010 cm3?mol?1?s?1 was determined. For CH? formation, best agreement was achieved when incorporating reactions C2+OH=CH?+CO (k 5=2.0×1014 cm3?mol?1?s?1) and C2H+O=CH?+CO (k 6=3.6×1012exp(?10.9 kJ?mol?1/RT) cm3?mol?1?s?1) and neglecting the C2H+O2=CH?+CO2 reaction.  相似文献   

8.
In this paper, the importance of fluctuations in flow field parameters is studied under MILD combustion conditions. In this way, a turbulent non-premixed CH4+H2 jet flame issuing into a hot and deficient co-flow air is modeled using the RANS Axisymmetric equations. The modeling is carried out using the EDC model to describe the turbulence-chemistry interaction. The DRM-22 reduced mechanism and the GRI2.11 full mechanism are used to represent the chemical reactions of H2/methane jet flame. Results illustrate that although the fluctuations in temperature field are small and the reaction zone volume are large in the MILD regime, the fluctuations in temperature and species concentrations are still effective on the flow field. Also, inappropriate dealing with the turbulence effect on chemistry leads to errors in prediction of temperature up to 15% in the present flame. By decreasing of O2 concentration of hot co-flow air, the effect of fluctuations in flow field parameters on flame characteristics are still significant and its effect on species reaction rates does not decrease. On the other hand, although decreasing of jet inlet Reynolds number at constant inlet turbulence intensity addresses to smaller fluctuations in flow filed, it does not lead to lower the effect of turbulence on species distribution and temperature field under MILD combustion conditions.  相似文献   

9.
Laminar flame speeds of ammonia with oxygen-enriched air (oxygen content varying from 21 to 30 vol.%) and ammonia-hydrogen-air mixtures (fuel hydrogen content varying from 0 to 30 vol.%) at elevated pressure (1–10 bar) and temperature (298–473 K) were determined experimentally using a constant volume combustion chamber. Moreover, ammonia laminar flame speeds with helium as an inert were measured for the first time. Using these experimental data along with published ones, we have developed a newly compiled kinetic model for the prediction of the oxidation of ammonia and ammonia-hydrogen blends in freely propagating and burner stabilized premixed flames, as well as in shock tubes, rapid compression machines and a jet-stirred reactor. The reaction mechanism also considers the formation of nitrogen oxides, as well as the reduction of nitrogen oxides depending on the conditions of the surrounding gas phase. The experimental results from the present work and the literature are interpreted with the help of the kinetic model derived here. The experiments show that increasing the initial temperature, fuel hydrogen content, or oxidizer oxygen content causes the laminar flame speed to increase, while it decreases when increasing the initial pressure. The proposed kinetic model predicts the same trends than experiments and a good agreement is found with measurements for a wide range of conditions. The model suggests that under rich conditions the N2H2 formation path is favored compared to stoichiometric condition. The most important reactions under rich conditions are: NH2+NH=N2H2+H, NH2+NH2=N2H2+H2, N2H2+H=NNH+H2 and N2H2+M=NNH+H+M. These reactions were also found to be among the most sensitive reactions for predicting the laminar flame speed for all the cases investigated.  相似文献   

10.
Understanding the ion chemistry in flames is crucial for developing ion sensitive technologies for controlling combustion processes. In this work, we measured the spatial distributions of positive ions in atmospheric-pressure burner-stabilized premixed flames of ethylene/oxygen/argon mixtures in a wide range of equivalence ratios ϕ = 0.4÷1.5. A flame sampling molecular beam system coupled with a quadrupole mass spectrometer was used to obtain the spatial distributions of cations in the flames, and a high mass resolution time-of-flight mass spectrometer was utilized for the identification of the cations having similar m/z ratios. The measured profiles of the flame ions were corrected for the contribution of hydrates formed during sampling in the flames slightly upstream the flame reaction zone. We also proposed an updated ion chemistry model and verified it against the experimental profiles of the most abundant cations in the flames. Our model is based on the kinetic mechanism available in the literature extended with the reactions for C3H5+ cation. Highly accurate W2-F12 quantum chemical calculations were used to obtain a reliable formation enthalpy of C3H5+. The model was found to reproduce properly the measured relative abundance of the key oxygenated cations (viz., CH5O+, C2H3O+) in the whole range of equivalence ratios employed, and the C3H5+ cation abundance in the richest flame with ϕ=1.5, but significantly underpredicts the relative mole fraction of C3H3+, which becomes a key species under fuel-rich conditions. Apart from this, several aromatic and cyclic CxHy cations dominating under fuel-rich conditions were identified. We also considered the most important directions for the further refinement of the mechanism.  相似文献   

11.
Having a better understanding of polycyclic aromatic hydrocarbon (PAH) formation under flame conditions contributes to optimizing the fuel reforming process, where soot poisons the downstream catalyst. In this work, the phenyl + 1,3-Butadiyne reaction is systematically investigated to examine its contribution to naphthalene formation. The reaction potential energy surfaces were calculated using DFT/M06–2X/cc-pvtz and G4 methods. The temperature- and pressure-dependent reaction rate constants were calculated using RRKM theory with solving master equation. The results revealed that 2-naphthyl could be directly formed by phenyl + 1,3-Butadiyne reaction. With H assistance, naphthalene could be formed by the pathway of phenyl + 1,3-Butadiyne → C6H5CHCCCH (+H) → C6H5CHCHCCH (+H) →naphthalene +H. The proposed pathway is kinetically favorable, and featured by relatively low energy barrier. The importance of the proposed pathway reaction was confirmed in a premixed and a diffusion C2H4/O2/Ar flame simulations, where the enhancement of naphthalene by the investigated reactions is notable. The mole fraction of A2 is promoted by a factor of 10% in premix C2H4/O2/Ar flame and 30% in C2H4/O2/Ar counterflow flame, bringing the prediction results closer to the experimental results. The relative contribution of different reaction route to A2 formation is evaluated for HACA, cyclopentadienyl radical-cyclopentadienyl radical, phenyl-vinylacetylene[1], benzyl radical-propargyl radical, indene-CH2 and phenyl-1,3-Butadiyne routes in premixed and diffusion C2H4/O2/Ar flames. This work suggests that the PAH growth by 1,3-Butadiyne addition reaction is an effective pathway for A2 formation, which should be considered in future PAH mechanism.  相似文献   

12.
Co-firing methane (CH4) and ammonia (NH3) has attracted growing concerns as a feasible greenhouse gas reduction strategy in gas turbine-based power generation, which raises the need to better understand the interaction of methane and nitric oxide (NO) under flame conditions. In this work, laminar flame propagation of CH4/NO mixtures at initial pressure (Pu) of 1 atm, initial temperature (Tu) of 298 K and equivalence ratios of 0.4–1.8 was experimentally investigated using a constant-volume combustion vessel. Laminar burning velocities (LBVs) and Markstein lengths were experimentally determined. A kinetic model of CH4/NO combustion was developed with rate constants of several important reactions updated, presenting reasonable predictions on the measured LBVs of CH4/NO mixtures. The modeling analyses reveal that the reduction of NO can proceed through two mechanisms, i.e. the hydrocarbon NO reduction mechanism and non-hydrocarbon NO reduction mechanism. Among the two mechanisms, the non-hydrocarbon NO reduction mechanism which includes reactions NO+H = N+OH, NO+O = N + O2 and NO+N = N2+O has a higher contribution to NO reduction at the equivalence ratio of 0.6, while the hydrocarbon NO reduction mechanism with hydrocyanic acid (HCN) as the key intermediate plays a more important role at the equivalence ratio of 1.8. NO+H = N+OH and CH3+NOHCN+H2O are found to be the two most sensitive reactions to promote the flame propagation, while the LBVs measured in this work are demonstrated to provide strong constraint for these reactions. Furthermore, previous CH4/O2/NO oxidation data measured in flow reactor and rapid compression machine were also simulated, which provides extended validation of the present model over wider conditions.  相似文献   

13.
On the basis of a multi-step kinetic mechanism for flame inhibition by organophosphorus compounds including more than 200 reactions, a skeletal mechanism for flame inhibition by trimethylphosphate was developed. The mechanism consists of 22 irreversible elementary reactions, involving nine phosphorus-containing species. Selection of the crucial steps was performed by analysing P-element fluxes from species to species and by calculating net reaction rates of phosphorus-involving reactions versus the flames zone. The developed mechanism was validated by comparing the modelling results with the measured and simulated (using the starting initial mechanism) speed and the chemical structure of H2/O2, CH4/O2 and syngas/air flames doped with trimethylphosphate. The mechanism was shown to satisfactorily predict the speed of H2/O2/N2 flames with various dilution ratios, CH4/air and syngas/air flames doped with trimethylphosphate. The skeletal mechanism was also shown to satisfactorily predict the spatial variation of H and OH radicals and the final phosphorus-containing products of the inhibitor combustion. Further reduction of the skeletal mechanism without modification of the rate constants recommended in the starting mechanism was shown to result in noticeable disagreement of the flame speed and structure.  相似文献   

14.
A theoretical model is proposed for the chemical and vibrational kinetics of hydrogen oxidation based on consistent accounting of the vibrational non-equilibrium of the HO2 radical that forms as a result of the bimolecular recombination H+O2 → HO2. In the proposed model, the chain branching H+O2 = O+OH and inhibiting H+O2+M = HO2+M formal reactions are treated (in the terms of elementary processes) as a single multi-channel process of forming, intramolecular energy redistribution between modes, relaxation, and unimolecular decay of the comparatively long-lived vibrationally excited HO2 radical, which is able to react and exchange energy with the other components of the mixture. The model takes into account the vibrational non-equilibrium of the starting (primary) H2 and O2 molecules, as well as the most important molecular intermediates HO2, OH, O2(1Δ), and the main reaction product H2O. It is shown that the hydrogen–oxygen reaction proceeds in the absence of vibrational equilibrium, and the vibrationally excited HO2(v) radical acts as a key intermediate in a fundamentally important chain branching process and in the generation of electronically excited species O2(1Δ), O(1D), and OH(2Σ+). The calculated results are compared with the shock tube experimental data for strongly diluted H2–O2 mixtures at 1000 < T < 2500 K, 0.5 < p < 4 atm. It is demonstrated that this approach is promising from the standpoint of reconciling the predictions of the theoretical model with experimental data obtained by different authors for various compositions and conditions using different methods. For T < 1500 K, the nature of the hydrogen–oxygen reaction is especially non-equilibrium, and the vibrational non-equilibrium of the HO2 radical is the essence of this process. The quantitative estimation of the vibrational relaxation characteristic time of the HO2 radical in its collisions with H2 molecules has been obtained as a result of the comparison of different experimental data on induction time measurements with the relevant calculations.  相似文献   

15.
In order to study the combustion chemistry of carboxyl functionality, the laminar burning velocity of acetic acid/air and propanoic acid/air mixtures was investigated in a high-pressure constant-volume cylindrical combustion vessel at 423 K, 1 atm and equivalence ratios of 0.7–1.4. Experimental results reveal that the flame propagation of propanoic acid flame is much faster than that of acetic acid flame, especially under rich conditions, and the laminar burning velocity of propanoic acid/air mixtures peaks at richer conditions than that of acetic acid. The present theoretical calculations for the isomerization and decomposition of propanoic acid radicals indicate that the primary radical products are HOCO, H and C2H5, while those in acetic acid flame are CH3 and OH based on previous studies. A kinetic model of the two acids was developed mainly based on previous and the present theoretical calculation results. It could reasonably capture the measured laminar burning velocities of acetic acid/air and propanoic acid/air mixtures in this work, as well as the previous experimental data in literature. Based on the present model, CH3- and ketene-related pathways play an important role in acetic acid flames. Under rich conditions, ketene is mostly converted to CH3 via CH2CO+HCH3+CO, and the chain-termination reaction of CH3+H(+M)=CH4(+M) is enhanced, which strongly inhibits the propagation of rich acetic acid flames. In contrast, C2H5 and ethylene chemistry play an important role in propanoic acid flames. Rich conditions promote the decomposition of C2H5, yielding ethylene and H, which can facilitate the flame propagation. This can explain the shift of the peak laminar burning velocity of propanoic acid/air mixtures towards a slightly richer condition compared with that of acetic acid/air mixtures.  相似文献   

16.
The dissociation CO2(X1Σ) + M → CO(X1Σ) + O(3P) + M and recombination CO(X1Σ) + O(3P) + M → CO2(X1Σ) + M processes are considered with the spin—orbit coupling taken into account in the ground and several excited states of the CO2 molecule. Because of the specific features of mutual position of potential energy surfaces of the CO2 molecule in the ground and several excited states and the large values of spin—orbit interaction matrix elements, which causes the quantum nonadiabatic transition of the molecule from one state to another, these processes become effectively spin-allowed and the rate constants for the nonadiabatic reactions have large values. The proposed dissociation and recombination mechanisms include reactions involving singlet—triplet crossings.  相似文献   

17.
A detailed comparison has been conducted between chemiluminescence (CL) species profiles of OH?, CH?, and C2 ?, obtained experimentally and from detailed flame kinetics modeling, respectively, of atmospheric pressure non-premixed flames formed in the forward stagnation region of a fuel flow ejected from a porous cylinder and an air counterflow. Both pure methane and mixtures of methane with hydrogen (between 10 and 30 % by volume) were used as fuels. By varying the air-flow velocities methane flames were operated at strain rates between 100 and 350 s?1, while for methane/hydrogen flames the strain rate was fixed at 200 s?1. Spatial profiles perpendicular to the flame front were extracted from spectrograms recorded with a spectrometer/CCD camera system and evaluating each spectral band individually. Flame kinetics modeling was accomplished with an in-house chemical mechanism including C1–C4 chemistry, as well as elementary steps for the formation, removal, and electronic quenching of all measured active species. In the CH4/air flames, experiments and model results agree with respect to trends in profile peak intensity and position. For the CH4/H2/air flames, with increasing H2 content in the fuel the experimental CL peak intensities decrease slightly and their peak positions shift towards the fuel side, while for the model the drop in mole fraction is much stronger and the peak positions move closer to the fuel side. For both fuel compositions the modeled profiles peak closer to the fuel side than in the experiments. The discrepancies can only partly be attributed to the limited attainable spatial resolution but may also necessitate revised reaction mechanisms for predicting CL species in this type of flame.  相似文献   

18.
The laminar burning velocities (LBVs) and cellular instability of 2-methyltetrahydrofuran (2-MTHF) were investigated at the unburned temperature of 423 K and pressures from 1 to 10 atm in a cylindrical constant-volume vessel. The LBVs of 2-MTHF/air flame exhibit a notably dropping with increasing pressure. The cellular instability analysis indicates that the critical flame radius of 2-MTHF/air mixture monotonically increases with increasing pressure and the flame surface suffers more badly cellularity under higher pressures. The critical flame radius exhibits non-monotonic variation versus ? and the most unstable flames appear at ? ≈ 1.3. It is observed that the measured Markstein length of 2-MTHF/air mixture decreases with increasing ? and Pu, leading to an earlier formation of wrinkling and cracks with respect to preferential-diffusional instability. Further investigation found that by using a mixture of 14.2% oxygen with 85.8% helium in place of air as bath gas at 10 atm can effectively suppress the cellular instability. Two recently developed models were used to simulate the experimental results and explore the chemical kinetic effects on LBV. Reaction path analysis reveals that the most consumption of 2-MTHF/air at stoichiometric conditions is through the abstraction of H-atom to form radical C5H9O-5. While the competitiveness of decomposition by CC scission yielding CH3 and tetrahydrofuran radical is relatively weak. Sensitivity analysis illustrates that small-species reactions show a controlling effect on LBV. The increasing pressure leads to an evident increase in the sensitivity coefficient of the recombination reaction H + O2 (+M)=HO2 (+M). The reduction of H atom concentration will cause competition to the initiation reaction H + O2 = O+OH. This could lower the overall oxidation rate and reduce the burning velocity.  相似文献   

19.
We have applied density functional theory calculations to study the reactions of NH2 + CnHn (n = 20, 40, 50, 60, 70 and 80). Due to the hard curvature in C20 cage, the NH2? + C20H20 → NH3 + C20H19? reaction is nearly thermoneutral with a high potential barrier height. For the CnHn fulleranes with n > 20 the transition states appear earlier on the reaction paths, as can be anticipated for exothermic reactions. Using the spherical excess parameter, we distinguished different curvatures on the surfaces of fullerane cages. The reaction enthalpies ΔH°298 and potential barrier heights ΔETS of the considered reactions indicate good correlation with the values of ?i parameter, showing an upward trend with the curvature increasing at carbon sites. We have also investigated the H-abstraction of the chemical derivatives of the C20H20 cage (C20H19–CH3, C20H19–CH2CH3 and C20H19–CH2CH2CH3) in comparison to the corresponding isolated alkanes (CH4, C2H6 and C3H8). Overall, it could be inferred that the H-abstraction from the primary and secondary C–H bonds of isolated alkanes could occur more easily than fullarane derivatives.  相似文献   

20.
NH is a key short-lived radical involved in the prompt-NO formation. Quantification of NH is thus particularly important for testing the NO kinetic mechanisms. However, quantitative measurements of native NH in hydrocarbon/oxygen/nitrogen flames remain very scarce. Therefore, in this work, the mole fractions of native NH were obtained using a combination of laser-based diagnostics; Laser Induced Fluorescence (LIF) and Cavity Ring-Down Spectroscopy (CRDS). The NH species was probed after exciting the transition R1(6) in the A3Π-X3Σ? (0-0) system at 333.9?nm. The mole fraction profiles of NH were successfully obtained in premixed low-pressure flames of CH4/O2/N2 and C2H2/O2/N2 at two equivalence ratios of 1.00 and 1.25. The estimated detection limit for the NH radical was around 4.5?×?108 molecule cm?3 (i.e. 2 ppb in mole fraction at 1600?K), which is nearly 2 orders of magnitude lower than previous values reported in the literature. These new experimental results were compared with predictions by a recently developed NO model (namely NOMecha2.0). In the case of the CH4 flames, a satisfying agreement between the experiment and model was observed. However, in the case of the C2H2 flames, some discrepancies were observed. Model analysis has highlighted the importance of the HCCO radicals in the NH formation through the HCNO→HNCO→NH2 reactions pathway. Modification of the rate constant values of the reactions C2H2+?O and HCCO?+?O2, which are key reactions for both the acetylene laminar flame speed and the HCCO predictions, has enabled the model to satisfactorily predict the experimental NH and NO profiles also in the C2H2 flames.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号