首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Usually premixed flame propagation and laminar burning velocity are studied for mixtures at normal or elevated temperatures and pressures, under which the ignition delay time of the premixture is much larger than the flame resistance time. However, in spark-ignition engines and spark-assisted compression ignition engines, the end-gas in the front of premixed flame is at the state that autoignition might happen before the mixture is consumed by the premixed flame. In this study, laminar premixed flames propagating into an autoigniting dimethyl ether/air mixture are simulated considering detailed chemistry and transport. The emphasis is on the laminar burning velocity of autoigniting mixtures under engine-relevant conditions. Two types of premixed flames are considered: one is the premixed planar flame propagating into an autoigniting DME/air without confinement; and the other is premixed spherical flame propagating inside a closed chamber, for which four stages are identified. Due to the confinement, the unburned mixture is compressed to high temperature and pressure close to or under engine-relevant conditions. The laminar burning velocity is determined from the constant-volume propagating spherical flame method as well as PREMIX. The laminar burning velocities of autoigniting DME/air mixture at different temperatures, pressures, and autoignition progresses are obtained. It is shown that the first-stage and second-stage autoignition can significantly accelerate the flame propagation and thereby greatly increase the laminar burning velocity. When the first-stage autoignition occurs in the unburned mixture, the isentropic compression assumption does not hold and thereby the traditional method cannot be used to calculate the laminar burning velocity. A modified method without using the isentropic compression assumption is proposed. It is shown to work well for autoigniting mixtures. Besides, a power law correlation is obtained based on all the laminar burning velocity data. It works well for mixtures before autoignition while improvement is still needed for mixtures after autoignition.  相似文献   

2.
The outwardly propagating spherical flame (OPF) method is popularly used to measure the laminar flame speed (LFS). Recently, great efforts have been devoted to improving the accuracy of the LFS measurement from OPF. In the OPF method, several assumptions are made. For examples, the burned gas is assumed to be static and in chemical equilibrium. However, these assumptions may not be satisfied under certain conditions. Here we consider low-pressure and super-adiabatic propagating spherical flames, for which chemical non-equilibrium exists and the burned gas may not be static. The objective is to assess the chemical non-equilibrium effects on the accuracy of LFS measurement from the OPF method. Numerical simulations considering detailed chemistry and transport are conducted. Stoichiometric methane/air flames at sub-atmospheric pressures and methane/oxygen flames at different equivalence ratios are considered. At low pressures, broad heat release zone is observed and the burned gas cannot quickly reach the adiabatic flame temperature, indicating the existence of chemical non-equilibrium of burned gas. Positive flow in the burned gas is identified and it is shown to become stronger at lower initial pressure. Consequently, the LFS measurement from OPF at low pressures is not accurate if the burned gas is assumed to be static and at chemical equilibrium. For super-adiabatic spherical flames, the burned gas speed is found to be negative due to the local temperature overshoot at the flame front. Such negative speed of burned gas can also reduce the accuracy of LFS measurement. It is recommended that the direct method measuring both flame propagation speed and flow speed of unburned gas should be used to determine the LFS at low pressures or for mixtures with super-adiabatic flame temperature.  相似文献   

3.
In this paper, an experimental and numerical investigation of premixed methane/air flame dynamics in a closed combustion vessel with a thin obstacle is described. In the experiment, high-speed video photography and a pressure transducer are used to study the flame shape changes and pressure dynamics. In the numerical simulation, four sub-grid scale viscosity models and three sub-grid scale combustion models are evaluated for their individual prediction compared with the experimental data. High-speed photographs show that the flame propagation process can be divided into five stages: spherical flame, finger-shaped flame, jet flame, mushroom-shaped flame and bidirectional propagation flame. Compared with the other sub-grid scale viscosity models and sub-grid scale combustion models, the dynamic Smagorinsky–Lilly model and the power-law flame wrinkling model are better able to predict the flame behaviour, respectively. Thus, coupling the dynamic Smagorinsky–Lilly model and the power-law flame wrinkling model, the numerical results demonstrate that flame shape change is a purely hydrodynamic phenomenon, and the mushroom-shaped flame and bidirectional propagation flame are the result of flame–vortex interaction. In addition, the transition from “corrugated flamelets” to “thin reaction zones” is observed in the simulation.  相似文献   

