首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
An experiment is carried out here to investigate the characteristics of the evaporation heat transfer for refrigerants R-134a and R-407C flowing in horizontal small tubes having the same inside diameter of 0.83 or 2.0 mm. In the experiment for the 2.0-mm tubes, the refrigerant mass flux G is varied from 200 to 400 kg/m2 s, imposed heat flux q from 5 to 15 kW/m2, inlet vapor quality xin from 0.2 to 0.8 and refrigerant saturation temperature Tsat from 5 to 15 °C. While for the 0.83-mm tubes, G is varied from 800 to 1500 kg/m2 s with the other parameters varied in the same ranges as those for Di = 2.0 mm. In the study the effects of the refrigerant vapor quality, mass flux, saturation temperature and imposed heat flux on the measured evaporation heat transfer coefficient hr are examined in detail. The experimental data clearly show that both the R-134a and R-407C evaporation heat transfer coefficients increase almost linearly and significantly with the vapor quality of the refrigerant, except at low mass flux and high heat flux. Besides, the evaporation heat transfer coefficients also increase substantially with the rises in the imposed heat flux, refrigerant mass flux and saturation temperature. At low R-134a mass flux and high imposed heat flux the evaporation heat transfer coefficient in the smaller tubes (Di = 0.83 mm) may decline at increasing vapor quality when the quality is high, due to the partial dryout of the refrigerant flow in the smaller tubes at these conditions. We also note that under the same xin, Tsat, G, q and Di, refrigerant R-407C has a higher hr when compared with that for R-134a. Finally, an empirical correlation for the R-134a and R-407C evaporation heat transfer coefficients in the small tubes is proposed.  相似文献   

2.
The current paper presents experimental investigation of nucleate pool boiling of R-134a and R-123 on enhanced and smooth tubes. The enhanced tubes used were TBIIHP and TBIILP for R-134a and R-123, respectively. Pool boiling data were taken for smooth and enhanced tubes in a single tube test section. Data were taken at a saturation temperature of 4.44 °C. Each test tube had an outside diameter of 19.05 mm and a length of 1 m. The test section was water heated with an insert in the water passage. The insert allowed measurement of local water temperatures down the length of the test tube. Utilizing this instrumentation, local heat transfer coefficients were determined at five locations along the test tube. The heat flux range was 2.5–157.5 kW/m2 for the TBIIHP tube and 3.1–73.2 kW/m2 for the TBIILP tube. The resulting heat transfer coefficient range was 4146–23255 W/m2. °C and 5331–25950 W/m2. °C for both tubes, respectively. For smooth tube testing, the heat flux ranges were 7.3–130.7 kW/m2 and 7.5–60.7 kW/m2 for R-134a and R-123, respectively; with resulting heat transfer coefficient ranges of 1798.9–11,379 W/m2. °C and 535.4–3181.8 W/m2. °C. The study provided one of the widest heat flux ranges ever examined for these types of tubes and showed significant structure to the pool boiling curve that had not been traditionally observed. Additionally, this paper presented an investigation of enhanced tubes pool boiling models.  相似文献   

3.
An experimental study on in-tube flow boiling heat transfer of R-134a/R-290/R-600a refrigerant mixture has been carried out under varied heat flux test conditions. The heat transfer coefficients are experimentally measured at temperatures between ?8 and 5 °C for mass flow rates of 3–5 g s?1. Acetone is used as a hot fluid which flows in the outer tube of diameter 28.57 mm while the refrigerant mixture flows in the inner tube of diameters 9.52 and 12.7 mm. By regulating the acetone flow conditions, the heat flux is maintained between 2 and 8 kW/m2 and the pressure of the refrigerant is maintained between 3.2 and 5 bar. The comparison of experimental results with the familiar correlations shows that the correlations over predict the heat transfer coefficients for this mixture when stratified and stratified-wavy flow prevail. Multiple regression technique is used to evolve and modify existing correlations to predict the heat transfer coefficient of the refrigerant mixture. It is found that the modified version of Lavin–Young correlation (1965) predicts the heat transfer coefficient of the considered mixture within an average deviation of ±20.5 %.  相似文献   

