首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
The accretion of hot slowly rotating gas onto a supermassive black hole is considered. The important case where the velocities of turbulent pulsations at the Bondi radius r B are low, compared to the speed of sound c s, is studied. Turbulence is probably responsible for the appearance of random average rotation. Although the angular momentum at r B is low, it gives rise to the centrifugal barrier at a depth r c = l 2 /GM BH ? r B, that hinders supersonic accretion. The numerical solution of the problem of hot gas accretion with finite angular momentum is found taking into account electron thermal conductivity and bremsstrahlung energy losses of two temperature plasma for density and temperature near Bondi radius similar to observed in M87 galaxy. The saturation of the Spitzer thermal conductivity was also taken into account. The parameters of the saturated electron thermal conductivity were chosen similar to the parameters used in the numerical simulations of interaction of the strong laser beam radiation with plasma targets. These parameters are confirmed in the experiments. It is shown that joint action of electron thermal conductivity and free-free radiation leads to the effective cooling of accreting plasma and formation of the subsonic settling of accreting gas above the zone of a centrifugal barrier. A toroidal condensation and a hollow funnel that separates the torus from the black hole emerge near the barrier. The barrier divides the flow into two regions: (1) the settling zone with slow subKeplerian rotation and (2) the zone with rapid supersonic nearly Keplerian rotation. Existence of the centrifugal barrier leads to significant decrease of the accretion rate ? in comparison with the critical Bondi solution for γ = 5/3 for the same values of density and temperature of the hot gas near Bondi radius. Shear instabilities in the torus and related friction cause the gas to spread slowly along spirals in the equatorial plane in two directions.As a result, outer (r > r c) and inner (r < r c) disks are formed. The gas enters the immediate neighborhood of the black hole or the zone of the internal ADAF flow along the accretion disk (r < r c). Since the angular momentum is conserved, the outer disk removes outward an excess of angular momentum along with part of the matter falling into the torus. It is possible, that such outer Keplerian disk was observed by Hubble Space Telescope around the nucleus of the M87 galaxy in the optical emission lines. We discuss shortly the characteristic times during which the accretion of the gas with developed turbulence should lead to the changes in the orientation of the torus, accretion disk and, possibly, of the jet.  相似文献   

2.
We analyze the spectra of DR Tau in the wavelength range 1200 to 3100 Å obtained with the GHRS and STIS spectrographs from the Hubble Space Telescope. The profiles for the C IV 1550 and He II 1640 emission lines and for the absorption features of some lines indicate that matter falls to the star at a velocity ~300 km s?1. At the same time, absorption features were detected in the blue wings of the N I, Mg I, Fe II, Mg II, C II, and Si II lines, suggesting mass outflow at a velocity up to 400 km s?1. The C II, Si II, and Al II intercombination lines exhibit symmetric profiles whose peaks have the same radial velocity as the star. This is also true for the emission features of the Fe II and H2 lines. We believe that stellar activity is attributable to disk accretion of circumstellar matter, with matter reaching the star mainly through the disk and the boundary layer. At the time of observations, the accretion luminosity was Lac ? 2L at an accretion rate ?10?7M yr?1. Concurrently, a small (<10%) fraction of matter falls to the star along magnetospheric magnetic field lines from a height ~R*. Within a region of size ?3.5R*, the disk atmosphere has a thickness ~0.1R* and a temperature ?1.5 × 104 K. We assume that disk rotation in this region significantly differs from Keplerian rotation. The molecular hydrogen lines are formed in the disk at a distance <1.4 AU from the star. Accretion is accompanied by mass outflow from the accretion-disk surface. In a region of size <10R*, the wind gas has a temperature ~7000 K, but at the same time, almost all iron is singly ionized by H I L α photons from inner disk regions. Where the warm-wind velocity reaches ?400 km s?1, the gas moves at an angle of no less than 30° to the disk plane. We found no evidence of regions with a temperature above 104 K in the wind and leave open the question of whether there is outflow in the H2 line formation region. According to our estimate, the star has the following set of parameters: M* ? 0.9M, R* ? 1.8R, L* ? 0.9L, and \(A_V \simeq 0\mathop .\limits^m 9\). The inclination i of the disk axis to the line of sight cannot be very small; however, i≤60°.  相似文献   

