首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Abstract

It is demonstrated in laboratory experiments with rotating shallow water that large scale Rossby vortices, greater than the Rossby-Obukhov radius in size, have dispersive and non-linear properties that are fundamentally different for the two possible polarities. We call this “cyclonic-anticyclonic asymmetry”. This asymmetry manifests itself in the following way: first, anticylones, unlike cyclones, do not undergo the dispersive spreading inherent in a linear wave packet. and therefore, having a considerably longer natural lifetime, are obvious candidates for Rossby solitons; second, dipolar vortices are, because of the comparatively rapid decay of a cyclone, transformed into anticyclonic solitons; third, anticyclones are much more readily generated by zonal flows of the type existing in planetary atmospheres. The evident dominance of anticyclones amongst the long-lived vortices in the atmospheres of giant planets strongly suggests that the cyclonic-anticyclonic symmetry plays a decisive role in the atmospheric cyclogenesis of large planets.

According to our concept, the Rossby soliton is a “real” vortex; unlike a wave, it retains some fluid particles within it throughout its lifetime. Two similar solitons can merge by mutual collisions. This picture of a “vortical” soliton differs in an essential way from the earlier idea due to Maxworthy and Redekopp (1976) of purely “wave-like” Rossby solitons that can freely pass through one another.

Laboratory experiments were performed by us to simulate the new Rossby solition, with special reference to naturally-occurring vortices of the same general type as Jupiter's great red spot. The experimental data presented contradict the idea of “pure wave solitons” but confirm our concept of “vortical solutions”.  相似文献   

2.
Large, long-lived vortices are abundant in the atmospheres of the giant planets. Some of them survive a few orders of magnitude longer than the dispersive linear Rossby wave packets, e.g. the Great Red Spot (GRS), Little Red Spot (LRS) and White Ovals (WO) of Jupiter, Big Bertha, Brown Spot and Anne's Spot of Saturn, the Great Dark Spot (GDS) of Neptune, etc. Nonlinear effects which prevent their dispersion spreading are the main subject of our consideration. Particular emphasis is placed on determining the dynamical processes which may explain the remarkable properties of observed vortices such as anticyclonic rotation in preference to cyclonic one and the uniqueness of the GRS, the largest coherent vortex, along the perimeter of Jupiter at corresponding latitude.We review recent experimental and theoretical studies of steadily translating solitary Rossby vortices (anticyclones) in a rotating shallow fluid. Two-dimensional monopolar solitary vortices trap fluid which is transported westward. These dualistic structures appear to be vortices, on the one hand, and solitary waves, on the other hand. Owing to the presence of the trapped fluid, such solitary structures collide inelastically and have a memory of the initial disturbance which is responsible for the formation of the structure. As a consequence, they have no definite relationship between the amplitude and characteristic size. Their vortical properties are connected with geostrophic advection of local vorticity. Their solitary properties (nonspreading and stationary translation) are due to a balance between Rossby wave dispersion and nonlinear effects which allow the anticyclones, with an elevation of a free surface, to propagate faster than the linear waves, without a resonance with linear waves, i.e. without wave radiation. On the other hand, cyclones, with a depression of a free surface, are dispersive and nonstationary features. This asymmetry in dispersion-nonlinear properties of cyclones and anticyclones is thought to be one of the essential reasons for the observed predominance of anticyclones among the long-lived vortices in the atmospheres of the giant planets and also among the intrathermocline oceanic eddies.The effects of shear flows and differences between the properties of monopolar vortices in planetary flows and various laboratory experiments are discussed. General geostrophic (GG) theory of Rossby vortices is presented. It differs essentially from the traditional quasi-geostrophic (QG) and intermediate-geostrophic (IG) approximations by the account of (i) all scales between the deformation radius and the planetary scale and (ii) the arbitrary amplitudes of vortices. It is shown that, unlike QG- and IG-models, the GG-model allows for explaining the mentioned cyclonic-anticyclonic asymmetry not only in planetary flows, but also in laboratory modeling with vessels of near paraboloidal form.  相似文献   

3.

