首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 62 毫秒
1.
Effect of spanwise jet-to-jet spacing on local heat transfer distribution due to an in-line rectangular array of confined multiple circular air jets impinging on a surface parallel to the jet plate are studied experimentally. Length-to-diameter ratio of nozzles of the jet plate is 1.0. The flow, after impingement, is constrained to exit in two opposite directions from the confined passage formed between jet plate and target plate. Mean jet Reynolds numbers based on the nozzle exit diameter (d) covered are 3000, 5000, 7500 and 10,000 and jet-to-plate spacings studied are d, 2d and 3d. Spanwise pitches considered are 2d, 4d and 6d in steps of 2d keeping the streamwise pitch at 5d. For all the configurations, the jet-plates have ten spanwise rows in streamwise direction and six jets in each spanwise row. Flat heat transfer surface is made of thin stainless steel metal foil. Local temperature distribution on a target plate is measured using thermal infrared camera. Wall static pressure on the target plate is measured in the streamwise direction to estimate crossflow velocities and individual jet velocities. Heat transfer characteristics are explained on the basis of the flow distribution. A simple correlation to predict streamwise distribution of heat transfer coefficients averaged over each spanwise strip resolved to one jet hole is developed.  相似文献   

2.
An experimental investigation is performed to study the effect of jet to plate spacing and low Reynolds number on the local heat transfer distribution to normally impinging submerged circular air jet on a smooth and flat surface. A single jet from a straight circular nozzle of length-to-diameter ratio (l/d) of 83 is tested. Reynolds number based on nozzle exit condition is varied between 500 and 8,000 and jet-to-plate spacing between 0.5 and 8 nozzle diameters. The local heat transfer characteristics are obtained using thermal images from infrared thermal imaging technique. It was observed that at lower Reynolds numbers, the effect of jet to plate distances covered during the study on the stagnation point Nusselt numbers is minimal. At all jet to plate distances, the stagnation point Nusselt numbers decrease monotonically with the maximum occurring at a z/d of 0.5 as opposed to the stagnation point Nusselt numbers at high Reynolds numbers which occur around a z/d of 6.  相似文献   

3.
The rate of heat transfer from a plate due to impinging of an array of jets was investigated. The effect of jet-to-plate spacing in a confined array of impinging laminar square jets was investigated numerically through the solution of Navier Stokes and energy equations. The simulation is carried out for the jet-to-plate spacing between 2 B and 20 B and for jet-to-jet spacing of 4 B, where B is the jet width. Five in-line jets subjected to across-flow were used in this investigation. Also, six different ratios of jet to cross-flow velocity are simulated (0.5, 1.0, 2.5, 5, 7.5 and 10) for the jet Reynolds number of 200. The predicted results show a formation of one or two ground horseshoe vortices between the jets. In addition, a horseshoe vortex forms at different position between the orifice and impinging plates due to the interaction of two jets before they combine. The number of the ground horseshoe vortex and its size are strongly affected by the jet-to-plate spacing and by jet to cross-flow velocity ratio. The effect of jet-to-plate spacing and jet to cross-flow velocity ratio on heat transfer is presented and discussed.  相似文献   

4.
A computational investigation is carried out to study the flow and heat transfer from a row of circular jets impinging on a concave surface. The computational domain simulates the impingement cooling zone of a gas turbine nozzle guide vane. The parameters, which are varied in the study include jet Reynolds number (Re d = 5000–67800), inter-jet distance to jet diameter ratio (c/d = 3.33 and 4.67) and target plate distance to jet diameter ratio (H/d = 1, 3 and 4). The flow field, predicted with K-ω turbulence model and using Fluent 6.2.16, is characterized with the presence of a pair of counter rotating vortices, an upwash fountain flow and entrainment. The local pressure coefficient and Nusselt number variations along the concave plate are presented and these values are found to under predict the available experimental data by about 12%.  相似文献   

