Theoretical investigation of the effect of an internal heat source on turbulent liquid metal heat transfer in a channel with uniform heat flux

Summary The effect of an internal heat source on the heat transfer characteristics for turbulent liquid metal flow between parallel plates is studied analytically. The analysis is carried out for the conditions of uniform internal heat generation, uniform wall heat flux, and fully established temperature and velocity profiles. Consideration is given both to the uniform or slug flow approximation and the power law approximation for the turbulent velocity profile. Allowance is made for turbulent eddying within the liquid metal through the use of an idealized eddy diffusivity function. It is found that the Nusselt number is unaffected by the heat source strength when the velocity profile is assumed to be uniform over the channel cross section. In the case of a 1/7-power velocity expression, the Nusselt numbers are lower than those in the absence of internal heat generation, and decrease with diminishing eddy conduction. Nusselt numbers, in the absence of an internal heat source, are compared with existing calculations, and indications are that the present results are adequate for preliminary design purposes.Nomenclature A hydrodynamic parameter - a half height of channel - a ₁ a constant, 1+0.01 Pr Re ^0.9 - a ₂ a constant, 0.01 Pr Re ^0.9 - C _p specific heat at constant pressure - D _h hydraulic diameter of channel, 4a - h heat transfer coefficient, q _w/(t _w–t _b) - I ₁ integral defined by (17) - I ₂ integral defined by (18) - k diffusivity parameter, (1+0.01 Pr Re ^0.9)^1/2 - m exponent in power velocity expression - Nu Nusselt number, hD _h/ - Nu ₀ Nusselt number in absence of internal heat generation - Pr Prandtl number, / - Q heat generation rate per volume - q _w wall heat flux - Re Reynolds number for channel, 2/ - s ratio of heat generation rate to wall heat flux, Qa/q _w - T dimensionless temperature, (t _w–t)/(t _w–t _b) - t fluid temperature, t _w wall temperature, t _b fluid bulk temperature - u fluid velocity in x direction, , fluid mean velocity - x longitudinal coordinate measured from channel entrance - x ⁺ dimensionless longitudinal coordinate, 2(x/a)/Pr Re - y transverse coordinate measured from channel centerline - z transverse coordinate measured from channel wall, a–y - molecular diffusivity of heat, /C _p - dummy variable of integration - dummy variable of integration - _H eddy diffusivity of heat - _M eddy diffusivity of momentum - dummy variable of integration - fluid thermal conductivity - _T dimensionless diffusivity, Pr ( _H/) - fluid kinematic viscosity - dummy variable of integration - fluid density - dummy variable of integration - ratio of eddy diffusivity for heat transfer to that for momentum transfer, _H/ _M - average value of - dimensionless velocity distribution, u/     

Peristaltic pumping of a second-order fluid in a planar channel

The peristaltic motion of a non-Newtonian fluid represented by the constitutive equation for a second-order fluid was studied for the case of a planar channel with harmonically undulating extensible walls. A perturbation series for the parameter ( half-width of channel/wave length) obtained explicit terms of 0(²), 0(²Re²) and 0(₁Re²) respectively representing curvature, inertia and the non-Newtonian character of the fluid. Numerical computations were performed and compared to the perturbation analysis in order to determine the range of validity of the terms.Presented at the second conference Recent Developments in Structured Continua, May 23–25, 1990, in Sherbrooke, Québec, Canada     

Spectral and correlation characteristics of the turbulent boundary layer on an extended flexible cylinder

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 Re_X=(2–4)·10⁷, a dimensionless diameter of the pressure fluctuation sensors d _p ⁺ =d_pu^*/=500, where d_p 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.     

An approximate theoretical analysis and experimental verification of turbulent entrance region flow of drag reducing fluids

Summary Entry lengths for pipe flows of moderately drag reducing fluids are determined using momentum integral technique. It is shown theoretically that the entry lengths for drag reducing fluids could be significantly larger than the Newtonian fluids flowing turbulently under otherwise identical conditions. The experimental data from the literature bear out the theoretical calculations.

