首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
A technique of high-image-density particle image velocimetry is employed to characterize the instantaneous and averaged patterns of velocity, vorticity and Reynolds stress due to flow past two cylinders in tandem. These features of the flow patterns are characterized in the gap region as a function of the distance between the cylinders. In turn, they are related to the patterns in the near-wake of the two-cylinder system. Along the gap between the cylinders, small-scale concentrations of vorticity are formed in the separated shear layers. These concentrations buffet the surface boundary layer on the downstream cylinder, and thereby influence the eventual shedding of large-scale vortices. Within the gap, the instantaneous structure of the recirculation zones can exhibit both symmetrical and asymmetrical patterns. In the near-wake of the downstream cylinder, the form of the vortex shedding, as well as the averaged patterns of the flow structure, are substantially altered, relative to the case of a single cylinder. The width of the near-wake, as represented by averaged patterns of vorticity, is substantially narrower and the magnitudes of the peak Reynolds stress are significantly attenuated. On the other hand, if the gap region is sufficiently large such that Kármán-like vortices form between the cylinders, the near-wake of the downstream cylinder shows distinctive patterns, and both the wake width and the magnitude of the Reynolds stresses become larger, relative to those at smaller gap width.  相似文献   

2.
 The instantaneous, quantitative patterns of vortices arising from sinusoidal oscillation of a cylinder in quiescent fluid are experimentally characterized for the first time using high-image-density particle image velocimetry. The near-wake does not indicate a separated layer of distributed vorticity leading to a single, large-scale vortex. Rather, for sufficiently high Reynolds number, a sequence of small-scale vorticity concentrations is formed. Agglomeration of only a fraction of the adjacent concentrations forms a larger-scale vortex. Simultaneously, vorticity concentrations of opposite sense are formed along the base (rear) of the cylinder. Streamline patterns typically indicate, however, only the larger-scale vortex; it has a circulation smaller than the total circulation of all vorticity concentrations that are not revealed by the streamlines. These observations are interpreted in the context of the effective resolution of the flow images. Received: 27 October 1995 / Accepted: 27 August 1996  相似文献   

3.
The division of flow regimes in a square cylinder wake at various angles of attack (α) is studied. This study provides evidence of the existence of modes A and B instabilities in the wake of an inclined square cylinder. The critical Reynolds numbers for the inception of these instability modes were identified through the determination of discontinuities in the Strouhal number versus Reynolds number curves. The spectra and time traces of wake streamwise velocity were observed to display three distinct patterns in different flow regimes. Streamwise vortices with different wavelengths at various Reynolds numbers were visualized. A PIV technique was employed to quantitatively measure the parameters of wake vortices. The wavelengths of the streamwise vortices in the modes A and B regimes were measured by using the auto-correlation method. From the present investigation, the square cylinder wake at various angles of attack undergoes a similar transition path to that of a circular cylinder, although various quantitative parameters measured which include the critical Reynolds numbers, spanwise wavelength of secondary vortices, and the circulation and vorticity of wake vortices all show an α dependence.  相似文献   

4.
Three-dimensional vortical structures have been measured in a circular-cylinder wake using particle imaging velocimetry (PIV) for the Reynolds number range of 2×103 to 1×104. The PIV was modified, compared with the conventional one, in terms of its light sheet arrangement to capture reliably streamwise vortices. While in agreement with previous reports, the presently measured spanwise structures complement the data in the literature in the streamwise evolution of the near-wake spanwise vortex in size, strength, streamwise and lateral convection velocities, shedding new light upon vigorous interactions between oppositely signed spanwise structures. The longitudinal vortices display mushroom patterns in the (x, z)-plane in the immediate proximity to the cylinder. Their most likely inclination in the (x, y)-plane is inferred from the measurements in different (x, z)-planes. The longitudinal vortices in the (y, z)-plane show alternate change in sign, though not discernible at x/d > 15. They decay in the maximum vorticity and circulation rapidly from x/d = 5 to 10 and slowly for x/d > 10, and are further compared with the spanwise vortices in size, strength and rate of decay.  相似文献   