5.
We report flow visualisations and laser Doppler anemometry (LDA) velocity measurements in the near field of two swirling jets. The Reynolds number based on jet diameter and bulk velocity at the nozzle exit is 1.4 × 105. In the first jet, a small recirculation region is formed around the jet axis, while, in the second, the streamwise velocity remains positive and overshoots near the jet centre. In both cases, flow visualisations show that the vortex core of the jets is depleted of seeding particles. By using time-averaged distributions of the streamwise and tangential velocities measured at the nozzle outlet, the dynamics of the particles is simulated, by integrating their simplified equations of motion. The particles trajectory thus computed agrees well with that observed in the flow visualisations. Although the turbulence intensity is substantially different in the core of the two jets, its effect on the seeding concentration is localised near the edge of the core.  相似文献   

6.
This paper reports a numerical investigation on the prediction of the thermal and hydrodynamic flow fields of multi micro jet impingement cooling of three dimensional turbine vanes. A three dimensional vane is modeled with an in-line array of impinging jets of diameters 0.5 and 0.25 mm. The numerical model consists of the steady, Reynolds-Averaged Navier–Stokes equations and the Kω SST Turbulence model. The governing equations are solved using a finite volume method. The crossflow mass velocity (G c ) to jet mass velocity (G j ) ratio, and the average and local heat transfer distributions are analyzed with varying mass velocity and jet-to-target spacing. It is found out that a significant decrease in crossflow ratio occurs with the smaller diameters. Due to the lower crossflow and higher exit velocities of the smaller jets, the penetration into the crossflow is much higher. Moreover, at a constant mass flow, the use of micro-jets enhanced the overall average heat transfer coefficient by 63%, while at a fixed pressure drop across the vane instead of the mass flow, the smaller diameters will still yield an enhancement of 34.3% in the overall average heat transfer coefficient.  相似文献   

7.
Large eddy simulations are performed for an unsteady flow and heat transfer in the region of interaction of a circular turbulent jet with a normally positioned flat obstacle (target). Space-filtered Navier-Stokes equations are closed by the RNG model of eddy viscosity, which takes into account the curvature of streamlines in the region of flow turning. The computations are performed for different dimensionless distances between the nozzle exit and the target and for different Reynolds numbers. The dependence between the Nusselt number distribution over the target surface and the vortex structure of the jet is analyzed. The local and integral characteristics of the flow are compared with the data of a physical experiment. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 55–67, January–February, 2007.  相似文献   

8.
The results of measurements of all three components of the mean velocity vector, the Reynolds normal and primary shear stresses and the mean static pressure in a turbulent free jet, issuing from a sharp-edged cruciform orifice, are presented in this paper. The measurements were made with an x-array hot-wire probe and a pitot-static tube in the near flow field of the jet. The Reynolds number, based upon the equivalent diameter of the orifice, was 1.70 × 105. In addition to the quantities measured directly, the mean streamwise centreline velocity decay, the jet half-velocity widths, the jet spreading rate, the mean streamwise vorticity, the mass entrainment rate, the integral momentum flux and the one-dimensional energy spectra have been derived from the measured data. The results show that the mean streamwise centreline velocity decay rate of the cruciform jet is higher than that of a round jet issuing from an orifice with the same exit area as that of the cruciform orifice. The mean streamwise velocity field changed shape continuously from a cruciform close to the orifice exit plane to circular at 12 and half equivalent diameters downstream. The mean streamwise vorticity field, up to about three equivalent diameters downstream of the orifice exit plane, consists of four pairs of counter-rotating cells, which are aligned with the four edges in the centre of the cruciform orifice.  相似文献   

9.
An experimental study has been carried out of the low speed Coanda wall jet with both streamwise and axisymmetric curvature. A single component laser Doppler technique was used, and by taking several orientations at a given point, values of the three mean velocities and five of the six Reynolds stresses were obtained. The lateral divergence and convex streamwise curvature both enhanced the turbulence in the outer part of the jet compared with a plane two-dimensional wall jet. The inner layer exhibited a large separation of the positions of maximum velocity and zero shear stress. It was found that the streamwise mean velocity profile became established very rapidly downstream of the slot exit. The profile appeared fairly similar at later downstream positions, but the mean radial velocity and turbulence parameters showed the expected nonself preservation of the flow. Removal of the streamwise curvature resulted in a general return of the jet conditions toward those expected of a plane wall jet. The range and accuracy of the data may be used for developing turbulence models and computational techniques for this type of flow.  相似文献   