4.
A complete solution for boiling phenomena in smooth tubes has been giving as a procedure regarding with the calculation of convective heat transfer coefficient and pressure drop using accurate experimental data validated by flow regime maps and sight glasses on the experimental facility. The experimental study is conducted in order to investigate the effect of operating parameters on flow boiling convective heat transfer coefficient and pressure drop of R134a. The smooth tube having 8.62 mm inner diameter and 1100 mm length is used in the experiments. The effect of mass flux, saturation temperature and heat flux is researched in the range of 290–381 kg/m2 s, 15–22 °C and 10–15 kW/m2, respectively. The experiments revealed that the heat transfer coefficient and pressure drop are significantly affected by mass flux for all tested conditions. Moreover, the experimental results are compared with well-known heat transfer coefficient and frictional pressure drop correlations given in the literature. In addition, 122 number of heat transfer and pressure drop raw experimental data is given for researchers to validate their theoretical models.  相似文献   

5.
An experimental study is carried out to investigate the characteristics of the evaporation heat transfer for different fluids. Namely, pure refrigerants fluids (R22 and R134a), azeotropic and quasi-azeotropic mixtures (R404A, R410A, R507) and zeotropic mixtures (R407C and R417A).The test section is a smooth, horizontal, stainless steel tube (6 mm ID, 6 m length) uniformly heated by the Joule effect. The flow boiling characteristics of the refrigerant fluids are evaluated in 250 different operating conditions. Thus, a data-base of more than 2000 data points is produced.The experimental tests are carried out varying: (i) the refrigerant mass fluxes within the range 200–1100 kg/m2 s; (ii) the heat fluxes within the range 3.50–47.0 kW/m2; (iii) the evaporating pressures within the range 3.00–12.0 bar.In this study, the effect on measured heat transfer coefficient of vapour quality, mass flux, saturation temperature, imposed heat flux, thermo-physical properties are examined in detail.  相似文献   

6.
Flow boiling heat transfer of R-134a refrigerant in a circular mini-channel, 600 mm long with a diameter of 1.75 mm, is investigated experimentally in this study. The test section is a stainless steel tube placed horizontally. Flow pattern and heat transfer coefficient data are obtained for a mass flux range of 200–1000 kg/m2 s, a heat flux range of 1–83 kW/m2 and saturation pressures of 8, 10, and 13 bar. Five different flow patterns including slug flow, throat-annular flow, churn flow, annular flow and annular-rivulet flow are observed and the heat transfer coefficient data for different flow patterns are presented. The heat transfer coefficient increases with increasing heat flux but is mostly independent of mass flux and vapour quality. In addition, it is indicated from the experiments that the higher the saturation pressure, the lower is the heat transfer coefficient. Comparisons of the present data with the existing correlations are also presented.  相似文献   

7.
The boiling heat transfer of refrigerant R-134a flow in horizontal small-diameter tubes with inner diameter of 0.51, 1.12, and 3.1 mm was experimentally investigated. Local heat transfer coefficient and pressure drop were measured for a heat flux ranging from 5 to 39 kW/m2, mass flux from 150 to 450 kg/m2 s, evaporating temperature from 278.15 to 288.15 K, and inlet vapor quality from 0 to 0.2. Flow patterns were observed by using a high-speed video camera through a sight glass at the entrance of an evaporator. Results showed that with decreasing tube diameter, the local heat transfer coefficient starts decreasing at lower vapor quality. Although the effect of mass flux on the local heat transfer coefficient decreased with decreasing tube diameter, the effect of heat flux was strong in all three tubes. The measured pressure drop for the 3.1-mm-ID tube agreed well with that predicted by the Lockhart–Martinelli correlation, but when the inner tube diameter was 0.51 mm, the measured pressure drop agreed well with that predicted by the homogenous pressure drop model. With decreasing tube diameter, the flow inside a tube approached homogeneous flow. The contribution of forced convective evaporation to the boiling heat transfer decreases with decreasing the inner tube diameter.  相似文献   

