首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
Mass conservation and linear momentum balance relations for a porous body and any fluid therein, valid at any given length scale in excess of nearest-neighbour molecular separations, are established in terms of local weighted averages of molecular quantities. The mass density field for the porous body at a given scale is used to identify its boundary at this scale, and a porosity field is defined for any pair of distinct length scales. Specific care is paid to the interpretation of the stress tensor associated with each of the body and fluid at macroscopic scales, and of the force per unit volume each exerts on the other. Consequences for the usual microscopic and macroscopic viewpoints are explored.Nomenclature material system; Section 2.1. - porous body (example of a material system); Sections 2.1, 3.1, 4.1 - fluid body (example of a material system); Sections 2.1, 3.1, 4.1 - weighting function; Sections 2.1, 2.3 - ,h weighting function corresponding to spherical averaging regions of radius and boundary mollifying layer of thicknessh; Section 3.2 - Euclidean space; Section 2.1 - V space of all displacements between pairs of points in; Section 2.1 - mass density field corresponding to; (2.3)1 - P , f mass density fields for , ; (4.1) - P momentum density field corresponding to; (2.3)2 - v velocity field corresponding to; (2.4) - S r (X) interior of sphere of radiusr with centre at pointx; (3.3) - boundary ofany region - region in which p > 0 with = ,h; (3.1) - subset of whose points lie at least+h from boundary of ; (3.4) - abbreviated versions of ; Section 3.2, Remark 4 - strict interior of ; (3.7) - analogues of for fluid system ; Section 3.2 - general version of corresponding to any choice of weighting function; (4.6) - interfacial region at scale; (3.8) - 0 scale of nearest-neighbour separations in ; Section 3.2. Remark 1 - porosity field at scales ( 1; 2); (3.9) - pore space at scales ( 1; 2); (3.12)  相似文献   

2.
Two-phase flow in stratified porous media is a problem of central importance in the study of oil recovery processes. In general, these flows are parallel to the stratifications, and it is this type of flow that we have investigated experimentally and theoretically in this study. The experiments were performed with a two-layer model of a stratified porous medium. The individual strata were composed of Aerolith-10, an artificial: sintered porous medium, and Berea sandstone, a natural porous medium reputed to be relatively homogeneous. Waterflooding experiments were performed in which the saturation field was measured by gamma-ray absorption. Data were obtained at 150 points distributed evenly over a flow domain of 0.1 × 0.6 m. The slabs of Aerolith-10 and Berea sandstone were of equal thickness, i.e. 5 centimeters thick. An intensive experimental study was carried out in order to accurately characterize the individual strata; however, this effort was hampered by both local heterogeneities and large-scale heterogeneities.The theoretical analysis of the waterflooding experiments was based on the method of large-scale averaging and the large-scale closure problem. The latter provides a precise method of discussing the crossflow phenomena, and it illustrates exactly how the crossflow influences the theoretical prediction of the large-scale permeability tensor. The theoretical analysis was restricted to the quasi-static theory of Quintard and Whitaker (1988), however, the dynamic effects described in Part I (Quintard and Whitaker 1990a) are discussed in terms of their influence on the crossflow.Roman Letters A interfacial area between the -region and the -region contained within V, m2 - a vector that maps onto , m - b vector that maps onto , m - b vector that maps onto , m - B second order tensor that maps onto , m2 - C second order tensor that maps onto , m2 - E energy of the gamma emitter, keV - f fractional flow of the -phase - g gravitational vector, m/s2 - h characteristic length of the large-scale averaging volume, m - H height of the stratified porous medium , m - i unit base vector in the x-direction - K local volume-averaged single-phase permeability, m2 - K - {K}, large-scale spatial deviation permeability - { K} large-scale volume-averaged single-phase permeability, m2 - K * large-scale single-phase permeability, m2 - K ** equivalent large-scale single-phase permeability, m2 - K local volume-averaged -phase permeability in the -region, m2 - K local volume-averaged -phase permeability in the -region, m2 - K - {K } , large-scale spatial deviation for the -phase permeability, m2 - K * large-scale permeability for the -phase, m2 - l thickness of the porous medium, m - l characteristic length for the -region, m - l characteristic length for the -region, m - L length of the experimental porous medium, m - characteristic length for large-scale averaged quantities, m - n outward unit normal vector for the -region - n outward unit normal vector for the -region - n unit normal vector pointing from the -region toward the -region (n = - n ) - N number of photons - p pressure in the -phase, N/m2 - p 0 reference pressure in the -phase, N/m2 - local volume-averaged intrinsic phase average pressure in the -phase, N/m2 - large-scale volume-averaged pressure of the -phase, N/m2 - large-scale intrinsic phase average pressure in the capillary region of the -phase, N/m2 - - , large-scale spatial deviation for the -phase pressure, N/m2 - pc , capillary pressure, N/m2 - p c capillary pressure in the -region, N/m2 - p capillary pressure in the -region, N/m2 - {p c } c large-scale capillary pressure, N/m2 - q -phase velocity at the entrance of the porous medium, m/s - q -phase velocity at the entrance of the porous medium, m/s - Swi irreducible water saturation - S /, local volume-averaged saturation for the -phase - S i initial saturation for the -phase - S r residual saturation for the -phase - S * { }*/}*, large-scale average saturation for the -phase - S saturation for the -phase in the -region - S saturation for the -phase in the -region - t time, s - v -phase velocity vector, m/s - v local volume-averaged phase average velocity for the -phase, m/s - {v } large-scale averaged velocity for the -phase, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - V large-scale averaging volume, m3 - y position vector relative to the centroid of the large-scale averaging volume, m - {y}c large-scale average of y over the capillary region, m Greek Letters local porosity - local porosity in the -region - local porosity in the -region - local volume fraction for the -phase - local volume fraction for the -phase in the -region - local volume fraction for the -phase in the -region - {}* { }*+{ }*, large-scale spatial average volume fraction - { }* large-scale spatial average volume fraction for the -phase - mass density of the -phase, kg/m3 - mass density of the -phase, kg/m3 - viscosity of the -phase, N s/m2 - viscosity of the -phase, Ns/m2 - V /V , volume fraction of the -region ( + =1) - V /V , volume fraction of the -region ( + =1) - attenuation coefficient to gamma-rays, m-1 - -   相似文献   