5.
On the study of vortex-induced vibration of a cylinder with helical strakes   总被引:1,自引:0,他引:1  
While the effect of helical strakes on suppression of Vortex-Induced Vibrations (VIV) has been studied extensively, the mechanism of VIV mitigation using helical strakes is much less well documented in the literature. In the present study, a rigid circular cylinder of diameter d=80 mm attached with three-strand helical strakes of dimensions of 10d in pitch and 0.12d in height was tested in a wind tunnel. It was found that the helical strakes can reduce VIV by about 98%. Unlike the bare cylinder, which experiences lock-in over the reduced velocity in the range of 5-8.5, the straked cylinder does not show any lock-in region. In exploring the mechanism of VIV reduction by helical strakes, measurements in stationary bare and straked cylinder wakes using both a single X-probe at four different Reynolds numbers, i.e. Re=10 240, 20 430, 30 610 and 40 800, and two X-probes with variable separations in the spanwise direction at Re=20 430 were conducted. It was found that vortices shed from the straked cylinder are weakened significantly. The dominate frequency varies by about 30% over the range of x/d=10-40 in the streamwise direction while that differs by about 37.2% of the averaged peak frequency over a length of 3.125d in the spanwise direction. The latter is supported by the phase difference between the velocity signals measured at two locations separated in the spanwise direction. The correlation length of the vortex structures in the bare cylinder wake is much larger than that obtained in the straked cylinder wake. As a result, the straked cylinder wake agrees more closely with isotropy than the bare cylinder wake. Flow visualization on the plane perpendicular to the cylinder axis at Reynolds number of about 300 reveals small-scale vortices in the shear layers of the straked cylinder wake. However, these vortices do not roll up and interact with each other to form the well-organized Karman-type vortices. Flow visualization on the plane parallel to the cylinder axis shows vortex dislocation and swirling flow, which should be responsible for the variations of the peak frequency in the streamwise as well as spanwise directions.  相似文献   

6.
PIV measurements of the near-wake behind a sinusoidal cylinder   总被引:2,自引:0,他引:2  
The three-dimensional near-wake structures behind a sinusoidal cylinder have been investigated using a particle image velocimetry (PIV) measurement technique at Re=3,000. The mean velocity fields and spatial distributions of ensemble-averaged turbulence statistics for flows around the sinusoidal and corresponding smooth cylinders were compared. The near-wake behind the sinusoidal cylinder exhibited pronounced spanwise periodic variations in the flow structure. Well-organized streamwise vortices with alternating positive and negative vorticity were observed along the span of the sinusoidal cylinder. They suppress the formation of the large-scale spanwise vortices and decrease the overall turbulent kinetic energy in the near-wake of the sinusoidal cylinder. The sinusoidal surface geometry significantly modifies the near-wake structure and strongly controls the three-dimensional vortices formed in the near-wake.  相似文献   

7.
Vortex formation from a vertical cylinder in shallow water is controlled by placement of a narrow transverse strip of roughness elements on the bed (bottom surface). A technique of high-image-density particle image velocimetry is employed to obtain global, instantaneous representations of the flow patterns, which lead to phase- and time-averaged patterns of streamline topology and Reynolds stress on planes at and above the bed. Near the bed, the overall form of the streamline topology is maintained, even at larger heights of the roughness elements. With increasing height of the elements, the downstream saddle point is further displaced in the streamwise direction. Correspondingly, the streamwise extent of the negative pocket of the streamwise velocity component, i.e., the region of reverse flow along the bed surface, increases substantially in the streamwise direction. The Reynolds stress in the very near-wake, at locations upstream of the roughness elements, is significantly attenuated, even for small height of roughness. This attenuation occurs not only near the bed surface, but also at the midplane of the shallow water wake, and thereby indicates that the consequence of localized roughness is to exert a global influence. In fact, corresponding patterns of instantaneous velocity and vorticity indicate that consistent formation of large-scale vortices in the very near-wake region is attenuated with relatively small surface roughness on the bed. Downstream of the roughness elements, the patterns of Reynolds stress near the bed surface, as well as at the midplane of the water layer, are significantly altered relative to the case of no roughness. Near the bed, highly concentrated patterns of positive and negative Reynolds stress in the absence of roughness give way to lower-level regions of Reynolds stress in the form of alternating concentrations; the particular pattern depends on the height of the roughness elements. At the midplane of the water layer, the Reynolds stress patterns maintain their same overall form, but the extrema of the Reynolds stress concentrations are attenuated in magnitude and are shifted in the downstream direction, with increasing height of the roughness elements. These observations are complemented by patterns of instantaneous velocity and vorticity.  相似文献   