8.
An experiment is carried out here to investigate the evaporation heat transfer and associated evaporating flow pattern for refrigerant R-134a flowing in a horizontal narrow annular duct. The gap of the duct is fixed at 1.0 and 2.0 mm. In the experiment, the effects of the duct gap, refrigerant vapor quality, mass flux and saturation temperature and imposed heat flux on the measured evaporation heat transfer coefficient hr are examined in detail. For the duct gap of 2.0 mm, the refrigerant mass flux G is varied from 300 to 500 kg/m2 s, imposed heat flux q from 5 to 15 kW/m2, vapor quality xm from 0.05 to 0.95, and refrigerant saturation temperature Tsat from 5 to 15 °C. While for the gap of 1.0 mm, G is varied from 500 to 700 kg/m2 s with the other parameters varied in the same ranges as that for δ = 2.0 mm. The experimental data clearly show that the evaporation heat transfer coefficient increases almost linearly with the vapor quality of the refrigerant and the increase is more significant at a higher G. Besides, the evaporation heat transfer coefficient also rises substantially at increasing q. Moreover, a significant increase in the evaporation heat transfer coefficient results for a rise in Tsat, but the effects are less pronounced in the narrower duct at a low imposed heat flux and a high refrigerant mass flux. Furthermore, the evaporation heat transfer coefficient increases substantially with the refrigerant mass flux except at low vapor quality. We also note that reducing the duct gap causes a significant increase in hr. In addition to the heat transfer data, photos of R-134a evaporating flow taken from the duct side show the change of the dominant two-phase flow pattern in the duct with the experimental parameters. Finally, an empirical correlation for the present measured heat transfer coefficient for the R-134a evaporation in the narrow annular ducts is proposed.  相似文献   

9.
The flow pattern based flow boiling heat transfer and two-phase pressure drop models for CO2, recently developed by Cheng et al. [L. Cheng, G. Ribatski, J. Moreno Quibén, J.R. Thome, New prediction methods for CO2 evaporation inside tubes: Part I – A two-phase flow pattern map and a flow pattern based phenomenological model for two-phase flow frictional pressure drops, Int. J. Heat Mass transfer 51 (2008) 111–124; L. Cheng, G. Ribatski, J.R. Thome, New prediction methods for CO2 evaporation inside tubes: Part II – An updated general flow boiling heat transfer model based on flow patterns, Int. J. Heat Mass transfer 51 (2008) 125–135], have been used to predict the thermal performance of CO2 in a silicon multi-microchannel evaporator (67 parallel channels with a width of 0.223 mm, a height of 0.68 mm and a length of 20 mm) for cooling of a microprocessor. First, some simulation results of CO2 flow boiling heat transfer and two-phase pressure drops in microscale channels are presented. The effects of channel diameter, mass flux, saturation temperature and heat flux on flow boiling heat transfer coefficients and two-phase pressure drops are next addressed. Then, simulations of the base temperatures of the silicon multi-microchannel evaporator using R236fa and CO2 were performed for the following conditions: base heat fluxes from 20 to 100 W/cm2, a mass flux of 987.6 kg/m2s and a saturation temperature of 25 °C. These show that the base temperatures using CO2 are much lower than those using R236fa. Compared to R236fa, CO2 has much higher heat transfer coefficients and lower pressure drops in the multi-microchannel evaporator. However, the operation pressure of CO2 is much higher than that of R236fa. Based on the analysis and comparison, CO2 appears to be a promising coolant for microprocessors at low operating temperatures but also presents a great technological challenge like other new cooling technologies.  相似文献   

10.
《Applied Thermal Engineering》2007,27(16):2713-2726
Enhanced heat transfer surfaces are used in heat exchangers to improve performance and to decrease system volume and cost. In-tube heat transfer enhancement usually takes the form of either micro-fin tubes (of the helical micro-fin or herringbone varieties), or of helical wire inserts. Despite a substantial increase in heat transfer, these devices also cause non-negligible pressure drops.By making use of well-proven flow pattern maps for smooth tubes and the new ones for smooth and enhanced tubes, it is shown from the refrigerant condensation data that flow patterns have a strong influence on heat transfer and pressure drop. This is done for data obtained from in-tube condensation experiments for mass fluxes ranging from 300 to 800 kg/m2 s at a saturation temperature of 40 °C, for refrigerants R-22, R-134a, and R-407C. The flow regimes, pressure drops, heat transfer coefficients, and the overall performance of three different tubes, namely a smooth-, 18° helical micro-fin-, and a herringbone micro-fin tube (each having a nominal diameter of 9.51 mm), are presented and compared to the performance of smooth tubes with helical wire inserts (with pitches of 5 mm, 7.77 mm and 11 mm corresponding to helical angles of 78.2°, 72°, and 65.3°, respectively).  相似文献   