4.
An initially laminar premixed flame front accelerates extremely fast and may even trigger a detonation when propagating in a semi-open obstructed channel (one end of the channel is closed; the flame is ignited at the closed end and moves towards the open one). However, industrial and laboratory conduits oftentimes have both ends open, or vented, with a flame ignited at one of these ends. The latter constitutes the focus of the present work. Specifically, premixed flame propagation through a comb-shaped array of obstacles, in-built in a channel with both ends open, is studied by means of computational simulation of the reacting flow equations with fully-compressible hydrodynamics and an Arrhenius chemical kinetics. The parametric study includes various blockage ratios and spacing as well as the thermal expansion ratios, with oscillations of the burning rate observed in the majority of the cases, which conceptually differs from fast flame acceleration in semi-open channels. Such a difference is devoted to the fact that while the entire flame-generated jet-flow is pushed towards a single exit in a semi-open channel, in a channel with two ends open, this jet-flow is distributed between the upstream and downstream flows, thereby moderating flame propagation. The flame oscillations are nonlinear in all cases where they are observed. The oscillation period grows with the blockage ratio but decreases with the thermal expansion. The present results also support the recent experiments, modeling and theory of flames in obstructed channels with both ends open, which all yielded steady or quasi-steady flame propagation prior to the onset of flame acceleration. Indeed, the present oscillations can be treated as the fluctuations around a quasi-steady solution.  相似文献   

5.
This paper presents an assessment of Large Eddy Simulations (LES) in calculating the structure of turbulent premixed flames propagating past solid obstacles. One objective of the present study is to evaluate the LES simulations and identify the drawbacks in accounting the chemical reaction rate. Another objective is to analyse the flame structure and to calculate flame speed, generated overpressure at different time intervals following ignition of a stoichiometric propane/air mixture. The combustion chamber has built-in repeated solid obstructions to enhance the turbulence level and hence increase the flame propagating speed. Various numerical tests have also been carried out to determine the regimes of combustion at different stages of the flame propagation. These have been identified from the calculated results for the flow and flame characteristic parameters. It is found that the flame lies within the ‘thin reaction zone’ regime which supports the use of the laminar flamelet approach for modelling turbulent premixed flames. A submodel to calculate the model coefficient in the algebraic flame surface density model is implemented and examined. It is found that the LES predictions are slightly improved owing to the calculation of model coefficient by using submodel. Results are presented and discussed in this paper are for the flame structure, position, speed, generated pressure and the regimes of combustion during all stages of flame propagation from ignition to venting. The calculated results are validated against available experimental data.  相似文献   

6.
The present study aims to clarify the effects of turbulence intensity and coal concentration on the spherical turbulent flame propagation of a pulverized coal particle cloud. A unique experimental apparatus was developed in which coal particles can be dispersed homogeneously in a turbulent flow field generated by two fans. Experiments on spherical turbulent flame propagation of pulverized coal particle clouds in a constant volume spherical chamber in various turbulence intensities and coal concentrations were conducted. A common bituminous coal was used in the present study. The flame propagation velocity was obtained from an analysis of flame propagation images taken using a high-speed camera. It was found that the flame propagation velocity increased with increasing flame radius. The flame propagation velocity increases as the turbulence intensity increases. Similar trends were observed in spherical flames using gaseous fuel. The coal concentration has a weak effect on the flame propagation velocity, which is unique to pulverized coal combustions in a turbulent field. These are the first reports of experimental results for the spherical turbulent flame propagation behavior of pulverized coal particle clouds. The results obtained in the present study are obviously different from those of previous pulverized coal combustion studies and any other results of gaseous fuel combustion research.  相似文献   

