首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Physical modeling of gas/liquid mass transfer in a gas stirred ladle   总被引:3,自引:0,他引:3  
The absorption of gas through the plume eye and of an injected gas in a steelmaking ladle process was investigated, using a physical model of CO2 absorption into a NaOH solution. The results show that the inert gas escaping through the plume eye is ineffective in protecting the bath from the atmosphere, and placing an oil layer (simulated slag) decreases the absorption rate significantly. Increasing the flow rate of the inert gas not only exposes more of the liquid surface to the CO2 atmosphere, but also increases the mass transfer coefficient at the surface. The overall mass transfer between an injected CO2 gas and NaOH solution includes the mass transfer through the surface of the bath as well as the mass transfer in the bubble dispersion zone. The difference between the mass transfer in the bubble dispersion zone and the overall mass transfer was found to be significant for relatively low gas flow rates. The mass transfer coefficient of CO2 in the bubble dispersion zone was estimated using available information regarding the bubble size and velocity. Mass transfer coefficient estimated for the constant bubble frequency regime shows a dependence on gas flow rate. However, if a constant characteristic size of bubbles is assumed as an alternative approach, the mass transfer coefficient is independent of the gas flow rate.  相似文献   

2.
Various forms of plumes have been identified following the injection of air at different rates through a porous plug into water contained in a ladle-shaped vessel. Discrete bubbles form at the plug and rise uniformly through the column of liquid at gas flow rates up to 14 cm3/s cm2 of plug surface; at higher flow rates, groups of bubbles increasingly coalesce into larger gas pockets, and beyond about 40 cm3/s cm2, the gas globes are large enough to cover the entire plug surface before detachment and gradual disintegration as they rise through the body of liquid. The gas fraction, as well as bubble frequency, bubble velocity, and bubble size, have been measured in the various dispersion regimes by means of an electroresistivity probe. The radial distributions of gas fraction and bubble frequency are approximately bell-shaped about the axis of flow, and the reduced values are close to Gaussian functions of the reduced radial distance from the axis. The gas fraction along the axis has been correlated to the reduced height of the plume; it increases with decreasing distance above the plug and with increasing gas flow rate. The axial bubble frequency shows a decrease in the vicinity of the plug with the onset of bubble coalescence, but the values of the frequencies at all gas injection rates converge to about 12 s−1 toward the surface of the bath. The mean bubble velocity increases with increasing flow rate but drops once coalescence is fully established. Conversely, there is a sudden increase in the mean bubble diameter with the onset of coalescence. The axial and radial components of the velocity of the liquid surrounding the plume have been measured by means of a Laser-Doppler Velocimeter (LDV), and the results show that the circulation patterns are identical, irrespective of the dispersion regime. The axial flow which is upward in the vicinity of the plume decreases in magnitude with increasing radial distance, ultimately reversing to an in-creasing downward flow beyond a certain distance from the plug axis. Similarly, the radial flow which is outward from the plume near the liquid surface decreases steadily with depth and eventually reverses to an inward flow at a depth independent of the gas injection rate. The profiles of the axial velocities are almost sigmoidal, except in the coalescence regime, where the effect of turbulence is profound at the upper liquid layers. The radial liquid velocities are generally small relative to the axial components, only about one-fifth as large, considering the maximum average values.  相似文献   

3.
Measurement of the velocities of bubbles and liquid with a two-element electroresistivity probe and laser-Doppler velocimeter, respectively, during bottom injection of air into a water bath, has confirmed the existence of a critical gas-injection rate. Above the critical flow rate, the change of axial bubble velocity in the air jet, and of liquid velocity with increasing volume flow rate, diminishes markedly. The existence of the critical flow rate is explicable from high-speed motion pictures of the vertical gas jets, which reveal four zones of gas dispersion axially distributed above the orifice: primary bubble at the orifice, free bubble, plume consisting of disintegrated bubbles, and spout at the bath surface. With increasing gas-injection rate, the free-bubble zone expands such that the point of bubble disintegration rises closer to the bath surface. Above the critical flow rate, the free bubbles rise with minimal breakup and erupt from the bath surface with maximum energy discharge. The combined Kelvin-Helmholtz, Rayleigh-Taylor instability theory has been applied to analyze the bubble breakup in the bath and the critical gas-injection rate in a gas-stirred ladle. The criterion for the critical diameter of bubble breakup has been found to depend primarily on the surface tension and density of the liquid. In the analysis, the propagation time of a disturbance on a bubble surface at the “most unstable” wave number has been compared with the bubble rising time in the bath in order to determine the critical gas-flow rate. The predicted critical values are in close agreement with the measured results. M. ZHOU formerly was Post Doctoral Fellow with the Centre for Metallurgical Process Engineering, University of British Columbia, Vancouver, BC, Canada V6T 1Z4 J.K. BRIMACOMBE holds the Alcan Chair in Materials Process Engineering  相似文献   

