首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
The study of heat transfer in turbulent flow over a flat plate is very important, not only because this situation frequently arises in practice, but also in that data for an isothermal flat plate are used to calculate heat transfer in more complex cases. In particular, such data are necessary when one uses the limiting relative laws which allow calculation of the effect of compressibility, pressure gradient, blowing, and other perturbing factors [1]. Most papers dealing with heat transfer for an isothermal flat plate refer to comparatively low Re values, when the velocity distribution in the boundary layer over almost its entire thickness can be described by the universal law of the wall. However, as Re increases there is an increasing layer adjacent to the outer boundary in which the velocity distribution cannot be described by the law of the wall, and therefore the results obtained for low Re are inapplicable. In the present paper coefficients of heat transfer from a turbulent flow to an isothermal flat plate have been obtained by numerical integration of the thermal boundary-layer equations over a wide range of the parameters 3 · 105 Re 2.5·1012, 102 Pr 103.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 94–100, July–August, 1976.  相似文献   

2.
The drag of flat plates with spherical recesses was measured using the direct balance method. The experiments were run in a low-turbulence wind tunnel with a cross-section measuring 1000 × 1000 mm and a length of 400 mm. Three surfaces with recesses 7.0, 3.9, and 1.3 mm in diameter and 0.5, 0.3, and 0.2 mm deep, respectively, were tested. It is shown that on the Reynolds and Mach number ranges Re= (3–9)· 106 and M 0.3 the spherical recesses add to the drag of a flat plate in turbulent flow. The recesses have almost no effect on the location of laminar-turbulent transition, which occurs at Ret 3· 106.  相似文献   

3.
The results of balance aerodynamic tests on model straight wings with smooth and ribbed surfaces at an angle of attack =–4°–12°, Mach number M=0.15–0.63, and Reynolds number Re=2.4·106–3.5·106 are discussed. The nondimensional riblet spacings +, which determines the effect of the riblets on the turbulent friction drag, and the effect of riblets on the upper and/or lower surface of a straight wing on its drag, lift, and moment characteristics are estimated.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 33–38, March–April, 1995.  相似文献   

4.
In marine geophysical seismological prospecting extensive use is made of towed receiving systems consisting of extended flexible cylinders containing acoustic sensors over which the water flows in the longitudinal direction. The boundary layer pressure fluctuations on the cylinder surface are picked up by the sensors as hydrodynamic noise. This paper is concerned with the study of the energy and spacetime characteristics of the pressure fluctuations in the turbulent boundary layer on an extended flexible cylinder in a longitudinal flow. The pressure fluctuations on the cylinder surface have been investigated experimentally for ReX=(2–4)·107, a dimensionless diameter of the pressure fluctuation sensors d p + =dpu*/=500, where dp is the sensor diameter, u* the dynamic viscosity, and the viscosity coefficient, and frequencies 0.02311.259 (=*/U). The spectral and correlation characteristics of the pressure fluctuations on the surface of the flexible cylinder are found to differ from the corresponding characteristics for a rigid cylinder [1–4].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i aza, No, 5, pp. 49–54, September–October, 1989.  相似文献   

5.
The possibility of simplifying the formulas obtained by the Maxwell-Loyalka method for the velocity u, temperature T and diffusion d slip coefficients and the temperature jump coefficient in a binary gas mixture with frozen internal degrees of freedom of the molecules is considered. Special attention is paid to gases not having sharply different physicochemical properties. The formulas are written in a form convenient for use without linearization in the thermal diffusion coefficient. They are systematically analyzed for mixtures of inert gases, N2, O2, CO2, and H2 at temperatures extending from room temperature to 2500°K. It is shown that for the molecular weight ratios m* = m2/m1 considered the expressions for u and can be radically simplified. With an error acceptable for practical purposes (up to 10%) it is possible to employ expressions of the same structural form as for a single-component gas: for u if 1 m* 6, and for if 1 m* 3. When 1 m* 2 the expression for T can be simplified with a maximum error of 5%. Within the limits of accuracy of the method the expression for t can be linearized in the thermal diffusion coefficient. An approximate expression convenient for practical calculations is proposed for d Finally, the , u, and T for a single-component polyatomic gas with easy excitation of the internal degrees of freedom of the molecules are similarly analyzed; it is shown that these expressions can be considerably simplified.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 152–159, November–December, 1990.  相似文献   