---

Zusammenfassung Mit Hilfe der Impuls-Methode wird die Einlauflänge in einer Rohrströmung für Flüssigkeiten mit mäßig starker Widerstandsverminderung berechnet. Es wird vorausgesagt, daß die Einlauflänge für derartige Flüssigkeiten erheblich größer sein kann als für newtonsche Flüssigkeiten unter sonst identischen Bedingungen. Aus der Literatur entnommene experimentelle Daten bestätigen diese theoretischen Berechnungen.

---

Nomenclature A ₁ Coefficient in eq. [7] - A Slope of logarithmic velocity profile - a Exponent in eq. [10] - B Intercept function for logarithmic velocity profile - De Deborah number, - f Friction factor - F Function, eq. [30] - G Function given in eq. [11] - Static pressure, dynes/cm² - q Index of power law velocity profile - R Pipe radius, cm - r Radial distance, cm - R Core radius, cm - Re Reynolds number - Axial velocity, cm/s - u _c Core velocity, cm/s - u ⁺ Dimensionless velocity, eq. [5] - u ^* Friction velocity, , cm/s - Radial velocity, cm/s - V Average velocity, cm/s - x Axial distance, cm - x _e Entry length, cm - y Distance from the wall, cm - y ⁺ Dimensionless distance, eq. [5] - y _I ⁺ Dimensionless viscous sublayer thickness - coefficient in eq. [17] - exponent of Reynolds number in eq. [17] - shear rate, s^–1 - boundary layer thickness, cm - _fl fluid relaxation time, s - µ fluid viscosity, gm/cm s - v kinematic viscosity, cm²/s - _l laminar sublayer thickness, dimensionless - fluid density, gm/cm³ - shear stress, dynes/cm² - _w shear stress at the wall, dynes/cm² - ₁, ₂, ₃, ₄ functions in eq. [27] - ~ time averaged quantities - — dimensionless quantity With 3 figures and 1 table     

Viscoelastic properties of semidilute poly(methyl methacrylate) solutions

Viscoelastic properties were examined for semidilute solutions of poly(methyl methacrylate) (PMMA) and polystyrene (PS) in chlorinated biphenyl. The number of entanglement per molecule, N, was evaluated from the plateau modulus, G _N . Two time constants, _s and ₁, respectively, characterizing the glass-to-rubber transition and terminal flow regions, were evaluated from the complex modulus and the relaxation modulus. A time constant _k supposedly characterizing the shrink of an extended chain, was evaluated from the relaxation modulus at finite strains. The ratios ₁/ _s and _k / _s were determined solely by N for each polymer species. The ratio ₁/ _s was proportional to N ^4.5 and N ^3.5 for PMMA and PS solutions, respectively. The ratio _k / _s was approximately proportional to N ^2.0 in accord with the prediction of the tube model theory, for either of the polymers. However, the values for PMMA were about four times as large as those for PS. The result is contrary to the expectation from the tube model theory that the viscoelasticity of a polymeric system, with given molecular weight and concentration, is determined if two material constants _s and G _N are known.     

An experimental study of the lateral migration of a droplet in a creeping flow

The distribution of droplets in a plane Hagen-Poiseuille flow of dilute suspensions has been measured by a special LDA technique. This method assumes a well defined relation between the velocity of the droplets and their lateral position in the channel. The measurements have shown that the droplet distribution is non-uniform and depends on the viscosity ratio between the droplets and the carrier liquid. The results have been compared with a theory by Chan and Leal describing the lateral migration of suspended droplets.List of symbols a particle radius, m - d half width of the channel, m - Re flow Reynolds number, = 2 _m · d · /µ - flow velocity, m/s - _m flow velocity at the channel axis, m/s - We Weber number, = 2 _m ^Emphasis>/2 · d · / - x distance from center line (x = 0) of the channel, m - non-dimensional distance from the channel center line, _x ^d - y distance along the channel (y = 0 at channel inlet), m - non-dimensional distance along the channel, = y/2d - non-dimensional, normalized distance along the channel, = · _m · µ/ - interfacial tension, N/m - viscosity ratio of dispersed (droplet) phase to viscosity of continuous phase - µ viscosity of continuous phase, Pa · s - density of continuous phase, kg/m³ - phase density difference, kg/m³ Experiments were performed at Max-Planck-Institut, Göttingen     