4.
Experimental measurements of the flow velocity were carried out with liquid Wood's metal in a ladle-shaped vessel with an inner diameter of 40 cm at centric gas blowing. By means of permanent magnet probes the liquid flow field was measured under various blowing conditions. The results show that a circulating flow field is established in the vessel. In the bubble plume zone an upwardly directed liquid flow is formed. The radial distribution of the flow velocity in this zone follows a Gaussian function. The axial flow velocity increases with growing gas flow rate and is nearly constant in vertical direction. The width of the upward flow becomes larger with increasing distance from the nozzle and its dependence from the gas flow rate is not considerable. At centric gas blowing the liquid in the upper part of the bath streams quickly, whereas in the lower part so-called “dead zones” with very low flow velocity are present. Besides the time-averaged value of the flow velocity, the turbulent behaviours of the liquid flow such as fluctuation velocity, the turbulent kinetic energy and its dissipation rate were investigated on the basis of measured data. It was found that the liquid flow is turbulent particularly in the region of bubble plume and of bath surface. The radial profiles of these parameters can also be described by a Gaussian function. Only a small part of the gas stirring energy is changed into the kinetic energy of the directed liquid flow. Most of the stirring energy is already dissipated in the bubble plume zone.  相似文献   

5.
It is shown that with the aid of digital simulation methods complex multiphase interrelated systems, such as gas-injection process can be analysed. Interdependencies can be revealed and quantitative evaluation of characteristic system quantities are provided. The method of digital system simulation is a very convenient tool for process analysis or system engineering. Results of the computer-aided process simulations (Caps) yield a better understanding of complex phenomena and better aimed engineering of gas dispersion techniques in metallurgical processes. A particular interest of this investigation is to reveal the effect of mass-transfer rate on the hydrodynamic behaviour of a gas-injection process. The combined effects of total flow rate of injected gas and mass-transfer rate on the system quantities such as mixing power, induced liquid flow rate, holdup, interfacial area and volumetric mass transfer coefficient are evaluated under steady state conditions of the investigated systems and illustrated in simulation plots. The liquid velocity has a minor effect on bubble size at some distance from the orifice but controls the location of bubble breakup. The frequency of bubble breakup and final bubble size depends on the intensity of mass transfer. Mixing power due to gas bubbles and circulation velocity of the steel bath increase appreaciably if there is a chance of bath reactions producing more gas. The integral mean values of mixing power, induced velocity of liquid and holdup in plume, specific interfacial area and volumetric mass-transfer coefficient increase with increasing total flow rate of injected gas and intensity of mass transfer.  相似文献   

6.
In a model converter containing oil and water the volumetric mass transfer of caprylic acid was determined in dependence of the gas flow rate. Several other parameters were varied such as the number and configuration of nozzles, gas velocity, oil viscosity and volume and model size. Partially the results could be described by dimensionless equations. Visual observations of the bath surface with special regard to the wave formation complete the results of the measurements. By knowledge of the drop size spectrum independently determined, the mass transfer coefficient under dispersion conditions could be evaluated. These values were compared with mass transfer coefficients which had been directly measured by the single rising drop technique as well as the bubble stirred interface without dispersion. It seems that turbulence in dispersed systems does not enlarge the mass transfer coefficient.  相似文献   