3.
Incoherent phase transitions are more difficult to treat than their coherent counterparts. The interface, which appears as a single surface in the deformed configuration, is represented in its undeformed state by a separate surface in each phase. This leads to a rich but detailed kinematics, one in which defects such as vacancies and dislocations are generated by the moving interface. In this paper we develop a complete theory of incoherent phase transitions in the presence of deformation and mass transport, with phase interface structured by energy and stress. The final results are a complete set of interface conditions for an evolving incoherent interface.Frequently used symbols Ai,Ci generic subsurface of St - Bi undeformed phase-i region - C configurational bulk stress, Eshelby tensor - F deformation gradient - G inverse deformation gradient - H relative deformation gradient - J bulk Jacobian of the deformation - ¯K, Ki total (twice the mean) curvature of and Si - Lin (U, V) linear transformations from U into V - Lin+ linear transformations of 3 with positive determinant - Orth+ rotations of 3 - Qa external bulk mass supply of species a - ¯S bulk Cauchy stress tensor - S bulk Piola-Kirchhoff stress tensor - Si undeformed phase i interface - Ui relative velocity of Si - Unim+ linear transformations of 3 with unit determinant - ¯V, Vi normal velocity of and Si - intrinsic edge velocity of S and A i S - Wi volume flow across the phase-i interface - X material point - b external body force - e internal bulk configurational force - fi external interfacial force (configurational) - ¯g external interfacial force (deformational) - grad, div spatial gradient and divergence - gradient and divergence on - h relative deformation - ha, diffusive mass flux of species a and list of mass fluxes - ¯m outward unit normal to a spatial control volume - ¯n, ni unit normal to and Si - n subspace of 3 orthogonal to n - ¯qa external interfacial mass supply of species a - s ......... - ¯v, vi compatible velocity fields of and Si - ¯w, wi compatible edge velocity fields for and Ai - x spatial point - yi deformation or motion of phase i - y. material velocity - generic subsurfaces of - , i deformed body and deformed phase-i region - () energy supplied to by mass transport - symmetry group of the lattice - i, surface jacobians - lattice - () power expended on - spatial control volume - S deformed phase interface - lattice point density - interfacial power density - , A total surface stress - C configurational surface stress for phase 1 (material) - ¯Ci configurational surface stress (spatial) - Fi tangential deformation gradient - Gi inverse tangential deformation gradient - H incoherency tensor - ¯1(x), 1i(X) inclusions of ¯n(x) and n i (X) into 3 - K configurational surface stress for phase 2 (material) - ¯L, li curvature tensor of and Si - ¯P(x), Pi(X) projections of 3 onto ¯n(x) and ni (X) - ¯S, S deformational surface stress (spatial and material) - ¯a, a normal part of total surface stress - c normal part of configurational surface stress for phase 1 (material) - ei internal interfacial configurational force - ¯v, vi unit normal to and A i - (x),i(X) projections of 3 onto ¯n(x) and n i (X) - i normal internal force (material) - bulk free energy - slip velocity - i=(–1)i i ......... - a, chemical potential of species a and list of potentials - a, bulk molar density of species a and list of molar densities - i normal internal force (spatial) - surface tension - , i effective shear - referential-to-spatial transform of field - interfacial energy - grand canonical potential - l unit tensor in 3 - x, vector and tensor product in 3 - (...)., t(...) material and spatial time derivative - , Div material gradient and divergence - gradient and divergence on Si - (...), (...) normal time derivative following and Si - (...) limit of a bulk field asx ,xi - [...],...> jump and average of a bulk field across the interface - (...)ext extension of a surface tensor to 3 - tangential part of a vector (tensor) on and Si  相似文献   

4.
In this paper we continue previous studies of the closure problem for two-phase flow in homogeneous porous media, and we show how the closure problem can be transformed to a pair of Stokes-like boundary-value problems in terms of pressures that have units of length and velocities that have units of length squared. These are essentially geometrical boundary value problems that are used to calculate the four permeability tensors that appear in the volume averaged Stokes' equations. To determine the geometry associated with the closure problem, one needs to solve the physical problem; however, the closure problem can be solved using the same algorithm used to solve the physical problem, thus the entire procedure can be accomplished with a single numerical code.Nomenclature a a vector that maps V onto , m-1. - A a tensor that maps V onto . - A area of the - interface contained within the macroscopic region, m2. - A area of the -phase entrances and exits contained within the macroscopic region, m2. - A area of the - interface contained within the averaging volume, m2. - A area of the -phase entrances and exits contained within the averaging volume, m2. - Bo Bond number (= (=(–)g2/). - Ca capillary number (= v/). - g gravitational acceleration, m/s2. - H mean curvature, m-1. - I unit tensor. - permeability tensor for the -phase, m2. - viscous drag tensor that maps V onto V. - * dominant permeability tensor that maps onto v , m2. - * coupling permeability tensor that maps onto v , m2. - characteristic length scale for the -phase, m. - l characteristic length scale representing both and , m. - L characteristic length scale for volume averaged quantities, m. - n unit normal vector directed from the -phase toward the -phase. - n unit normal vector representing both n and n . - n unit normal vector representing both n and n . - P pressure in the -phase, N/m2. - p superficial average pressure in the -phase, N/m2. - p intrinsic average pressure in the -phase, N/m2. - p p , spatial deviation pressure for the -phase, N/m2. - r 0 radius of the averaging volume, m. - r position vector, m. - t time, s. - v fluid velocity in the -phase, m/s. - v superficial average velocity in the -phase, m/s. - v intrinsic average velocity in the -phase, m/s. - v v , spatial deviation velocity in the -phase, m/s. - V volume of the -phase contained within the averaging volmue, m3. - averaging volume, m3. Greek Symbols V /, volume fraction of the -phase. - viscosity of the -phase, Ns/m2. - density of the -phase, kg/m3. - surface tension, N/m. - (v +v T ), viscous stress tensor for the -phase, N/m2.  相似文献   