6.
The purpose of this study is the construction of interpolation formulas for the dependence of Maxwell viscosity, a quantity which is the reciprocal of shear-strain relaxation time , on shear-strain intensity and temperature for several metals: iron, aluminum, copper, and lead. This function was interpolated in various temperature and deformation velocity ranges in accordance with available experimental data for iron (0 107 sec–1, 200 ° T 1500 °); aluminum (0 107 sec–1, 300 ° T 900 °); copper (0 105 sec–1, 300 ° T 1300 °); lead (0 106 sec–1, 90 ° T 400 °); temperatures in °K.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 114–118, July–August, 1974.  相似文献   

7.
Two new schemes of incompressible inviscid steady-state symmetric (circulationless) cavitation flow are considered. In these schemes the cavity is closed on a flow domain bounded by a closed streamline encircling either a dipole or twin point vortices with circulations equal in magnitude and opposite in sign. A comparative analysis of the efficiency of applying these schemes is performed for the case of transverse flow past a flat plate. On the cavitation number range 0.01 0.15 the drag coefficient and the relative cavity diameter and length obtained in this study are compared with their values calculated from the Riabouchinsky and Efros schemes. For a particular value of the cavitation number, the shapes of the cavities are also compared. The solutions obtained possess the properties of univalence and closure and satisfy the Brillouin condition.  相似文献   

8.
A new method for describing the rheological properties of reactive polymer melts, which was presented in an earlier paper, is developed in more detail. In particular, a detailed derivation of the equation of a first-order rheometrical flow surface is given and a procedure for determining parameters and functions occurring in this equation is proposed. The experimental verification of the presented approach was carried out using our data for polyamide-6.Notation E Dimensionless reduced viscosity, eq. (34) - E 0 Newtonian asymptote of the function (36) - E power-law asymptote of the function (36) - E = 1 the value ofE at = 1 - k degradation reaction rate constant, s–1 - k 1 rate constant of function (t), eq. (26), s–1 - k 2 rate constant of function (t), eq. (29), s–1 - K(t) residence-time-dependent consistency factor, eq. (22) - M w weight-average molecular weight - M x x-th moment of the molecular weight distribution - R gas constant - S x M x /M w - t residence time in molten state, s - t j thej-th value oft, s - T temperature, K - % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xd9vqpe0x% c9q8qqaqFn0dXdir-xcvk9pIe9q8qqaq-xir-f0-yqaqVeLsFr0-vr% 0-vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaaieGaceWFZo% Gbaiaaaaa!3B4E!\[\dot \gamma \] shear rate, s–1 - i thei-th value of , s–1 - r =1 the value of at = 1, s–1 - * reduced shear rate, eq. (44), s–1 - dimensionless reduced shear rate, eq. (35) - viscosity, Pa · s - shear-rate and residence-time dependent viscosity, Pa · s - zero-shear-rate degradation curve - degradation curve at - t0 (t) zero-residence-time flow curve - Newtonian asymptote of the RFS - instantaneous flow curve - power-law asymptote of the RFS - 0,0 zero-shear-rate and zero-residence-time viscosity, Pa · s - E=1 value of viscosity atE=1, Pa · s - * reduced viscosity, eq. (43), Pa · s - zero-residence-time rheological time constant, s - density, kg/m3 - (t),(t) residence time functions  相似文献   