7.
The structural development of air-water bubble plumes during upward injection into a ladle-shaped vessel has been measured under different conditions of air flow rate, orifice diameter, and bath depth. The measured radial profiles of gas fraction at different axial positions in the plume were found to exhibit good similarity, and the distribution of the phases in the plume was correlated to the modified Froude number. Different regions of flow behavior in the plume were identified by changes in bubble frequency, bubble velocity, and bubble pierced length which occur as bubbles rise in the plume. Measurement of bubble velocity indicates that close to the nozzle the motion of the gas phase is strongly affected by the injection velocity; at injection velocities below 41 m/s, the velocity of the bubbles along the centerline exhibits an increase with height, while above, the tendency reverses. High-speed film observations suggest that this effect is related to the nature of gas discharge,i.e., whether the gas discharge produces single bubbles or short jets. In this region of developing flow, measurement of bubble frequency and pierced length indicates that break-up of the discharging bubbles occurs until a nearly constant bubble-size distribution is established in a region of fully developed flow. In this largest zone of the plume the bubbles influence the flow only through buoyancy, and the spectra of bubble pierced length and diameter can be fitted to a log-normal distribution. Close to the bath surface, a third zone of bubble motion behavior is characterized by a faster decrease in bubble velocity as liquid flows radially outward from the plume.  相似文献   

8.
Based on the Euler–Euler approach, a mathematical model is established to describe gas and liquid two-phase flow in the gas-stirred system for steelmaking, and the influences of the interphase force including turbulent dispersion force, drag force, and lift force are investigated. The modified k–ε model with extra source terms to account for the bubble-induced turbulence is adopted to model the turbulence in the system, and the simulation results of gas volume fraction, liquid velocity, and turbulent kinetic energy are compared with the measured data. The results show that the turbulent dispersion force dominates the bubbly plume shape and is responsible for successful prediction of the gas volume fraction. The bubble-induced turbulence has a significant influence on the liquid turbulence, and the conversion coefficient C b, which denotes the fraction of bubble-induced energy converted into liquid turbulence, should be in the range of 0.8 and 0.9. The drag force also strongly influences the bubbly plume dynamics, and the coefficient model proposed by Kolev performs the best for determining the drag force; however, the lift force and bubble diameter do not have much effect on the current bubbly plume system. For different gas flow rates, the current Euler–Euler approach predictions are more consistent with the measured data than the Euler–Lagrange approach and the early Euler–Euler model.  相似文献   

9.
A computation fluid dynamics–population balance model (CFD–PBM) coupled model has been proposed to investigate the bubbly plume flow and inclusion behavior including growth, size distribution, and removal in gas-stirred ladles, and some new and important phenomena and mechanisms were presented. For the bubbly plume flow, a modified k-ε model with extra source terms to account for the bubble-induced turbulence was adopted to model the turbulence, and the bubble turbulent dispersion force was taken into account to predict gas volume fraction distribution in the turbulent gas-stirred system. For inclusion behavior, the phenomena of inclusions turbulent random motion, bubbles wake, and slag eye forming on the molten steel surface were considered. In addition, the multiple mechanisms both that promote inclusion growth due to inclusion–inclusion collision caused by turbulent random motion, shear rate in turbulent eddy, and difference inclusion Stokes velocities, and the mechanisms that promote inclusion removal due to bubble-inclusion turbulence random collision, bubble-inclusion turbulent shear collision, bubble-inclusion buoyancy collision, inclusion own floatation near slag–metal interface, bubble wake capture, and wall adhesion were investigated. The importance of different mechanisms and total inclusion removal ratio under different conditions, and the distribution of inclusion number densities in ladle, were discussed and clarified. The results show that at a low gas flow rate, the inclusion growth is mainly attributed to both turbulent shear collision and Stokes collision, which is notably affected by the Stokes collision efficiency, and the inclusion removal is mainly attributed to the bubble-inclusion buoyancy collision and inclusion own floatation near slag–metal interface. At a higher gas flow rate, the inclusions appear as turbulence random motion in bubbly plume zone, and both the inclusion–inclusion and inclusion-bubble turbulent random collisions become important for inclusion growth and removal. With the increase of the gas flow rate, the total removal ratio increases, but when the gas flow rate exceeds 200 NL/min in 150-ton ladle, the total removal ration almost does not change. For the larger size inclusions, the number density in bubbly plume zone is less than that in the sidewall recirculation zones, but for the small size inclusions, the distribution of number density shows the opposite trend.  相似文献   