5.
The effective permeability of a heterogeneous porous medium   总被引:5,自引:0,他引:5  
The effective (single-phase) permeability of an (infinite) heterogeneous porous medium is studied using a formalism of Green's functions. We give formal expressions for it in the form of a series expansion involving the microscopic random-permeability field many-body correlation functions of higher and higher order.The particular case of a log-normal medium of infinite extent is studied using field-theoretical methods. Using partial series resummation techniques, we derivea formula up to all orders in the local correlations which was first reckoned by many authors by means of a first-order calculation. The formula — which remains an approximation — works whatever the dimensionality of the space, and gives the following simple estimate for the effective permeability in 3 D:K eff=k 1/33. The method is general and the approximations can be systematically improved on when more complex situations are studied.Roman Letters D number of dimensions of the space in which the flow takes place - f(r) body force field,N - f(q) Fourier-transformed body-force field, Nm3 - G 0(r, r) Green's function of the Laplace operator, m–1 - g(k,r, r) velocity propagator before averaging, m–1 - G(r, r) velocity propagator after averaging, m–1 - j(r) a scalar dimensionless field - k(r) local value of the permeability at point r, m2 - K eff effective permeability - K g geometric average of the local permeability, m2 - l typical size of the averaging volume, m - L characteristic length of the porous medium or of the reservoir, m - L(r, r) projection operator, m–2 - M(r, r) scattering operator, m–3 - p(r) local value of the pressure, Nm–2 - p(k,r, r) pressure propagator before averaging, m–1 - P(r, r) pressure propagator after averaging, m–1 - r position vector, m - r modulus of vectorr, m - unit vector pointing in the direction ofr - q Fourier wave vector, m–1 - q modulus of the Fourier wave-vectorq, m–1 - unit vector pointing in the direction ofq - projector over vector - 1 unit tensor - X(r) a local random variable - ¯X(r) volume averaged local random variable - X (r) ensemble averaged local random variable - V large-scale averaging volume, m3 - Z(j) generating functional of a random field - Z(r,j) modified generating functional of a random field - Z normalization factor Greek Letters 0 average value of the logarithm of the permeability - (r) fluctuation of the logarithm of permeability at pointr - viscosity of the fluid, Nt/m2 - (r–r) two-point correlation function of the fluctuations of the logarithm of the permeability - k correlation length of the permeability correlation function, m - u correlation length of the velocity correlation function, m  相似文献   

6.
The paper reports the outcome of a numerical study of fully developed flow through a plane channel composed of ribleted surfaces adopting a two-equation turbulence model to describe turbulent mixing. Three families of riblets have been examined: idealized blade-type, V-groove and a novel U-form that, according to computations, achieves a superior performance to that of the commercial V-groove configuration. The maximum drag reduction attained for any particular geometry is broadly in accord with experiment though this optimum occurs for considerably larger riblet heights than measurements indicate. Further explorations bring out a substantial sensitivity in the level of drag reduction to the channel Reynolds number below values of 15 000 as well as to the thickness of the blade riblet. The latter is in accord with the trends of very recent, independent experimental studies.Possible shortcomings in the model of turbulence are discussed particularly with reference to the absence of any turbulence-driven secondary motions when an isotropic turbulent viscosity is adopted. For illustration, results are obtained for the case where a stress transport turbulence model is adopted above the riblet crests, an elaboration that leads to the formation of a plausible secondary motion sweeping high momentum fluid towards the wall close to the riblet and thereby raising momentum transport.Nomenclature c f Skin friction coefficient - c f Skin friction coefficient in smooth channel at the same Reynolds number - k Turbulent kinetic energy - K + k/ w - h Riblet height - S Riblet width - H Half height of channel - Re Reynolds number = volume flow/unit width/ - Modified turbulent Reynolds number - R t turbulent Reynolds numberk 2/ - P k Shear production rate ofk, t (U i /x j + U j /x i ) U i /x j - dP/dz Streamwise static pressure gradient - U i Mean velocity vector (tensor notation) - U Friction velocity, w/ where w=–H dP/dz - W Mean velocity - W b Bulk mean velocity through channel - y + yU /v. Unless otherwise stated, origin is at wall on trough plane of symmetry - Kinematic viscosity - t Turbulent kinematic viscosity - Turbulence energy dissipation rate - Modified dissipation rate – 2(k 1/2/x j )2 - Density - k , Effective turbulent Prandtl numbers for diffusion ofk and   相似文献   

7.
Summary A three-parameter model is introduced to describe the shear rate — shear stress relation for dilute aqueous solutions of polyacrylamide (Separan AP-30) or polyethylenoxide (Polyox WSR-301) in the concentration range 50 wppm – 10,000 wppm. Solutions of both polymers show for a similar rheological behaviour. This behaviour can be described by an equation having three parameters i.e. zero-shear viscosity 0, infinite-shear viscosity , and yield stress 0, each depending on the polymer concentration. A good agreement is found between the values calculated with this three-parameter model and the experimental results obtained with a cone-and-plate rheogoniometer and those determined with a capillary-tube rheometer.
Zusammenfassung Der Zusammenhang zwischen Schubspannung und Schergeschwindigkeit von strukturviskosen Flüssigkeiten wird durch ein Modell mit drei Parametern beschrieben. Mit verdünnten wäßrigen Polyacrylamid-(Separan AP-30) sowie Polyäthylenoxidlösungen (Polyox WSR-301) wird das Modell experimentell geprüft. Beide Polymerlösungen zeigen im untersuchten Schergeschwindigkeitsbereich von ein ähnliches rheologisches Verhalten. Dieses Verhalten kann mit drei konzentrationsabhängigen Größen, nämlich einer Null-Viskosität 0, einer Grenz-Viskosität und einer Fließgrenze 0 beschrieben werden. Die Ergebnisse von Experimenten mit einem Kegel-Platte-Rheogoniometer sowie einem Kapillarviskosimeter sind in guter Übereinstimmung mit den Werten, die mit dem Drei-Parameter-Modell berechnet worden sind.

