The Thermodynamic Significance of the Local Volume Averaged Temperature

Since the temperature is not an additive function, the traditional thermodynamic point of view suggests that the volume integral of the temperature has no precise physical meaning. This observation conflicts with the customary analysis of non-isothermal catalytic reactors, heat pipes, driers, geothermal processes, etc., in which the volume averaged temperature plays a crucial role. In this paper we identify the thermodynamic significance of the volume averaged temperature in terms of a simple two-phase heat transfer process. Given the internal energy as a function of the point temperature and the density

Flow of a simple non-newtonian fluid past a sphere

Summary Creeping flow past a sphere is solved for a limiting case of fluid behaviour: an abrupt change in viscosity.List of Symbols d _ij Component of rate-of-deformation tensor - F _d Drag force exerted on sphere by fluid - G ^(d) Coefficients in expression for _ij in terms of d _ij - G _YOJK ^(d) Coefficients in power series representing G ^(d) - R Radius of sphere - r Spherical coordinate - V Velocity of fluid very far from sphere - v _i Component of the velocity vector - x Dimensionless radial distance, r/R - x _i Rectangular Cartesian coordinate - Dimensionless quantity defined by (26) - ^(d) Potential defined by (7) - Value of x denoting border between Regions 1 and 2 as a function of - ₁, ₂ Lower and upper limiting viscosities defined by (10) - Spherical coordinate - * Value of for which =1 - Value of denoting border between regions 1 and 2 as a function of x - Newtonian viscosity - _ij Component of the stress tensor - Spherical coordinate - ₁, ₂ Stream functions defined by (12) and (14) - Second and third invariants of the stress tensor and of the rate-of-deformation tensor, defined by (3)     

Einfluß der Beheizungsart auf den Wärmeübergang im horizontalen Verdampferrohr

Zusammenfassung Durch die Phasenverteilung von Flüssigkeit und Dampf ergeben sich im horizontalen Verdampferrohr i.a. unterschiedliche lokale Wärmeübergangskoeffizienten (WÜK) am Rohrumfang. Die beiden technisch vorkommenden Randbedingungen — konstante Wärmestromdichte und konstante WandemperaturT _W am Rohrinnendurchmesser — können unterschiedliche Auswirkungen auf umfangsgemittelte WÜK haben. Während bei konstanter aufgeprägter Wärmestromdichte ein interner Transport der zu übertragenden Wärmemenge durch die Wärmeleitung der Wand erfolgen kann, ist dies beiT _W =konst. nicht möglich. Die Ergebnisse zeigen, daß sich der umfangsgemittelte WÜK bei konstanter Wärmestromdichte gegenüber dem Wert bei konstanter Wandtemperatur verringert. Dies ist jedoch nur dann der Fall, wenn der Rohrumfang unvollständig benetzt ist. Dabei hängt die Verminderung des umfangsgemittelten WÜK stark vom Wärmeleitvermögen der Wand, s, ab ( Wärmeleitfähigkeit,s Wandstärke). Es wurde festgestellt, daß die Verminderung sowohol aus der höheren Wandtemperatur am Rohrscheitel als auch aus der Änderung des benetzten Umfangs resultiert, da dann der Anteil des Rohrumfangs mit schlechtem Wärmeübergang vergrößert wird.

---

Effect of thermal boundary conditions on heat transfer in a horizontal evaporator tube

As a result of the non uniform phase distribution of liquid and vapour local heat transfer coefficients vary considerably along the circumference in a horizontal evaporator tube. It was expected that both modes of heating, uniform heat flux ( = const) and uniform inside wall temperature (T _W =const), lead to different circumferential averaged heat transfer coefficients (HTC). The experiments show that the circumferential averaged HTC is lower in the case of uniform heat flux compared to the value of uniform wall temperature if the perimeter of the tube is only partly wetted. However, the reduction is affected by circumferential heat conduction. This reduction is a result of both, the higher wall temperature on the top of the tube and the reduction of the wetted perimeter.

---

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

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.     

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     

First order phase equilibrium for classical bodies