7.
A premixed propane–air flame stabilised on a triangular bluff body in a model jet-engine afterburner configuration is investigated using large-eddy simulation (LES). The reaction rate source term for turbulent premixed combustion is closed using the transported flame surface density (TFSD) model. In this approach, there is no need to assume local equilibrium between the generation and destruction of subgrid FSD, as commonly done in simple algebraic closure models. Instead, the key processes that create and destroy FSD are accounted for explicitly. This allows the model to capture large-scale unsteady flame propagation in the presence of combustion instabilities, or in situations where the flame encounters progressive wrinkling with time. In this study, comprehensive validation of the numerical method is carried out. For the non-reacting flow, good agreement for both the time-averaged and root-mean-square velocity fields are obtained, and the Karman type vortex shedding behaviour seen in the experiment is well represented. For the reacting flow, two mesh configurations are used to investigate the sensitivity of the LES results to the numerical resolution. Profiles for the velocity and temperature fields exhibit good agreement with the experimental data for both the coarse and dense mesh. This demonstrates the capability of LES coupled with the TFSD approach in representing the highly unsteady premixed combustion observed in this configuration. The instantaneous flow pattern and turbulent flame behaviour are discussed, and the differences between the non-reacting and reacting flow are described through visualisation of vortical structures and their interaction with the flame. Lastly, the generation and destruction of FSD are evaluated by examining the individual terms in the FSD transport equation. Localised regions where straining, curvature and propagation are each dominant are observed, highlighting the importance of non-equilibrium effects of FSD generation and destruction in the model afterburner.  相似文献   

8.
A tribrachial (or triple) flame is one kind of edge flame that can be encountered in nonpremixed mixing layers, consisting of a lean and a rich premixed flame wing together with a trailing diffusion flame all extending from a single point. The flame could play an important role on the characteristics of various flame behaviors including lifted flames in jets, flame propagation in two-dimensional mixing layers, and autoignition fronts. The structure of tribrachial flame suggests that the edge is located along the stoichiometric contour in a mixing layer due to the coexistence of all three different types of flames. Since the edge has a premixed nature, it has unique propagation characteristics. In this review, the propagation speed of tribrachial flames will be discussed for flames propagating in mixing layers, including the effects of concentration gradient, velocity gradient, and burnt gas expansion. Based on the tribrachial edge structure observed experimentally in laminar lifted flames in jets, the flame stabilization characteristics including liftoff height, reattachment, and blowout behaviors and their buoyancy-induced instability will be explained. Various effects on liftoff heights in both free and coflow jets including jet velocity, the Schmidt number of fuel, nozzle diameter, partial premixing of air to fuel, and inert dilution to fuel are discussed. Implications of edge flames in the modeling of turbulent nonpremixed flames and the stabilization of turbulent lifted flames in jets are covered.  相似文献   

9.
Laminar premixed cool flames, induced by the coupling of low-temperature chemistry and convective-diffusive transport process, have recently attracted extensive interest in combustion and engine research. In this work, numerical simulations have been conducted using a recently developed open-source reacting flow platform reactingFOAM-SCT, to investigate the minimum ignition energy (MIE) and propagation dynamics of premixed cool flames in a 1D spherical coordinate. Results have shown that when ignition energy is below the MIE of regular hot flames, a class of cool flames could be initiated, which allow much wider flammability limits, both lean and rich, compared to hot flames. Furthermore, the overall cool flame propagation dynamics exhibit intrinsic similarity to those of hot flames, in that, they begin with an ignition kernel propagation regime, followed by two transition regimes, and eventually reach a normal flame propagation regime. However, a spherical expanding cool flame responds completely differently to stretch. Specifically, a regular outwardly propagating hot spherical flame accelerates with increasing stretch rate when the mixture Le < 1 and decelerates when Le > 1. However, it is found that a cool flame always tends to decelerate with increasing stretch rate regardless of mixture composition, exhibiting unique flame aerodynamic characteristic. This research discovers novel features of premixed cool flame initiation and propagation dynamics and sheds light on flame transition, spark-ignition system design, and advanced engine combustion control.  相似文献   

10.
Three-dimensional (3D) unsteady Reynolds-averaged Navier–Stokes simulations of a spark-ignited turbulent methane/air jet flame evolving from ignition to stabilisation are conducted for different jet velocities. A partially premixed combustion model is used involving a correlated joint probability density function and both premixed and non-premixed combustion mode contributions. The 3D simulation results for the temporal evolution of the flame's leading edge are compared with previous two-dimensional (2D) results and experimental data. The comparison shows that the final stabilised flame lift-off height is well predicted by both 2D and 3D computations. However, the transient evolution of the flame's leading edge computed from 3D simulation agrees reasonably well with experiment, whereas evident discrepancies were found in the previous 2D study. This difference suggests that the third physical dimension plays an important role during the flame transient evolution process. The flame brush's leading edge displacement speed resulting from reaction, normal and tangential diffusion processes are studied at different typical stages after ignition in order to understand the effect of the third physical dimension further. Substantial differences are found for the reaction and normal diffusion components between 2D and 3D simulations especially in the initial propagation stage. The evolution of reaction progress variable scalar gradients and its interaction with the flow and mixing field in the 3D physical space have an important effect on the flame's leading edge propagation.  相似文献   