3.
The following conclusions about the kinematics and parameters of the gas in the vicinity of TW Hya have been drawn from an analysis of optical and ultraviolet line profiles and intensities. The accreting matter rises in the magnetosphere to a distance z>R* above the disk plane and falls to the star near its equator almost perpendicular to its plane. The matter outflows from a disk region with an outer radius of ≤0.5 AU. The [OI], [SII], and H2 lines originate in the disk atmosphere outside the outflow region, where the turbulent gas velocity is close to the local speed of sound. In the formation region of the forbidden lines, T?8500 K and Ne?5×106 cm?3, and the hydrogen is almost neutral: xe<0.03. The absorption features observed in the blue wings of some of the ultraviolet lines originate in the part of the wind that moves almost perpendicular to the disk plane, i.e., in the jet of TW Hya. The V z gas velocity component in the jet decreases with increasing distance from the jet axis from 200 to 30 km s?1. The matter outflowing from the inner disk boundary, moves perpendicular to the disk plane in the formation region of blue absorption line components, at a distance of ~0.5 AU from the axis of symmetry of the disk. This region of the wind is collimated into the jet at a distance of <3 AU from the disk plane. The gas temperature in the formation region of absorption components is ?2×104 K, and the gas density is <3×106 cm?3. This region of the jet is on the order of several AU away from the disk plane, while free recombination in the jet begins even farther from the disk. The mass-loss rate for TW Hya is \(\dot M_w < 7 \times 10^{ - 10} M_ \odot yr^{ - 1}\), which is a factor of 3lower than the mean accretion rate. The relative abundance of silicon and aluminum in the jet gas is at least an order of magnitude lower than its standard value.  相似文献   

4.
We analyze the superfine structure of the supermaser H2O emission region in Orion KL over the period 1979–1999. The angular resolution reached 0.1 mas, which corresponds to 0.045 AU at a distance to Orion KL of 450 pc. We determined the velocity of the local standard of rest, VLSR = 7.65 km s?1. The formation of a protostar is accompanied by a structure that consists of an accretion disk, a bipolar outflow, and a surrounding envelope. The disk is at the stage of separation into protoplanetary rings. The disk plane is warped like the brim of a hat. The disk is 27 AU in diameter and ~0.3 AU in thickness. The rings contain ice granules. Radiation and stellar wind sublimate and blow away the water molecules to form halos around the rings, maser rings. The radiation from the rings is concentrated in the azimuthal plane, and its directivity reaches 10?3. The relative velocities of the rings located in the central part of the disk 15 AU in diameter correspond to rigid-body rotation, Vrot = ΩR. The rotation period is T ≈ 170 yr. The injector is surrounded by a toroidal structure 1.2 AU in diameter. The diameter of the injected flow does not exceed 0.05 AU. A highly collimated bipolar outflow with a diameter of ~0.1 AU is observed at a distance as large as 3 AU. Precession of the injector axis with a period of ~10 yr forms a spiral flow structure. The flow velocity is ~10 km s?1. The kinetic energy of the accreting matter and the disk is assumed to be transferred to the bipolar outflow, causing the rotation velocity distribution of the rings to deviate from the Keplerian velocity. The surrounding envelope amplifies the emission from the structure at a velocity of 7.65 km s?1 in a band of ~0.5 km s?1 by more than two orders of magnitude, which determines the supermaser emission.  相似文献   

5.
We consider the problem of dust grain survival in the disk winds from T Tauri and Herbig Ae stars. For our analysis, we have chosen a disk wind model in which the gas component of the wind is heated through ambipolar diffusion to a temperature of ~104 K. We show that the heating of dust grains through their collisions with gas atoms is inefficient compared to their heating by stellar radiation and, hence, the grains survive even in the hot wind component. As a result, the disk wind can be opaque to the ultraviolet and optical stellar radiation and is capable of absorbing an appreciable fraction of it. Calculations show that the fraction of the wind-absorbed radiation for T Tauri stars can be from 20 to 40% of the total stellar luminosity at an accretion rate ? a = 10?8-10?6 M yr?1. This means that the disk winds from T Tauri stars can play the same role as the puffed-up inner rim in current accretion disk models. In Herbig Ae stars, the inner layers of the disk wind (r ≤ 0.5 AU) are dust-free, since the dust in this region sublimates under the effect of stellar radiation. Therefore, the fraction of the radiation absorbed by the disk wind in this case is considerably smaller and can be comparable to the effect from the puffed-up inner rim only at an accretion rate of the order of or higher than 10?6 M yr?1. Since the disk wind is structurally inhomogeneous, its optical depth toward the observer can be variable, which should be reflected in the photometric activity of young stars. For the same reason, moving shadows from gas and dust streams with a spiral-like shape can be observed in high-angular-resolution circumstellar disk images.  相似文献   