8.
Different flow models have been proposed for the flow around surface-mounted finite-height square prisms, but there is still a lack of consensus about the origin and connection of the streamwise tip vortices with the other elements of the wake. This numerical study was performed to address this gap, in addition to clarifying the relationship of the near-wake structures with the far wake and the near-wall flow, which is associated with the fluid forces. A large-eddy simulation approach was adopted to solve the flow around a surface-mounted finite-height square prism with an aspect ratio of AR = 3 and a Reynolds number Re = 500. The mean drag and normal forces and the bending moment for the prism were quantitatively compared in terms of skin-friction and pressure contributions, and related to the near-wall flow. Both three-dimensional visualizations and planar projections of the time-averaged flow field were used to identify, qualitatively, the main structures of the wake, including the horseshoe vortex, corner vortices and regions of high streamwise vorticity in the upper part of the wake. These features showed the same qualitative behavior as reported in high Reynolds number studies. It was found that some regions of high streamwise vorticity magnitude, like the tip vortices, are associated with the three-dimensional bending of the flow, and the tip vortices did not continuously extend to the free end of the prism. The three-dimensional flow analysis, which integrated different observations of the flow field around surface-mounted finite-height square prisms, also revealed that the mean near-wake structure is composed of two sections of different origin and location of dominance.  相似文献   

9.
The near-wake behind a circular cylinder undergoing rotational oscillatory motion with a relatively high forcing frequency has been investigated experimentally. Experiments were carried out varying the ratio of the forcing frequency ff to the natural vortex shedding frequency fn in the range of 0.0 (stationary) to 1.6 at an oscillation amplitude of θA=30° and Reynolds number of Re=4.14×103. Depending on the frequency ratio (FR=ff /fn), the near-wake flow could be divided into three regimes—non-lock-on (FR=0.4), transition (FR=0.8, 1.6) and lock-on (FR=1.0) regimes—with markedly different flow structures. When the frequency ratio was less than 1.0 (FR⩽1.0), the rotational oscillatory motion of the cylinder decreased the length of the vortex formation region and enhanced the mutual interaction between large-scale vortices across the wake centerline. The entrainment of ambient fluid seemed to play an important role in controlling the near-wake flow and shear-layer instability. In addition, strong vortex motion was observed throughout the near-wake region. The flow characteristics changed markedly beyond the lock-on flow regime (FR=1.0) due to the high frequency forcing. At FR=1.6, the high frequency forcing decreased the size of the large-scale vortices by suppressing the lateral extent of the wake. In addition, the interactions between the vortices shed from both sides of the cylinder were not so strong at this forcing frequency. As a consequence, the flow entrainment and momentum transfer into the wake center region were reduced. The turbulent kinetic energy was large in the region near the edge of the recirculation region, where the vortices shed from both sides of the cylinder cross the wake centerline for all frequency ratios except for the case of FR=1.6. The temporally resolved quantitative flow information extracted in the present work is useful for understanding the effects of open-loop active flow control on the near-wake flow structure.  相似文献   

10.
In order to find the intrinsic physical mechanism of the original Kármán vortex wavily distorted across the span due to the introduction of three-dimensional (3-D) geometric disturbances, a flow past a peak-perforated conic shroud is numerically simulated at a Reynolds number of 100. Based on previous work by Meiburg and Lasheras (1988), the streamwise and vertical interactions with spanwise vortices are introduced and analyzed. Then vortex-shedding patterns in the near wake for different flow regimes are reinspected and illustrated from the view of these two interactions. Generally, in regime I, spanwise vortices are a little distorted due to the weak interaction. Then in regime II, spanwise vortices, even though curved obviously, are still shed synchronously with moderate streamwise and vertical interactions. But in regime III, violently wavy spanwise vortices in some vortex-shedding patterns, typically an \(\Omega \)-type vortex, are mainly attributed to the strong vertical interactions, while other cases, such as multiple vortex-shedding patterns in sub-regime III-D, are resulted from complex streamwise and vertical interactions. A special phenomenon, spacial distribution of streamwise and vertical components of vorticity with specific signs in the near wake, is analyzed based on two models of streamwise and vertical vortices in explaining physical reasons of top and bottom shear layers wavily varied across the span. Then these two models and above two interactions are unified. Finally two sign laws are summarized: the first sign law for streamwise and vertical components of vorticity is positive in the upper shear layer, but negative in the lower shear layer, while the second sign law for three vorticity components is always negative in the wake.  相似文献   