10.
The mass transfer characteristics in a steel bath during the AOD refining process with the conditions of combined side and top blowing were investigated. The experiments were conducted on a water model unit of 1/4 linear scale for a 120‐t combined side and top blowing AOD converter. Sodium chloride powder of analytical purity was employed as the flux for blowing, and the mass transfer coefficient of solute (NaCI) in the bath was determined under the conditions of the AOD process. The effects of the gas flow rates of side and top blowing processes, the position arrangement and number of side tuyeres, the powdered flux particle (bubble) size and others on the characteristics were examined. The results indicated that, under the conditions of the present work, the mass transfer coefficient of solute in the bath liquid is in the range of (7.31×10?5‐3.84×10?4) m/s. The coefficient increases non‐linearly with increasing angle between each tuyere, for the simple side blowing process at a given side tuyere number and gas side blowing rate. The gas flow rate of the main tuyere has a governing influence on the characteristics, and the gas jet from the top lance decreases the mass transfer rate, the relevant coefficient being smaller than that for a simple side blowing. Also, in the range of particle (bubble) size used in the present work and with all other factors being constant, raising particle (bubble) size increases the coefficient. Excessively fine powder particle (bubble) sizes are not advantageous to strengthening the mass transfer. With the oxygen top blowing rate practiced in the industrial technology, the side tuyere arrangements of 7 and 6 tuyeres with an angular separation of 22.5° and 27° between each tuyere, as well as 5 tuyeres with an angle of 22.5° between each tuyere can provide a larger mass transfer rate in the bath. Considering the relative velocity of the particles to the liquid, the energy dissipation caused by the fluctuation in the velocity of the liquid in turbulent flow and regarding the mass transfer as that between a rigid bubble and molten steel, the related dimensionless relationships for the coefficient were obtained.  相似文献   

11.
Experimental investigation of mixing phenomena in a gas stirred liquid bath   总被引:1,自引:0,他引:1  
Mixing phenomena in a room temperature water bath, agitated by injecting air through a straight circular nozzle fitted axially at the bottom of the vessel, were characterized by experimentally measuring mixing time(t mix) by electrical conductivity technique. It was found thatt mix defined at 99.5 pct homogenization did not depend on location and size of conductivity probe, location of tracer injection, and the amount of tracer injected. tpet decreased with increasing gas flow rate and bath height, but decreasing nozzle diameter. Visual observations of the two-phase plume and flow conditions in the bath revealed that the plume swirled above a certain gas flow rate which enhanced the mixing rates in the bath. The transitions in Int mix vs In εb curves were found to correspond to onset of swirling; εb is the rate of buoyancy energy input per unit bath volume. Systematic analysis of experimental data revealed that a fraction of gas kinetic energy contributed to mixing in the bath. It was a function of bath height, being negligible at lower bath heights and almost 1 at larger bath heights. Further, it was experimentally found thatt mix decreased with increasing bath height only up to a certain value, beyond which it started increasing. Visual observations of the bath revealed that the height at whicht mix started increasing corresponded to a transition in which the bath was converted into a bubble column. The experimental data, for a particular bath height, were fitted into two separate straight lines of the formt mix = −n wherec andn are empirical constants and ε is the rate of energy input per unit bath volume. Formerly Graduate Student in the Department of Metallurgical Engineering at the Indian Institute of Technology, Kanpur, India  相似文献   

12.
Profiles of the rate of dissipation of turbulent kinetic energy were inferred from temperature microstructure measurements near a bubble plume at the center of a tank with diameter of 13.7 m and maximum depth of 8.3 m. Six sets of between 18 and 51 profiles were collected at airflow rates of 0.1–0.6 L/s, measured at atmospheric pressure, and ensemble-averaged dissipation profiles were calculated. The dissipation in all cases was between 10?8 and 10?6?m2/s3 in most of the profile, but it increased sharply near the water surface. Energy considerations are used to discuss the experimental results in terms of previous numerical models of bubble plume turbulence. Two previous numerical studies show that the turbulence dissipates between 15 and 30% of the available power. In the experiments, the fraction is less than 1% because some of the energy of the plume is used to generate waves on the water surface and the profiles used to compute the volume-averaged dissipation were relatively far from the bubble plume.  相似文献   