Rheologie pathologischer Gelenkflüssigkeiten

Zusammenfassung Zur Charakterisierung der Lösungsstruktur pathologischer Synovialflüssigkeiten wurden parallel rheologische Untersuchungen sowohl im Bereich konzentrierter als auch verdünnter Lösungen durchgeführt. Durch die Berücksichtigung der Scherabhängigkeit der Grenzviskositätszahl [] und Messung bei niedrigen Schergeschwindigkeiten (D 2,9 10^–3 – 4,6 s^–1) errechneten sich höhere Molekulargewichte der Hyaluronsäure (HA) als bisher angegeben. Für entzündliche Gelenkergüsse betragen sie (2,4–12,0) 10⁶, für nicht entzündliche (3,1–11,4) 10⁶ und sind damit gegenüber der normalen Synovia mit 16,2 10⁶ erniedrigt. Unterschiedlichek _SB-Werte sprechen zusätzlich für eine stärkere Verknäuelung der HA-Makromoleküle in der pathologischen Synovia.Durch Ermittlung der kritischen Polymerkonzentration der HA sowie der Verhängungszahlenv aus Anlaufmessungen konnten quantitative Aussagen zum makromolekularen Netzwerk getroffen werden. Die normale Synovia entspricht einem sehr engmaschigen und dichten Netz, aber auch in den pathologischen Gelenkflüssigkeiten ist die HA nicht im Sinne einer Partikellösung isoliert, sondern bildet ein, wenn auch lockeres Netzwerk mit wesentlich niedrigerer Verhängungszahl. Eine Erniedrigung der Viskosität und Elastizität isolierter HA weist darauf hin, daß molekulare Wechselwirkungen zwischen HA und anderen Makromolekülen in der Synovialflüssigkeit für das rheologische Verhalten von wesentlicher Bedeutung sind. Relaxationszeitmessungen in Ruhe und unter Scherung zeigten eine vermehrte Segmentbeweglichkeit des Netzwerkes der pathologischen Synovia als Folge der HA-Konzentrationsabnahme und/oder Molekulargewichtsverminderung.

---

Summary Rheological investigations were done in concentrated as well as diluted solutions to characterize the solution structure of pathological synovial fluids. By looking at the shear rate dependence of [] at low shear rates (D 2,9 10^–3 – 4,6 s^–1) higher molecular weights were determined than reported before. The so determined molecular weights of hyaluronic acid (HA) are in the range ofM = 2,4 10⁶–12 10⁶ g/mol in inflammatory joint diseases and in non-inflammatory in a range ofM = 3,1 10⁶–11,4 10⁶ g/mol and therefore both are lower than in normal synovial fluids,M = 16,2 10⁶ g/mol. Additional distinctk _SB-values indicate a compacter coil of HA macromolecules in the pathological synovial fluid.Quantitative statements can be made to the macromolecular network structure by determination of the critical HA polymer concentration and the number of entanglements resulting from prestationary measurements. The normal synovial fluid is comparable with a very narrow-meshed and dense network, but also in the pathological synovial fluid the HA is not isolated like in a particle solution but forms a much looser entanglement network structure. The reduction of viscosity and elasticity of isolated HA indicate that molecular interactions between HA and other macromolecules in synovial fluid are of essential significance for the rheological behaviour. Measurements of relaxation time at rest as well as shear conditions show a higher segment-flexibility of macromolecular network structure resulting from a reduction of HA-concentration and/or molecular weight in the pathological synovial fluid.

---