11.
The study herein focuses on the vortex shedding characteristics and near-wake vorticity patterns of a square cylinder having self-issuing jets through holes along its span. Three different values of spacing between the consecutive holes λ with respect to the cylinder diameter D, i.e., λ/D = 1.5, 3 and 4 are studied experimentally via Digital Particle Image Velocimetry for the Reynolds number range extending from 200 to 1,000. It has been observed that the three-dimensionality of the wake flow depends on the spacing between the holes and Re number. For sufficiently low Reynolds numbers, the jet flows issuing from the holes yield a non-uniform distribution of mean flow characteristics like the shedding frequency and the formation length of vortices along the span of the cylinder when the spacing between jets along centerline is close to wavelength of the naturally existing three-dimensional wake instability. Additionally, for Re number up to 500, the self-issuing jets emanating from the holes show an indirect interaction with shear layers originating from upper and lower separation lines of the cylinder. However, for higher Re numbers of 750 and 1,000, they directly interact with and modify the vortices forming from the cylinder.  相似文献   

12.
This paper investigates flow past a rotating circular cylinder at 3600?Re?5000 and α?2.5. The flow parameter α is the circumferential speed at the cylinder surface normalized by the free-stream velocity of the uniform cross-flow. With particle image velocimetry (PIV), vortex shedding from the cylinder is clearly observed at α<1.9. The vortex pattern is very similar to the vortex street behind a stationary circular cylinder; but with increasing cylinder rotation speed, the wake is observed to become increasing narrower and deflected sideways. Properties of large-scale vortices developed from the shear layers and shed into the wake are investigated with the vorticity field derived from the PIV data. The vortex formation length is found to decrease with increasing α. This leads to a slow increase in vortex shedding frequency with α. At α=0.65, vortex shedding is found to synchronize with cylinder rotation, with one vortex being shed every rotation cycle of the cylinder. Vortex dynamics are studied at this value of α with the phase-locked eduction technique. It is found that although the shear layers at two different sides of the cylinder possess unequal vorticity levels, alternating vortices subsequently shed from the cylinder to join the two trains of vortices in the vortex street pattern exhibit very little difference in vortex strength.  相似文献   

13.
This work aims to understand the difference in the correlations between the fluctuating temperature and the vorticity from that between the fluctuating temperature and the velocity in a turbulent cylinder near wake. Measurements are made at x/d = 10, 20 and 40, where x is the streamwise distance from the cylinder axis and d is the cylinder diameter, with a Reynolds number of 2.5×103 based on d and the free-stream velocity. The three components of the fluctuating velocity vector ui(i = 1, 2 and 3), vorticity vector ωi (i = 1, 2 and 3), and temperature θ in the plane of the mean shear are measured simultaneously with a multi-wire probe consisting of four X-hotwires and four cold wires. It is found that at x/d = 10, both correlations between uiand θ and between ωi and θ predominantly take place at St = 0.21, due to the concentric distribution of the Kármán vortices and the heat. With increasing x/d, the correlation between ωi (i = 1, 2 and 3) and θ drops rapidly, as a result of the weakened Kármán vortices; in contrast, the correlation between u1 and θ increases appreciably, largely due to an enhanced correlation between u1 and θ at low frequencies or scales of motions larger than the Kármán vortex. The slowly decreasing (along x) two-point autocorrelations of u1 and θ suggest that the very-large-scale motions (VLSMs) found in wall flows occur also in the turbulent wake and are responsible for the high correlation between u1 and θ at low frequencies.  相似文献   