9.
Various aspects of the problem of intense blowing through the surface of bodies have, been theoretically studied by a number of authors, within the framework of inviscid flow theory. A detailed bibliography on this topic is given, e.g., in [1, 2]. The well-known approaches to solution of this problem have a limited area of application. For example, asymptotic methods can be used for hypersonic flow regimes only at relatively low levels of the blown gas momentum ( = 2 = ovo 2/ V 2 1). The same limitation applies to the numerical method of straight lines [2]. The forward Eulerian calculation schemes [3, 4] smear the contact discontinuity severely, and cannot handle the case where the blown gas and the gas in the incident flow have different thermodynamic properties (o ). This paper presents results of a numerical investigation of supersonic flow over two-dimensional and axisymmetric bodies with intense blowing on the forward surface, performed using a time-dependent finite-difference method [5] with an explicit definition of the contact interface between the two cases. The calculations encompass a family of elliptic cylinders with semiaxis ratio 0.5 4, a flat-face cylinder, and a flat plate with rounding near the midsection, with variations in the blowing law, the incident flow Mach number M (3 M 10), the adiabatic indices, and the blowing parameter 0 0.5.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 117–124, January–February, 1977.In conclusion, the authors thank T. S. Novikov and I. D. Sandomirskii, who took part In the present calculations.  相似文献   

10.
Two jet methods for saturating the fluid boundary layer with microbubbles for drag reduction in contrust with gas injection through porous materials are considered. The first method is the gas injection through the slot under a special fluid wall jet. The second method is the saturation of boundary layer by microbubbles via the gas-water mixture injection through the slot. Experimental data, reflecting the skin friction drag reduction on the flat plate and total drag reduction of axisymmetric bodies, are presented. The comparison between a jet methods of gas injection and gas injection through porous materials is made.Nomenclature v free-stream velocity - v j mean velocity of a water through slot - v g mean velocity of a gas through slot - h width of slot for realizing water jet - h 1 width of slot for gas injection - incidence angle - Q volume airflow rate - C Q airflow rate coefficient (v g/v ) - C f skin friction coefficient - v j/v - C f0 C f ifQ=0 andv j=0 - f C f/C f 0 - d diameter of an axisymmetric body - L length of body - C Q 4 · ·Q/d 2 v - C D 4 ·D/1/2v 2 ·d 2 - C Q 4 ·Q/d 2 v - Q j volume flow rate of water jet - C 8 ·Q jvj/d 2 v 2 - 1 fluid density of main flow - 2 fluid density of wall jet - B 1 main stream total pressure - B 2 wall jet total pressure - v 1 main stream velocity - Be (B 2B 1)/1/21 v 1 2 = Bernoulli number - 2 v 2/1 v 1 - p st static pressure - p at atmospheric pressure - p st/p at - D hydrodynamic drag of body  相似文献   

11.
In water flows with velocities of up to 9 m/s the friction drag of a body of revolution in axial flow was investigated for dependence on the body surface structure. This was done for different types of riblet film fixed on the surface with the riblet direction aligned with the flow. The lateral spacing between the triangular shaped riblets varied between 0.033 mm and 0.152 mm. In all cases the riblet spacing was equal to the riblet height. For comparison a smooth reference film was used.Depending on the Reynolds number and the non-dimensional riblet spacings +, a turbulent drag reduction of up to 9% could be achieved with riblets in comparison with the flow over a smooth surface.In the region of transition to turbulent flow and with non-dimensional riblet spacings ofs +10–15 drag reductions of up to 13% were obtained. It is therefore conjectured, that in addition to hampering the near wall momentum exchange, the riblets can delay the development of initial turbulent structures in time and space.  相似文献   