11.
This paper presents the experimental results of condensing heat transfer coefficients and pressure gradients of HC refrigerants (e.g. R-1270, R-290 and R-600a) and R-22 in horizontal double pipe heat exchangers, having two different internal diameters of 12.70 mm and 9.52 mm (OD), respectively. Both the local condensing heat transfer coefficients and pressure drops (inside the tube) of hydrocarbon refrigerants were higher than R-22. The average condensing heat transfer coefficient increased with the mass flux. The experimental heat transfer coefficients agreed with the correlations of Shah, Travis and Cavallini–Zecchin’s to within ±20%. These results can be useful in the design of new heat exchangers involving hydrocarbon refrigerants for future air-conditioning systems.  相似文献   

12.
A physical model of gas–liquid two-phase annular flow regime is presented for predicting the enhanced evaporation heat transfer characteristics in horizontal microfin tubes. The model is based on the equivalence of a periodical distortion of the disturbance wave in the substrate layer. Corresponding to the stratified flow model proposed previously by authors, the dimensionless quantity Fr0 = G/[gdeρv(ρl ? ρv)]0.5 may be used as a measure for determining the applicability of the present theoretical model, which was used to restrict the transition boundary between the stratified-wavy flow and the annular/intermittent flows. Comparison of the prediction of the circumferential average heat transfer coefficient with available experimental data for four tubes and three refrigerants reveals that a good agreement is obtained or the trend is better than that of the previously developed stratified flow model for Fr0 > 4.0 as long as the partial dry out of tube does not occur. Obviously, the developed annular model is applicable and reliable for evaporation in horizontal microfin tubes under the case of high heat flux and high mass flux.  相似文献   

13.
A new flow boiling heat transfer model and a new flow pattern map based on the flow boiling heat transfer mechanisms for horizontal tubes have been developed specifically for CO2. Firstly, a nucleate boiling heat transfer correlation incorporating the effects of reduced pressure and heat flux at low vapor qualities has been proposed for CO2. Secondly, a nucleate boiling heat transfer suppression factor correlation incorporating liquid film thickness and tube diameters has been proposed based on the flow boiling heat transfer mechanisms so as to capture the trends in the flow boiling heat transfer data. In addition, a dryout inception correlation has been developed. Accordingly, the heat transfer correlation in the dryout region has been modified. In the new flow pattern map, an intermittent flow to annular flow transition criterion and an annular flow to dryout region transition criterion have been proposed based on the changes in the flow boiling heat transfer trends. The flow boiling heat transfer model predicts 75.5% of all the CO2 database within ±30%. The flow boiling heat transfer model and the flow pattern map are applicable to a wide range of conditions: tube diameters (equivalent diameters for non-circular channels) from 0.8 to 10 mm, mass velocities from 170 to 570 kg/m2 s, heat fluxes from 5 to 32 kW/m2 and saturation temperatures from −28 to 25 °C (reduced pressures from 0.21 to 0.87).  相似文献   

14.
This paper presents an experimental study on the convective boiling heat transfer and the critical heat flux (CHF) of ethanol–water mixtures in a diverging microchannel with artificial cavities. The results show that the boiling heat transfer and the CHF are significantly influenced by the molar fraction (xm) as well as the mass flux. For the single-phase convection region except for the region near the onset of nucleate boiling with temperature overshoot, the single-phase heat transfer coefficient is independent of the wall superheat and increases with a decrease in the molar fraction. After boiling incipience, the two-phase heat transfer coefficient is much higher than that of single-phase convection. The two-phase heat transfer coefficient shows a maximum in the region of bubbly-elongated slug flow and deceases with a further increase in the wall superheat until approaching a condition of CHF, indicating that the heat transfer is mainly dominated by convective boiling. A flow-pattern-based empirical correlation for the two-phase heat transfer coefficient of the flow boiling of ethanol–water mixtures is developed. The overall mean absolute error of the proposed correlation is 15.5%, and more than 82.5% of the experimental data were predicted within a ±25% error band. The CHF increases from xm = 0–0.1, and then decreases rapidly from xm = 0.1–1 at a given mass flux of 175 kg/m2 s. The maximum CHF is reached at xm = 0.1 due to the Marangoni effect, indicating that small additions of ethanol into water could significantly increase the CHF. On the other hand, the CHF increases with increasing the mass flux at a given molar fraction of 0.1. Moreover, the experimental CHF results are compared with existing CHF correlations of flow boiling of the mixtures in a microchannel.  相似文献   