14.
The hairpin packet's structure and its statistical scale in the later stage of bypass transition induced by a cylinder wake are investigated by time-resolved particle image velocimetry from the side and top view, respectively. Linear stochastic estimation is used to achieve the conditionally averaged velocity fields. For the side view case, the conditionally averaged structure consists of a series of swirling motions located along a line inclining at a large angle (18°) from the wall and a low-speed region occupied by the cylinder wake appearing right above them. In the (x, z)-plane at the wall-normal height y/???=?0.2, the dominant structures are shown to be the large-scale regions of low momentum elongated almost over 3?? along the streamwise. The low-speed regions are consistently bordered by small-scale counter-rotating vortice pairs organized in the streamwise with a statistical spanwise width of 0.55??. The results suggest that in the later part of the transitional zone, the upward induction of the cylinder wake enhances both the wall-normal and spanwise extent of the hairpin packets.  相似文献   

15.
The interaction between the wake of a transverse circular cylinder and the underlying flat-plate boundary layer with a moderate gap ratio G/D=1.0 is investigated using both hydrogen-bubble-based and PIV-based visualization techniques. The spanwise rollers in the cylinder wake are found to be capable of inducing secondary vortices in the near-wall region. The mutual induction from the counter-clockwise rollers, which are closer to the wall, plays a primary role, so that these secondary vortices present linear lift-up motion at first. Their subsequent evolution dominantly determines the characteristics of the wake/boundary-layer interaction. Two different vortex interaction scenarios are observed: the secondary vortices can be either entrained into the rollers or pushed down towards the wall. This leads to a rapid three-dimensional destabilization process, through which streamwise vortices are generated. And it is suggested that these streamwise vortices are the dominant structures to promote the following boundary layer transition.  相似文献   

16.
The wake of a surface-mounted finite-height circular cylinder and the associated vortex patterns are strongly dependent on the cylinder aspect ratio and the thickness of the boundary layer on the ground plane relative to the dimensions of the cylinder. Above a critical aspect ratio, the mean wake is characterized by streamwise tip vortex structures and Kármán vortex shedding from the sides of the cylinder. Below a critical aspect ratio, a unique mean wake structure is observed. Recent experimental studies in the literature that used phase-averaged techniques, as well as recent numerical simulations, have led to an improved physical understanding of the near-wake vortex flow patterns. However, the flow above the free end of the finite circular cylinder, and its relationship to the near wake, has not been systematically studied. The effects of aspect ratio and boundary layer thickness on the free-end flow field are also not completely understood, nor has the influence of Reynolds number on the free-end flow field been fully explored. Common features associated with the free end include separation from the leading edge, a mean recirculation zone containing a prominent cross-stream arch (or mushroom) vortex, and reattachment onto the free-surface. Other flow features that remain to be clarified include a separation bubble near the leading edge, one or two cross-stream vortices within this separation bubble, the origins of the streamwise tip or trailing vortices, and various critical points in the near-surface flow topology. This paper reviews the current understanding of the flow above the free end of a surface-mounted finite-height circular cylinder, with a focus on models of the flow field, surface oil flow visualization studies, pressure and heat flux distributions on the free-end surface, measurements of the local velocity field, and numerical simulations, found in the literature.  相似文献   

17.
A comparative study of the wakes behind cylinders with grooved and smooth surfaces was performed with a view to understand the wake characteristics associated with the adult Saguaro cacti. A low-speed recirculation water channel was established for the experiment; the Reynolds number, based on the free-stream velocity and cylinder diameter (D), was kept at ReD=1500. State-of-the-art time-resolved particle image velocimetry (TR-PIV) was employed to measure a total of 20 480 realizations of the wake field at a frame rate of 250 Hz, enabling a comprehensive view of the time- and phase-averaged wake pattern. In comparison to the wake behind the smooth cylinder, the length of the recirculation zone behind the grooved cylinder was extended by nearly 18.2%, yet the longitudinal velocity fluctuation intensity was considerably weakened. A global view of the peaked spectrum of the longitudinal velocity component revealed that the intermediate region for the grooved cylinder, which approximately corresponds to the transition region where the shear layer vortices interact, merge and shed before the formation of the Karman-like vortex street, was much wider than that for the smooth one. The unsteady events near St=0.3-0.4 were detected in the intermediate region behind the grooved cylinder, but no such events were found in the smooth cylinder system. Although the formation of the Karman-like vortex street was delayed by about 0.6D downstream for the grooved cylinder, no prominent difference in the vortex street region was found in the far wake for both cylinders. The Proper Orthogonal Decomposition (POD) method was used extensively to decompose the vector and swirling strength fields, which gave a close-up view of the vortices in the near wake. The first two POD modes of the swirling strength clarified the spatio-temporal characteristics of the shear layer vortices behind the grooved cylinder. The small-scale vortices superimposed on the shear layers behind the grooved cylinder were found to be generated and convected downstream in the same phase, which would significantly reduce the fluctuating force on the cylinder surface.  相似文献   