We examine the three-dimensional, nonlinear evolution of columnar vortices in a rotating environment. As the initial vorticity distribution, a wavetrain of finite amplitude Kelvin-Helmholtz vortices in shear is employed. Through direct numerical simulation of the Navier-Stokes equations we seek to better understand the process of maturation of the various three-dimensional modes of instability to which such vortical flows are subject, especially those which exist as a consequence of the action of the Coriolis force. In the absence of rotational influence, we thereby demonstrate that the nonlinear evolution of columnar vortices is most strongly controlled by one or the other of two mechanisms. One mechanism of instability is identifiable as a so-called elliptical instability, which promotes the initial bending of vortex tubes in a sinusoidal fashion, while the other is a hyperbolic mode, which is responsible for the development of streamwise vortex streaks in the "braids" between adjacent vortex cores. In the rotating case, anticyclonic vortices are strongly destabilized by weak background rotation, while rapid rotation stabilizes both the cyclones and anticyclones. The strong anticyclones are subject to two distinct forms of instability, namely a Coriolis force modified elliptical instability and an inertial (centrifugal) instability. The former instability is very similar to the nonrotating form of the elliptical instability as it promotes bending of vortex tubes, while the latter instability grows on the edge of the vortex core and generates streaks of vorticity, which surround the vortex core itself. These results of direct numerical simulation fully verify the results of previous linear stability analyses. Taken together, they provide a simple explanation for the broken symmetry that is often observed to be characteristic of the von Karman vortex streets that develop in the atmospheric lee of oceanic islands.  相似文献   

4.
The dynamics of a single vortex on a beta-plane is discussed in this paper. A barotropic, an equivalent barotropic, one-and-a half and two-layer models are considered. The momentum and energy balances are used to describe the evolution of a vortex. A quasi-stationary balance of the Rossby, Zhukovsky-Kutta forces and the force induced by Rossby-wave radiation, describes the dynamics of the barotropic vortex. A net Coriolis force occurs if the fluid is stratified. The difference between the dynamics of cyclones and anticyclones results directly from the Coriolis force acting on a single vortex in a stratified fluid.All vortices radiate Rossby waves in the quasigeostrophic approximation but intense anticyclones propagate steadily in a one-and-a half layer model. A critical amplitude that bounds radiating and steadily propagating anticyclones is found. Steady propagation of anticyclones in general is impossible in a two-layer fluid due to the radiation of a barotropic Rossby-wave. Some solutions of solitary wave type which are known for a two-layer model, survive owing to wave interference.A single vortex can extract energy from a Rossby wave if synchronism conditions are satisfied. The wave interference again plays a crucial role in this case. The wave interference also determines the energy exchange of vortices located at larger distances. If the distance between the vortices is shorter than the length of the radiated waves, modon may be formed due to a small energy loss.The unbounded monotonic variation of the planetary vorticity is a characteristic feature of a beta-plane approximation. As a result, a single vortex propagates up to a 'rest latitude' where it disappears. The evolution of a single barotropic vortex over bottom topography provides another example of a background vorticity distribution with a local extremum above hills (valleys) or ridges (troughs). Physics of its movement differs from a beta-plane case, but if a vortex lies over broad topography, equations are similar and the evolution of a vortex manifests the same typical features. Particularly, a cyclonic vortex tends to drift to the top of a hill or a ridge. An anticyclonic vortex, on the contrary, slides to the bottom of a valley or a trough.An interaction of a barotropic vortex with a broad mean flow is tractable qualitatively on the basis of previous results. Numerical examples illustrating absorption of a small vortex by a larger one and a vortex movement across the flow, are direct analogies of the vortex evolution over a hill and a ridge, respectively. At the same time, strong influence of strain drastically changes the vortex structure.  相似文献   

5.

We present results from a new series of experiments on the geophysically important issue of the instability of anticyclonic columnar vortices in a rotating fluid in circumstances such that the Rossby number exceeds unity. The vortex pair consisting of a cyclonic and an anticyclonic vortex is induced by a rotating flap in a fluid which is itself initially in a state of solid-body rotation. The anticyclonic vortex is then subject to either centrifugal or elliptical instability, depending on whether its initial ellipticity is small or large, while the cyclone always remains stable. The experimental results demonstrate that the perturbations due to centrifugal instability have a typical form of toroidal vortices of alternating sign (rib vortices). The perturbations due to elliptical instability are of the form of sinuous deformation of the vortex filament in the plane of maximal stretching which corresponds to the plane of symmetry for the vortex pair. The initial perturbations in both cases are characterized by a definite wave number in the vertical direction. The characteristics of the unstable anticyclone are determined by the main nondimensional parameter of the flow - the Rossby number. The appearance of both centrifugal and elliptical instabilities are in accord with the predictions of theoretical criteria for these cases.  相似文献   