10.
An experimental investigation of flow-induced vibrations of gates with multiple degrees of freedom is presented. An underflown vertical gate plate with submerged discharge was allowed to oscillate both in the cross-flow (z -) and in the streamwise (x -) direction. The two purposes of the investigation were to further the insight into the hydrodynamic coupling mechanisms of the two vibration modes and to determine the interaction of the unsteady lift and drag forces. Self-excited vibration tests were run with reduced velocities VrzandVrx from 0·8 to 14, covering a range in which the instability-induced excitation (IIE) due to impinging-leading-edge vortices (ILEV) as well as the transition to galloping (MIE) occurred. The ratio of the natural frequencies of the two vibration modes fx 0/fz 0, the gate opening ratios /d, and the submergence of the gate plate were varied. Depending on the ranges of reduced velocities and frequency ratios, a complex interaction of two different kinds of instability-induced excitation was detected. Furthermore, it was found that streamwise IIE-excitation and cross-flow galloping coexist. To assess the relevant fluid dynamic amplification and attenuation mechanisms, simultaneous body response and flow velocity measurements were carried out.  相似文献   

11.
Heat transfer from an obliquely impinging circular, air jet to a flat plate   总被引:5,自引:0,他引:5  
A series of experiments was conducted for the measurement of local convective heat transfer coefficients for an obliquely impinging circular air jet to a flat plate. In the experiments, the oblique angles selected were 90°, 75°, 60° and 45°, with 90° being a vertical jet. Two different Reynolds numbers of 10,000 and 23,000 were considered for the purpose of comparison with previous data available in the literature. Another parameter varied in the measurements was the dimensionless jet-to-plate distance, L/D. Four values of L/D(2, 4, 7, and 10) were considered in the experiments. The experiments were conducted using the preheated wall transient liquid-crystal technique. Liquid-crystal color changes were recorded with a video system. Local convective heat transfer coefficients were obtained through the surface transient temperatures that were related to the recorded color information. Detailed local heat transfer coefficients were presented and discussed in relation to the asymmetric wall jet upon impingement of the jet flow. Results of experiments show that, for a given flow situation, the point of maximum heat transfer shifts away from the geometrical impingement point toward the compression side of the wall jet on the axis of symmetry. The shift is more pronounced with a smaller oblique angle (larger jet inclination) and a smaller jet-to-plate distance. Comparisons of experimental results with existing heat transfer data for both obliquely impinging jets and vertical impinging jets are made. The effect of oblique angles on heat transfer was assessed.  相似文献   

12.
A comparative study of the length scales and morphology of dissipation fields in turbulent jet flames and non-reacting jets provides a quantitative analysis of the effects of heat release on the fine-scale structure of turbulent mixing. Planar laser Rayleigh scattering is used for highly resolved measurements of the thermal and scalar dissipation in the near fields of CH4/H2/N2 jet flames (Re d  = 15,200 and 22,800) and non-reacting propane jets (Re d  = 7,200–21,700), respectively. Heat release increases the dissipation cutoff length scales in the reaction zone of the flames such that they are significantly larger than the cutoff scales of non-reacting jets with comparable jet exit Reynolds numbers. Fine-scale anisotropy is enhanced in the reaction zone. At x/d = 10, the peaks of the dissipation angle PDFs in the Re d  = 15,200 and 22,800 jet flames exceed those of non-reacting jets with corresponding jet exit Reynolds numbers by factors of 2.3 and 1.8, respectively. Heat release significantly reduces the dissipation layer curvature in the reaction zone and in the low-temperature periphery of the jet flames. These results suggest that the reaction zone shields the outer regions of the jet flame from the highly turbulent flow closer to the jet axis.  相似文献   