15.
This paper presents the heat transfer coefficients and the pressure drop measured during HFC refrigerants 236fa, 134a and 410A saturated vapour condensation inside a brazed plate heat exchanger: the effects of saturation temperature (pressure), refrigerant mass flux and fluid properties are investigated. The heat transfer coefficients show weak sensitivity to saturation temperature (pressure) and great sensitivity to refrigerant mass flux and fluid properties. A transition point between gravity controlled and forced convection condensation has been found for a refrigerant mass flux around 20 kg/m2s that corresponds to an equivalent Reynolds number around 1600–1700. At low refrigerant mass flux (Gr < 20 kg/m2s) the heat transfer coefficients are not dependent on mass flux and are well predicted by the Nusselt [20] analysis for vertical surface: the condensation process is gravity controlled. For higher refrigerant mass flux (Gr > 20 kg/m2s) the heat transfer coefficients depend on mass flux and are well predicted by Akers et al. [21] equation: forced convection condensation occurs. In the forced convection condensation region the heat transfer coefficients show a 25–30% increase for a doubling of the refrigerant mass flux.The frictional pressure drop shows a linear dependence on the kinetic energy per unit volume of the refrigerant flow and therefore a quadratic dependence on mass flux.HFC-410A shows heat transfer coefficients similar to HFC-134a and 10% higher than HFC-236fa together with frictional pressure drops 40-50% lower than HFC-134a and 50–60% lower than HFC-236fa.  相似文献   

16.
An experimental test rig for study on the pooling-boiling heat transfer performance of pure and mixed refrigerants were designed and established. With this test system, the heat transfer coefficients (HTCs) of the nucleate boiling on a smooth flat surface were measured for pure fluids of HFC134a, HC290, HC600a and their binary and ternary mixtures. Extensive experimental measures were made for those pure and mixed refrigerants at different heat fluxes from 10 kW m−2 to 300 kW m−2 and different pressures from 0.2 to 0.6 MPa. Comprehensive measured data are presented in this paper. From experimental results, these binary mixtures and ternary mixtures show different heat transfer features according to their vapor–liquid phase equilibria behaviors. New heat transfer correlations were regressed from the measured data with average deviations within ±15% for pure refrigerants and within ±20% for mixtures.  相似文献   

17.
Flow boiling heat transfer experiments using R134a were carried out for jet impingement on smooth and enhanced surfaces. The enhanced surfaces were circular micro pin fins, hydrofoil micro pin fins, and square micro pin fins. The effects of saturation pressure, heat flux, Reynolds number, pin fin geometry, pin fin array configuration, and surface aging on flow boiling heat transfer characteristics were investigated. Flow boiling experiments were carried out for two different saturation pressures, 820 kPa and 1090 kPa. Four jet exit velocities ranging from 1.1–4.05 m/s were investigated. Flow boiling jet impingement on smooth surfaces was characterized by large temperature overshoots, exhibiting boiling hysteresis. Flow boiling jet impingement on micro pin fins displayed large heat transfer coefficients. Heat transfer coefficients as high as 150,000 W/m2 K were observed at a relatively low velocity of 2.2 m/s with the large (D = 125 μm) circular micro pin fins. Jet velocity, surface aging, and saturation pressure were found to have significant effects on the two-phase heat transfer characteristics. Subcooled nucleate boiling was found to be the dominant heat transfer mechanism.  相似文献   