6.
Vortex merger is a phenomenon characterizing the whole class of geophysical vortices, from atmospheric storms and large oceanic eddies up to small scale turbulence. Here we focus on the merger of subsurface oceanic anticyclones in an idealized primitive equations model. This study has been motivated by past and recent observations of colliding lens-like anticyclones off of Gibraltar Strait. The critical conditions for merger (critical merger distance and time needed for merger) are determined. We will show that the predictions of classical two-dimensional merger are not verified for subsurface isolated vortices. For instance, critical merger distances will be reduced because of the vortex potential vorticity (PV) structure. The post-merger characteristics of the vortex (radius, extension and PV), are also determined. Merger-related effects, like production of peripheral filaments and small-scale eddies are also investigated and suggest the contribution of merger in both direct and inverse energy cascades.  相似文献   

7.
In a previous paper, Caillol [Geophys. Astrophys. Fluid Dyn., 2014, 108] investigated the steady nonlinear vortical structure of a singular vortex Rossby mode that has survived to a strong critical-layer-like interaction with a linearly stable, columnar, axisymmetric and dry vortex. We presented a general theory for this wave/mean flow interaction through the nonlinear critical layer theory and calculated the mean azimuthal and axial winds induced at the critical radius at the end of this interaction in the final stage. We here apply that theory to rapidly rotating geophysical vortices: tropical cyclones, cold-air mesocyclones and tornadoes. We find that the numerous assumptions invoked in that paper agree well with the reality of those intense vortices. We also find that in spite of a lack of moist-convection modelling, this dry vortex is fairly well accelerated at the critical radius by such a shear wave with a magnitude of order the square root of the damped-wave amplitude. The intensification level strongly depends on the aspect ratio, height of the system: rapid vortex and parent vortex, over core radius. The thinner the vortex is, the sharper the intensification is. This result is in sharp contrast to the numerous numerical simulations on VR wave/vortex interactions that yield a much smaller intensification of order the square of the wave amplitude. This weakly nonlinear approach nevertheless fails to model small vertical wavelength VR wave/vortex interactions for their related asymptotic expansions are divergent and for they yield strongly nonlinear VR waves coupled with evolving critical layers whose extent can no longer be considered as thin.  相似文献   

8.
The interaction of meddies with a complex distribution of seamounts is studied in a three-layer quasi-geostrophic model on the f-plane. This study aims at understanding if and how this seamount chain can represent a barrier to the propagation of these eddies and how it can be involved in their decay. The eddies are idealized as vortex patches in the middle layer, interacting with a regional cyclonic current and with ten idealized seamounts. The numerical code is based on the contour surgery technique. The initial position, radius, shape, number and polarity of the eddies are varied. The main results are the following: (1) Though they do not describe the unsteady flow, the streamlines of the regional and topographic flow provide a useful estimate of the vortex trajectories, in particular towards the major seamounts, where stronger velocity shears are expected. (2) The tallest and widest seamounts which have the largest vorticity reservoir are able to considerably erode the vortices, but also to draw anticyclones towards the seamount top. The ability of narrower seamounts to erode vortices is related to their multiplicity. (3) Only 1/3 of the anticyclones with about 30-km radius reach the southern boundary of the seamount chain, and their erosion is larger than 50 %. The other anticyclones are either completely eroded or trapped over a wide seamount top. Cyclones are less affected by seamounts because they oppose the topographic draft towards the seamount top and they drift along the side of the seamount. (4) Large vortices resist topographic erosion more efficiently. The rate of erosion grows from a few percent to about 35–50 % as the vortex radius decreases from about 60 to 30 km. Small cyclones are not eroded, contrary to small anticyclones (which completely decay), in relation with the different trajectories of these eddies in the vicinity of the seamounts. (5) The detailed vortex shape does not appear critical for their evolution, if they are close enough to the seamount chain initially. The interaction between a group of vortices initially north of the seamount chain can modify their trajectory to such an extent that they finally avoid collision with seamounts. (6) Finally, meddy trajectories across the Horseshoe Seamounts (data from the AMUSE experiment) show qualitative similarity with the vortex paths in the model. Several events of vortex decay also occur at comparable locations (in particular over the wide and tall seamounts) in the model and observations.  相似文献   