13.
Simultaneous velocity and concentration fields in a confined liquid-phase rectangular jet with a Reynolds number based on the hydraulic diameter of 50,000 (or 10,000 based on the velocity difference between streams and the jet exit dimension) and a Schmidt number of 1,250 were obtained by means of a combined particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF) system. Data were collected at the jet exit and six further downstream locations. The velocity and concentration field data were analyzed for flow statistics such as turbulent fluxes, turbulent viscosity and diffusivity, and turbulent Schmidt number (Sc T ). The streamwise turbulent flux was found to be larger than the transverse turbulent flux, and the mean concentration gradient was not aligned with the turbulent flux vector. The average Sc T was found to vary both in streamwise and in cross stream directions and had a mean value around 0.8, a value consistent with the literature. Spatial correlation fields of turbulent fluxes and concentration were then determined. The R u′ϕ′ correlation was elliptical in shape with a major axis tilted downward with respect to the streamwise axis, whereas the R v′ϕ′ correlation was an ellipse with a major axis aligned with the cross-stream direction. Negative regions of R u′ϕ′ were observed in the outer streams, and these negatively correlated regions decayed with downstream distance and finally disappeared altogether. The R ϕ′ϕ′ correlation field was found to be an ellipse with the major axis inclined at about 45° with respect to the streamwise direction. Linear stochastic estimation was used to interpret spatial correlation data and to determine conditional flow structures. It is believed that a vortex street formed near the splitter plate is responsible for the negatively correlated region observed in the R u′ϕ′ spatial correlations of turbulent fluxes. A positive concentration fluctuation event was observed to correspond to a finger of nearly uniform concentration fluid reaching out into the outer stream, whereas a negative event corresponds to a pocket of nearly uniform fluid being entrained from the outer stream into the center jet region. Large-scale vortical structures were observed in the conditional velocity fields with an elliptical shape and a streamwise major axis. The growth of the structure size increased linearly initially but then grew more slowly as the flow transitioned toward channel flow. Support of this work was provided by the National Science Foundation through grants CTS-9985678 and CTS-0336435 and by the Dow Chemical Company. The author greatly acknowledge Charles Lipp at Dow Chemical and Ken Junk at Emerson Fisher for their valuable assistance in the design and construction of the flow system.  相似文献   

14.
Pairs of counter-rotating streamwise vortical structures have been observed using the LIF (Laser Induced Fluorescence) flow visualization technique by means of regular, video and high-speed photography in the near-field of an axisymmetric water jet. The temporal evolution of these structures at fixedx/d=1.0, 2.0, 4.0, and 6.0 in planes perpendicular to the flow direction for Reynolds numbersRe=cd/v=6000 has been the focus of attention in this investigation (d = nozzle diameter,c = mean exit velocity). These streamwise vortical structures also appear when the formation of Taylor-Görtler vortices in the nozzle is suppressed. So far basic questions pertaining to the generation and development of these secondary structures remain unanswered. Understanding of instability and interaction mechanisms can be gained by the analysis of these structures. These also contribute to entrainment and mixing due to their dynamics and large interfacial area.
Sommario Coppie di strutture tubolari controrotanti sono state osservate usando la tecnica di visualizzazione del flusso LIF (Laser Induced Fluorescence), per mezzo di una telecamera ad alta e bassa velocità, nel campo vicino ad un getto d'acqua assialsimmetrico. Lo studio dell'evoluzione temporale di queste strutture, fissatox/d=1.0, 2.0, 4.0 e 6.0 in piani perpendicolari alla direzione del flusso ed in sezioni longitudinali perx/d=0÷6.0, per numeri di Reynolds Re=cd/v=6000, è stato il punto principale della presente indagine. Queste strutture vorticose appaiono persino quando la formazione dei vortici di Taylor-Görtler è soppressa attivamente e passivamente. Si osserva che rimangono senza risposta questioni basilari concernenti la generazione e lo sviluppo di queste strutture secondarie. Attraverso l'analisi di queste strutture è possibile far progressi nella conoscenza dell'instabilità e dei meccanismi d'interazione. Inoltre, tali strutture contribuiscono all'intrappolamento ed al mescolamento dovuti alla loro dinamica ed alla larga area interfacciale.
  相似文献   