11.
A numerical study of one-dimensional n-heptane/air spray flames is presented. The objective is to evaluate the flame propagation speed in the case where droplets evaporate inside the reaction zone with possibly non-zero relative velocity. A Direct Numerical Simulation approach for the gaseous phase is coupled to a discrete particle Lagrangian formalism for the dispersed phase. A global two-step n-heptane/air chemical mechanism is used. The effects of initial droplet diameter, overall equivalence ratio, liquid loading and relative velocity between gaseous and liquid phases on the laminar spray flame speed and structure are studied. For lean premixed cases, it is found that the laminar flame speed decreases with increasing initial droplet diameter and relative velocity. On the contrary, rich premixed cases show a range of diameters for which the flame speed is enhanced compared to the corresponding purely gaseous flame. Finally, spray flames controlled by evaporation always have lower flame speeds. To highlight the controlling parameters of spray flame speed, approximate analytical expressions are proposed, which give the correct trends of the spray flame propagation speed behavior for both lean and rich mixtures.  相似文献   

12.
As a carbon-free fuel, hydrogen has received significant attention recently since it can help enable low-carbon-economy. Hydrogen has very broad flammability range and very low minimum ignition energy, and thereby there are severe safety concerns for hydrogen transportation and utilization. Cryo-compressed hydrogen is popularly used in practice. Therefore, it is necessary to investigate the combustion properties of hydrogen at extremely low or cryogenic temperatures. This study aims to assess and interpret the effects of cryogenic temperature on premixed hydrogen/air flame propagation and acceleration in a thin closed channel. Different initial temperatures ranging from normal temperature (T0 = 300 K) to cryogenic temperature (T0 = 100 K) are considered. Both one- and two-dimensional hydrogen/air flames are investigated through transient simulations considering detailed chemistry and transport. It is found that when the initial temperature decreases from T0 = 300 K to T0 = 100 K, the expansion ratio and equilibrium pressure both increase substantially while the laminar flame speeds relative to unburned and burned gasses decrease moderately. The one-dimensional flame propagation is determined by laminar flame speed and thereby the combustion duration increases as the initial temperature decreases. However, the opposite trend is found to happen to two-dimensional flame propagation, which is mainly controlled by the flame surface area increase due to the no-slip side wall constraint and flame instability. Based on the change in flame surface area, three stages including the initial acceleration, steady burning and rapid acceleration are identified and investigated. It is demonstrated that the large expansion ratio and high pressure rise at cryogenic temperatures can significantly increase the flame surface area in early stage and promote both Darrieus-Landau instability (hydrodynamic instability) and Rayleigh-Taylor instability in later stage. These two instabilities can substantially increase the flame surface area and thereby accelerate flame propagation in hydrogen/air mixtures at cryogenic temperatures. The present study provides useful insights into the fundamental physics of hydrogen flames at extremely low temperatures, and is closely related to hydrogen safety.  相似文献   

13.
This study investigates the influence of large-scale flow features, including flow structure and velocity magnitude, on the early-burn period variability in a homogenous-charge spark-ignited engine fueled with premixed propane-air mixture. Particle image velocimetry and in-cylinder pressure measurement data from a previous study - were processed to enable simultaneous flow characterization and flame-front tracking as well as apparent heat-release analysis. By combining probability analysis of flame development with conditional sampling of fast and slow early-burn cycles using 10% fuel mass fraction burned, it is shown that an undesirable flow structure produces an asymmetric flame development at the initial flame growth period. This asymmetric flame structure persists through the whole initial-to-turbulent transition period until the flame becomes fully turbulent. The undesirable flow condition is characterized by large-scale convective flows near spark plug rather than flows that lead to increased flame spread in multiple directions. The simultaneous flow and flame characterization enables the quantifications of flame-front propagation speed, unburned fuel-air mixture velocity ahead of flame front and local burning velocity at flame surface. Here the local burning velocity is referred to as laminar or turbulent flame speed. A simplified approach is introduced to derive integrated values for these quantities per crank-angle-degree, enabling the quantitative comparison of the trend-wise difference in these integrated metrics between fast and slow early-burn cycles. It is revealed that for the transition period, the CCV in the velocity magnitude of unburned fuel-air mixture ahead of the flame front accounts for nearly 50% to the variability of flame propagation speed. The burning velocity provides the remaining source of the flame propagation variability in this period. The flame propagation variations in the initial flame growth and fully turbulent periods are smaller than those in the transition period and are primarily dependent on the variability of large-scale flow features.  相似文献   