D Schergeschwindigkeit - D _K D-Koordinate des Tangentenschnittpunktes der Fließkurve - t _k Knickzeit (Relaxationszeit) - scheinbare Viskosität (Scherviskosität) - ₀ Anfangs-, Nullviskosität - _rel relative Viskosität - _sp spezifische Viskosität - _red = _sp/c reduzierte Viskosität - [] Grenzviskositätszahl (Staudinger-Index) - A Flächenmaß der master-curve - G ₀ Ruheschermodul - G _s Schermodul der strömenden Lösung - _N = ₁₁ – ₂₂ erste Normalspannungsdifferenz - ₁₂ Schubspannung - ₀ scheinbare Relaxationszeit der ruhenden Lösung - scheinbare Relaxationszeit der strömenden Lösung - v Verhängungszahl - M Molekulargewiht - s ₂₀ Sedimentationskoeffizient bezogen auf 20°C - s ₂₀ ⁰ Sedimentationskoeffizient bezogen auf 20°C bei unendlicher Verdünnung (c 0) - c Hyaluronsäure-Konzentration - c _krit kritische Hyaluronsäure-Konzentration - c _p Protein-Konzentration - _äqu Dichte äquivalenter Kugeln - Dichte des Mediums im Kegel/Platte-Spalt - Winkelgeschwindigkeit - R Radius von Kegel und Platte des Meßsystems - r Korrelationskoeffizient - p Signifikanzniveau Herrn Prof. Dr. med.Fritz Hartmann zu seinem 60. Geburtstag gewidmet.Mit 14 Abbildungen und 5 Tabellen     

Evaluation of the dynamic-mechanical properties of solids from a cooperative model for stress relaxation

For many solid materials the stress relaxation process obeys the universal relationF = – (d/d lnt)_max = (0.1 ± 0.01) ( ₀ — _i ), regardless of the structure of the material. Here denotes the stress,t the time, ₀ the initial stress of the experiment and _i the internal stress. A cooperative model accounting for the similarity in relaxation behaviour between different materials was developed earlier. Since this model has a spectral character, the concepts of linear viscoelasticity are used here to evaluate the corresponding prediction of the dynamic mechanical properties, i.e. the frequency dependence of the storageE () and lossE () moduli. Useful numerical approximations ofE () andE () are also evaluated. It is noted that the universal relation in stress relaxation had a counterpart in the frequency dependence ofE (). The theoretical prediction of the loss factor for high-density polyethylene is compared with experimental results. The agreement is good.     

New concepts in relative permeabilities at high capillary numbers for surfactant flooding

