Instationäre Messung der Wärmeleitfähigkeit mit optischer Registrierung

Zusammenfassung Die Wärmeleitfähigkeit von Wasser wird im Temperaturbereich von 20 bis 90 °C und bei 1 bar mit einem neuen instationären Absolutverfahren bestimmt. Zur Aufzeichnung des instationären Feldes des Brechungsindex und der Temperatur werden zwei interferometrische Anordnungen benutzt: Die Methode vonMach-Zehnder und das Biprismaverfahren. Die Ergebnisse stehen in guter Übereinstimmung zu den Messungen anderer Autoren, die stationäre Methoden benutzten. Die Unsicherheit der 374 Einzelmessungen wird auf höchstens±1% geschätzt. Damit ist nachgewiesen, daß ein instationäres Meßverfahren mit optischer Registrierung mit den klassischen stationären Verfahren hinsichtlich der Meßunsicherheit konkurrieren kann. Das instationäre Verfahren kommt ohne kalorische Messungen aus und besteht bei optischer Registrierung im wesentlichen aus Längenmessungen.

---

Unsteady-state measurements of the thermal conductivity with optical recording

The thermal conductivity of water in the temperature region from 20 to 90°C and at 1 bar was measured by means of a new unsteady-state absolute method. To record the unsteady-state field of the index of refraction and of the temperature, two interferometric arrangements were used: TheMach-Zehnder and the biprisma methods. The results are in good agreement with measurements of other authors, who had used steady-state methods. The maximum degree of uncertainty of the 374 measurements is estimated to be±1%. Thus it is shown that unsteady-state methods with optical recording can well be compared with classical steady-state methods regarding uncertainties. The method does not require caloric measurements and uses primarly determinations of lengths.

---

Bezeichnungen A=2 a in den Auswerteverfahren gebrauchte Abkürzung, m - B=q _x=0/ in den Auswerteverfahren gebrauchte Abkürzung, grd/m - F Fläche, m² - Fo=a/x ² Fourierzahl - I Strom, A - R Widerstand, - U Spannung, V - T Kelvintemperatur, °K - a=/c _p Temperaturleitfähigkeit, m²/s - Wärmeeindringzahl, Ws^1/2/m²grd - c _p isobare spezifische Wärmekapazität, kJ/kggrd - l Modellänge, m - n Brechungsindex - q Wärmestromdichte, W/m² - t Celsiustemperatur, °C - =1/ spezifisches Volumen, m³/kg - x Wandabstand, Ortskoordinate, m - z Ordnungszahl der Interferenzstreifen - =/_x=0 dimensionslose Übertemperatur - n=n -n _x=0 Brechungsindexdifferenz - Verhältnis der Wärmeeindringzahlen - Übertemperatur, grd - Wärmeleitfähigkeit, W/mgrd - Lichtwellenlänge, m - =1/ Dichte, kg/msu3 - Zeit, s Indices Zustand des Bades, Umgebung - x=0 Wand - z Stelle der z-ten Ordnung - i Laufparameter - err errechneter Wert - mess Meßwert - Bez Bezug Auszug aus der von der Fakultät für Maschinenwesen, und Elektrotechnik der Technischen Hochschule München genehmigten Dissertation von J.Bach.     

Transport in ordered and disordered porous media I: The cellular average and the use of weighting functions

In this work we consider transport in ordered and disordered porous media using singlephase flow in rigid porous mediaas an example. We defineorder anddisorder in terms of geometrical integrals that arise naturally in the method of volume averaging, and we show that dependent variables for ordered media must generally be defined in terms of thecellular average. The cellular average can be constructed by means of a weighting function, thus transport processes in both ordered and disordered media can be treated with a single theory based on weighted averages. Part I provides some basic ideas associated with ordered and disordered media, weighted averages, and the theory of distributions. In Part II a generalized averaging procedure is presented and in Part III the closure problem is developed and the theory is compared with experiment. Parts IV and V provide some geometrical results for computer generated porous media.Roman Letters A interfacial area of the- interface contained within the macroscopic region, m² - A_e area of entrances and exits for the-phase contained within the macroscopic system, m² - g gravity vector, m/s² - I unit tensor - K traditional Darcy's law permeability tensor, m² - L general characteristic length for volume averaged quantities, m - characteristic length (pore scale) for the-phase - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - N unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p0 reference pressure in the-phase, N/m² - p traditional intrinsic volume averaged pressure, N/m² - r0 radius of a spherical averaging volume, m - r position vector, m - r position vector locating points in the-phase, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - v velocity vector in the-phase, m/s - v traditional superficial volume averaged velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V/V, volume average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m²     

Flow visualization of compressible vortex structures using density gradient techniques