9.
Abstract

Supercritically unstable density fronts near a vertical wall in a rotating, two-layer fluid were created on a laboratory turntable by withdrawing the outer wall of an annulus with a narrow gap, and allowing buoyant fluid from within the annulus to collapse toward a state of quasi-geostrophic balance. The resulting “coastal” current has a nearly uniform potential vorticity and is bounded by a front on which ageostrophic, wave-like disturbances grow. If the current width is comparable to the Rossby radius of deformation, the dominant length scale of disturbances is proportional to the width of the current. On the other hand, if the upper layer is much wider than the Rossby radius, then the observed length scale is a constant multiple of the Rossby radius. If the vertical boundary is omitted in the experiments, so that we are left with a circular anticyclonic vortex, the observed length scales and large-amplitude behaviour of disturbances are identical to those for the boundary currents, indicating that the wall has no significant influence on the flow.

At very large amplitude the growing waves lead to the formation of cyclone-anticyclone vortex pairs. For very wide currents, both the mean flow and the disturbances are first confined to a region within a few Rossby radii of the front. However, both the mean flow and the turbulent eddy motions slowly propagate into the previously stationary upper layer until, eventually, the whole of the upper layer is turbulent.  相似文献   

10.
We examine the equilibrium form, properties, stability and nonlinear evolution of steadily-rotating simply-connected vortex patches in the single-layer quasi-geostrophic model of geophysical fluid dynamics. This model, valid for rotating shallow-water flow in the limit of small Rossby and Froude numbers, has an intrinsic length scale L D called the “Rossby deformation length” relating the strength of the stratification to that of the background rotation. Here, we generate steadily-rotating vortex equilibria for a wide range of γ?=?L/L D , where L is the typical horizontal length scale of the vortex. We vary both γ (over the range 0.02?≤?γ?≤?10) and the vortex aspect ratio λ (over the range 0?<?λ?<?1). We find two modes of instability arising at sufficiently small aspect ratio λ?<?λ c (γ): an asymmetric (dominantly wave 3) mode at small γ (or large L D ) and a symmetric (dominantly wave 4) mode at large γ (or small L D ). At marginal stability, the asymmetric mode dominates for γ???3, while the symmetric mode dominates for γ???3. The nonlinear evolution of weakly-perturbed unstable equilibria results in major structural changes, in most cases producing two dominant vortex patches and thin, quasi-passive filaments. Overall, the nonlinear evolution can be classified into three principal types: (1) vacillations for a limited range of aspect ratios λ when 5?≤?γ?≤?6, (2) filamentation and a single-dominant vortex for γ???1, and (3) vortex splitting – asymmetric for 1???γ???4 and symmetric for γ???4.  相似文献   

11.
Abstract

Various interactions between small numbers (two and four) of baroclinic, geostrophic point vortices in a two-layer system are studied with attention to the qualitative changes in behavior which occur as size of the deformation radius is varied.

A particularly interesting interaction, which illustrates the richness of baroclinic vortex dynamics, is a collision between two hetons. (A heton is a vortex pair in which the constituent vortices have opposite signs and are in opposite layers. The “breadth” of a heton is the distance between its constituent vortices. A translating heton transports heat.) When two hetons, which initially have different breadths, collide, the result is either an exchange of partners, or a “slip-through” collision in which the initial structures are preserved. It is shown here that the outcome is always an exchange, provided the deformation radius is sufficiently small. This strongly contrasts with a collision between pairs of classical, one-layer vortices in which no exchange occurs if the initial ratio of the breadths is sufficiently extreme.

Finally the transport of passive fluid by a translating baroclinic pair is investigated. A pair of vortices in the top layer transports no lower layer fluid if the distance between the vortices is less than 1.72 deformation radii. By contrast, the size of the region trapped by a heton increases without bound as the spacing between the vortices increases.  相似文献   