12.
Zusammenfassung An längsangeströmten dicken Platten löst sich die Strömung im Bereich des Anströmprofils ab und legt sich nach einer gewissen Entfernung wieder an die Platte an. Dies führt gegenüber der dünnen Platte zu signifikanten Veränderungen von Umströmung, Wärmeübergang und Druckverlust.Die Aussagen zur Umströmung werden aus örtlichen Stoffübergangsmessungen (Sc=0,616) mit Hilfe einer remissionsfotometrischen Stoffübergangsmeßmethode auf der Basis von Absorption, chemischer Reaktion und gekoppelter Farbreaktion gewonnen. Als wesentliche Einflußgrößen auf Umströmung und zugeordneten Stoffübergang sind das Anströmprofil, die Anström-Reynolds-Zahl ResB sowie insbesondere die Turbulenzintensität Tu zu nennen. Die Untersuchungen erstrecken sich auf Platten mit Dicken 0,8 mms91, 3 mm mit stumpfem, halbrundem und keilförmigen Anströmprofilen sowie auf längsüberströmte Kreiszylinder mit stumpfem, halbkugelförmigem und kegelförmigen Anströmprofilen. Die mit der Plattendicke s gebildete ResB -Zahl wurde in weiten Grenzen 102 < ResB < 2 · 105, die Turbulenzintensität zwischen 0,8 Tu 6% variiert.Die Ergebnisse zeigen, daß sich die Umströmung dicker Platten generell in drei Hauptströmungsformen Plattengrenzschicht, Ablöseblase und Querwirbelablösung untergliedern läßt. Der Übergang von einer Strömungsform zur anderen wird durch kritische ResB -Zahlen erfaßt und die Abhängigkeit vom Anströmprofil über einen Profilfaktor beschrieben. Bei Ablöseblasen konnten erstmals an Platten Längswirbel nachgewiesen werden mit ihren charakteristischen Auswirkungen auf Umströmung und Stoffübergang.Dicke Platten mit Strömungsablösung lassen sich in drei Abschnitte unterteilen: 1. Profilbereich, 2 Bereich abgelöster Strömung, 3. Plattengrenzschicht ab Wiederanlegen der Strömung. Im Hinblick auf eine Berechnung des örtlichen Stoffübergangs wird die Lage des Stoffübergangsmaximums beim Wiederanlegen der Strömung — als Trennlinie der beiden letzten Plattenabschnitte — unter Einbeziehung der wesentlichen Parameter Anströmprofil, ResB-Zahl und Turbulenzintensität erfaßt und als Berechnungsgleichung angegeben.
The influence of nose section and turbulence intensity on the flow around thick plates in parallel flow
The flow parallel to thick plates separates in the nose section and reattaches on the flat plate after a certain distance. Compared with the thin flat plate the flow separation with reattachement causes significant changes of flow, heat transfer and pressure drop.The informations about the flow are recieved from measurements of local mass transfer coefficients (Sc=0,616), which are fotometrically determined by light-remission. The fundamentals of the measuring technique are based on absorbed ammonia and subsequent chemical reaction with immediate color reaction. The main parameters of flow and related mass transfer are the nose of the plate, the Reynolds-number ResB and especially the turbulence intensity Tu. The measurements were accomplished with plates of a thickness 0.8 mm s 91.3 mm of truncated, hemicylindrical and wedge-shaped noses and with cylinders in parallel flow of truncated, hemispherical and conical noses. The ResB -number with the plate thickness as characteristic length was varied in a wide range between 102 < ResB < 2 · 105, the turbulence intensity Tu between 0,8% Tu 6%.The measured results indicate, that the flow around thick plates may generally be subdivided into three main forms of the flow: boundary layer of flat plate, separation bubble and vortex shedding with reattachement. The critical Reynolds-numbers for the transition of one form to another were determined as a function of the different shapes of noses. For the first time longitudinal vortices can be observed for separation bubbles with significant influence on the flow and on the mass transfer.Thick plates with flow separation and reattachement may be subdivided into three sections: 1. the nose section, 2. the section of separated flow, 3. the section of boundary layer of flat plate downstream of reattachement. For the computation of local mass transfer rates, the position of the maximum mass transfer at the point of reattachement — the dividing line of the 2nd and 3rd section — are determined as a function of the main parameters plate nose, ResB -number and turbulence intensity.