6.
We have analyzed the Hubble Space Telescope spectrum of the young star FU Ori in the range 2300–3100 Å. The long-wavelength part of the spectrum is similar to the spectrum of a supergiant with T eff ? 5000–6000 K, but the range of wavelengths shorter than ?2600 Å is dominated by radiation from a region with T eff ? 9000 K. We discuss the possibility of explaining these peculiarities of the spectrum, the Al II] 2669.2 emission line profile, and the results of X-ray observations for FU Ori in terms of an accretion disk model whose thickness increases as the star is approached starting from distances ?1012 cm. Near the star, the disk has the shape of a cone in which only the part of its surface on the far (from the observer) side is visible. The suggested model is a kind of a compromise between the models of a thin α-disk and a supergiant: basically, this is an accretion model, but it resembles a supergiant in observational manifestations. Numerous absorption lines originating in the disk wind are superimposed on the disk spectrum. The wind is a cold (T ? 5000 K), dense (N e ? 1011 cm?3) gas. The number of wind absorption lines in the ultraviolet spectrum of FU Ori increases with decreasing wavelength. This causes a rapid decline in intensity in the short-wavelength part of the spectrum. As a result, the maximum temperature in the disk estimated from low-resolution IUE spectra has been underestimated.  相似文献   

7.
During the period 1979–1999, we investigated the hyperfine structure of the H2O supermaser region located in the core of the molecular cloud OMC-1 in Orion KL. The angular resolution is 0.1 mas, which corresponds to 0.045 AU. The detected structure, which consists of a central object, an accretion disk, a bipolar outflow, and an envelope, corresponds to the initial formation stage of a low-mass star. The accretion disk is at the stage of separation into groups of concentric rings. The bipolar outflow is a neutral, highly collimated jet of accreted material that includes H2O molecules and dust grains in the icy envelope. The injector is a bright compact source with a size <0.05 AU and a brightness temperature Tb≈1017 K. The velocity of the bipolar outflow is v≈10 km s?1. The rotation velocity of the jet is vrot≈1.5 km s?1. The jet has the shape of a conical helix due to the precession of the rotation axis. Occasionally, dense blobs (comet-shaped bullets) are ejected. The envelope amplifies the radio emission from the structures in a ~0.5 km s?1 maser window band with velocities v≈7.65 km s?1 by more than two orders of magnitude.  相似文献   

8.
We present the optical spectroscopic study of two classical Be stars, 59 Cyg and OT Gem obtained over a period of few months in 2009. We detected a rare triple-peak H α emission phase in 59 Cyg and a rapid decrease in the emission strength of H α in OT Gem, which are used to understand their circumstellar disks. We find that 59 Cyg is likely to be rapid rotator, rotating at a fractional critical rotation of ~0.80. The radius of the H α emission region for 59 Cyg is estimated to be R d/R ? ~ 10.0, assuming a Keplerian disk, suggesting that it has a large disk. We classify stars which have shown triple-peaks into two groups and find that the triple-peak emission in 59 Cyg is similar to ζ Tau. OT Gem is found to have a fractional critical rotation of ~0.30, suggesting that it is either a slow rotator or viewed in low inclination. In OT Gem, we observed a large reduction in the radius of the H α emission region from ~6.9 to ~1.7 in a period of three months, along with the reduction in the emission strength. Our observations suggest that the disk is lost from outside to inside during this disk loss phase in OT Gem.  相似文献   