12.
刘爽  钟玮  刘宇迪 《地球物理学报》2018,61(6):2207-2219
本文基于正压浅水模型,分析基态位涡(Potential Vorticity:PV)结构对热带气旋(Tropical Cyclone:TC)类涡旋系统稳定性及其波动特征的影响.通过引入基态PV结构参数:宽度δ(眼墙内外边界涡度发生陡变的半径长度之比)和中空度γ(眼心相对涡度与内核区域平均相对涡度之比),设计具有相同基流最大切向风速和最大风速半径的170组不同基态PV环结构的敏感性试验,并讨论了不同基态PV结构下涡旋系统最不稳定波数(the most unstable wavenumber:MUWN)和系统最不稳定模态(the most unstable mode of System:MUMS)的特征频率及其不稳定增长率的大小.结果指出:当PV环较宽,系统表现为低波数最不稳定,相应的MUMS为低频波且增长率小;当PV环较窄,系统表现为高波数不稳定,且PV环越实最不稳定波数越高;当PV环窄且空时,MUMS均为中高频波动,且不稳定增长率随PV环的宽度变窄和中空度变空而明显增大.分析典型PV结构下系统演变特征可知,当PV环较宽,MUMS表现为具有平衡约束的低频波动的线性不稳定特征;当PV环趋向窄且空时,MUMS的平衡性约束趋向弱化,同时不稳定增长表现为明显的指数型增长.进一步讨论系统内部非对称结构的形成和传播机制发现,对于弱不稳定的PV环来说,低波数波最不稳定的特征波动具有典型涡旋Rossby波特征;而对于强不稳定的PV环来说,高波数不稳定的特征波动混合波性质明显.  相似文献   

13.
Abstract

spin-up and spin-down in a circular tank with a uniformly sloping bottom are studied experimentally and numerically for small values of the relative change in the angular velocity of the tank. Generally, the initial single-cell flow evolves into a number of smaller vortices. The evolution is compared with an analytical model based on an expansion of the flow field in linear Rossby waves (Pedlosky and Greenspan, 1967). Although it is possible to tune the experimental parameters in such a way that agreement with the theory is found, in most cases the experiments show shedding of vortices in the initial stage of the spin-up or spin-down, a phenomenon not described by the analytical model. Nonetheless, in such cases the analytical model still accounts for other observations: the alternating generation of cyclonic and anticyclonic vortices in the eastern part of the tank and their subsequent westward motion.  相似文献   

14.
Based on an eddy-permitting numerical model, the mesoscale variability in the East-Sakhalin Current is investigated during the winter-spring period. Analysis of necessary conditions for the development of baroclinic instability showed that the nearshore component of the East-Sakhalin Current is potentially baroclinic unstable in the first half-year. The simulated circulation uncovered a generation of anticyclonic eddies on the eastern Sakhalin shelf. It was established that a spatial scale of these eddies and the first baroclinic Rossby radius of deformation are values of the same order; a lifetime of these eddies varies from 4 to 6 weeks, given the Rossby number varies from 0.05 to 0.2. Analysis of the rate of eddy energy conversion on the eastern Sakhalin shelf showed that the generation of the revealed mesoscale eddies results from, mainly, baroclinic instability, whereas barotropic instability can be both favoring and preventing to the generation of these eddies.  相似文献   

15.
Abstract

An investigation is made of the evolution of small-scale, axisymmetric vortices in a stratified fluid with spatially variable Coriolis parameter. The criteria for smadness are a horizontal scale less than or equal to the first internal radius of deformation and a vertical scale less than or equal to that of the ambient stratification. These circumstances match those of Sub-mesoscale, Coherent Vortices frequently observed in the oceans. The dynamical model is the balance equations, which include various effects of finite Rossby number. The principal topics are the regime of nearly uniform propagation, the development of an equilibrium ratio of vertical and horizontal scales (i.e., Burger number selection), and the occurrence of various types of instability for vortices with extremes in amplitude or shape.  相似文献   

16.
17.
Abstract