a Pa–1 physical quantity defined by:a = {1 – ( / 0)}/ 0 - c l concentration (wppm) - D m capillary diameter - L m length of capillary tube - P Pa pressure drop - R m radius of capillary tube - u m s–1 average velocity - v r m s–1 local axial velocity at a distancer from the axis of the tube - shear rate (–dv r /dr) - local shear rate in capillary flow - s–1 wall shear rate in capillary flow - Pa s dynamic viscosity - a Pa s apparent viscosity defined by eq. [2] - ( a ) Pa s apparent viscosity in capillary tube at a distanceR from the axis - 0 Pa s zero-shear viscosity defined by eq. [4] - Pa s infinite-shear viscosity defined by eq. [5] - l ratior/R - kg m density - Pa shear stress - 0 Pa yield stress - r Pa local shear stress in capillary flow - R Pa wall shear stress in capillary flow R = (PR/2L) - v m3 s–1 volume rate of flow With 8 figures and 1 table  相似文献   

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.
The theory of a vibrating-rod densimeter   总被引:1,自引:0,他引:1  
The paper presents a theory of a device for the accurate determination of the density of fluids over a wide range of thermodynamic states. The instrument is based upon the measurement of the characteristics of the resonance of a circular section tube, or rod, performing steady, transverse oscillations in the fluid. The theory developed accounts for the fluid motion external to the rod as well as the mechanical motion of the rod and is valid over a defined range of conditions. A complete set of working equations and corrections is obtained for the instrument which, together with the limits of the validity of the theory, prescribe the parameters of a practical design capable of high accuracy.Nomenclature A, B, C, D constants in equation (60) - A j , B j constants in equation (18) - a j + , a j wavenumbers given by equation (19) - C f drag coefficient defined in equation (64) - C f /0 , C f /1 components of C f in series expansion in powers of - c speed of sound - D b drag force of fluid b - D 0 coefficient of internal damping - E extensional modulus - force per unit length - F j + , F j constants in equation (24) - f, g functions of defined in equations (56) - G modulus of rigidity - I second moment of area - K constant in equation (90) - k, k constants defined in equations (9) - L half-length of oscillator - Ma Mach number - m a mass per unit length of fluid a - m b added mass per unit length of fluid b - m s mass per unit length of solid - n j eigenvalue defined in equation (17) - P power (energy per cycle) - P a , P b power in fluids a and b - p pressure - R radius of rod or outer radius of tube - R c radius of container - R i inner radius of tube - r radial coordinate - T tension - T visc temperature rise due to heat generation by viscous dissipation - t time - v r , v radial and angular velocity components - y lateral displacement - z axial coordinate - dimensionless tension - a dimensionless mass of fluid a - b dimensionless added mass of fluid b - b dimensionless drag of fluid b - dimensionless parameter associated with - 0 dimensionless coefficient of internal damping - dimensionless half-width of resonance curve - dimensionless frequency difference defined in equation (87) - spatial resolution of amplitude - R, , , s , increments in R, , , s , - dimensionless amplitude of oscillation - dimensionless axial coordinate - ratio of to - a , b ratios of to for fluids a and b - angular coordinate - parameter arising from distortion of initially plane cross-sections - f thermal conductivity of fluid - dimensionless parameter associated with - viscosity of fluid - a , b viscosity of fluids a and b - dimensionless displacement - j jth component of - density of fluid - a , b density of fluids a and b - s density of tube or rod material - density of fluid calculated on assumption that * - dimensionless radial coordinate - * dimensionless radius of container - dimensionless times - rr rr, r radial normal and shear stress components - spatial component of defined in equation (13) - j jth component of - dimensionless streamfunction - 0, 1 components of in series expansion in powers of - phase angle - r phase difference - ra , rb phase difference for fluids a and b - streamfunction - j jth component defined in equation (22) - dimensionless frequency (based on ) - a , b dimensionless frequency in fluids a and b - s dimensionless frequency (based on s ) - angular frequency - 0 resonant frequency in absence of fluid and internal damping - r resonant frequency in absence of internal fluid - ra , rb resonant frequencies in fluids a and b - dimensionless frequency - dimensionless frequency when a vanishes - dimensionless frequencies when a vanishes in fluids a and b - dimensionless resonant frequency when a , b, b and 0 vanish - dimensionless resonant frequency when a , b and b vanish - dimensionless resonant frequency when b and b vanish - dimensionless frequencies at which amplitude is half that at resonance  相似文献   

10.
The mixed convection flow in a vertical duct is analysed under the assumption that , the ratio of the duct width to the length over which the wall is heated, is small. It is assumed that a fully developed Poiseuille flow has already been set up in the duct before heat from the wall causes this to be changed by the action of the buoyancy forces, as measured by a buoyancy parameter . An analytical solution is derived for the case when the Reynolds numberRe, based on the duct width, is of 0 (1). This is extended to the case whenRe is 0 (–1) by numerical integrations of the governing equations for a range of values of representing both aiding and opposing flows. The limiting cases, || 1 andR=Re of 0 (1), andR and both large, with of 0 (R 1/3) are considered further. Finally, the free convection limit, large with R of 0 (1), is discussed.
Mischkonvektion in engen senkrechten Rohren
Zusammenfassung Mischkonvektion in einem senkrechten Rohr wird unter der Voraussetzung untersucht, daß das Verhältnis der Rohrbreite zur Länge, über welche die Wand beheizt wird, klein ist. Es wird angenommen, daß sich bereits eine voll entwickelte Poiseuille-Strömung in dem Rohr eingestellt hat, bevor Antriebskräfte, gemessen mit dem Auftriebsparameter , aufgrund der Wandbeheizung die Strömung verändern. Es wird eine analytische Lösung für den Fall erhalten, daß die mit der Rohrbreite als charakteristische Länge gebildete Reynolds-ZahlRe konstant ist. Dies wird mittels einer numerischen Integration der wichtigsten Gleichungen auf den FallRe =f (–1) sowohl für Gleich- als auch für Gegenstrom ausgedehnt. Weiterhin werden die beiden Grenzfälle betrachtet, wenn || 1 undR=Re konstant ist, sowieR und beide groß mit proportionalR 1/3. Schließlich wird der Grenzfall der freien Konvektion, großes mit konstantem R, diskutiert.