Summary In this note the Author gives a general and unified treatment of first order phase equilibria for classical bodies like those considered by Truesdell and Toupin in[3].The Author reaches a system of partial differential equations (generalized Clapeyron equations) the conditions of whose solution are shown always to be satisfied.In particular, the Author derives the equations governing the polarized phase equilibrium for a fluid.Besides the equations ruling the phase equilibrium for a two phase n-component fluid mixture are given and the equivalence with the statical Gibbs-Duhem relation is shown.

---

Sommario In questa nota l'Autore presenta una trattazione generale ed unificata degli equilibri di fase del primo ordine per corpi classici come quelli definiti da Truesdell e Toupin in[3].L'Autore perviene ad un sistema di equazioni alle derivate parziali (equazioni di Clapeyron generalizzate) del quale si dimostra la integrabilità.In particolare, si deducono le equazioni che governano gli equilibri di fase polarizzati.Inoltre si ottengono le equazioni che regolano l'equilibrio di fase per una miscela fluida a n componenti; in questo caso si dimostra l'equivalenza delle equazioni con la relazione statica di Gibbs-Duhem.

---

---

This work was supported by the Gruppo Nazionale per la Fisica Matematica of C.N.R.     

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     

A three-parameter model describing the experimental relation between shear stress and shear rate for laminar flow of aqueous polymer solutions

Summary A three-parameter model is introduced to describe the shear rate — shear stress relation for dilute aqueous solutions of polyacrylamide (Separan AP-30) or polyethylenoxide (Polyox WSR-301) in the concentration range 50 wppm – 10,000 wppm. Solutions of both polymers show for a similar rheological behaviour. This behaviour can be described by an equation having three parameters i.e. zero-shear viscosity ₀, infinite-shear viscosity , and yield stress ₀, each depending on the polymer concentration. A good agreement is found between the values calculated with this three-parameter model and the experimental results obtained with a cone-and-plate rheogoniometer and those determined with a capillary-tube rheometer.

---

Zusammenfassung Der Zusammenhang zwischen Schubspannung und Schergeschwindigkeit von strukturviskosen Flüssigkeiten wird durch ein Modell mit drei Parametern beschrieben. Mit verdünnten wäßrigen Polyacrylamid-(Separan AP-30) sowie Polyäthylenoxidlösungen (Polyox WSR-301) wird das Modell experimentell geprüft. Beide Polymerlösungen zeigen im untersuchten Schergeschwindigkeitsbereich von ein ähnliches rheologisches Verhalten. Dieses Verhalten kann mit drei konzentrationsabhängigen Größen, nämlich einer Null-Viskosität ₀, einer Grenz-Viskosität und einer Fließgrenze ₀ beschrieben werden. Die Ergebnisse von Experimenten mit einem Kegel-Platte-Rheogoniometer sowie einem Kapillarviskosimeter sind in guter Übereinstimmung mit den Werten, die mit dem Drei-Parameter-Modell berechnet worden sind.

---

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

An application of nonsymmetric Lax-Milgram lemma to nonassociated plasticity

Based on a general assumption for plastic potential and yield surface, some properties of the nonassociated plasticity are studied, and the existence and uniqueness of the distribution of incremental stress and displacement for work-hardening materials are proved by using nonsymmetric Lax-Milgram lemma, when the work-hardening parameter A>F/Q/–F/, Q/.     

The theory of a vibrating-rod densimeter

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

Transport in ordered and disordered porous media III: Closure and comparison between theory and experiment

In this paper we examine the closure problem associated with the volume averaged form of the Stokes equations presented in Part II. For both ordered and disordered porous media, we make use of a spatially periodic model of a porous medium. Under these circumstances the closure problem, in terms of theclosure variables, is independent of the weighting functions used in the spatial smoothing process. Comparison between theory and experiment suggests that the geometrical characteristics of the unit cell dominate the calculated value of the Darcy's law permeability tensor, whereas the periodic conditions required for thelocal form of the closure problem play only a minor role.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² - A interfacial area of the- interface associated with the local closure problem, m² - A _p surface area of a particle, m² - b vector used to represent the pressure deviation, m^–1 - B ⁰ B+I, a second order tensor that maps v _m ontov - B second-order tensor used to represent the velocity deviation - d _p 6V _p/A_p, effective particle diameter, m - d a vector related to the pressure, m - D a second-order tensor related to the velocity, m² - g gravity vector, m/s² - I unit tensor - K traditional Darcy's law permeability tensor calculated on the basis of a spatially periodic model, m² - 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 characteristic length for the volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - characteristic length (pore scale) for the-phase - _i i=1, 2, 3 lattice vectors, m - weighting function - m(-y) , convolution product weighting function - 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 - 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 _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function - r position vector, m - r position vector locating points in the-phase, m. - V averaging volume, m³ - B 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 _m superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - v traditional superficial volume averaged velocity, m/s - v – v _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²     