Bezeichnungen d Zylinderdurchmesser - DA Diffusionskoeffizient von Ammoniak in Luft - KP Koeffizient in Gl. (7) - KTu Koeffizient in Gl. (6) - RedB=UB · d/ Reynolds-Zahl - ResB=UB · s/ Reynolds-@#@ Zahl - Reskrit1=UB · s/ Reynolds für den Übergang von Plattengrenzschicht zu Strömungsablösung - Reskrit2=UB · s/ Reynolds-Zahl für den Übergang von Ablöseblase zu Querwirbelablösung - Rexkrit=UB · xkrit/ Reynolds-Zahl für den Übergang von laminarer zu turbulenter Plattengrenzschicht - s Plattendicke einschließlich Trägerfolie - Shs=A · s/DA Sherwood-Zahl - Tu Turbulenzintensität in % - UB Hauptstromgeschwindigkeit im verengten Kanalquerschnitt - U Hauptstromgeschwindigkeit im freien Kanal querschnitt - mittlere turbulente Geschwindigkeitsschwankung in x-Richtung - x Koordinate in Strömungsrichtung, tangential zur Oberfläche - XP Profillänge des Anströmprofils - x Koordinate in x-Richtung ab Profilende (Plattenbeginn) - xmax' Entfernung von Plattenbeginn bis zum Wiederanlegen der Strömung nach Strömungsablösung (Stoffübergangs-maximum) - xmin Entfernung von Plattenbeginn bis zur Stelle minimalen Stoffübergangs - x* Koordinate in x-Richtung ab Wiederanlegen der Strömung nach Strömungsablösung - A Stoffübergangskoeffizient (Ammoniak) - Am mittlerer Stoffübergangskoeffizient - kinematische Zähigkeit Unter dem Titel Einfluß des Anströmprofils auf die Umströmung von Platten endlicher Dicke gekürzt vorgetragen auf der Sitzung des GVC-Fachausschusses Wärme- und Stoffübertragung am 5./6. April 1976 in Schliersee.Herrn Prof. Dr. phil. Dr.-Ing. E.h. Peter Grassmann zum 70. Geburtstag am 13. August 1977.  相似文献   

13.
The supersonic perfect-gas flow past a circular cylinder is studied on the basis of a numerical analysis of the time-dependent two-dimensional Reynolds equations using a differential q– turbulence model with reference to the experimental conditions. The calculations are carried out at Reynolds and Mach numbers Re=2× 105 and M=1.1, 1.3, and 1.7 and the experimental investigations at Re=1.62×105–2×105 and Mach numbers on the interval 0.7 M 1.7. The calculated and experimental data on the pressure coefficient distribution over the cylinder surface, the location of the separation point on the surface, and the pressure drag coefficient are compared.  相似文献   

14.
Internal waves generated by the turbulent wake of a sphere   总被引:1,自引:0,他引:1  
Internal waves generated by the turbulent wake of a sphere travelling horizontally through a linearly stratified fluid were studied using shadowgraph and particle-streak photography. The Reynolds and internal Froude number ranges considered were 2,000 Re 12,900 and 2.0 Fi 28.0, respectively. Two quite distinct flow regimes based on the structure of the turbulent wake were identified. In one, the wake is characterized by large-scale coherent structures. In the other, the wake, as viewed on a side-view shadowgraph, grows in a roughly symmetric fashion to a maximum height and then collapses slowly; such flows are termed the smallscale structures regime.Wave lengths and maximum wave heights of the internal waves were measured as functions of Nt and Fi, where N is the Brunt-Väisälä frequency and t the time. It was found that the wave lengths scale well with the streamwise dimension of the spiralling coherent structures. The maximum amplitude of the internal waves were found to scale with the vertical dimension of the turbulent wake, upon varying the internal Froude number.  相似文献   

15.
The problem of the three-dimensional incompressible turbulent boundary layer developing ahead of a circular cylinder mounted at right angles on a flat plate is considered. The direction of the uniform approach stream is normal to the leading edge of the plate. The turbulence is simulated by means of five different isotropic algebraic models of eddy viscosity. The boundary layer equations are solved numerically by means of a second-order-accurate implicit finite-difference method. The principal characteristics of the flow obtained on the basis of the turbulence models selected are compared for a free-stream Reynolds number Re = 107.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 36–43, November–December, 1991.  相似文献   