13.
A model study was carried out to elucidate bubble and liquid flow characteristics in the reactor of metals refining processes stirred by gas injection. Wood’s metal with a melting temperature of 70 °C was used as the model of molten metal. Helium gas was injected into the bath through a centered single-hole bottom nozzle to form a vertical bubbling jet along the centerline of the bath. The bubble characteristics specified by gas holdup, bubble frequency, and so on were measured using a two-needle electroresistivity probe, and the liquid flow characteristics, such as the axial and radial mean velocity components, were measured with a magnet probe. In the axial region far from the nozzle exit, where the disintegration of rising bubbles takes place and the radial distribution of gas holdup follows a Gaussian distribution, the axial mean velocity and turbulence components of liquid flow in the vertical direction are predicted approximately by empirical correlations derived originally for a water-air system, although the physical properties of the two systems are significantly different from each other. Under these same conditions, those turbulent parameters in high-temperature metals refining processes should thus be accurately predicted by the same empirical correlations.  相似文献   

14.
A macroscopic, steady state energy balance model has been formulated to describe mixing phenom-ena in a liquid bath stirred by injecting gas through a straight nozzle fitted axially at the bottom of the vessel. This, along with experimental data on a water model previously reported, was employed to make predictions. Input energy terms considered in the model consist of buoyancy energy and empirically determined fraction of gas kinetic energy. Dissipation of energy was attributed to liquid circulation and bubble slip. The two-phase plume was assumed to be a truncated cone whose dimen-sions depended upon operating conditions. Numerical solution of model equations gave liquid velocity and gas hold-up inside the plume as well as liquid circulation rate and liquid velocity in the region outside the plume. Influence of process variables, e.g., gas flow rate, bath height, and nozzle diameter, have been predicted. Validity of the model has been established by comparing some pre-dicted entrainment ratios with those experimentally measured by other investigators. Empirical cor-relations to predict circulation time and circulation number have been proposed. Circulation number was found to vary between 2 and 12 in contrast to the existing assumption in the literature of a con-stant value of 3. Usefulness of these correlations in predicting mixing time for industrial vessels has been demonstrated. Formerly a Graduate Student in the De-partment of Metallurgical Engineering at the Indian Institute of Technol-ogy, Kanpur  相似文献   

15.
Basic fluid dynamic features of combined electromagnetic stirring, EMS, and gas stirring (EMGAS) have been studied in the present work. A transient and turbulent multiphase numerical flow model was built. Simulations of a real size ladle furnace were conducted for 7 cases, operating with and without combined stirring and varying the argon gas inlet plug position. The results of these simulations are compared considering melt velocity, melt turbulence, melt/slag‐interface turbulence and dispersion of gas bubbles. An experimental water model was also built to simulate the effects of combined stirring. The water model was numerically simulated and visual comparison of the gas plume shape and flow pattern in the numerical and in the experimental model was also done for 3 flow situations. The results show that EMGAS has a strong flexibility regarding the flow velocity, gas plume, stirring energy, mixing time, slag layer, etc.  相似文献   

16.
The water model experiments were carried out to study the bubble morphology in the tundish and mold with the process of annular argon blowing at tundish upper nozzle. The effects of the position of gas permeable brick, the casting speed and the argon flow rate on the bubble size distribution, the bubble migration behavior and the flow behavior of liquid steel near the liquid level in tundish were further investigated, coupled with the numerical simulation. The results show that with the process of annular argon blowing at tundish upper nozzle, a frustum cone shaped bubble plume can be formed around the stopper rod. The concentration of argon bubbles gradually decreases outward along the radial direction of the stopper rod. Owing to the wall attached effect, the bubble plumes float upward along the stopper rod, which can increase the collision probability between bubbles and the velocity of bubble plumes, causing a larger impact strength on the liquid level in tundish. In addition, a part of small bubbles are wrapped into the nozzle and the mold due to the drag force of liquid steel. With increasing argon flow rate, the number of bubbles in annular bubble plumes and the vertical velocity of liquid steel near the liquid level in tundish increase significantly. With increasing casting speed, the width and the bubble number of annular bubble plumes gradually decrease, leading to a decrease of the vertical velocity of liquid steel near the liquid level in tundish. Increasing the distance between the annular gas permeable brick and the center of tundish upper nozzle, the dispersion of bubbles and the width of bubble plumes increase, and the impact strength of bubbles acting on the liquid level in tundish becomes weaker. As the argon flow rate and the casting speed increase, and the distance between the gas permeable brick and the center of tundish upper nozzle decreases, the gas volume and bubble size in the mold increase. Under the experimental conditions, when the inner and outer diameters of the annular gas permeable brick are 110mm and 140mm, respectively, and the casting speed is 1.2m/min, the appropriate argon flow rate is 4L/min.  相似文献   