9.
We present the results of our infrared JHK photometry for the unusual UX Ori star V1184 Tau. Comparison with previous observations performed before the catastrophic decline in its optical brightness in 2004 (when the star faded approximately by a factor of 100) has shown the following: the star faded approximately by 2 m and 1 m in the J and H bands, respectively, while its K brightness remained almost constant. This pattern of infrared variability seems incompatible with the mechanism of variable circumstellar extinction responsible for the dramatic decline in the star’s optical brightness. However, if this mechanism is considered in the context of an accretion disk model with a puffed-up inner rim in the dust sublimation zone and with a disk wind producing an expanding gas-dust atmosphere above the disk surface, then the paradox can be resolved. In this model, the photometric activity of V1184 Tau in both visible and near-infrared spectral ranges, including the sharp brightness decline in 2004, can be explained by an increase in the geometric thickness of the disk in the dust sublimation zone caused by enhanced accretion of circumstellar matter onto the star. There is reason to believe that such events occur periodically and result from the presence of a companion to V1184 Tau moving in a highly eccentric orbit. The offered interpretation of the photometric activity of V1184 Tau allows this object to be classified as an UX Ori star based on the observed photometric effect and, at the same time, as a FU Ori star based on the pattern of the physical process that produced this effect.  相似文献   

10.
11.
Available velocity dispersion estimates for the old stellar population of galactic disks at galactocentric distances r?2L (where L is the photometric radial scale length of the disk) are used to determine the threshold local surface density of disks that are stable against gravitational perturbations. The mass of the disk Md calculated under the assumption of its marginal stability is compared with the total mass Mt and luminosity L B of the galaxy within r=4L. We corroborate the conclusion that a substantial fraction of the mass in galaxies is probably located in their dark halos. The ratio of the radial velocity dispersion to the circular velocity increases along the sequence of galactic color indices and decreases from the early to late morphological types. For most of the galaxies with large color indices (B–V)0>0.75, which mainly belong to the S0 type, the velocity dispersion exceeds significantly the threshold value required for the disk to be stable. The reverse situation is true for spiral galaxies: the ratios Md/LB for these agree well with those expected for evolving stellar systems with the observed color indices. This suggests that the disks of spiral galaxies underwent no significant dynamical heating after they reached a quasi-equilibrium stable state.  相似文献   

12.
13.
We examine the dependence of the total hydrogen mass M HI in late-type star-forming galaxies on rotation velocity V rot and optical size D 25 or radial scale length R 0 of the disk for two samples of galaxies: (i) isolated galaxies (AMIGA) and (ii) galaxies with edge-on disks (flat galaxies according to Karachentsev et al.). M HI given in the HYPERLEDA database for flat galaxies have turned out to be, on average, overestimated by ~0.2 dex compared to isolated galaxies with similar V rot or D 25, which is apparently due to an overestimation of the self-absorption in the HI line. The hydrogen mass in the galaxies of both samples closely correlates with the total specific angular momentum of the galactic disk J, which is proportional to V rot D 25 or V rot R 0, with the low-surface-brightness galaxies lying along the common V rot R 0 sequence. We discuss the possibility of explaining the relationship between M HI and V rot D 25 by assuming that the gas mass in the disk is regulated by the marginal gravitational stability condition for the gas layer. Comparison of the observed and theoretically expected dependences leads us to conclude that either the gravitational stability corresponds to higher values of the Toomre parameter than is usually assumed, or the threshold stability condition formost galaxies was fulfilled only in the past, when the gasmass in the disks was a factor of 2–4 higher than that at present (except for the galaxies with an anomalously high observed HI content). The latter condition requires that for most galaxies the conversion of gas into stars be not compensated by the external accretion of gas onto the disk.  相似文献   