14.
The peak flame surface density within the turbulent flame brush is central to turbulent premixed combustion models in the flamelet regime. This work investigates the evolution of the peak surface density in spherically expanding turbulent premixed flames with the help of direct numerical simulations at various values of the Reynolds and Karlovitz number. The flames propagate in decaying isotropic turbulence inside a closed vessel. The effects of turbulent transport, transport due to mean velocity gradient, and flame stretch on the peak surface density are identified and characterized with an analysis based on the transport equation for the flame surface density function. The three mechanisms are governed by distinct flow time scales; turbulent transport by the eddy turnover time, mean transport by a time scale related to the pressure rise in the closed chamber, and flame stretch by the Kolmogorov time scale. Appropriate scaling of the terms is proposed and shown to collapse the data despite variations in the dimensionless groups. Overall, the transport terms lead to a reduction in the peak value of the surface density, while flame stretch has the opposite effect. In the present configuration, a small imbalance between the two leads to an exponential decay of the peak surface density in time. The dimensionless decay rate is found to be consistent with the evolution of the wrinkling scale as defined in the Bray-Moss-Libby model.  相似文献   

15.
The influence of Soret diffusion on lean premixed flames propagating in hydrogen/air mixtures is numerically investigated with a detailed chemical and transport models at normal and elevated pressure and temperature. The Soret diffusion influence on the one-dimensional (1D) flame mass burning rate and two-dimensional (2D) flame propagating characteristics is analysed, revealing a strong dependency on flame stretch rate, pressure and temperature. For 1D flames, at normal pressure and temperature, with an increase of Karlovitz number from 0 to 0.4, the mass burning rate is first reduced and then enhanced by Soret diffusion of H2 while it is reduced by Soret diffusion of H. The influence of Soret diffusion of H2 is enhanced by pressure and reduced by temperature. On the contrary, the influence of Soret diffusion of H is reduced by pressure and enhanced by temperature. For 2D flames, at normal pressure and temperature, during the early phase of flame evolution, flames with Soret diffusion display more curved flame cells. Pressure enhances this effect, while temperature reduces it. The influence of Soret diffusion of H2 on the global consumption speed is enhanced at elevated pressure. The influence of Soret diffusion of H on the global consumption speed is enhanced at elevated temperature. The flame evolution is more affected by Soret diffusion in the early phase of propagation than in the long run due to the local enrichment of H2 caused by flame curvature effects. The present study provides new insights into the Soret diffusion effect on the characteristics of lean hydrogen/air flames at conditions that are relevant to practical applications, e.g. gas engines and turbines.  相似文献   

16.
We have developed a general theory of non-adiabatic premixed flames that is valid for flames of arbitrary shape that fully accounts for the hydrodynamic and diffusive-thermal processes, and incorporates the effects of volumetric heat losses. The model is used to describe aspects of experimentally observed phenomena of self-extinguishing (SEFs) and self-wrinkling flames (SWFs), in which radiative heat losses play an important role. SEFs are spherical flames that propagate considerable distances in sub-limit conditions before suddenly extinguishing. Our results capture many aspects of this phenomenon including an explicit determination of flame size and propagation speed at quenching. SWFs are hydrodynamically unstable flames in which cells spontaneously appear on the flame surface once the flame reaches a critical size. Our results yield expressions of the critical flame size at the onset of wrinkling and expected cell size beyond the stability threshold. The various possible burning regimes are mapped out in parameter space.  相似文献   