Transient laminar free convection in horizontal cylinders

A numerical study of laminar natural convection inside uniformly heated, partially or fully filled horizontal cylinders is made. A coordinate transformation which simplifies the discretization of the equations of motion and energy is utilized. The resulting system of partial differential equations with their boundary conditions is solved using central differences for various Prandtl and Grashof numbers for two different grid sizes. The flow in completely filled cylinders for which experimental data are available is predicted. Close agreement between steady-state predictions and experiments is obtained for temperature and velocity profiles as well as for the streamline contours and isotherms. The technique is further demonstrated by solving the transient natural convection flow inside a partially filled horizontal cylinder with an adiabatic free surface and subjected to uniform wall heating.

---

Laminare freie Konvektion in horizontalen Zylindern

Zusammenfassung Es wurde eine numerische Berechnung der laminaren, freien Konvektion in gleichmäßig beheizten, teilweise oder ganz gefüllten, horizontalen Zylindern durchgeführt. Dabei wird eine Koordinatentransformation benützt, welche die Diskretisierung der Bewegungs- und der Energiegleichung vereinfacht. Das so resultierende System von partiellen Differentialgleichungen wird, zusammen mit seinen Randbedingungen, unter Verwendung einer Differenzenmethode für verschiedene Prandtl und Grashof-Zahlen sowie für zwei verschiedene Gittergrößen gelöst. Für den vollständig gefüllten Zylinder, für den experimentelle Daten verfügbar sind, wird die Strömung vorhergesagt. Dabei wird für stationäre Zustände gute Übereinstimmung zwischen Rechnung und Experiment erzielt. Dies gilt sowohl für den Verlauf der Stromlinien als auch für den der Isothermen. Das Verfahren wird weiterhin am Beispiel der Berechnung instationärer, freier Konvektion in einem partiell gefüllten, horizontalen Zylinder demonstriert, wobei eine adiabate, freie Oberfläche und gleichmäßige Beheizung der Wand angenommen sind.

---

Nomenclature g acceleration due to gravity, m/s² - Gr _R ^* modified Grashof number =gqR⁴/kv² - Gr _R Grashof number =gTR³/v² - H heat function vector, dimensionless - k thermal conductivity, W/mK - L(Y) cord length associated with coordinateY, dimensionless - Pr Prandtl number=v/ - q wall heat flux, W/m² - R radius, m - r(X, Y,Z) distance of a boundary point from the reference axis, dimensionless - S vector derived from the flow field solution, dimensionless - T temperature, K - T _w wall temperature, K - T reference temperature, K - t time, s - u, v velocity components inx, y directions, m/s - U, V dimensionless velocity components inX- and Y-direction normalized withU - U reference velocity=gqR²/k or gTR, m/s - V velocity vector, dimensionless - W vorticity vector, dimensionless - W vorticity, dimensionless - x, y, z cartesian coordinates, m - X, Y, Z cartesian coordinates normalized with a reference length, dimensionless Greek letters thermal diffusivity, m²/s - coefficient of thermal expansion, K^–1 - ,,, non-dimensional coordinates in the transformed domain - non-dimensional temperature =(T–T)_k/q^R or T–T/T_w–T - v kinematic viscosity, m²/s - non-dimensional time=v/R² Gr_Rt or v/R² G _R ^* t - angle measured from the bottom of the cylinder, rads - ^* angle measured from the axis on (– ) plane, rads - heat potential, dimensionless - angle of incidence of the heat flux vector, rads - non-dimensional stream function - vector potential, dimensionless - grid size, dimensionless - ² Laplacian operator - gradient vector     

Hyperbolic phenomena in a strongly degenerate parabolic equation