17.
摘要:通过水模型实验研究了上水口环形吹氩工艺下中间包和结晶器内气泡形貌,并结合数值模拟分析了透气砖位置、拉坯速度和吹氩量对中间包和结晶器内气泡尺寸、气泡迁移和中间包近液面钢液流动的影响。结果表明:上水口环形吹氩形成以塞棒为中心的圆台状气泡羽流,气泡浓度沿径向向外逐渐减少;附壁效应使得气泡羽流偏向塞棒壁面流动,增大气泡的碰撞聚并概率和近塞棒壁面的羽流上升速度,对中间包液面产生较大冲击作用;同时,部分细小气泡会随钢液进入水口及结晶器内部;增大吹氩量,中间包内环形气泡羽流中气泡数目明显增多,中间包近液面钢液上升速度增大;增大拉坯速度,环形气泡羽流的宽度和气泡数量逐渐减小,近液面速度减小;增大透气环距水口中心距离,中间包内气泡弥散度增大,环形气泡羽流宽度也随之增大,气泡羽流对中间包液面冲击作用减弱;增大吹氩量和拉坯速度、减小透气环距水口中心距离,进入结晶器的气量和气泡尺寸逐渐增大。实验条件下,透气环内外径为110mm/140mm、拉坯速度为1.2m/min时,吹氩量为4L/min较为合适。  相似文献   

18.
Averaged and turbulent fluctuating liquid velocities in the gas/liquid plume zone of a gas-stirred water model ladle were measured with a combined laser Doppler anemometer (LDA) and elec-trical probe technique. The measured turbulence fields, void fraction distribution, and gas and liquid velocities in the plume zone were used for evaluation of various turbulence models. It was found that, among all of the turbulence models tested, only a modified k-ε model, with extra source terms to take into account the generation and dissipation resulting from the inter-action of the bubbles with the liquid, yielded good agreement with both the mean liquid flow field and the turbulent kinetic energy distribution. However, the values of the coefficients orig-inally proposed by their authors were found inapplicable to the bubbly plume situation; more appropriate values of the coefficients were determined based on comparison with experimental measurement.  相似文献   

19.
A mathematic model for an argon bubble containing magnesium vapour in hot metal was developed for interpreting the desulfurization process by injection of granulated magnesium. The change of rising velocity, pressure, size and magnesium vapour content of the bubble and the rates of desulfurization, deoxidation, magnesium dissolution by the bubble were investigated with the aid of the model. Most of the magnesium vapour in an argon bubble is consumed by the reaction with hot metal while the bubble is rising from 3.5 m depth in hot metal to the bath surface. The magnesium vapour content in the bubble decreases to a minimum halfway on the bubble's rising and slightly increases thereafter as the bubble approaches the bath surface. During the initial or middle stages of the injection process, magnesium vapour in a gas bubble is mainly consumed by the desulfurization and deoxidation of hot metal. During the final stage or within the active zone of the process, magnesium dissolution into hot metal and unreacted magnesium vapour in the bubble arriving at the bath surface become significant.  相似文献   

20.
Water model experiments were performed to study the effect of top slag on the mean flow and turbulence characteristics in a steel bath agitated by bottom gas injection. The slag was modeled by silicone oil with a density of 0.968 g/cm3 and a kinematic viscosity 100 times larger than that of water at 25 °C. Velocity measurements were made using a two-dimensional (2-D) laser Doppler velocimeter (LDV) in the absence of swirl motions. The output signals of the LDV system were processed on a personal computer to obtain the axial and radial mean velocity components, the root-mean-square (rms) values of the axial and radial turbulence fluctuations, the Reynolds shear stress, and the turbulence production for two cases with and without top slag. The bubbling jet (or the bubble dispersion) region was localized near the centerline of the bath by the presence of the top oil layer. The mean flow and turbulence motions in the recirculation region located outside the bubbling jet region were also suppressed significantly by the top layer. This result could be attributed to the entrainment of top slag into steel in a real system.  相似文献   

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

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