15.
Investigation of the mixing process is one of the main issues in chemical engineering and combustion and the configuration of a jet into a cross-flow (JCF) is often employed for this purpose. Experimental data are gained for the symmetry plane in a JCF-arrangement of an air flow using a combination of particle image velocimetry (PIV) with laser-induced fluorescence (LIF). The experimental data with thoroughly measured boundary conditions are complemented with direct numerical simulations, which are based on idealized boundary conditions. Two similar cases are studied with a fixed jet-to-cross-flow velocity ratio of 3.5 and variable cross-flow Reynolds numbers equal to 4,120 and 8,240; in both cases the jet issues from the pipe at laminar conditions. This leads to a laminar-to-turbulent transition, which depends on the Reynolds number and occurs quicker for the case with higher Reynolds number in both experiments and simulations as well. It was found that the Reynolds number only slightly affects the jet trajectory, which in the case with the higher Reynolds number is slightly deeper. It is attributed to the changed boundary layer shape of the cross-flow. Leeward streamlines bend toward the jet and are responsible for the strong entrainment of cross-flow fluid into the jet. Velocity components are compared for the two Reynolds numbers at the leeward side at positions where strongest entrainment is present and a pressure minimum near the jet trajectory is found. The numerical simulations showed that entrainment is higher for the case with the higher Reynolds number. The latter is attributed to the earlier transition in this case. Fluid entrainment of the jet in cross-flow is more than twice stronger than for a similar flow of a jet issuing into a co-flowing stream. This comparison is made along the trajectory of the two jets at a distance of 5.5 jet diameters downstream and is based on the results from the direct numerical simulations and recently published experiments of a straight jet into a co-flow. Mixing is further studied by means of second-order statistics of the passive scalar variance and the Reynolds fluxes. Windward and leeward sides of the jet exhibit different signs for the time-averaged streamwise Reynolds flux 〈v x c′〉. The large coherent structures which contribute to this effect are investigated by means of timely correlated instantaneous PIV-LIF camera snapshots and their contribution to the average statistics of 〈v x c′〉 are discussed. The discussion on mixing capabilities of the jet in cross-flow is supported by simulation results showing the instantaneous three-dimensional coherent structures defined in terms of the pressure fluctuations.  相似文献   

16.
The intrinsic characteristics of coherent structures in the near field of a plane jet are extensively studied by hot-wire measurements. The instability modes which are responsible for the dynamics of the coherent structures are found to exhibit distinct evolution characters at different transverse positions of the shear layer along downstream direction. The occurrence of multiple peaks in the energy spectra depicts the formation of the sideband instabilities in the early stage of the jet flow field. These sideband instabilities are investigated to be induced by the mechanisms of the nonlinear interactions between neighboring fundamental and subharmonic instabilities, and the feedback effects of the preferred mode near the end of the potential core. Also, from the spatial distributions of the instability modes over the jet flow field, Ho's subharmonic evolution model (1982) is further examined with more interpretations.List of symbols E (f) energy content of streamwise velocity fluctuation at spe cific frequency - f e excitation frequency - f 0 fundamental frequency - f p preferred frequency - f r response frequency in an excited jet - H height of the plane jet at the exit - U streamwise mean velocity - U 0 mean velocity at the nozzle exit - U c mean velocity at the jet center line - u streamwise RMS velocity fluctuation - u p peak streamwise velocity fluctuation alongY axis - u (f) amplitude of streamwise velocity fluctuation at specific frequency - X, Y streamwise and transverse coordinates - Y a transverse position whereU = aU c ,a = 0.99, 0.9,..., etc. - Y c transverse position at the jet center line - 0 initial instability wave length (=U 0/2f 0) - 0 momentum thickness - 0 initial boundary layer momentum thickness A version of this paper was presented at the 11th Symposium on Turbulence, University of Missouri-Rolla, Oct. 17–19, 1988  相似文献   

17.
An experimental investigation is performed to study the effect of the finned surfaces and surfaces with vortex generators on the local heat transfer coefficient between impinging circular air jet and flat plate. Reynolds number is varied between 7000 and 30,000 based on the nozzle exit condition and jet to plate spacing between 0.5 and 6 nozzle diameters. Thermal infrared imaging technique is used for the measurement of local temperature distribution on the flat plate. Fins used are in the form of cubes of 2 mm size spaced at a pitch of 5 mm on the target plate and hexagonal prism of side 2.04 mm and height of 2 mm spaced at a pitch of 7.5 mm. Vortex generators in the form of a equilateral triangle of side 4 mm are used. Effect of number of rows of vortex generators, radius of a row, number of vortex generators in a row and inclination angle (i.e., the angle between the plane of the target plate and the plane of the vortex generators) on Nusselt number is studied. It is observed that the heat transfer coefficient between the impinging jet and the target plate is sensitive to the shape of the fin. The increase in the heat transfer coefficient up to 77% depending on the shape of the fin, nozzle plate spacing and the Reynolds number is observed. The augmentation in the heat transfer for the surfaces vortex generators are higher than that of the finned surfaces. The heat transfer augmentation in case of vortex generator is as high as 110% for a single row of six vortex generators at a radius of 1 nozzle diameter as compared to the smooth surface at a given nozzle plate spacing of 1 nozzle diameter and a Reynolds number of 25,000 at extreme radial location.  相似文献   