We consider the equation u _t =((u) (u _x )) _x , where >0 and where is a strictly increasing function with lim _s = <. We solve the associated Cauchy problem for an increasing initial function, and discuss to what extent the solution behaves qualitatively like solutions of the first-order conservation law u _t = ((u)) _x . Equations of this type arise, for example, in the theory of phase transitions where the corresponding free-energy functional has a linear growth rate with respect to the gradient.     

Heat transfer from laminar flow in ducts with elliptic cross-section

Summary The cooling of a hot fluid in laminar Newtonian flow through cooled elliptic tubes has been calculated theoretically. Numerical data have been computed for the two values 1.25 and 4 of the axial ratio of the elliptic cross-section . For =1.25 the influence of non-zero thermal resistance between outmost fluid layer and isothermal surroundings has also been investigated. Special attention has been given to the distribution of heat flux around the perimeter; when increases the flux varies more with the position at the circumference. This positional dependence becomes less pronounced, however, as the (position-independent) thermal resistance of the wall increases.Flattening of the conduit, while maintaining its cross-sectional area constant, improves the cooling. Comparison with rectangular pipes shows that this improvement is not as marked with elliptic as with rectangular pipes.Nomenclature A _k =A _{m, n} coefficients of expansion (6) - a, b half-axes of ellipse, b<a - a _p =a _{r, s} coefficients of representation (V) - D hydraulic diameter, = 4S/P; S = cross-sectional area, P = perimeter - D _e equivalent diameter, according to (13) - n coordinate (outward) normal to the tube wall - T temperature of fluid - T _i temperature of fluid at the inlet - T _s temperature of surroundings - v ₀ mean velocity of fluid - v _z longitudinal velocity of fluid - x, y carthesian coordinates coinciding with axes of ellipse - z coordinate in flow direction - , dimensionless half-axes of ellipse, =a/D and =b/D - _t heat transfer coefficient from fluid at bulk temperature to surroundings; equation (11) - _w heat transfer coefficient at the wall; equation (3) - axial ratio of ellipse, = a/b = / - , , , dimensionless coordinates; =x/D, =y/D, =z/D, =n/D - dimensionless temperature, = (T–T _s)/(T _i–T _s) - ₀ cup-mixing mean value of ; equation (10) - thermal conductivity of fluid - _m,n = _k eigenvalue - c volumetric heat capacity of fluid - _{m, n} = _k = _k eigenfunction; equations (6) and (I) - Nu total Nusselt number, = _t D/ - Nusselt number at large distance from the inlet - Nu _w wall Nusselt number, = _w D/, based on _w - Pé Péclet number, = ₀ Dc/     

Instability of a plane axisymmetric flow with a shear discontinuity in the velocity profile

It is proposed to investigate the stability of a plane axisymmetric flow with an angular velocity profile (r) such that the angular velocity is constant when r < r_O – L and r > r_O + L but varies monotonically from ₁ to ₂ near the point r_O, the thickness of the transition zone being small L r_O, whereas the change in velocity is not small ¦₂ – ₁¦ ₂, ₁. Obviously, as L O short-wave disturbances with respect to the azimuthal coordinate (k=m/r_O 1/r_O) will be unstable with a growth rate-close to the Kelvin—Helmholtz growth rate. In the case L=O (i.e., for a profile with a shear-discontinuity) we find the instability growth rate _O and show that where the thickness of the discontinuity L is finite (but small) the growth rate does not differ from _O up to terms proportional to kL 1 and 1/m 1. Using this example it is possible to investigate the effect of rotation on the flow stability. It is important to note that stabilization (or destabilization) of the flow in question by rotation occurs only for three-dimensional or axisymmetric perturbations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 111–114, January–February, 1985.     

On the continuum modelling of porous media containing fluid: a molecular viewpoint with particular attention to scale

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

Numerische Berechnung und Simulation des instationären Stofftransports bei wandernden Phasengrenzflächen in Zweistoffsystemen