Nomenclature g acceleration due to gravity - Gr Grashof number - G modified Grashof number - h duct width - l length of the heated section of the duct wall - p pressure - Pr Prandtl number - Q flow rate through the duct - Q 0 heat transfer on the wally=0 - Q 1 heat transfer on the wally=1 - Re Reynolds number - R modified Reynolds number - T temperature of the fluid - T 0 ambient temperature - T applied temperature difference - u, velocity component in thex-direction - v, velocity component in they-direction - x, co-ordinate measuring distance along the duct - y, co-ordinate measuring distance across the duct - buoyancy parameter - 0 modified buoyancy parameter, 0=R –1/3 - coefficient of thermal expansion - ratio of duct width to heated length, =h/l - (non-dimensional) temperature - w applied temperature on the wally=0 - kinematic viscosity - density of the fluid - 0 shear stress on the wally=0 - 1 shear stress on the wally=1 - stream function  相似文献   

11.
The documentation and control of flow disturbances downstream of various open inlet contractions was the primary focus with which to evaluate a spatial sampling technique. An X-wire probe was rotated about the center of a cylindrical test section at a radius equal to one-half that of the test section. This provided quasi-instantaneous multi-point measurements of the streamwise and azimuthal components of the velocity to investigate the temporal and spatial characteristics of the flowfield downstream of various contractions. The extent to which a particular contraction is effective in controlling ingested flow disturbances was investigated by artificially introducing disturbances upstream of the contractions. Spatial as well as temporal mappings of various quantities are presented for the streamwise and azimuthal components of the velocity. It was found that the control of upstream disturbances is highly dependent on the inlet contraction; for example, reduction of blade passing frequency noise in the ground testing of jet engines should be achieved with the proper choice of inlet configurations.List of symbols K uv correlation coefficient= - P percentage of time that an azimuthal fluctuating velocity derivative dv/d is found - U streamwise velocity component U=U (, t) - V azimuthal or tangential velocity component due to flow and probe rotation V=V (, t) - mean value of streamwise velocity component - U m resultant velocity from and - mean value of azimuthal velocity component induced by rotation - u fluctuating streamwise component of velocity u=u(, t) - v fluctuating azimuthal component of velocity v = v (, t) - u phase-averaged fluctuating streamwise component of velocity u=u(0) - v phase-averaged fluctuating azimuthal component of velocity v=v() - û average of phase-averaged fluctuating streamwise component of velocity (u()) over cases I-1, II-1 and III-1 û = û() - average of phase-averaged fluctuating azimuthal component of velocity (v()) over cases I-1, II-1 and III-1 - u fluctuating streamwise component of velocity corrected for non-uniformity of probe rotation and/or phase-related vibration u = u(0, t) - v fluctuating azimuthal component of velocity corrected for non-uniformity or probe rotation and/or phase-related vibration v=v (, t) - u 2 rms value of corrected fluctuating streamwise component of velocity - rms value of corrected fluctuating azimuthal component of velocity - phase or azimuthal position of X-probe  相似文献   

12.
Zusammenfassung Zur Berechnung turbulenter Strömungen wird das k--Modell im Ansatz für die turbulente Scheinzähigkeit erweitert, so daß es den Querkrümmungs- und Dichteeinfluß auf den turbulenten Transportaustausch erfaßt. Die dabei zu bestimmenden Konstanten werden derart festgelegt, daß die bestmögliche Übereinstimmung zwischen Berechnung und Messung erzielt wird. Die numerische Integration der Grenzschichtgleichungen erfolgt unter Verwendung einer Transformation mit dem Differenzenverfahren vom Hermiteschen Typ. Das erweiterte Modell wird auf rotationssymmetrische Freistrahlen veränderlicher Dichte angewendet und zeigt Übereinstimmung zwischen Rechnung und Experiment.
On the influence of transvers-curvature and density in inhomogeneous turbulent free jets
The prediction of turbulent flows based on the k- model is extended to include the influence of transverse-curvature and density on the turbulent transport mechanisms. The empirical constants involved are adjusted such that the best agreement between predictions and experimental results is obtained. Using a transformation the boundary layer equations are solved numerically by means of a finite difference method of Hermitian type. The extended model is applied to predict the axisymmetric jet with variable density. The results of the calculations are in agreement with measurements.

Bezeichnungen Wirbelabsorptionskoeffizient - ci Massenkonzentration der Komponente i - cD, cL, c, c1, c2 Konstanten des Turbulenzmodells - d Düsendurchmesser - E bezogene Dissipationsrate - f bezogene Stromfunktion - f Korrekturfunktion für die turbulente Scheinzähigkeit - j turbulenter Diffusionsstrom - k Turbulenzenergie - ki Schrittweite in -Richtung - K dimensionslose Turbulenzenergie - L turbulentes Längenmaß - Mi Molmasse der Komponente i - p Druck - allgemeine Gaskonstante - r Querkoordinate - r0,5 Halbwertsbreite der Geschwindigkeit - r0,5c Halbwertsbreite der Konzentration - T Temperatur - u Geschwindigkeitskomponente in x-Richtung - v Geschwindigkeitskomponente in r-Richtung - x Längskoordinate - y allgemeine Funktion - Yi diskreter Wert der Funktion y - Relaxationsfaktor für Iteration - turbulente Dissipationsrate - transformierte r-Koordinate - kinematische Zähigkeit - Exponent - transformierte x-Koordinate - Dichte - k, Konstanten des Turbulenzmodells - Schubspannung - allgemeine Variable - Stromfunktion - Turbulente Transportgröße Indizes 0 Strahlanfang - m auf der Achse - r mit Berücksichtigung der Krümmung - t turbulent - mit Berücksichtigung der Dichte - im Unendlichen - Schwankungswert oder Ableitung einer Funktion - – Mittelwert Herrn Professor Dr.-Ing. R. Günther zum 70. Geburtstag gewidmet  相似文献   