14.
We have analyzed the optical (U BV) and ultraviolet (λ1000–2700 Å) observations of the nuclear variability of the Seyfert galaxy NGC 4151 in the period 1987–2001 (the second cycle of activity). The fast (tens of days) and slow (~10 years) components of the nuclear variability, F and S, respectively, are shown to be completely different, but thermal in nature. We associate the S component with the formation and evolution of an accretion disk and the F component (flares) with instabilities in the accretion disk and their propagation over the disk in the form of a shock wave. The S component is present not only in the optical, but also in the ultraviolet range, with its amplitude being comparable over the entire range λ1000–5500 Å under study. The amplitude of the average flare (the F component) doubles as the wavelength decreases from 5500 to 1000 Å, while the rise time of the brightness to its maximum Δt (the variability time scale) decreases from
to 6d ± 2d. The brightness decline (flare decay) time decreases by a factor of 16. The extinction in the ultraviolet is shown to have been grossly underestimated: beginning from the first IUE data, only the extinction in our Galaxy,
, has been taken into account. A proper allowance for the total extinction, i.e., for the extinction in the nucleus of NGC 4151 as well
leads to a large increase in the luminosity of the variable source in the nucleus of NGC 4151: L = (6–8) × 1046 erg s?1. The spectral energy distribution for the variable source (λ950–5500 Å) agrees well with two Planck distributions: Te = 65 000 (λmax = 450 Å) and 8000 K. The radiation with Te = 8000 K is the reprocessing of the bulk of the ultraviolet radiation by the accretion disk with a lag of 0.5–0.6 days in the V band. The lag in the U-B variability of the slow component revealed the existence of an extended broad line region (EBLR) at an effective distance of 1.5 lt-years, as confirmed by spectroscopic data obtained at the Crimean Astrophysical Observatory. This yields the following mass of the central object in NGC 4151: Mc = (1–3) × 109M. The luminosity of the variable source then accounts for 50–60% of LEdd rather than 1–2%, as has been thought previously. In general, the pattern of ultraviolet and optical variability in NGC 4151 agrees excellently with the theory of disk accretion instability for a supermassive black hole suggested by N. Shakura and R. Sunyaev 30 years ago: the energy release is at a maximumin the ultraviolet (in the case under consideration, at λ450 Å), the luminosity is ~1047 erg s?1 for Mc ~ 109M (several tens of percent of LEdd), and the variability time scale ranges from several days to many years.
  相似文献   

15.
We consider a spherically symmetric general relativistic perfect fluid in its comoving frame. It is found that, by integrating the local energy momentum conservation equation, a general form of g 00 can be obtained. During this study, we get a cue that an adiabatically evolving uniform density isolated sphere having ρ(r,t)=ρ 0(t), should comprise “dust” having p 0(t)=0; as recently suggested by Durgapal and Fuloria (J. Mod. Phys. 1:143, 2010) In fact, we offer here an independent proof to this effect. But much more importantly, we find that for the homogeneous and isotropic Friedmann-Robertson-Walker (FRW) metric having p(r,t)=p 0(t) and ρ(r,t)=ρ 0(t), \(g_{00} = e^{-2p_{0}/(p_{0} +\rho_{0})}\). But in general relativity (GR), one can choose an arbitrary tt ?=f(t) without any loss of generality, and thus set g 00(t ?)=1. And since pressure is a scalar, this implies that p 0(t ?)=p 0(t)=0 in the Big-Bang model based on the FRW metric. This result gets confirmed by the fact the homogeneous dust metric having p(r,t)=p 0(t)=0 and ρ(r,t)=ρ 0(t) and the FRW metric are exactly identical. In other words, both the cases correspond to the same Einstein tensor \(G^{a}_{b}\) because they intrinsically have the same energy momentum tensor \(T^{a}_{b}=\operatorname {diag}[\rho_{0}(t), 0,0, 0]\).  相似文献   

16.
An intense outburst of hard radiation (with a peak flux of ~50 mCrab) was detected from the X-ray transient AX J1749.1-2733 by the IBIS/ISGRI gamma-ray telescope onboard the INTEGRAL observatory when the Galactic center field was monitored on September 8–10, 2003. Previously, this source had never been observed in a bright X-ray state. During the outburst, the source’s radiation spectrum was gently sloping and hard (extended to ~100 keV), followed a power law in the standard X-ray energy range, and had an exponential cutoff above 40–50 keV. The spectral hardness decreased with increasing flux. These and other properties described here and the shortness of the outburst (<2 days) allow the source to be attributed to the group of fast X-ray transients many representatives of which have an early O-B supergiant as their optical counterpart. Possible causes of the outbursts of fast transients are discussed. We show that accretion from the supergiant’s stellar wind should have led to intense persistent radiation from transients. The absence of radiation can be explained by the ejection of accreting matter from the system (propeller effect) during its contact with the magnetosphere of a rapidly rotating neutron star. Transient outbursts could originate in sources of this type if the spin period of their neutron star P s differed only slightly from the critical period P s * ? 3 s at which the propeller effect is still possible. The outburst is triggered by an insignificant rise in the local stellar wind density, by a factor of (P s * /P s)7/3. The entire outburst profile cannot be explained by an individual inhomogeneity in the wind, but is the reflection of a long-term (~2 days for AX J1749.1-2733) change in the rate of wind outflow from the supergiant’s surface facing the compact source. The rate of wind outflow could be enhanced through X-ray heating of the supergiant’s surface.  相似文献   