Zusammenfassung Der instationäre Stofftransport vor wandernden Phasengrenzflächen flüssig/fest wurde für Zweistoffsysteme unter Berücksichtigung allgemeiner Kristallisationsbedingungen und einfacher Körpergeometrien abgeleitet. Dies ermöglichen spezielle Transformationsbeziehungen sowie die Randwertbestimmung für die Aufschmelzfront. Ein neuartiges, für Hybridrechenanlagen entwickeltes Zerlegungsverfahren lieferte die analytischen Hilfsmittel für numerisch außerordentlich stabile und konvergente Lösungen des allgemeinen, linearen Randwertproblems. Damit ließ sich ein Rechenprogramm erstellen, welches die Simulation weitgehend beliebiger Kristallisationsprozesse erlaubt. Am Beispiel des Zonenschmelzens und der Normalkristallisation werden einige, derart simulierte Kristallisationsprozesse diskutiert.

---

Numeric calculation and simulation of nonsteady masstransfer by moving phaseboundaries in binary alloys

Non-Steady mass transfer in front of moving phase boundaries liquid/solid in binary alloys was derived. This was allowed by means of transformation relations and boundary-values for the melting-front. A newly analyzing-method primarily developed for serial hybrid computer integration made it possible to obtain very steady and convergent solutions of the general linear boundary-value-problem. These solutions allow to create computer-programs for the simulation of arbitrary processes of crystallization. For example processes of non-steady zone melting and normal freezing were simulated and discussed.

---

Bezeichnungen f_e Erstarrungsfrontgeschwindigkeit - f_a Aufschmelzfrontgeschwindigkeit - s Schmelzzonenlänge - s zeitliche Änderung der Schmelzzonenlänge - K₀ Gleichgewicht-Verteilungskoeffizient - K_eff effektiver Verteilungskoeffizient - k Konzentration - k₀ Ausgangskonzentration in flüssiger Phase - k_0B Ausgangskonzentration in fester Phase - ¯k₀ Konzentration in fester Phase vor der Aufschmelzfront - ¯k_B Konzentration in fester Phase nach der Kristallisation - D Diffusionskoeffizient - m Masse - m Massenstromdichte - t Zeit - T absolute Temperatur - vektorielle Strömungsgeschwindigkeit der Legierungskomponente - v_x,v_r axiale und radiale Strömungsgeschwindigkeit (z.B. durch Schwerkraftseigerung) - x,r Ortskoordinaten - L Stablänge - L_F vorwärtsstabiler Operator - L_B rückwärtsstabiler Operator - S Speicherfunktion - R Anfangsradius - , mitwandernde Ortskoordinaten - Länge des kristallisierten Volumens - Gewichtsfaktor - Wurzel der charakteristischen Gleichung - a,b,c, d,e,,, } Hilfsfunktionen - Grenzschichtdicke Auszug aus der vom Fachbereich für Maschinenwesen der Technischen Universität München zur Erlangung des akademischen Grades eines Doktor-Ingenieurs genehmigten Dissertation Instationärer und quasistationärer Stoffransport in binär legierten Schmelzen vor wandernden Phasengrenzflächen des Dipl.-Ing. Richard Bung. Berichterstatter Privatdozent Dr.-Ing. H.M. Tensi undProf. Dr.-Ing. F. Fischer; Tag der Einreichung 2.2.1976, Tag der Annahme 20.2.1976, Tag der Promotion 2.3.1976.     

Experimental and numerical analysis of oscillatory tube flow of viscoelastic fluids represented at the example of human blood

Summary The 4 constant parameters of an Oldroyd type constitutive equation for normal human blood of 45% haematocrit were determined by means of the steady flow curve and the material functions () and (), measured at 2 Hz in an oscillatory capillary viscometer. It was found, as other authors did before, that the phenomenological behaviour of blood flow can be reproduced qualitatively thereby. The quantitative behaviour, however, cannot be described by thus developed parameters. The parameters of the constitutive equationµ ₀, ₁ and ₂ were therefore generalized to become dependent of shear rate and frequency respectively.In itself this is nothing but a transformation of the material functions ( ),() and (), but these can be used as parameters in a constitutive equation though having lost the property of constancy.In this way the linear region of viscoelasticity and the steady flow curve can be reproduced quantitatively. A computer simulation of oscillatory flow for large amplitudes shows another tendency for the phase shift between pressure and flow than the experiment in the oscillatory capillary rheometer does.The applicability of a constitutive equation modified in this manner for other than oscillatory flow should be further examined especially for pulsatile flow.

---