13.
An analysis is presented for laminar source flow between infinite parallel porous disks. The solution is in the form of a perturbation from the creeping flow solution. Expressions for the velocity, pressure, and shear stress are obtained and compared with the results based on the assumption of creeping flow.Nomenclature a half distance between disks - radial coordinate - r dimensionless radial coordinate, /a - axial coordinate - z dimensionless axial coordinate, /a - radial coordinate of a point in the flow - R dimensionless radial coordinate of a point in the flow, /a - velocity component in radial direction - u =a/, dimensionless velocity component in radial direction - velocity component in axial direction - v = a/}, dimensionless velocity component in axial direction - static pressure - p = (a 2/ 2), dimensionless static pressure - =p(r, z)–p(R, z), dimensionless pressure drop - V magnitude of suction or injection velocity - Q volumetric flow rate of the source - Re source Reynolds number, Q/4a - reduced Reynolds number, Re/r 2 - critical Reynolds number - R w wall Reynolds number, Va/ - viscosity - density - =/, kinematic viscosity - shear stress at upper disk - 0 = (a 2/ 2), dimensionless shear stress at upper disk - shear stress ratio, 0/( 0)inertialess - u = , dimensionless average radial velocity - u/u, ratio of radial velocity to average radial velocity - dimensionless stream function  相似文献   

14.
We study the modelling of purely conductive heat transfer between a porous medium and an external fluid within the framework of the volume averaging method. When the temperature field for such a system is classically determined by coupling the macroscopic heat conduction equation in the porous medium domain to the heat conduction equation in the external fluid domain, it is shown that the phase average temperature cannot be predicted without a generally negligible error due to the fact that the boundary conditions at the interface between the two media are specified at the macroscopic level.Afterwards, it is presented an alternative modelling by means of a single equation involving an effective thermal conductivity which is a function of point inside the interfacial region.The theoretical results are illustrated by means of some numerical simulations for a model porous medium. In particular, temperature fields at the microscopic level are presented.Roman Letters sf interfacial area of thes-f interface contained within the macroscopic system m2 - A sf interfacial area of thes-f interface contained within the averaging volume m2 - C p mass fraction weighted heat capacity, kcal/kg/K - g vector that maps to s , m - h vector that maps to f , m - K eff effective thermal conductivity tensor, kcal/m s K - l s,l f microscopic characteristic length m - L macroscopic characteristic length, m - n fs outwardly directed unit normal vector for thef-phase at thef-s interface - n outwardly directed unit normal vector at the dividing surface. - R 0 REV characteristic length, m - T i macroscopic temperature at the interface, K - error on the external fluid temperature due to the macroscopic boundary condition, K - T * macroscopic temperature field obtained by solving the macroscopic Equation (3), K - V averaging volume, m3 - V s,V f volume of the considered phase within the averaging volume, m3. - mp volume of the porous medium domain, m3 - ex volume of the external fluid domain, m3 - s , f volume of the considered phase within the volume of the macroscopic system, m3 - dividing surface, m2 - x, z spatial coordinates Greek Letters s, f volume fraction - ratio of the effective thermal conductivity to the external fluid thermal conductivity - * macroscopic thermal conductivity (single equation model) kcal/m s K - s, f microscopic thermal conductivities, kcal/m s K - spatial average density, kg/m3 - microscopic temperature, K - * microscopic temperature corresponding toT *, K - spatial deviation temperature K - error in the temperature due to the macroscopic boundary conditions, K - * i macroscopic temperature at the interface given by the single equation model, K - spatial average - s , f intrinsic phase average.  相似文献   

15.
Zusammenfassung Durch die Phasenverteilung von Flüssigkeit und Dampf ergeben sich im horizontalen Verdampferrohr i.a. unterschiedliche lokale Wärmeübergangskoeffizienten (WÜK) am Rohrumfang. Die beiden technisch vorkommenden Randbedingungen — konstante Wärmestromdichte und konstante WandemperaturT W am Rohrinnendurchmesser — können unterschiedliche Auswirkungen auf umfangsgemittelte WÜK haben. Während bei konstanter aufgeprägter Wärmestromdichte ein interner Transport der zu übertragenden Wärmemenge durch die Wärmeleitung der Wand erfolgen kann, ist dies beiT W =konst. nicht möglich. Die Ergebnisse zeigen, daß sich der umfangsgemittelte WÜK bei konstanter Wärmestromdichte gegenüber dem Wert bei konstanter Wandtemperatur verringert. Dies ist jedoch nur dann der Fall, wenn der Rohrumfang unvollständig benetzt ist. Dabei hängt die Verminderung des umfangsgemittelten WÜK stark vom Wärmeleitvermögen der Wand, s, ab ( Wärmeleitfähigkeit,s Wandstärke). Es wurde festgestellt, daß die Verminderung sowohol aus der höheren Wandtemperatur am Rohrscheitel als auch aus der Änderung des benetzten Umfangs resultiert, da dann der Anteil des Rohrumfangs mit schlechtem Wärmeübergang vergrößert wird.
Effect of thermal boundary conditions on heat transfer in a horizontal evaporator tube
As a result of the non uniform phase distribution of liquid and vapour local heat transfer coefficients vary considerably along the circumference in a horizontal evaporator tube. It was expected that both modes of heating, uniform heat flux ( = const) and uniform inside wall temperature (T W =const), lead to different circumferential averaged heat transfer coefficients (HTC). The experiments show that the circumferential averaged HTC is lower in the case of uniform heat flux compared to the value of uniform wall temperature if the perimeter of the tube is only partly wetted. However, the reduction is affected by circumferential heat conduction. This reduction is a result of both, the higher wall temperature on the top of the tube and the reduction of the wetted perimeter.