17.
The two-dimensional structure of a thin accretion disk in the vicinity of a Schwarzschild black hole after passing a marginally stable orbit (r< 3r g is discussed in terms of the Grad-Shafranov hydrodynamic equation. The accretion disk is shown to be sharply compressed as the sonic surface is approached, so the mass flow here is no longer radial. As a result, the dynamic forces ρ[(v ?)v] θ , which are equal in magnitude to the pressure gradient ? θ P on the sonic surface, become significant in vertical balance. Therefore, the disk thickness in the supersonic region (and, in particular, near the black-hole horizon) may be assumed to be determined not by the pressure gradient but by the shape of ballistic trajectories.  相似文献   

18.
We calculate the parameters of the two-point correlation function of quasars w(r) = (r c /r) γ on the basis of the SDSS DR3 data. The correlation functions are first determined from projected distances with the use of a special technique for compiling randomized catalogs. Next the parameters of the spatial correlation function are obtained with the assumption of local isotropy. For the quasars with redshifts z = 0.8–2.1, we obtained the estimates γ = 1.76 ± 0.14, r c = 6.60 ± 0.85 h ?1 Mpc in the comoving distance range 2–30 Mpc and γ = 1.90 ± 0.11, r c = 6.95±0.57 h ?1 Mpc in the range 2–50 Mpc. These estimates agree, within the limits of errors, with the estimates obtained for the redshifts 0.4 < z < 2.1. The original catalog shows some deficit of pairs with separations less than 1 Mpc.  相似文献   

19.
We determine the color excesses, photometric distances, ages, astrometric parallaxes and proper motions for 94 open clusters in the northern part of the Milky Way. We estimate the color excesses and photometric distances based on the data from IPHAS photometric survey of the northern Galactic plane using individual total-to-selective extinction ratios Rr = Ar/Er?i for each cluster computed via the color-difference method based on IPHAS r, i, and Hα-band, 2MASS J-, H-, and Ks-band, WISE W1-band, and Pan-STARRS i-, z-, and y-band data. The inferred Rr values vary significantly from cluster to cluster spanning the Rr = 3.1–5.2 interval with a mean and standard deviation equal to 〈Rr〉 = 3.99 and σRr = 0.34, respectively.We identified cluster members using (1) absolute proper motions determined from individual-epoch positions of stars retrieved from IPHAS, 2MASS,URAT1, ALLWISE,UCAC5, and Gaia DR1 catalogs and positions of stars on individual Palomar Sky Survey plates reconstructed based on the data provided in USNO-B1.0 catalog and (2) absolute proper motions provided in Gaia DR2 catalog, and computed the averageGaia DR2 trigonometric parallaxes and propermotions of the clusters. Themean formal error of the inferred astrometric parallaxes of clusters is of about 7 μas, however, a comparison of astrometric and photometric parallaxes of our cluster sample implies that Gaia DR2 parallaxes are, on the average, systematically underestimated by 45 ± 9 μas. This result agrees with estimates obtained by other authors using other objects. At the same time, we find our photometric distance scale to be correct within the quoted errors (the inferred correction factor is equal to unity to within a standard error of 0.025).  相似文献   

20.
We investigate the properties of an axisymmetric gas flow without angular momentum onto a small compact object, in particular, on a Schwarzschild black hole in the supersonic region; the velocity of the object itself is assumed to be low compared to the speed of sound at infinity. First of all, we show that the streamlines intersect (i.e., a caustic is formed) on the symmetry axis at a certain distance r x from the center on the front side if the pressure is ignored. The characteristic radial size of the region in which the streamlines emerging from the sonic surface at an angle no larger than θ0 to the axis intersect is Δr = r x θ 0 2 /3. To refine the flow structure in this region, we have numerically computed the system without ignoring the pressure in the adiabatic approximation. We have estimated the parameters of the inferred region with anomalously high matter temperature and density accompanied by anomalously high energy release.  相似文献   

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

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