Mathematical results are derived for the schlieren and shadowgraph contrast variation due to the refraction of light rays passing through two-dimensional compressible vortices with viscous cores. Both standard and small-disturbance solutions are obtained. It is shown that schlieren and shadowgraph produce substantially different contrast profiles. Further, the shadowgraph contrast variation is shown to be very sensitive to the vortex velocity profile and is also dependent on the location of the peak peripheral velocity (viscous core radius). The computed results are compared to actual contrast measurements made for rotor tip vortices using the shadowgraph flow visualization technique. The work helps to clarify the relationships between the observed contrast and the structure of vortical structures in density gradient based flow visualization experiments.Nomenclature a Unobstructed height of schlieren light source in cutoff plane, m - c Blade chord, m - f Focal length of schlieren focusing mirror, m - C _T Rotor thrust coefficient, T/( ² R ⁴) - I Image screen illumination, Lm/m ² - l Distance from vortex to shadowgraph screen, m - n _b Number of blades - p Pressure,N/m ² - p Ambient pressure, N/m ² - r, , z Cylindrical coordinate system - r _c Vortex core radius, m - Non-dimensional radial coordinate, (r/r _c ) - R Rotor radius, m - Tangential velocity, m/s - Specific heat ratio of air - Circulation (strength of vortex), m ²/s - Non-dimensional quantity, ² 8²p r _c ² - Refractive index of fluid medium - ₀ Refractive index of fluid medium at reference conditions - Gladstone-Dale constant, m ³/kg - Density, kg/m ³ - Density at ambient conditions, kg/m ³ - Non-dimensional density, (/ ) - Rotor solidity, (n _b c/ R) - Rotor rotational frequency, rad/s     

Transport in ordered and disordered porous media II: Generalized volume averaging

In this paper we develop the averaged form of the Stokes equations in terms of weighting functions. The analysis clearly indicates at what point one must choose a media-specific weighting function in order to achieve spatially smoothed transport equations. The form of the weighting function that produces the cellular average is derived, and some important geometrical theorems are presented.Roman Letters A interfacial area of the- interface associated with the local closure problem, m² - A _e area of entrances and exits for the-phase contained within the averaging system, m² - A _p surface area of a particle, m² - d _p 6V _p/A_p, effective particle diameter, m - g gravity vector, m/s² - I unit tensor - K _m permeability tensor for the weighted average form of Darcy's law, m² - L general characteristic length for volume averaged quantities, m - L _p general characteristic length for volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - l characteristic length (pore scale) for the-phase - l _i i=1, 2, 3 lattice vectors, m - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - m _v special convolution product weighting function associated with the traditional averaging volume - m _g general convolution product weighting function - m _V unit cell convolution product weighting function - m _C special convolution product weighting function for ordered media which produces the cellular average - m _D special convolution product weighting function for disordered media - m _M master convolution product weighting function for ordered and disordered media - n unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p_m superficial weighted average pressure, N/m² - p _m intrinsic weighted average pressure, N/m² - p traditional intrinsic volume averaged pressure, N/m² - p p p _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function, m - r position vector, m - r position vector locating points in the-phase, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - V velocity vector in the-phase, m/s - vm superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - V volume of the-phase contained in the averaging volume, m³ - V _p volume of a particle, m³ - v traditional superficial volume averaged velocity, m/s - v v p _m spatial deviation velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V /V, volume average porosity - _m m * . weighted average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m² - V /V, volume fraction of the-phase     

Laminar source flow between parallel porous disks with equal suction and injection

An analysis is presented for laminar source flow between parallel stationary porous disks with suction at one of the disks and equal injection at the other. The solution is in the form of an infinite series expansion about the solution at infinite radius, and is valid for all suction and injection rates. Expressions for the velocity, pressure, and shear stress are presented and the effect of the cross flow is discussed.Nomenclature a distance between disks - A, B, ..., J functions of R _w only - F static pressure - p dimensionless static pressure, p(a ²/ ²) - Q volumetric flow rate of the source - r radial coordinate - r dimensionless radial coordinate, r/a - R radial coordinate of a point in the flow region - R dimensionless radial coordinate of a point in the flow region, R - Re source Reynolds number, Q/2a - R _w wall Reynolds number, Va/ - reduced Reynolds number, Re/r ² - critical Reynolds number - velocity component in radial direction - u dimensionless velocity component in radial direction, a/ - average radial velocity, Q/2a - u dimensionless average radial velocity, Re/r - ratio of radial velocity to average radial velocity, u/u - velocity component in axial direction - v dimensionless velocity component in axial direction, v - V magnitude of suction or injection velocity - z axial coordinate - z dimensionless axial coordinate, z a - viscosity - density - kinematic viscosity, / - shear stress at lower disk - shear stress at upper disk - ₀ dimensionless shear stress at lower disk, - ₁ dimensionless shear stress at upper disk, - dimensionless stream function     

Mixed convection from a horizontal plate to fluids of any Prandtl number

Laminar mixed convection over a horizontal plate with uniform wall temperature or uniform wall heat flux is analyzed by introducing proper buoyancy parameters and transformation variables for fluids of any Prandtl number between 0.001 and 10,000. Both cases of buoyancy assisting and opposing flow conditions are investigated. For the buoyancy-assisting case, the obtained numerical results are very accurate over the entire range of mixed convection intensity from pure forced convection limit to pure free convection limit. For the buoyancy-opposing case, solutions are obtained from the forced convection limit to the point of breakdown.