18.
The streamwise evolution of an inclined circular cylinder wake was investigated by measuring all three velocity and vorticity components using an eight-hotwire vorticity probe in a wind tunnel at a Reynolds number Red of 7,200 based on free stream velocity (U ) and cylinder diameter (d). The measurements were conducted at four different inclination angles (α), namely 0°, 15°, 30°, and 45° and at three downstream locations, i.e., x/d = 10, 20, and 40 from the cylinder. At x/d = 10, the effects of α on the three coherent vorticity components are negligibly small for α ≤ 15°. When α increases further to 45°, the maximum of coherent spanwise vorticity reduces by about 50%, while that of the streamwise vorticity increases by about 70%. Similar results are found at x/d = 20, indicating the impaired spanwise vortices and the enhancement of the three-dimensionality of the wake with increasing α. The streamwise decay rate of the coherent spanwise vorticity is smaller for a larger α. This is because the streamwise spacing between the spanwise vortices is bigger for a larger α, resulting in a weak interaction between the vortices and hence slower decaying rate in the streamwise direction. For all tested α, the coherent contribution to [`(v2)] \overline{{v^{2}}} is remarkable at x/d = 10 and 20 and significantly larger than that to [`(u2)] \overline{{u^{2}}} and [`(w2)]. \overline{{w^{2}}}. This contribution to all three Reynolds normal stresses becomes negligibly small at x/d = 40. The coherent contribution to [`(u2)] \overline{{u^{2}}} and [`(v2)] \overline{{v^{2}}} decays slower as moving downstream for a larger α, consistent with the slow decay of the coherent spanwise vorticity for a larger α.  相似文献   

19.
Three-dimensional (3D) proper orthogonal decomposition (POD) analyses are conducted to investigate the near wake of sinusoidal wavy cylinders. For a wave amplitude a/Dm = 0.152, three typical spanwise wavelengths (λz) of the wavy cylinder are taken into account, i.e., λz/Dm = 1.89, 3.79 and 6.06, where Dm is the mean diameter of the wavy cylinder, among which λz/Dm = 1.89 and 6.06 are the optimum wavelengths corresponding to the largest reduction/suppression of fluid forces acting on the wavy cylinder. Time- and space-resolved three-component velocities of the near wake flow, obtained from large eddy simulation (LES) at a subcritical Reynolds number Re = 3 × 103, are used in the 3D POD analyses. Comparison is made among the wavy cylinders of the three λz/Dm values as well as between them and a smooth cylinder, in terms of POD modes, mode energy, mode coefficients, as well as reconstructed flow structures by lower modes. For the optimum λz/Dm = 1.89 and 6.06, energy associated with the first two POD modes is significantly reduced compared with that for λz/Dm = 3.79 and the smooth cylinder. Distinct characteristics are observed on the lower POD modes for the wavy cylinders. It is found that the first two POD modes for λz/Dm = 1.89 and 6.06 are linked to large-scale streamwise vortices that are additionally introduced into the near wake due to the wavy geometry. Meanwhile, POD mode 3 suggests that the wavy cylinder with the larger optimum λz/Dm (= 6.06) generates dominant hairpin-like and spanwise coherent structures (CSs) shedding from the saddle at a different frequency from those shedding from the node. Evolutionary development of these CSs is discussed based on reconstructed flows.  相似文献   

20.
The wake of a finite circular cylinder of small aspect ratio was studied with a seven-hole probe and thermal anemometry. The cylinder was mounted normal to a ground plane and was partially immersed in a turbulent boundary layer. The time-averaged velocity and streamwise vorticity fields showed the development of the tip vortex structures, the extent of the near-wake recirculation zone, the downwash phenomenon and base vortex structures within the boundary layer. The wake structure and power spectra were similar for cylinder aspect ratios of 5 to 9, but a distinctly different behaviour was observed for an aspect ratio of 3.  相似文献   

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

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