Flooding oil reservoirs with surfactant solutions can increase the amount of oil that can be recovered. Macroscopic modelling of the process requires relative permeabilities to be functions of saturation and capillary number. With only limited experimental data, relative permeabilities have usually been assumed to be linear functions of saturation at high capillary numbers. The experimental data is reviewed, some of which suggest that this assumption is not necessarily correct. The basis for the assumption is therefore reviewed and it is concluded that the linear model corresponds to microscopically segregated flow in the porous medium. Based on new but equally plausible complementary assumptions about the flow pattern, a mixed flow model is derived. These models are then shown to be limiting cases of a droplet model which represents the mixing scale within the porous medium and gives a physical basis for interpolating between the models. The models are based on physical concepts of flow in a porous medium and so the approach described here represents a significant improvement in the understanding of high capillary number flow. This is shown by the fact that fewer parameters are needed to describe experimental data.Notation A total cross-sectional area assigned to capillary bundle - A ⁽ⁱ⁾ physical cross-sectional area of tube i - c ⁽ⁱ⁾ ordered configurational label for droplets in tube i - c configuration label for tube i (order not considered) - D defined by Equation (26) - E(...) expectation value with respect to the trinomial distribution - S _r ⁽⁾ fractional flow of phase - k absolute permeability - k _r relative permeability of phase - k _r ⁰ endpoint relative permeability of phase - L capillary tube length in bundle model - m ⁽ⁱ⁾ number of droplets of phase a occupying tube i - n exponent for phase a in Equation (2) - N number of droplets in bundle model - N _c capillary number - p pressure - p(c') probability of configuration c - Q ⁽ⁱ⁾ total volume flow rate in tube i - S saturation of phase - S flowing saturation of phase - S _r residual saturation of phase - S _r ⁽⁾ saturations when fractional flow of phase is 1 in the case of varying residual saturations for three-phase flow ( ) - t _c residence time for droplet configuration c - v ⁽ⁱ⁾ total fluid velocity in bundle tube i - , phase label - p pressure differential across capillary bundle - ⁽ⁱ⁾ tube conductivity defined by Equation (7) - viscosity of phase - interfacial tension - gradient operator - ... average over tube droplet configurations     

An experimental study on the two-dimensional opposed wall jet in a turbulent boundary layer: Change in the flow pattern with velocity ratio

Results of the measurement of flow properties in a two-dimensional turbulent wall jet which is injected into the turbulent boundary layer in the direction opposite to that of the boundary layer flow are presented by varying the ratio of the jet issuing velocity to the mainstream velocity . This flow includes the flow separation and the recirculating flow, and so it requires the magnitude and direction of instantaneous velocity be measured. In the present work, a tandem hot-wire probe is manufactured and its characteristics are examined experimentally. With the use of this probe the change in the penetration length of the jet with the velocity ratio is investigated. It is clarified that two regimes of flow patterns exist: in the low velocity ratio the penetration length of the jet changes little with , and in the high velocity ratio it goes far from the nozzle with increasing . Streamlines, turbulence intensity contours and static pressure contours are presented in the two typical velocity ratios corresponding to each of two flow patterns, and they are compared.List of symbols b ₀ nozzle width (Fig. 1) - , e mean and fluctuating output voltages of hot-wire anemometer - p, p mean static pressure, p = p – p _o - p ₀ static pressure in undisturbed mainstream - p _w , p _w mean wall pressure, p _w = p _w – p _o - U ₀ mainstream velocity - U _j jet velocity at the nozzle exit - , u mean and fluctuating velocity components in x-direction - u _e effective cooling velocity - x distance along the wall from nozzle exit - x _pmax, x _pmin positions where the wall pressure has maximum and minimum values respectively - x _s penetration length of jet - y distance from the wall - forward flow fraction - ₁, ₂ yaw and pitch angles of flow direction (Fig. 4) - velocity ratio, = U _j /U _o - density of the fluid - nondimensional stream function The authors wish to express their appreciation to Prof. Toshio Tanaka of Gifu University for his advice in the course of the experiment. They also would like to thank the Research Foundation for the Electrotechnology of Chubu which partly supported this work.     

Dielectric properties of the composite system polyurethane — sodium chloride

Summary The dielectric properties of the composite system polyurethane-sodium chloride have been measured at frequencies between 10^–4 Hz and 3 · 10⁵ Hz in the temperature range from –150 °C up to +90 dgC. Three dielectric loss mechanisms have been found; they are indicated by ₁, ₂ and. The activation energy of the ₁-transition is 35 kcal/mole, that of the-transition 6.7 kcal/mole. The ₂-loss peak was only observed at frequencies of 10³ Hz and higher, forming one broad peak with the ₁-loss peak at lower frequencies. In the composite materials, the- and ₂-loss peaks measured at fixed frequencies were found at the same temperature. The ₂-loss peak, however, was shifted to a lower temperature, due to the sodium chloride filler. Comparison of experimental data of and tan with theoretical predictions concerning the dielectric properties of composite systems showed only partial agreement. The difference mainly consisted in. the temperature shift in the tan-peak of the ₁-transition.

---

Zusammenfassung Die dielektrischen Eigenschaften des Verbundssystems Kochsalz-Polyurethankautschuk wurden im Frequenzgebiet zwischen 10^–4 Hz und 3.10⁵ Hz und im Temperaturbereich von –150 °C bis +90 °C gemessen. Es wurden drei dielektrische Verlustmaxima gefunden, die mit ₁, ₂ und angedeutet werden. Die Aktivierungsenergie des ₁-Überganges beträgt 35 kcal/Mol, die des-Überganges 6.7 kcal/Mol. Das ₂-Maximum konnte nur bei Frequenzen höher als 10³Hz vom ₁-Maximum gesondert erfaßt werden. Die Lage der ₂- und-Maxima war vom Füllgrad unabhängig. Das ₁-Maximum verschiebt sich mit steigendem Füllgrad zu niedrigeren Temperaturen. Die gemessenen Werte von und tan stimmen nur teilweise mit den Aussagen einer Theorie der dielektrischen Eigenschaften von Mischkörpern überein.

     

Mixed convection flow in narrow vertical ducts

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 ₀ heat transfer on the wally=0 - Q ₁ heat transfer on the wally=1 - Re Reynolds number - R modified Reynolds number - T temperature of the fluid - T ₀ 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 - ₀ modified buoyancy parameter, ₀=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 - ₀ shear stress on the wally=0 - ₁ shear stress on the wally=1 - stream function     

Unsteady effects in a thin viscous shock layer near a three-dimensional stagnation point of a body moving with given acceleration and deceleration

An unsteady viscous shock layer near a stagnation point is studied. The Navier-Stokes equations are analyzed in the limit 1, Re₀ , df/dt = ^n-mF(t/^m). The Reynolds number Re₀ is defined in the paper by Eq. (1.3) (df/dt is the velocity of the body with respect to an inertial frame of reference moving with the original steady velocity –V'_t8, ₂ = ( – 1)/( + 1)). Various flow regimes in the case 1, l, n max(2m, m + 1), m 0, where ² = 1/Re₀ are analyzed. Equations are derived that generalize the asymptotic analysis to the case of a viscous unsteady flow of gas in a thin three-dimensional shock layer. The problem of a thin unsteady viscous shock layer near the stagnation point of a body with two curvatures is formulated. Examples of numerical solution are given for different ratios of the principal curvatures of the body, the wall temperature, the parameters of the original steady flow, and the acceleration and deceleration regimes.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 100–111, March–April, 1981.I thank Yu. D. Shevelev for a fruitful discussion of the work.     

Experimental analysis of scaling laws in low Re λgrid-generated turbulence

An experimental investigation of scaling laws in turbulence generated by a biplane grid for low Reynolds numbers (Re ) is presented. The present study covers a wide range of flow field conditions (from Re 2.5 to Re 36) that have not been analyzed from this point of view. It is shown that following the Kolmogorov theory a scaling range can not be observed since the separation between the energy production scales and the dissipative scales is too short. On the other hand, an extended form of scaling, the Extended Self-Similarity (ESS), permits the identification of a range of scales characterized by the same scaling exponent much wider than the one previously examined. Thus the scaling laws for the first six moments of the velocity structure function are accurately calculated by means of the ESS and an anomalous scaling with respect to the Kolmogorov theory is observed for Re down to the order of 10. As a matter of fact the scaling exponents are in good agreement with the ones that were determined at higher Re by previous experimental and numerical investigations. For Re 6 a regularization of the scaling exponents is observed as an effect of the dissipation. We also present an analysis of the universality properties of the ESS form of scaling by means of the form function and an analysis of the sensitivity of the scaling range to the Re .     

Investigation of supersonic viscous flow past a sphere in the presence of subsonic and sonic injection

Supersonic flow past a sphere with a given rate of gas injection along the generator is investigated numerically on the range Re=10²–10⁴. Calculations have been made on the interval 0 90°, where is the angle between the axis of symmetry and the normal to the surface. It is shown that for high subsonic and sonic injection rates it is possible to observe qualitatively new features in the flow structure and in the distribution of the local supersonic flow characteristics around the perimeter of the sphere not previously noted in [9]. In the case of sonic injection the changes in flow structure occur only in the supersonic zone. In the neighborhood of the transition from a subsonic to sonic injection velocity the heat flux has a local maximum, which in absolute value does not exceed the heat flux in the absence of injection. It is shown that there may be qualitative differences in the pressure distribution over the surface of the body with increase in the injection parameter depending on the distribution and value of the injected gas flow rate and, moreover, the number Re.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 83–89, January–February, 1988.     

Hypersonic flow over a narrow delta plate at large angles of attack

The problem of hypersonic flow over a flat delta plate with a high sweepback anglex at angles of attack close to /2 is solved using a numerical algorithm based on transition to the conical solution. The existence of conical flow at /2 with the velocity vector directed towards the apex of the plate is established. Values ofC _p/sin² and the thickness of the shock layer in the plane of symmetry of the plate are given as functions of the hypersonic similarity parameterk=tan tanx. A comparison of the calculated and experimental data shows that they are in good agreement.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.5, pp. 183–185, September–October, 1992.     

Stability characteristics of condensate films

The linear stability theory is used to study stability characteristics of laminar condensate film flow down an arbitrarily inclined wall. A critical Reynolds number exists above which disturbances will be amplified. The magnitude of the critical Reynolds number is in all practical situations so small that a laminar gravity-induced condensate film can be expected to be unstable. Several stabilizing effects are acting on the film flow; at an inclined wall these effects are due to surface tension, gravity and condensation mass transfer.

---

Zusammenfassung Mit Hilfe der linearen Stabilitätstheorie werden die Stabilitätseigenschaften laminarer Kondensatfilme an einer geneigten Wand untersucht. Es zeigt sich, daß Kondensatfilme in jedem praktischen Fall ein unstabiles Verhalten aufweisen. Der stabilisierende Einfluß von Oberflächenspannung, Schwerkraft und Stoffübertragung durch Kondensation bewkkt jedoch, daß Störungen in bestimmten Wellenlängenbereichen gedämpft werden.

---

Nomenclature c=c^*/u₀ complex wave velocity, celerity, dimensionless - c^*=c _r ^* + i c _i ^* complex wave velocity, celerity, dimensional - c_p specific heat at constant pressure - g gravitational acceleration - h_fg latent heat - k thermal conductivity of liquid - p^* pressure - p=p^*/u₀ ² dimensionless pressure - Pe=Pr Re^* Peclet number - Pr Prandtl number - Re^*=u₀ / Reynolds number (defined with surface velocity) - S temperature perturbation amplitude - t^* time - t=t^* u₀/ dimensionless time - T temperature - T_s saturation temperature - T_w wall temperature - T=T_s-T_w temperature drop across liquid film - u^*, v^* velocity components - u=u^*/u₀ dimensionless velocity components - v=v^*/u₀ dimensionless velocity components - u₀ surface velocity of undisturbed film flow - v _g ^* vapor velocity - x^*, y^* coordinates - x=x^*/ dimensionless coordinates - y=y^*/ dimensionless coordinates Greek Symbols =^* wave number, dimensionless - ^*=2 /^* wave number dimensional - ^* wave length, dimensional - =^*/ wave length, dimensionless - local thickness of undisturbed condensate film - kinematic viscosity, liquid - density, liquid - _g density vapor - surface tension - = (1 +) film thickness of disturbed film, Fig. 1 - stream function perturbation amplitude - angle of inclination Base flow quantities are denoted by^–, disturbance quantities are denoted by.     

Natural convection heat transfer through an enclosed horizontal layer of supercritical carbon dioxide

In natural convection heat transfer through a thin horizontal layer of carbon dioxide, maxima in the equivalent thermal conductivities are obtained in the vicinity of the respective pseudocritical temperatures at pressures of 75.8, 89.6 and 103.4 bar. The maxima are the more pronounced, the closer the critical point is approached.Comparison of experimental results with Nusselt equations shows good agreement except for the immediate vicinity of the pseudocritical temperature.In visual observations a distinct change in flow structure appears in the immediate vicinity of the pseudocritical temperature. A steady state polygon pattern and a boiling-like action could not be observed in this geometry.

---

Zusammenfassung Beim Wärmetransport durch freie Konvektion in einer dünnen waagerechten Schicht von Kohlendioxid ergaben sich Maxima der scheinbaren Wärmeleitfähigkeit in der Nähe der pseudokritischen Temperaturen bei Drükken von 75,8, 89,6 und 103,4 bar. Die Maxima sind um so ausgeprägter, je mehr man sich dem kritischen Punkt nähert.Ein Vergleich der Versuchsergebnisse mit Nusseltbeziehungen ergibt gute Übereinstimmung außer in unmittelbarer Umgebung der pseudokritischen Temperatur. Direkte Beobachtungen der Konvektionsmuster zeigen in unmittelbarer Umgebung der pseudokritischen Temperatur eine deutliche Strukturänderung. Ein stationäres Zellmuster und siedeähnliche Vorgänge konnten in dieser Anordnung nicht beobachtet werden.

---

Nomenclature A area of the heating or cooling plate - C constant in the correlation - g acceleration of gravity - h heat transfer coefficient - k thermal conductivity of fluid in the gap - k _e equivalent thermal conductivity - m, n exponents of dimensionless numbers - q heat flux - T _C,PC absolute temperature; critical C, pseudocritical PC - Gr Grashof numberg ( _h– _c) ³/ ² - Nu Nusselt numberh/k - Pr Prandtl number/ - thermal diffusivity - coefficient of volume expansion - width of gap - _c,h temperature of cooling (c)-, heating (h)-plate - _m arithmetic mean temperature ( _c+ _h)/2 - kinematic viscosity - _c,h fluid density at the temperature of the cooling (c)- or heating (h)-plate - heat flow rate through the gap     

Critique of the computational Preston tube method

The general validity of the Computational Preston Tube Method (CPM) proposed by Nitsche et al. is disputed. The method involves the determination of both the skin-friction coefficient c _f and the von Karman constant from fitting two near-wall mean velocity measurements to a generalized van Driest family of velocity profiles. It is demonstrated by a detailed examination of measured mean velocity profiles that for at least one flow situation for which success with the CPM has been claimed, the laminar-turbulent transition of a flat-plate boundary layer, the method must either fail completely or lead to non-unique results. It is concluded that the applicability of the CPM is restricted to flows for which the normal law-of-the-wall is valid.     

Rayleigh problem in slip flow with transverse magnetic field

An analytical continuum solution of the Rayleigh problem in slip flow with applied magnetic field is obtained using a modified initial condition and slip boundary conditions. The results are uniformly valid for all times and show that the velocity slip and the local skin friction coefficient remain almost unaffected by the imposition of the magnetic field for small times. They increase however with the magnetic field for large times. The present results reduce to the corresponding results of the hydrodynamic case when there is no magnetic field.Nomenclature A constant - b characteristic length - B magnetic field vector - B ₀ magntidue of the applied magnetic field normale to the plate - B _x magnitude of the induced magnetic field parallel to the plate - C slip coefficient, (2–f)/f - C _f skin friction coefficient, - C _D average drag coefficient - erfc(_x) complementary error function, - E electric field vector - f Maxwell's reflection coefficient - H _a Hartmann number, (B ₀ ² b ²/)^1/2 - nondimensional magnetic parameter - J current vector - Kn=L/b Knudsen number - L mean free path - M Mach number - p constant parameter - P _m magnetic Prandtl number, Re _m/Re= ₀ - q velocity vector - Re Reynolds number, Ub/ - Re _m magnetic Reynolds number, ₀ Ub - t time - nondimensional time, tU/b - u velocity of the fluid parallel to the plate - nondimensional velocity, u/U - U velocity of the plate - Laplace transform of - x, y coordinates along and normal to the plate respectively - y nondimensional distance, y/b - Z nondimensional parameter, 1/Re ^1/2 Kn - ratio of specific heats - boundary layer thickness - velocity slip - viscosity - ₀ magnetic permeability - kinematic viscosity - nondimensional time parameter, ( /Re)^1/2/Kn - density - electrical conductivity