We study the formation of lenses of the ocean's intermediate water using a 2.5-layerβ-plane primitive equation model with localized injection of water mass. For the injecting rate of 1.0 Sv, we have observed that strong vortices are shed regularly. These vortices propagate westward much faster than the second baroclinic long Rossby wave. They are totally isolated from each other and show strong baroclinicity as well. Moreover, they remain stable over a sufficiently long period of time. Regular formation of such strong vortices in the intermediate layer has not been reported previously. The translation speed is explained using the Euler's momentum integral theorem for the nonlinear baroclinic vortex on the β-plane. We have demonstrated that coupling between the primary motion in the intermediate layer and the secondary motion in the upper layer with a meridional shift is crucial to the fast westward translation of the intense vortices. A simple dispersion formula relating the zonal translation speed with the vortex radius is also derived under the assumption of quasi-geostrophy. It has turned out that the analytical relation explains the numerical results surprisingly well despite the limitation of its derivation.  相似文献   

18.
Abstract

In this paper we use the CASL method to explore the role of boundary conditions in determining the long-time behaviour of rotating, stratified, quasi-geostrophic turbulence. We find that initially two-dimensional (sufficiently tall) columns of potential vorticity (PV) break down through three-dimensional instability to give a fully three-dimensional flow consisting of ellipsoidal structures. This is the case both for rigid-lid (isothermal) vertical boundary conditions and for vertically periodic boundaries. However, the rigid boundary case gives rise to semi-ellipsoids at both the top and bottom boundaries, and, for sufficient domain depths, preferred depths for the formation of ellipsoids in the interior. By contrast, in the vertically periodic case, the distribution of ellipsoids is homogeneous in depth.

The role of the horizontal boundaries is indirect, but still significant. In all cases doubly periodic horizontal boundary conditions are imposed. We consider a range of initial conditions where in each case equal numbers of two-dimensional columns of positive and negative vorticity are used, taking up a fixed, but relatively small fraction of the domain (approximately 5%). Thus when there is only a small number of vortices, they have larger radius. When the initial number of vortices is small enough (i.e., when the radius is not small compared with the horizontal domain width), at long time there is a two-dimensionalisation giving rise to a single column of positive PV and a single column of negative PV, as has been observed in some previous simulations. We find the same phenomenon for both vertically periodic and rigid lid boundary conditions, but it occurs over a broader range of initial conditions in the vertically periodic case. However, in all cases fully three-dimensional final states are regained when the number of vortices is increased while keeping the fraction of the domain occupied by vortices fixed, i.e., when the vortex radius/domain width ratio is sufficiently small.  相似文献   

19.
The problem of zonal jet formation and cyclone–anticyclone asymmetry in decaying rotating turbulence is addressed using both laboratory experiments and numerical simulations with a high-resolution shallow water model in a spherical geometry. Experiments are performed at different Rossby and Froude numbers and applying a rigid wall as meridional boundary in the numerical scheme to mimic the experimental apparatus. The formation of a zonally banded flow pattern, i.e. meridionally confined easterly/westerly jets, has observed; both experimental and numerical results confirmed that this tendency is favoured by high-planetary vorticity gradients. Also, in the experiments characterized by large rotation speeds and small Rossby deformation radius, an initial symmetric distribution of relative vorticity is found to evolve towards a dominance of anticyclonic structures, indicating a breaking of the cyclone–anticyclone symmetry. This aspect has deepened by numerically analysing the sensitivity of the temporal variations of the asymmetry index with respect to the position of the meridional confinement as well as the effect of relaxing the divergence of the fluid (i.e. non-divergent case) to zero. Results suggested that experiments characterized by the higher rotation speed and the lower fluid thickness are better reproduced by a divergent model with a high-latitude meridional boundary.  相似文献   

20.
The interaction of warm core rings with a western wall and shelf/slope is examined with a three-dimensional primitive-equation model. The model ring is initialized with an axisymmetric Gaussian-type anticyclonic eddy placing far from the coastal boundary to allow the ring to freely propagate towards the wall and shelf/slope. The ring initially propagates steadily to the southwest at about 3 km/day under the combined planetary β and nonlinear effects. When colliding with a wall, the ring adjusts into a ‘D’ shape and moves poleward under primarily the image effect. When colliding with a shelf and slope, the ring however becomes stalled and bounces on and off the shelf/slope with little net movement. Small cyclones marked by strong upwelling are generated near the shelfbreak. Cyclones and anticyclones also are spawned at the periphery of the ring. Satellite SST images and concurrent ADCP transects are used to illustrate the strong interaction of a Gulf Stream warm ring (99B) with the Middle Atlantic Bight.  相似文献   

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

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