Zusammenfassung Mit Hilfe der Fließkurve für die stationäre Strömung und Messungen der Materialfunktion () im Oszillations-Kapillarrheometer bei 2 Hz wurden die Konstanten der 4-Konstanten-Oldroyd-Stoffgleichung für gesundes Humanblut von 45% Hämatokrit bestimmt. Es zeigte sich, daß sich — wie auch schon von anderen Autoren mitgeteilt wurde — mit den so ermittelten Modellkonstanten die Phänomenologie des Fließverhaltens des Blutes qualitativ gut beschreiben läßt. Zur quantitativen Beschreibung reicht dieses Modell jedoch nicht aus, wie man an der Wiedergabe der stationären Fließkurve und der Materialfunktionen der linearen Viskoelastizität erkennt. Aus diesem Grund wurden die Parameterµ ₀, ₁ und ₂ in Abhängigkeit der Schergeschwindigkeit bzw. der Frequenz angesetzt.Dies bedeutet zunächst nichts anderes als eine Transformation der Materialfunktionen, jedoch in einer Art, daß letztere als Parameter in einer Stoffgleichung verwendet werden können, was allerdings mit dem Verlust der Konstanz der Parameter verbunden ist. Mit einer derart modifizierten Stoffgleichung lassen sich der Bereich der linearen Viskoelastizität und die stationäre Fließkurve quantitativ beschreiben. Eine Computersimulation der oszillierenden Rohrströmung zeigt für große Amplituden eine andere Tendenz für die Phasenverschiebung zwischen Druck und Volumenstrom, als sie sich bei Messungen am Oszillations-Kapillarrheometer ergibt.Die Anwendbarkeit der modifizierten Stoffgleichung für andere Strömungsformen, wie z.B. pulsierende Rohrströmungen, muß noch geprüft werden.

---

---

Paper, presented at the Annual Conference of the Deutsche Rheologische Gesellschaft in Berlin, May 8–10, 1978.

---

With 8 figures     

LDV measurements of a turbulent air-solid two-phase flow in a 90° bend

Simultaneous measurements of the mean streamwise and radial velocities and the associated Reynolds stresses were made in an air-solid two-phase flow in a square sectioned (10×10 cm) 90° vertical to horizontal bend using laser Doppler velocimetry. The gas phase measurements were performed in the absence of solid particles. The radius ratio of the bend was 1.76. The results are presented for two different Reynolds numbers, 2.2×10⁵ and 3.47×10⁵, corresponding to mass ratios of 1.5×10^–4 and 9.5×10^–5, respectively. Glass spheres 50 and 100 m in diameter were employed to represent the solid phase. The measurements of the gas and solid phase were performed separately. The streamwise velocity profiles for the gas and the solids crossed over near the outer wall with the solids having the higher speed near the wall. The solid velocity profiles were quite flat. Higher negative slip velocities are observed for the 100 m particles than those for the 50 gm particles. At angular displacement =0°, the radial velocity is directed towards the inner wall for both the 50 and 100 m particles. At =30° and 45°, particle wall collisions cause a clear change in the radial velocity of the solids in the region close to the outer wall. The 100 m particle trajectories are very close to being straight lines. Most of the particle wall collisions occur between the =30° and 60° stations. The level of turbulence of the solids was higher than that of the air.List of symbols D hydraulic diameter (100 mm) - De Dean number,De = - mass flow rate - number of particles per second (detected by the probe volume) - r radial coordinate direction - r _i radius of curvature of the inner wall - r ₀ radius of curvature of the outer wall - r ^* normalized radial coordinate, - R mean radius of curvature - Re Reynolds number, - R _r radius ratio, - U ,U _z mean streamwise velocity - U _r ,U _y mean radial velocity - U _b bulk velocity - , _z rms fluctuating streamwise velocity - _r , _y rms fluctuating radial velocity - -r shear stress component - z-y shear stress component - x spanwise coordinate direction - x ^* normalized spanwise coordinate, - y radial coordinate direction in straight ducts - y ^* normalized radial coordinate in straight ducts, - z streamwise coordinate direction in straight ducts - z ^* normalized streamwise coordinate in straight ducts, Greek symbols streamwise coordinate direction - kinematic viscosity of air