Formelzeichen a Verhältnis von -Werten - c p spezifische Wärmekapazität - C F Faktor - d Innendurchmesser - g Fallbeschleunigung - h v spezifische Verdampfungsenthalpie - Massenstromdichte - M Rippenkenngröße - n Exponent der Wärmestromdichte - Nu Nusselt-Zahl - p Druck - p k Kritischer Druck - p* reduzierter Druck (p/p k ) - Pr L Prandtl-Zahl Flüssigkeit ( c p /) - Pr G Prandtl-Zahl Dampf ( c p /) - Wärmestromdichte - r Radius - Re Reynolds-Zahl - R p Glättungstiefe nach DIN 4762 - s Wandstärke - T Temperatur - Strömungsdampfgehalt - z Koordinate in Strömungsrichtung - Wärmeübergangskoeffizient - Differenz - Dampfvolumenanteil - dynamische Viskosität - Faktor - Wärmeleitfähigkeit - Druckverlustbeiwert - Dichte - Oberflächenspannung - normierter unbenetzter Bogen - Zentriwinkel - unbenetzter Bogen - Funktion - Faktor Indizes B Blasensieden - G Dampf - G0 gesamte Massenstromdichte als Dampf - k konvektives Strömungssieden - kr kritische Größe - L Flüssigkeit - Lb mit Flüssigkeit benetzt - L0 gesamte Massenstromdichte als Flüssigkeit - s Siedezustand - W Wand - gesättigte Flüssigkeit - gesättigter Dampf - – Mittelwert Herrn Prof. Dr.-Ing. Dr. h.c./INPL E.-U. Schlünder zum 60. Geburtstag gewidmet  相似文献   

16.
The problem of laminar forced convection heat transfer in short elliptical ducts with (i) uniform wall temperature and (ii) prescribed wall heat flux is examined in detail with the well known Lévêque theory of linear velocity profile near the wall. Moreover, consideration is given to the variation of the slope of the linear velocity profile with the position on the duct wall. A correction factor for the temperature dependent viscosity is included. Expressions for the local and average Nusselt numbers and wall temperatures are obtained. For the case of constant heat flux the Nusselt numbers are higher than for constant wall temperature.The results corresponding to the classical Graetz and Purday problems are deduced as special cases.Nomenclature a, b semiaxes of ellipse, b Graetz number (average), Re Pr D e/Z - h i o local heat transfer coefficient - J n(x) Bessel function of order n - K thermal conductivity of the fluid - [X] Laplace transform of X - N u o local Nusselt number, h i o D e/K - perimeter average Nusselt number - overall average Nusselt number - Nu w wall Nusselt number - Nu Nusselt number at large distance from the inlet - p Laplace transform parameter - Pr Prandtl number, C a/K - Re Reynolds number, D e / a - T temperature of the fluid - T 1, T W inlet and wall temperatures, respectively - u z local isothermal velocity along the axis of the duct - average fluid velocity - x, y, z Cartesian coordinates, z-axis parallel to the axis of the duct (z=0 at duct inlet) - Z length of the duct - thermal diffusivity, K/C - * correction factor for the temperature dependent viscosity - (x) gamma function - coordinate measured normal to the wall of the duct - a, w viscosity of fluid at average and wall temperatures - , , z elliptic cylindrical coordinates - density of fluid - (z) heat flux  相似文献   

17.
In this paper, a method using the mean velocity profiles for the buffer layer was developed for the estimation of the virtual origin over a riblets surface in an open channel flow. First, the standardized profiles of the mixing length were estimated from the velocity measurement in the inner layer, and the location of the edge of the viscous layer was obtained. Then, the virtual origins were estimated by the best match between the measured velocity profile and the equations of the velocity profile derived from the mixing length profiles. It was made clear that the virtual origin and the thickness of the viscous layer are the function of the roughness Reynolds number. The drag variation coincided well with other results.Nomenclature f r skin friction coefficient - f ro skin friction coefficient in smooth channel at the same flow quantity and the same energy slope - g gravity acceleration - H water depth from virtual origin to water surface - H + u*H/ - H false water depth from top of riblets to water surface - H + u*H/ - I e streamwise energy slope - I b bed slope - k riblet height - k + u*k/ - l mixing length - l s standardized mixing length - Q flow quantity - Re Reynolds number volume flow/unit width/v - s riblet spacing - u mean velocity - u* friction velocity = - u* false friction velocity = - y distance from virtual origin - y distance from top of riblet - y 0 distance from top of riblet to virtual origin - y v distance from top of riblet to edge of viscous layer - y + u*y/ - y + u*y/ - y 0 + u*y 0/ - u + u*y/ - shifting coefficient for standardization - thickness of viscous layer=y 0+y - + u*/ - + u*/ - eddy viscosity - ridge angle - v kinematic viscosity - density - shear stress  相似文献   

18.
On the boundary conditions at the macroscopic level   总被引:2,自引:0,他引:2  
We study the problem of the boundary conditions specified at the boundary of a porous domain in order to solve the macroscopic transfer equations obtained by means of the volume-averaging method. The analysis is limited to the case of conductive transport but the method can be extended to other cases. A numerical study enables us to illustrate the theoretical results in the case of a model porous medium. Roman Letters sf interfacial area of the s-f interface contained within the macroscopic system m2 - A sf interfacial area of the s-f interface contained within the averaging volume m2 - C p mass fraction weighted heat capacity, kcal/kg/K - d s , d f microscopic characteristic length m - g vector that maps to s, m - h vector that maps to f , m - K eff effective thermal conductivity tensor, kcal/m s K - l REV characteristic length, m - L macroscopic characteristic length, m - n fs outwardly directed unit normal vector for the f-phase at the f-s interface - n e outwardly directed unit normal vector at the dividing surface - T * macroscopic temperature field obtained by solving the macroscopic equation (3), K - V averaging volume, m3 - V s , V f volume of the considered phase within the averaging volume, m3 - volume of the macroscopic system, m3 - s , f volume of the considered phase within the volume of the macroscopic system, m3 - dividing surface, m2 Greek Letters s , f volume fraction - ratio of thermal conductivities - s , f thermal conductivities, kcal/m s K - spatial average density, kg/m3 - microscopic temperature, K - * microscopic temperature corresponding to T * , K - spatial deviation temperature K - error on the temperature due to the macroscopic boundary conditions, K - spatial average - s , f intrinsic phase average  相似文献   