18.
Velocity and passive scalar (temperature) measurements have been made in the near field of a round jet with and without obstructing grids placed at the jet exit. The Reynolds number Re D (based on the exit centreline velocity and nozzle diameter) is 4.9 × 104 and the flow is incompressible, while the temperature rise does not affect the velocity behaviour. The streamwise development and radial spreading of the passive scalar are attenuated, relative to the unobstructed jet. Close to the jet outlet, the spatial similarity of the moments (up to the third-order) of velocity fluctuations is improved, when the jet is perturbed. An explanation, based on the reduced effect of the large coherent structures in the developing region, is provided.  相似文献   

19.
 Experiments have been performed to assess the impact of an extended surface on the heat transfer enhancement for axisymmetric, turbulent liquid jet impingement on a heated round disk. The disk, with an array of integral radial fins mounted on its surface, is placed at the bottom of an open vertical circular cavity. Hydrodynamic and heat transfer data were obtained for a dielectric fluorocarbon liquid FC-77. For a fixed circular heater of diameter D=22.23 mm, several geometric parameters were tested: the nozzle diameter (4.42≤d≤9.27 mm), the confining wall diameter of the vertical cavity (22.23≤D c≤30.16 mm), and the nozzle-to-heater spacing (0.5≤S/d≤5.0). The FC-77 flow rates varied from =0.2 to 11.0 l/min producing Reynolds numbers in the wide interval 700≤Re d ≤44,000. For d=4.42 mm, the heat transfer response to the separation distance S/d was small but increased gradually with increasing nozzle diameter up to d=9.27 mm. The thermal resistance R th increased with the confining wall diameter D c and also with the nozzle diameter d. A minimum value of the thermal resistance of R th,min=0.4 cm2 K/W was attained for a combination of d=4.42 mm, D c=22.23 mm, S/d=1, and =7.5 l/min. Based on a simplified heat transfer model, reasonable agreement was obtained between measured values of the thermal resistance and the R th-predictions. The total fin effectiveness ɛf was shown to increase with increasing nozzle diameter, but was invariant with the flow rate (or the jet exit velocity). More than a three-fold heat transfer enhancement was realized through the addition of the array of integral radial fins on the heated round disk. Received on 30 August 2000 / Published online: 29 November 2001  相似文献   

20.
We investigate the gas-particle dynamics of a device designed for biological pre-clinical experiments. The device uses transonic/supersonic gas flow to accelerate microparticles such that they penetrate the outer skin layers. By using a shock tube coupled to a correctly expanded nozzle, a quasi-one-dimensional, quasi-steady flow (QSF) is produced to uniformly accelerate the microparticles. The system utilises a microparticle “cassette” (a diaphragm sealed container) that incorporates a jet mixing mechanism to stir the particles prior to diaphragm rupture. Pressure measurements reveal that a QSF exit period – suitable for uniformly accelerating microparticles – exists between 155 and 220 mus after diaphragm rupture. Immediately preceding the QSF period, a starting process secondary shock was shown to form with its (x,t) trajectory comparing well to theoretical estimates. To characterise the microparticle, flow particle image velocimetry experiments were conducted at the nozzle exit, using particle payloads with varying diameter (2.7–48 μm), density (600–16,800 kg/m3) and mass (0.25–10 mg). The resultant microparticle velocities were temporally uniform. The experiments also show that the starting process does not significantly influence the microparticle nozzle exit velocities. The velocity distribution across the nozzle exit was also uniform for the majority of microparticle types tested. For payload masses typically used in pre-clinical drug and vaccine applications (≤ 1 mg), it was demonstrated that payload scaling does not affect the microparticle exit velocities. These characteristics show that the microparticle exit conditions are well controlled and are in agreement with ideal theory. These features combined with an attention to the practical requirements of a pre-clinical system make the device suitable for investigating microparticle penetration into the skin for drug delivery.  相似文献   

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

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