---

Mischkonvektion an einer horizontalen Platte für Fluide mit beliebiger Prandtl-Zahl

Zusammenfassung Es wurde laminare Mischkonvektion an einer horizontalen Platte mit einheitlicher Wandtemperatur oder einheitlicher Wandwärmestromdichte bei Einführung zweckmäßiger Auftriebsparameter und Transformationsvariablen für Fluide mit beliebiger Prandtl-Zahl zwischen 0,001 und 10 000 untersucht. Es wurden die Fälle der Strömung entgegen und in Richtung der Auftriebskraft untersucht. Für den Fall der Strömung in Richtung der Auftriebskraft wurden sehr genaue numerische Ergebnisse für den gesamten Bereich der gemischten Konvektion von rein erzwungener Konvektion bis zu rein freier Konvektion erhalten. Für den Fall der Strömung entgegen der Auftriebsrichtung wurden Lösungen für erzwungene Konvektion bis zum Umkehrpunkt erhalten.

---

Nomenclature C _f local friction coefficient - f reduced stream function - g gravitational acceleration - Gr local Grashof number for UWT,g (T _w –T )x ³/ ² - Gr* local Grashof number for UHF,g q _w x ⁴/k ² - m =10 for UWT; and =6 for UHF - n =5 for UWT; and =3 for UHF - Nu local Nusselt number - p pressure - Pr Prandtl number,/ - q _w wall heat flux - Ra local Rayleigh number for UWT,Gr Pr - Ra* local Rayleigh number for UHF,Gr*Pr - Re local Reynolds number,u x/ - T fluid temperature - T _w wall temperature - T free-stream temperature - u velocity component inx-direction - u free-stream velocity - v velocity component iny-direction - x coordinate parallel to the plate - y coordinate normal to the plate Greek symbols thermal diffusivity - thermal expansion coefficient - =0 for UWT; and =1 for UHF - buoyancy parameter, =( Ra)^1/5/( Re)^1/2 for UWT; and =( Ra*)^1/6/( Re)^1/2 for UHF - pseudo-similarity variable, (y/x) - dimensionless temperature, =(T–T )/(T _w –T ) for UWT; and =(T–T )/(q _w x/k) for UHF - =[( Re)^1/2+( Ra)^1/5] for UWT; and =[( Re)^1/2+( Ra*)^1/6] for UHF - dynamic viscosity - kinematic viscosity - /(1+) - dimensionless pressure - density - Pr/(1+Pr) - _w wall shear stress,(u/y) _y=0 - stream function - Pr/(1+Pr)^1/3     

Wärmeübergang bei erzwungener turbulenter Strömung in Ringspalten und örtlich beliebig veränderlichen Randbedingungen zweiter Art

Zusammenfassung Das Problem des Wärmeübergangs bei turbulenter Strömung in konzentrischen Ringspalten wird für den dreidimensionalen Fall theoretisch gelöst, wobei die Wandwärmestromdichte sowohl in azimutaler als auch in axialer Richtung beliebig variiert. Die Lösung der Energiegleichung erfolgt mit der klassischen Methode der Superposition und Trennung der Variablen, wobei das dabei auftretende Sturm-Liouvillesche Eigenwertproblem numerisch gelöst wird. Zur Lösung werden Verteilungen für die Geschwindigkeit und anisotropen turbulenten Austauschgrößen verwendet, die mit dem phänomenlogischen Turbulenzmodell von Ramm berechnet wurden. Ergebnisse werden über einen weiten Bereich der Reynolds-Zahl (10⁴ Re 10⁶), der Prandtl-Zahl (0 Pr 100) und für verschiedene Radienverhältnisse diskutiert.

---

Turbulent forced convection heat transfer in annuli with arbitrarily varying boundary conditions of second kind

The problem of turbulent flow heat transfer in concentric annuli is analysed for the general threedimensional case in which the wall heat flux varies arbitrarily in both the circumferential and axial directions. The energy equation is solved using the classical method of superposition and separating variables, where the resulting Sturm-Liouville problem are evaluated numerically. The solution is based on velocity profiles and anisotropic thermal turbulent transport properties evaluated by Ramm's phenomenological turbulence model. Results are discussed over a wide range of Reynolds number (10⁴ Re 10⁶), Prandtl number (0 Pr 100) and radius ratio.

---