17.
The effects of flow compression and flame stretch on the accurate determination of laminar flame speeds at normal and elevated pressures using propagating spherical flames at constant pressure or constant volume are studied theoretically and numerically. The results show that both the compression-induced flow motion and flame stretch have significant impacts on the accuracy of flame speed determination. For the constant pressure method, a new method to obtain a compression-corrected flame speed (CCFS) for nearly constant pressure spherical bomb experiments is presented. Likewise, for the constant volume method, a technique to obtain a stretch-corrected flame speed (SCFS) at elevated pressures and temperatures is developed. The validity of theoretical results for both constant pressure and constant volume methods is demonstrated by numerical simulations using detailed chemistry for hydrogen/air, methane/air, and propane/air mixtures. It is shown that the present CCFS and SCFS methods not only improve the accuracy of the flame speed measurements significantly but also extend the parameter range of experimental conditions. The results can be used directly in experimental measurements of laminar flame speeds.  相似文献   

18.
In this study, the influence of the negative velocity field formed ahead of an abruptly deformed flame tip on the propagation behaviour of a laminar premixed flame is numerically investigated. A strong deformation in the flame front is induced by imposing a very narrow, in-line pre-heating zone in the unburned region. The simulation is performed under low Mach number approximation by using a multi-scale multi-physics Computational Fluid Dynamics (CFD) solver FrontFlow/Red with one-step finite rate chemistry in order to track the time-dependent flame dynamics. The computed results unveil that the flame front is deformed significantly within a short time due to the narrow in-line pre-heating effect. The flame deformation induces a strong negative velocity field ahead of the deformed flame tip, acting in the direction of propagation, which gives rise to a strong pair vortex. This strong pair vortex interacts with the flame tip and then slides down along the flame surface in the upstream direction during propagation. This flame-vortex interaction causes further deformation in the flame surface in the upstream direction, and consequently, the flame exhibits a wave-like surface, which enhances the flame propagation speed. The auto-generation of a strong pair vortex ahead of the flame front due to the localised thermal input could be applied as one of the methods to control the combustion externally. It is also expected that the results obtained in this study could have a significant impact on the detailed understanding of the local thermo-fluid dynamical interaction process of turbulent combustion in practical combustors.  相似文献   

19.
The structure and dynamics of a turbulent partially premixed methane/air flame in a conical burner were investigated using laser diagnostics and large-eddy simulations (LES). The flame structure inside the cone was characterized in detail using LES based on a two-scalar flamelet model, with the mixture fraction for the mixing field and level-set G-function for the partially premixed flame front propagation. In addition, planar laser induced florescence (PLIF) of CH and chemiluminescence imaging with high speed video were performed through a glass cone. CH and CH2O PLIF were also used to examine the flame structures above the cone. It is shown that in the entire flame the CH layer remains very thin, whereas the CH2O layer is rather thick. The flame is stabilized inside the cone a short distance above the nozzle. The stabilization of the flame can be simulated by the triple-flame model but not the flamelet-quenching model. The results show that flame stabilization in the cone is a result of premixed flame front propagation and flow reversal near the wall of the cone which is deemed to be dependent on the cone angle. Flamelet based LES is shown to capture the measured CH structures whereas the predicted CH2O structure is somewhat thinner than the experiments.  相似文献   

20.
为解决瓦斯输送过程中的爆炸安全问题,探索寻找绿色环保且阻火性能优越的新型抑爆剂,开展了当量比下甲烷-空气预混气体爆炸传播过程中的七氟丙烷抑爆效果研究。实验采用长径比L/D=108的水平管道爆炸特性测试系统,研究了在强点火作用下不同体积分数的七氟丙烷对9.5%甲烷-空气预混气体最大爆炸压力、最大压力上升速率和火焰传播速度的影响。实验结果显示:将2.5 m长的管段作为七氟丙烷抑爆区时,七氟丙烷阻断9.5%甲烷-空气预混气体爆炸火焰传播的最小体积分数为5%;当七氟丙烷的体积分数为1%~4%时,不仅无法阻断爆炸火焰的传播,而且与对照组相比,会使火焰传播速度加快;当七氟丙烷的体积分数为1%~6%时,爆炸源及管道末端处的爆炸压力峰值随着七氟丙烷体积分数的增加而逐渐减小;当七氟丙烷的体积分数为3%时,抑爆区处的爆炸压力峰值与对照组相比增幅为10.9%。  相似文献   

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

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