16.
The natural convection of fresh water in a square cell is considered at a temperature close to the density inversion temperature for Grashof numbers 2.9 · 104 Gr 106. As a result of the numerical investigation, one steady-state and three self-oscillating regimes are found in addition to the three steady-state flows previously detected earlier and described for low Grashof numbers ( 0 Gr 2 · 105). The basic characteristics of the unsteady flows are analyzed by means of the Fourier method, the fundamental oscillation frequencies are found, and the flow evolution and the variation of the oscillation characteristics with increase in the Gr number are considered.  相似文献   

17.
Results are presented of a study of the gasdynamic parameters and the geometric characteristics of the mixing zone of axisymmetric jets of gases of differing density (Freon-12, air, and helium) propagating in a parallel air stream, within the limits of the initial segment (0x/R3–30). Experimental data are presented on the effect of different densities (0. 27 n8.2) and velocities (0m1.7) of the gas jet and the parallel stream on the mixing process.  相似文献   

18.
Several theoretical [1–4] and experimental [5–7] studies have been devoted to the study of the effect of distributed injection of a gaseous substance on the characteristics of the turbulent boundary layer. The primary study has been made of flow past a flat plate with gas injection. The theoretical methods are based primarily on the semiempirical theories of Prandtl [1] and Karman [2].In contrast with the previous studies, the present paper proposes a power law for the mixing length; this makes it possible to obtain velocity profiles which degenerate to the known power profiles [8] in the case of flow without blowing and heat transfer. This approach yields analytic results for flows with moderate pressure gradient.Notation x, y coordinates - U, V velocity components - density - T temperature - h enthalpy - H total enthalpy - c mass concentration - , , D coefficients of molecular viscosity, thermal conductivity, diffusion - cp specific heat - adiabatic exponent - r distance from axis of symmetry to surface - boundary layer thickness - U velocity in stream core - friction - cf friction coefficient - P Prandtl number - S Schmidt number - St Stanton number - M Mach number - j=0 plane case - j=1 axisymmetric case The indices 1 injected gas - 2 mainstream gas - w quantities at the wall - core of boundary layer - 0 flow of incompressible gas without injection - v=0 flow of compressible gas without injection - * quantities at the edge of the laminar sublayer - quantities at the initial section - turbulent transport coefficients  相似文献   

19.
The distribution of pressure, velocity, and electrical potential has been investigated for a mercury flow in insulated rectangular ducts with a large side ratio (Hartmann-type flow). The ranges of variation of the Reynolds, Hartmann, and Stewart numbers were 7·102R5·105, 0H490, and 0N24, respectively. Special attention is given to the sections of the channel where the flow enters and leaves the magnetic field. In these zones the pressure is sharply nonuniform and the velocity profiles in a plane perpendicular to the field acquire an M shape. A relation is established between the length of the entrance section, where the flow is three-dimensional, and the MHD similarity criteria. It is shown that ducts which are hydraulically smooth in the absence of a magnetic field become increasingly rough as the field grows stronger. Data are obtained on the resistance coefficient for a stabilized flow measured in a magnetic field and on the dependence of the critical Reynolds number on the Hartmann number.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 10–21, July–August, 1971.  相似文献   

20.
Steady flow of supersonic air over a sphere is examined, allowing for viscosity, heat conduction, and actual physical and chemical processes. Flow in the shock layer at flight speeds in the range 3 km/sec V10 km/sec (104R106) is investigated, under the assumption of local thermodynamic equilibrium. The flow is described by simplified Navier-Stokes equations, which are solved by a finite difference method. The case of a cooled surface is examined. The distribution of gasdynamic parameters is obtained in different flow regimes. The distribution of heat flux and friction coefficient is investigated as a function of the oncoming-stream parameters and the sphere radius. The shape and position of the shock wave are determined, and the stream lines and sonic lines are constructed.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 150–153, July–August, 1970.The authors thank Yu. P. Lun'kin and F. D. Popov for their help in formulating the problem and their constant interest.  相似文献   

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

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