Bezeichnungen a,b Fourierkoeffizienten - B geometrische Funktion, [s(1-r) + r]/(1–s) - C Koeffizienten - D hydraulischer Durchmesser, 2(r₂ – r₁) - E Energietransportfunktion - f axiale Wärmestromdichteverteilung - F azimutale Wärmestromdichteverteilung - g radiale Temperaturfunktion - l Kanallänge - L dimensionslose Kanallänge, 1/D - M axialer Temperaturgradient im thermisch ausgebildeten Bereich - n harmonischer Parameter - Nu Nusselt-Zahl - Pe Péclet-Zahl - Pr Prandtl-Zahl - q Wärmestromdichte - Q dimensionslose Wärmestromdichte, q/q₀ - r dimensionslose radiale Koordinate, (R-r₁)/(r₂-r₁) - r₁,r₂ innerer und äußerer Ringspaltradius - R radiale Koordinate - Re Reynolds-Zahl - s Ringspaltverhältnis, r₁/r₂ - T dimensionslose Temperatur, 2· · (-_E/(D· q₀ - u dimensionslose Geschwindigkeit, U/U_m - U Geschwindigkeit - x dimensionslose axiale Koordinate, X/D - X axiale Koordinate - Wärmeübergangskoeffizient - _un modifizierter Eigenwert - halber Segmentwinkel - turbulente Austauschgröe - Temperatur - dimensionslose Temperaturdifferenz, T - T_m - Wärmeleitfähigkeit - _un Eigenwerte - kinematische Viskosität - azimutale Koordinate - Eigenfunktionen Indizes e thermischer Einlauf - E Eintritt bei x=0 - H Wärme - i Bedingung an der i-ten benetzten Oberfläche (i=1 – Innenrohr, i=2 - Außenrohr) - j Bedingung, wenn nur an der j-ten Oberfläche des Ringspaltes die Wärme übertragen wird (j=1,2) - ij Bedingung an der i-ten Oberfläche, wenn nur an der j-ten Oberfläche des Ringspaltes die Wärme übertragen wird (ij=11, 12, 22, 21) - m mittel - n Ordnung der Harmonischen - r radiale Richtung - u Ordnung des Eigenwertproblems - azimutale Richtung - 0 umfangskonstant - thermisch ausgebildet     

Flow through a helically coiled annulus

In this paper the flow is studied of an incompressible viscous fluid through a helically coiled annulus, the torsion of its centre line taken into account. It has been shown that the torsion affects the secondary flow and contributes to the azimuthal component of velocity around the centre line. The symmetry of the secondary flow streamlines in the absence of torsion, is destroyed in its presence. Some stream lines penetrate from the upper half to the lower half, and if is further increased, a complete circulation around the centre line is obtained at low values of for all Reynolds numbers for which the analysis of this paper is valid, being the ratio of the torsion of the centre line to its curvature.Nomenclature A =constant - a outer radius of the annulus - b unit binormal vector to C - C helical centre line of the pipe - D rL - g 1000 - K Dean number=Re² - L 1+r sin - M (L ²+ ² r ²)^1/2 - n unit normal vector to C - P, P pressure and nondimensional pressure - p ₀, p pressures of O(1) and O() - Re Reynolds number=aW ₀/ - (r, , s), (r, , s) coordinates and nondimensional coordinates - nonorthogonal unit vectors along the coordinate directions - r ₀ radius of the projection of C - t unit tangent vector to C - V _r, V , V _s velocity components along the nonorthogonal directions - V_r, V, V _s nondimensional velocity components along - W ₀ average velocity in a straight annulus Greek symbols , curvature and nondimensional curvature of C - U, V, W lowest order terms for small in the velocity components along the orthogonal directions t - _r, , _s first approximations to V _r , V, V _s for small - =/=/ - kinematic viscosity - density of the fluid - , torsion and nondimensional torsion of C - , stream function and nondimensional stream function - nondimensional streamfunction for U, V - a inner radius of the annulus After this paper was accepted for publication, a paper entitled On the low-Reynolds number flow in a helical pipe, by C.Y. Wang, has appeared in J. Fluid. Mech., Vol 108, 1981, pp. 185–194. The results in Wangs paper are particular cases of this paper for =0, and are also contained in [9].     

Convective heat transfer within fibrous insulation slabs

A numerical study of convective heat flow within a fibrous insulating slab is presented. The material is treated as an anisotropic porous medium and the variation of properties with temperature is taken into account. Good agreement is obtained with available experimental data for the same geometry.

---

Zusammenfassung Für den konvektiven Wärmestrom in einem faserförmigen Isolierstoff wird eine numerische Berechnung angegeben. Der Stoff wird als anisotropes poröses Medium mit temperaturabhängigen Stoffwerten angesehen. Die Übereinstimmung mit verfügbaren Versuchswerten ist gut.

---

Nomenclature C_p specific heat of the gas at the mean temperature - Da Darcy number=k_y/H² - Gr^* modified Grashof number=gTHk_y/²= (Grashof number) × (Darcy number) - H thickness of the specimen - P gas pressure - Pr^* modified Prandtl number= C_p/_x - Ra^* modified Rayleigh number=Gr^* Pr^* - R_p ratio of permeabilities=k_y/k_x - R_k ratio of conductivities= _y/_x - T absolute temperature of the gas - t₁ absolute temperature of the hot face - T₂ absolute temperature of the cold face - T_m mean temperature of the gas=(T₁+T₂)/2 - k_x specific permeability of the porous medium along the x-direction - k_y specific permeability of the porous medium along the y-direction - p T/T_m - q exponent - r exponent - u gas velocity along the x-direction - v gas velocity along the y-direction - X_* distance along the x-direction - y_* distance along the y-direction - T temperature difference=t₁–T₂ - thermal coefficient of expansion of the gas - _m thermal coefficient of expansion of the gas at the mean temperature - _* T–T_m - dimensionless temperature= ^*/T - _a apparent thermal conductivity of the porous medium along the x-direction - _al local apparent thermal conductivity of the porous medium along the x-direction - _x thermal conductivity of the porous medium along the x-direction in the absence of convection - _y thermal conductivity of the porous medium along the y-direction in the absence of convection - dynamic viscosity of the gas - _m dynamic viscosity of the gas at the mean temperature - kinematic viscosity of the gas - _m kinematic viscosity of the gas at the mean temperature - density of the gas - _m density of the gas at the mean temperature - ^* stream function at any point - dimensionless stream function= ^*/( _m/_m)     

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     

Zur Variationsformulierung nichtlinearer Randwertprobleme

Übersicht Es werden verschiedene Bedingungen aufgestellt, die es erlauben, die durch die beiden (Systeme von) nichtlinearen DifferentialgleichungenA (u, ) = q, B (u, ) = und Randbedingungen zusammen mit den nichtlinearen algebraischen Relationenq = C(u, ), = D(u, ) beschriebene Aufgabe durch äquivalente Variationsprobleme zu ersetzen. Dabei zeigt sich ein enger Zusammenhang mit den in der Festkörpermechanik wohlbekannten Prinzipien der virtuellen Verschiebungen und der virtuellen Kräfte. Die auf systematischem Weg konstruierten Variationsfunktionale enthalten viele in der Physik bekannte Funktionale als Sonderfälle, insbesondere jene, die in der Elastomechanik nach Green, Castigliano, Hellinger, Reißner, Hu und Washizu benannt werden.

---

Summary In this paper there are established various conditions which allow a variational formulation of the problem described by the two (systems of) nonlinear differential equationsA(u, ) = q, B(u, ) = and boundary conditions together with the nonlinear algebraic relationsq = C(u, ), = D(u, ). Besides a close relationship is revealed to the principles of virtual displacements and virtual forces which are wellknown in solid mechanics. The systematically constructed variational functional contain many functionals in physics as special cases, mainly those of Green, Castigliano, Hellinger, Reißner, Hu and Washizu in elastomechanics.

     

Untersuchung des lokalen Stoffüberganges am laminaren und turbulenten Rieselfilm

Zusammenfassung Der lokale Stoffübergang wurde in Abhängigkeit von der Meßlänge, dem Startort und der Zulaufhöhe gemessen. Der Gültigkeitsbereich der Theorie von Nusselt wird ermittelt. Die Reynolds-Zahl nahm Werte zwischen 3,86 und 2496 an. Die örtlich wirkende Hydrodynamik ist entscheidend für das Anwachsen der örtlichen Sherwood-Zahl. Die Genauigkeit aller Versuchsergebnisse kann auf ± 5% abgeschätzt werden.

---

Investigation of the local mass transfer of a laminar and turbulent falling liquid film

The local mass transfer was measured as a function of the measuring length, the starting point and the liquid height above the ring-slot. The range of the Reynolds number was 3,86 Re 2496. The validity of the Nusselt theory and the range of it is shown. The local hydrodynamic is the most important factor of the increase of the local Sherwood number. The accuracy of the measurements is ± 5%.

---

Bezeichnungen a Temperaturleitfähigkeit m²/s=/(c_p) - c Konzentration, c=¯c + c kmol/m³ - c_i0 Konzentration im Flüssigkeitskern kmol/m³ - D Diffusionskoeffizient m²/s - EL-NR Elektrodennummer - Fa Faraday-Konstante A s/kgäq=96,5·10⁶ - g Erdbeschleunigung m/s² - i_G Grenzstromdichte A/m² - u Geschwindigkeit in x-Richtung, u= + u - U Umfang des Rohres m - v Geschwindigkeit in y-Rich- m/stung, v=¯v + v - V^* Volumenstrom m³/s - x Lauflänge, Koordinate in m Strömungsrichtung - x_M Meßlänge für den Stoff-Übergang m - x_ST Startort für den Stoff-Übergang m - y Wegkoordinate senkrecht zur Rohroberfläche m - z Wertigkeit der Elektro-denreaktion kgäq/kmol - ZH Zulaufhöhe m - Wärmeübergangskoeffizient W/m²C - Stoffübergangskoeffizient m/s - Filmdicke m - Wärmeleitfähigkeit W/(mC) - kinematische Viskosität m²/s - Re=u/=V^*/U Reynolds-Zahl - Pr=/a=c_p/ Prandtl-Zahl - Sc=/D Schmidt-Zahl - Nu= / Nusselt-Zahl - Sh= /D Sherwood-Zahl - SHL lokale Sherwood-Zahl - SHM mittlere Sherwood-Zahl - - zeitlich gemittelt - örtlich gemittelt Die Durchführung der Arbeit am Institut für Verfahrens — und Kältetechnik der ETH Zürich bei Prof. Dr. P. Grassmann wurde ermöglicht durch Zuschüsse der Kommission zur Förderung der wissenschaftlichen Forschung und meiner Eltern.     

Mean and turbulent velocity measurements of supersonic mixing layers

The behavior of supersonic mixing layers under three conditions has been examined by schlieren photography and laser Doppler velocimetry. In the schlieren photographs, some large-scale, repetitive patterns were observed within the mixing layer; however, these structures do not appear to dominate the mixing layer character under the present flow conditions. It was found that higher levels of secondary freestream turbulence did not increase the peak turbulence intensity observed within the mixing layer, but slightly increased the growth rate. Higher levels of freestream turbulence also reduced the axial distance required for development of the mean velocity. At higher convective Mach numbers, the mixing layer growth rate was found to be smaller than that of an incompressible mixing layer at the same velocity and freestream density ratio. The increase in convective Mach number also caused a decrease in the turbulence intensity ( _u/U).List of symbols a speed of sound - b total mixing layer thickness between U ₁ – 0.1 U and U ₂ + 0.1 U - f normalized third moment of u-velocity, f u³/(U)³ - g normalized triple product of u² , g u²/(U)³ - h normalized triple product of u ², h u ²/(U)³ - l _u axial distance for similarity in the mean velocity - l _u axial distance for similarity in the turbulence intensity - M Mach number - M _c convective Mach number (for ₁ = ₂), M _c (U ₁ – U ₂)/(a ₁ + a ₂) - P static pressure - r freestream velocity ratio, r U ₂/U ₁ - Re unit Reynolds number, Re U/ - s freestream density ratio, s ₂/₁ - T _t total temperature - u instantaneous streamwise velocity - u deviation of u-velocity, uu – U - U local mean streamwise velocity - U ₁ primary freestream velocity - U ₂ secondary freestream velocity - average of freestream velocities, (U ₁ + U ₂)/2 - U freestream velocity difference, U U ₁ – U ₂ - instantaneous transverse velocity - v deviation of -velocity, – V - V local mean transverse velocity - x streamwise coordinate - y transverse coordinate - y ₀ transverse location of the mixing layer centerline - ensemble average - ratio of specific heats - boundary layer thickness (y-location at 99.5% of free-stream velocity) - similarity coordinate, (y – y ₀)/b - compressible boundary layer momentum thickness - viscosity - density - standard deviation - dimensionless velocity, (U – U ₂)/U - 1 primary stream - 2 secondary stream A version of this paper was presented at the 11th Symposium on Turbulence, October 17–19, 1988, University of Missouri-Rolla     

Free convection heat transfer from an isothermal plate with arbitrary inclination

This paper presents a new formulation for the laminar free convection from an arbitrarily inclined isothermal plate to fluids of any Prandtl number between 0.001 and infinity. A novel inclination parameter is proposed such that all cases of the horizontal, inclined and vertical plates can be described by a single set of transformed equations. Moreover, the self-similar equations for the limiting cases of the horizontal and vertical plates are recovered from the transformed equations by setting=0 and=1, respectively. Heated upward-facing plates with positive and negative inclination angles are investigated. A very accurate correlation equation of the local Nusselt number is developed for arbitrary inclination angle and for 0.001 Pr .

---

Wärmeübertragung bei freier Konvektion an einer isothermen Platte mit beliebiger Neigung

Zusammenfasssung Diese Untersuchung stellt eine neue Formulierung der laminaren freien Konvektion von Flüssigkeiten mit einer Prandtl-Zahl zwischen 0,001 und unendlich an einer beliebig schräggestellten isothermen Platte dar. Ein neuer Neigungsparameter wird eingeführt, so daß alle Fälle der horizontalen, geneigten oder vertikalen Platte von einem einzigen Satz transformierter Gleichungen beschrieben werden können. Die unabhängigen Gleichungen für die beiden Fälle der horizontalen and vertikalen Platte wurden für=0 und=1 aus den transformierten Gleichungen wieder abgeleitet. Es wurden erwärmte aufwärtsgerichtete Platten mit positiven und negativen Neigungswinkeln untersucht. Eine sehr genaue Gleichung wurde für die lokale Nusselt-Zahl bei beliebigen Neigungswinkeln und für 0,001 Pr entwickelt.

---

Nomenclature C _p specific heat - f reduced stream function - g gravitational acceleration - Gr local Grashof number,g(T _w –T _w ) x³/v² - h local heat transfer coefficient - k thermal conductivity - n constant exponent - Nu local Nusselt number,hx/k - p pressure - Pr Prandtl number, v/ - Ra local Rayleigh number,g(T _w –T )J x³/v - T fluid temperature - T _w wall temperature - T temperature of ambient fluid - u velocity component in x-direction - v velocity component in y-direction - x coordinate parallel to the plate - y coordinate normal to the plate Greek symbols thermal diffusivity - thermal expansion coefficient - (Ra¦sin¦)^1/4/( Ra cos()^1/5 - pseudo-similarity variable, (y/) - dimensionless temperature, (T–T )/(T _w–T ) - ( Ra cos)^1/5+(Rasin)^1/4 - v kinematic viscosity - 1/[1 +(Ra cos)^1/5/( Ra¦sin)^1/4] - density of fluid - Pr/(1+Pr) - _w wall shear stress - angle of plate inclination measured from the horizontal - stream function - dimensionless dynamic pressure     

The effect of non-steady wall shear on pressure surges in conduits carrying viscous liquids

Summary A study is made of the attenuation of pressure surges in a two-dimension a channel carrying a viscous liquid when a valve at the downstream end is suddenly closed. The analysis differs from previous work in this area by taking into account the transient nature of the wall shear, which in the past has been assumed as equivalent to that existing in steady flow. For large values of the frictional resistance parameter the transient wall shear analysis results in less attenuation than given by the steady wall shear assumption.Nomenclature c /, ft/sec - e base of natural logarithms - F(x, y) integration function, equation (38) - (x) mean value of F(x, y) - g local acceleration of gravity, ft/sec² - h width of conduit, ft - k (2m–1)^{2 2} L/h ² c, equation (35) - k* 12L/h ² c, frictional resistance parameter, equation (46) - L length of conduit, ft - m positive integer - n positive integer - p pressure, lb/ft² - p ₀ constant pressure at inlet of conduit, lb/ft² - P pressure plus elevation head, p+gz, equation (4) - mean value of P over the conduit width h - P ₀ p ₀+gz ₀, lbs/ft² - R frictional resistance coefficient for steady state wall shear, lb sec/ft⁴ - s positive integer; also, condensation, equation (6) - t time, sec - t ct/L, dimensionless time - u x component of fluid velocity, ft/sec - u _m mean velocity in conduit, equation (12), ft/sec - u ₀(y) velocity profile in Poiseuille flow, equation (19), ft/sec - transformed velocity - U initial mean velocity in conduit - x distance along conduit, measured from valve (fig. 1), ft - x x/L, dimensionless distance - y distance normal to conduit wall (fig. 1), ft - y y/h, equation (25) - z elevation, measured from arbitrary datum, ft - z ₀ elevation of constant pressure source, ft - isothermal bulk compression modulus, lbs/ft² - _n , equation (37) - _n (2n–1)/2, equation (36) - viscosity, slugs/ft sec - / = kinematic viscosity, ft²/sec - density of fluid, slugs/ft³ - ₀ density of undisturbed fluid, slugs/ft³ - ø angle between conduit and vertical (fig. 1) The research upon which this paper is based was supported by a grant from the National Science Foundation.     

Reservoir transport equations by compositional approach

The objective of this paper is to present an overview of the fundamental equations governing transport phenomena in compressible reservoirs. A general mathematical model is presented for important thermo-mechanical processes operative in a reservoir. Such a formulation includes equations governing multiphase fluid (gas-water-hydrocarbon) flow, energy transport, and reservoir skeleton deformation. The model allows phase changes due to gas solubility. Furthermore, Terzaghi's concept of effective stress and stress-strain relations are incorporated into the general model. The functional relations among various model parameters which cause the nonlinearity of the system of equations are explained within the context of reservoir engineering principles. Simplified equations and appropriate boundary conditions have also been presented for various cases. It has been demonstrated that various well-known equations such as Jacob, Terzaghi, Buckley-Leverett, Richards, solute transport, black-oil, and Biot equations are simplifications of the compositional model.Notation List B reservoir thickness - B formation volume factor of phase - Cⁱ mass of component i dissolved per total volume of solution - C ⁱ mass fraction of component i in phase - C heat capacity of phase at constant volume - C_p heat capacity of phase at constant pressure - D ⁱ hydrodynamic dispersion coefficient of component i in phase - D_MTf thermal liquid diffusivity for fluid f - F = F(x, y, z, t) defines the boundary surface - f_p fractional flow of phase - g gravitational acceleration - H_p enthalpy per unit mass of phase - J_p volumetric flux of phase - k_rf relative permeability to fluid f - k₀ absolute permeability of the medium - M_p ⁱ mass of component i in phase - n porosity - N rate of accretion - P_f pressure in fluid f - p_ca capillary pressure between phases and =p-p - Rⁱ rate of mass transfer of component i from phase to phase - Rⁱ _source source rate of component i within phase - S saturation of phase - s gas solubility - T temperature - t time - U displacement vector - u velocity in the x-direction - v velocity in the y-direction - V volume of phase - V_s velocity of soil solids - W_i body force in coordinate direction i - x horizontal coordinate - z vertical coordinate Greek Letters _p volumetric coefficient of compressibility - _T volumetric coefficient of thermal expansion - _ij Kronecker delta - volumetric strain - _m thermal conductivity of the whole matrix - internal energy per unit mass of phase - gf suction head - density of phase - _ij tensor of total stresses - _ij tensor of effective stresses - volumetric content of phase - f viscosity of fluid f     

Uniqueness and the maximum principle for quasilinear elliptic equations with quadratic growth conditions

In this paper, we show that the maximum principle holds for quasilinear elliptic equations with quadratic growth under general structure conditions.Two typical particular cases of our results are the following. On one hand, we prove that the equation (1) {ie77-01} where {ie77-02} and {ie77-03} satisfies the maximum principle for solutions in H ¹()L(), i.e., that two solutions u ₁, u ₂H¹() L() of (1) such that u ₁u₂ on , satisfy u ₁u₂ in . This implies in particular the uniqueness of the solution of (1) in H ₀ ¹ ()L().On the other hand, we prove that the equation (2) {ie77-04} where fH^–1() and g(u)>0, g(0)=0, satisfies the maximum principle for solutions uH¹() such that g(u)¦Du|{²L¹(). Again this implies the uniqueness of the solution of (2) in the class uH ₀ ¹ () with g(u)¦Du|{²L¹().In both cases, the method of proof consists in making a certain change of function u=(v) in equation (1) or (2), and in proving that the transformed equation, which is of the form (3) {ie77-05}satisfies a certain structure condition, which using ((v₁ -v ₂)⁺)ⁿ for some n>0 as a test function, allows us to prove the maximum principle.     

Über den Einfluß von Querkrümmung und Dichte bei inhomogenen turbulenten Freistrahlen

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 - c_i Massenkonzentration der Komponente i - c_D, c_L, c, c₁, c₂ 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 - k_i Schrittweite in -Richtung - K dimensionslose Turbulenzenergie - L turbulentes Längenmaß - M_i Molmasse der Komponente i - p Druck - allgemeine Gaskonstante - r Querkoordinate - r_0,5 Halbwertsbreite der Geschwindigkeit - r_0,5c Halbwertsbreite der Konzentration - T Temperatur - u Geschwindigkeitskomponente in x-Richtung - v Geschwindigkeitskomponente in r-Richtung - x Längskoordinate - y allgemeine Funktion - Y_i 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     

Shear viscosity behavior of highly concentrated suspensions at low and high shear-rates

A method is proposed for determining the shear viscosity behavior of highly concentrated suspensions at low and high shear-rates through the use of a formulation that is a function of three parameters signifying the effects of particle size distribution. These parameters are the intrinsic viscosity [], a parametern that reflects the level of particle association at the initiation of motion and the maximum packing concentration _m. The formulation reduces to the modified Eilers equation withn = 2 for high shear rates. An analytical method was used for the calculation of maximum packing concentration which was subsequently correlated with the experimental values to account for the surface induced interaction of particles with the fluid. The calculated values of viscosities at low and high shear-rates were found to be in good agreement with various experimental data reported in literature. A brief discussion is also offered on the reliability of the methods of measuring the maximum packing concentration. _r = /₀ relative viscosity of the suspension - volumetric concentration of solids - k _n coefficient which characterizes a specific effect of particle interactions - _m maximum packing concentration - _r,0 relative viscosity at low shear-rates - [] intrinsic viscosity - n, n parameter that reflects the level of particle interactions at low and high shear-rates, respectively - _r, relative viscosity at high shear-rates - (_m)_s, (_m)_i, (_m)_l packing factors for small, intermediate and large diameter classes - v _s, v_i, v_l volume fractions of small, intermediate and large diameter classes, respectively - _si, _sl coefficient to be used in relating a smaller to an intermediate and larger particle group, respectively - _is, _il coefficient to be used in relating an intermediate to a smaller and larger particle group, respectively - _ls, _li coefficient to be used in relating a larger to a smaller and intermediate particle group, respectively - _m0 maximum packing concentration for binary mixtures - _m,e measured maximum packing concentration - _m,c calculated maximum packing concentration     

On laminar flow through a uniformly porous pipe

Numerous investigations ([1] and [4–9]) have been made of laminar flow in a uniformly porous circular pipe with constant suction or injection applied at the wall. The object of this paper is to give a complete analysis of the numerical and theoretical solutions of this problem. It is shown that two solutions exist for all values of injection as well as the dual solutions for suction which had been noted by previous investigators. Analytical solutions are derived for large suction and injection; for large suction a viscous layer occurs at the wall while for large injection one solution has a viscous layer at the centre of the channel and the other has no viscous layer anywhere. Approximate analytic solutions are also given for small values of suction and injection.

Nomenclature

General r distance measured radially - z distance measured along axis of pipe - u velocity component in direction of z increasing - v velocity component in direction of r increasing - p pressure - density - coefficient of kinematic viscosity - a radius of pipe - V velocity of suction at the wall - r ²/a ² - R wall or suction Reynolds number, Va/ - f() similarity function defined in (6) - u ₀() eigensolution - U(0) a velocity at z=0 - K an arbitrary constant - B _K Bernoulli numbers Particular Section 5 perturbation parameter, –2/R - ² a constant, –K - x / - g(x) f()/ Section 6 perturbation parameter, –R/2 - ² a constant, –K - g() f() - g _c ()=g() near centre of pipe - * point where g()=0 Section 7 2/R - ² K - t (1–)/ - w(t, ) [1–f(t)]/ - ₀, ₁ constants - g() f()– ₀ - ₀/ - ₀ a constant - * point where f()=0