19.
An in depth study into the development and decay of distorted turbulent pipe flows in incompressible flow has yielded a vast quantity of experimental data covering a wide range of initial conditions. Sufficient detail on the development of both mean flow and turbulence structure in these flows has been obtained to allow an implied radial static pressure distribution to be calculated. The static pressure distributions determined compare well both qualitatively and quantitatively with earlier experimental work. A strong correlation between static pressure coefficient Cp and axial turbulence intensity is demonstrated.List of symbols C p static pressure coefficient = (pw-p)/1/2 - D pipe diameter - K turbulent kinetic energy - (r, , z) cylindrical polar co-ordinates. / 0 - R, y pipe radius, distance measured from the pipe wall - U, V axial and radial time mean velocity components - mean value of u - u, u/ , / - u, , w fluctuating velocity components - axial, radial turbulence intensity - turbulent shear stress - u friction velocity, (u 2 = 0/p) - 0 wall shear stress - * boundary layer thickness A version of this paper was presented at the Ninth Symposium on Turbulence, University of Missouri-Rolla, October 1–3, 1984  相似文献   

20.
H. Potente 《Rheologica Acta》1978,17(4):406-414
Zusammenfassung In dieser Arbeit wird der Versuch unternommen, eine umfassende Modelltheorie für Kunststoff-Einschneckenmaschinen zu entwickeln. Ausgegangen wird vom analytischen Ähnlichkeitsprinzip.Wesentliche Merkmale der vorgestellten Theorie sind: a) das Energieprinzip, d. h. die dimensionslosen Kennzahlen werden aus Energiebilanzen hergeleitet, b) die strömungstechnische Ähnlichkeit wird durch die Konstanz des Drosselquotienten beschrieben, c) für die schergeschwindigkeitsabhängige Viskosität wird ein Potenzgesetz zugrunde gelegt, d) abgesehen von zwei Ausnahmen wird von einem konstantenL/D-Verhältnis ausgegangen.Es stellte sich heraus, daß man den bisher empirisch ermittelten Exponenten des Gangtiefenverhältnisses theoretisch ermitteln kann und daß dieser Exponent als dimensionslose rheologische Kennzahl deutbar ist. Mit dieser stoffabhängigen Kennzahl, die aufgrund bestimmter Voraussetzungen modifiziert werden kann, und dem dazugehörigen Drehzahlmodellgesetz lassen sich innerhalb bestimmter Grenzen für die übrigen Maschinen- und Verfahrensparameter Modellgesetze herleiten. Hier kann man zwischen einem universellen und zwei speziellen Modellgesetztypen, sowie deren Modifikationen aufgrund unterschiedlicher geometrischer Bedingungen unterscheiden.
Summary In this paper it is attempted to develop a comprehensive model theory for polymer single-screw machines proceeding from the analytical similarity principle.The essential points of the presented theory are: a) the energy principle, i.e. the dimensionless parameters are derived from energy balances, b) that the flow-characteristical similarity is described using a constant value of the throttle quotient, c) that a power-law approximation is taken as a basis for the shearrate dependent viscosity, d) a procedure using a constantL/D-ratio except for two special cases.It appeared that the until now empirically determined exponents of the channel depth ratios may be determined theoretically and that this exponent is a meaningful dimensionless rheological parameter. With this material-dependent parameter which may be modified on the basis of definite assumptions, and the related speed of rotation model law, model laws may be derived in some range of validity of the remaining machine and process parameters. Here it is distinguished between a universal and two special modellaw types including their modifications owing to different geometrical conditions.

Symbole A Betriebspunkt - a Drosselquotient - C Konstante - c p spezifische Wärme - D Schneckendurchmesser - D z Zylinderdurchmesser - H Heizleistung bzw. Kühlleistung der Ausstoßzone - H Heizleistung bzw. Kühlleistung - h 1 Gangtiefe der Einzugszone - h=h 2 Gangtiefe der Ausstoßzone - i(Index) Bezeichnung für die einzelnen Schneckenzonen - L Schneckenlänge - L Länge der Ausstoßzone - L w wirksame Schneckenlänge, Länge der Druckaufbauzone - Massedurchsatz - M d Drehmoment - n Drehzahl - P am Schneckenschaft abgegebene Leistung - P der Ausstoßzone zugeführte Leistung - p Druck in der Trennebene zwischen Werkzeug bzw. Düse und Maschine - Vom Zylinder pro Flächen- und Zeiteinheit zugeführte Wärme - Volumendurchsatz - Volumendurchsatz infolge Druckströmung - Volumendurchsatz infolge Schleppströmung - v Geschwindigkeit - x frei wählbarer Exponent im Drehzahlmodellgesetz - Wärmeübergangszahl - Schergeschwindigkeit - Viskosität - 1 Extrudattemperatur - Mittlere Massetemperatur der Ausstoßzone - z Zylinderwandtemperatur - h Temperaturdifferenz zwischen Masse und Zylinderwand - ( 10) Gesamtmassetemperaturdifferenz - L=(1 – ) Massetemperaturdifferenz der Ausstoßzone - Steigung der Viskositätskurve - Wärmeleitfähigkeit - v Steigung der Fließkurve - Steigung der Massetemperatur — Schergeschwindigkeitskurve bei konstanter Viskosität - Dichte - Schubspannung - Gangsteigungswinkel - dimensionslose rheologische Kennzahl Gekürzte Fassung eines Plenarvortrags, gehalten während des 9. Kunststofftechnischen Kolloquiums des IKV in Aachen vom 8.–10. März 1978.Mit 8 Abbildungen und 5 Tabellen  相似文献   

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

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