18.
《Applied Thermal Engineering》2007,27(14-15):2353-2362
Flow boiling experiments for refrigerant R-600a, and R-290 mixed with the lubricating oil (EMKARATE) in the serpentine small-diameter (2.46 mm) U-tubes are reported. The tests were conducted at the nominal inlet pressure of 186.2 kPa, vapor qualities (0–0.76), mass flux of 100–320 (kg/m2s) and inlet oil concentrations from 0 to 5 mass% oil. It was noted that a significant degradation of heat transfer coefficients presented at high qualities and high oil concentrations. The present study investigated that whether the average heat transfer coefficient increased as mass fluxes and the numbers of the U-bend increased. In addition, the ratios (the enhanced heat transfer factor, EF) for both R-600a and R-290 refrigerants increased as the oil concentration increased up to 1% then decreased with each increase in oil concentration. Moreover, pressure dropped during evaporation increased with the addition of a lubricant, mass fluxes and the numbers of the U-turn. The values of the pressure drop penalty factor PF were generally larger than l and increased rapidly as the oil concentration and the numbers of the U-turn increased.  相似文献   

19.
Corresponding to the updated flow pattern map presented in Part I of this study, an updated general flow pattern based flow boiling heat transfer model was developed for CO2 using the Cheng–Ribatski–Wojtan–Thome [L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside horizontal tubes, Int. J. Heat Mass Transfer 49 (2006) 4082–4094; L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, Erratum to: “New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside tubes” [Heat Mass Transfer 49 (21–22) (2006) 4082–4094], Int. J. Heat Mass Transfer 50 (2007) 391] flow boiling heat transfer model as the starting basis. The flow boiling heat transfer correlation in the dryout region was updated. In addition, a new mist flow heat transfer correlation for CO2 was developed based on the CO2 data and a heat transfer method for bubbly flow was proposed for completeness sake. The updated general flow boiling heat transfer model for CO2 covers all flow regimes and is applicable to a wider range of conditions for horizontal tubes: tube diameters from 0.6 to 10 mm, mass velocities from 50 to 1500 kg/m2 s, heat fluxes from 1.8 to 46 kW/m2 and saturation temperatures from ?28 to 25 °C (reduced pressures from 0.21 to 0.87). The updated general flow boiling heat transfer model was compared to a new experimental database which contains 1124 data points (790 more than that in the previous model [Cheng et al., 2006, 2007]) in this study. Good agreement between the predicted and experimental data was found in general with 71.4% of the entire database and 83.2% of the database without the dryout and mist flow data predicted within ±30%. However, the predictions for the dryout and mist flow regions were less satisfactory due to the limited number of data points, the higher inaccuracy in such data, scatter in some data sets ranging up to 40%, significant discrepancies from one experimental study to another and the difficulties associated with predicting the inception and completion of dryout around the perimeter of the horizontal tubes.  相似文献   

20.
Experiments of flow boiling heat transfer were conducted in four horizontal flattened smooth copper tubes of two different heights of 2 and 3 mm. The equivalent diameters of the flattened tubes are 8.6, 7.17, 6.25, and 5.3 mm. The working fluids were R22 and R410A. The test conditions were: mass velocities from 150 to 500 kg/m2 s, heat fluxes from 6 to 40 kW/m2 and saturation temperature of 5 °C. The experimental heat transfer results are presented and the effects of mass flux, heat flux, and tube diameter on heat transfer are analyzed. Furthermore, the flow pattern based flow boiling heat transfer model of Wojtan et al. [L. Wojtan, T. Ursenbacher, J.R. Thome, Investigation of flow boiling in horizontal tubes: Part I – A new diabatic two-phase flow pattern map, Int. J. Heat Mass Transfer 48 (2005) 2955–2969; L. Wojtan, T. Ursenbacker, J.R. Thome, Investigation of flow boiling in horizontal tubes: Part II – Development of a new heat transfer model for stratified-wavy, dryout and mist flow regimes, Int. J. Heat Mass Transfer 48 (2005) 2970–2985], using the equivalent diameters, were compared to the experimental data. The model predicts 71% of the entire database of R22 and R410A ±30% overall. The model predicts well the flattened tube heat transfer coefficients for R22 while it does not predicts well those for R410A. Based on several physical considerations, a modified flow boiling heat transfer model was proposed for the flattened tubes on the basis of the Wojtan et al. model and it predicts the flattened tube heat transfer database of R22 and R410A by 85.8% within ±30%. The modified model is applied to the reduced pressures up to 0.19.  相似文献   

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

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