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²     

Thermal effects in the capillary rheometer

Summary The effect of viscous heating in a capillary rheometer is analysed for a power-law fluid by means of a perturbation expansion based upon a boundary-layer-core structure. This expansion is found to complement the eigenfunction series solution obtained by earlier investigators. A similar analysis is presented for the work-of-expansion effect. These two thermal effects are superimposed together with a third perturbation effect due to the pressure dependence of viscosity.On the basis of the present theory, earlier work in this area is discussed and, in some cases, apparent inaccuracies or inconsistencies are pointed out. A means is indicated for correcting data on the basis of the present theory.

---

Zusammenfassung Es wird der Effekt der Erwärmung einer Potenzflüssigkeit infolge viskoser Reibung in einem Kapillar-Rheometer mittels einer Störungsrechnung untersucht, die auf der Unterteilung der Strömung in eine Grenzschicht und einen Kern basiert. Diese Störungsentwicklung ergänzt eine früher von anderen Autoren gefundene Reihenentwicklung mit Hilfe von Eigenfunktionen. Eine ähnliche Untersuchung wird für die thermische Ausdehnungsarbeit durchgeführt. Diese beiden thermischen Effekte sind zusammen einem dritten Störeffekt superponiert, der von der Druckabhängigkeit der Viskosität herrührt.Aufgrund der vorgelegten Theorie werden verschiedene auf diesem Gebiet früher durchgeführte Arbeiten diskutiert, und es werden in einigen Fällen offensichtliche Ungenauigkeiten und Folgewidrigkeiten aufgedeckt. Schließlich wird eine Methode zur Korrektur von Meßdaten mit Hilfe der vorliegenden Theorie angegeben.

---

Nomenclature a tube radius - b ; evaluated atT ₀ andp = 0 when used in perturbation expansion - C _p specific heat - f - f ^* - h defined by eq. [15] - k thermal conductivity - L tube length - m defined by eq. [8] - m ₀ m(T₀, 0) - n power-law index - p pressure - Pe C _p W a/k Peclet number - Pr C _pa/k Prandtl number - Q volumetric flow rate - Q ₀ unperturbed value ofQ in specified-p formulation - r radial coordinate - Re W a/ _a Reynolds number - T temperature - T ₀ inlet temperature - u radial velocity component - u ₀ 0 unperturbed radial velocity - w axial velocity component - w ₀ /W(1 – ) unperturbed axial velocity - W Q/(a ²) average axial velocity - W ₀ Q ₀/(a ²) - z axial coordinate - (3n + 1)/n - ^* ; evaluated atT ₀ andp = 0 when used in perturbation expansion - 4^1-n - ^* - (n + 1)/n - ... shear rate - 4W/a apparent shear rate - p total pressure drop - T _a W ²/k characteristic temperature difference - T _b total bulk-temperature rise - ^* T - r/a - shear viscosity - _a m₀ - (1 –)/ ^1/3 - —p/z - ₀ ... unperturbed value of - z-averaged value of - µ ^{n + 1}/n - z/(a Pe) - _L L/(a Pe) - mass density - _w shear stress at wall - streamfunction - ^*T₀ (absolute temperature scale) - ( )₁ leading-order effect due to viscous heating - ( ) ₁ ^* leading-order effect due to work-of-expansion Note: in specified-p formulation,W gets replaced byW ₀ in definition of Pe, Re, and. With 7 figures and 7 tables     

Flow in porous media I: A theoretical derivation of Darcy's law

Stokes flow through a rigid porous medium is analyzed in terms of the method of volume averaging. The traditional averaging procedure leads to an equation of motion and a continuity equation expressed in terms of the volume-averaged pressure and velocity. The equation of motion contains integrals involving spatial deviations of the pressure and velocity, the Brinkman correction, and other lower-order terms. The analysis clearly indicates why the Brinkman correction should not be used to accommodate ano slip condition at an interface between a porous medium and a bounding solid surface.The presence of spatial deviations of the pressure and velocity in the volume-averaged equations of motion gives rise to aclosure problem, and representations for the spatial deviations are derived that lead to Darcy's law. The theoretical development is not restricted to either homogeneous or spatially periodic porous media; however, the problem ofabrupt changes in the structure of a porous medium is not considered.Roman Letters A interfacial area of the - interface contained within the macroscopic system, m² - A _e area of entrances and exits for the -phase contained within the macroscopic system, m² - A interfacial area of the - interface contained within the averaging volume, m² - A ^* interfacial area of the - interface contained within a unit cell, m² - Ae area of entrances and exits for the -phase contained within a unit cell, m² - B second order tensor used to represent the velocity deviation (see Equation (3.30)) - b vector used to represent the pressure deviation (see Equation (3.31)), m^–1 - d distance between two points at which the pressure is measured, m - g gravity vector, m/s² - K Darcy's law permeability tensor, m² - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the -phase (see Figure 2), m - characteristic length scale for the -phase (see Figure 2), m - n unit normal vector pointing from the -phase toward the -phase (n =–n ) - n _e unit normal vector for the entrances and exits of the -phase contained within a unit cell - p pressure in the -phase, N/m² - p intrinsic phase average pressure for the -phase, N/m² - p –p , spatial deviation of the pressure in the -phase, N/m² - r ₀ radius of the averaging volume and radius of a capillary tube, m - v velocity vector for the -phase, m/s - v phase average velocity vector for the -phase, m/s - v intrinsic phase average velocity vector for the -phase, m/s - v –v , spatial deviation of the velocity vector for the -phase, m/s - V averaging volume, m³ - V volume of the -phase contained within the averaging volume, m³ Greek Letters V/V, volume fraction of the -phase - mass density of the -phase, kg/m³ - viscosity of the -phase, Nt/m² - arbitrary function used in the representation of the velocity deviation (see Equations (3.11) and (B1)), m/s - arbitrary function used in the representation of the pressure deviation (see Equations (3.12) and (B2)), s^–1     

LDA measurements in low Reynolds number turbulent boundary layer

LDA measurements of the mean velocity in a low Reynolds number turbulent boundary layer allow a direct estimate of the friction velocity U from the value of /y at the wall. The trend of the Reynolds number dependence of / is similar to the direct numerical simulations of Spalart (1988).     

Analytical evaluation of stress-strain test data

An analytical model for deducing the actual stress-strain properties from laboratory test results is discussed. As an illustration, an elastic bilinear material is used for unconfined cylindrical compression test conditions, as simulated with a finite element analysis. The results obtained are applicable for assisting in evaluating measured strength and stiffness properties of some clay soils, concrete test cylinders, concrete cores, and rock cores.The quantitative results of this study can be used for interpreting measured stress-strain data for unconfined compression test conditions. The error in measured results is shown to be influenced by Poisson's ratio, length-to-diameter ratio of the specimen, end condition, and ratio of inelastic modulus to initial elastic modulus. Curves for adjusting the measured results to the theoretical results are presented.Nomenclature D specimen diameter - E _i initial elastic stiffness modulus - E _y elastic stiffness modulus beyond the yield stress, plastic or inelastic modulus - L specimen length - axial strain - _av average strain - _g gage length strain - _y yield strain - Poisson's ratio - compressive stress - _av average stress - _t theoretical compressive stress - _y yield stress - _ym measured stress at the yield strain     

The effect of background particulates on the delayed transition of a heated 9:1 ellipsoid

An experimental study was done to quantify the effects of a variety of background particulates on the delayed laminar-turbulent transition of a thermally stabilized boundary layer in water. A Laser-Doppler Velocimeter system was used to measure the location of boundary layer transition on a 50 mm diameter, 9:1 fineness ratio ellipsoid. The ellipsoid had a 0.15 m RMS surface finish. Boundary layer transition locations were determined for length Reynolds numbers ranging from 3.0 × 10⁶ to 7.5 × 106. The ellipsoid was tested in three different heating conditions in water seeded with particles of four distinct size ranges. For each level of boundary layer heating, measurements of transition were made for clean water and subsequently, water seeded with 12.5 m, 38.9 m, 85.5 m and 123.2 m particles, alternately. The three surface heating conditions tested were no heating, T = 10°C and T = 15°C where T is the difference between the inlet model heating water temperature, T _i, and free stream water temperature, T . The effects of particle concentration were studied for 85.5 m and 123.2 m particulates.The results of the study can be summarized as follows. The 12.5 m and 38.9 m particles has no measurable effect on transition for any of the test conditions. However, transition was significantly affected by the 85.5 m and 123.2 m particles. Above a length Reynolds number of 4 × 10⁶ the boundary layer transition location moved forward on the body due to the effect of the 85.5 m particles for all heating conditions. The largest percentage changes in transition location from clean water, were observed for 85.5 m particles seeded water.Transition measurements made with varied concentrations of background particulates indicated that the effect of the 85.5 m particles on the transition of the model reached a plateau between 2.65 particulates/ml concentration and 4.2 particles/ml. Measurements made with 123.3 m particles at concentrations up to 0.3 part/ml indicated no similar plateau.     

The problems of nonlinear bending for orthotropic rectangular plate with four clamped edges

（黄家寅）（秦圣立）ＴＨＥＰＲＯＢＬＥＭＳＯＦＮＯＮＬＩＮＥＡＲＢＥＮＤＩＮＧＦＯＲＯＲＴＨＯＴＲＯＰＩＣＲＥＣＴＡＮＧＵＬＡＲＰＬＡＴＥＷＩＴＨＦＯＵＲＣＬＡＭＰＥＤＥＤＧＥＳ￥ＨｕａｎｇＪｉａｙｉｎ;ＱｉｎＳｈｅｎｇｌｉ（ＱｕｆｕＮｏｒｍａｌＵｎ...     

The rheology of polymeric liquid crystals

The molecular theory of Doi has been used as a framework to characterize the rheological behavior of polymeric liquid crystals at the low deformation rates for which it was derived, and an appropriate extension for high deformation rates is presented. The essential physics behind the Doi formulation has, however, been retained in its entirety. The resulting four-parameter equation enables prediction of the shearing behavior at low and high deformation rates, of the stress in extensional flows, of the isotropic-anisotropic phase transition and of the molecular orientation. Extensional data over nearly three decades of elongation rate (10^–2–10¹) and shearing data over six decades of shear rate (10^–2–10⁴) have been correlated using this analysis. Experimental data are presented for both homogeneous and inhomogeneous shearing stress fields. For the latter, a 20-fold range of capillary tube diameters has been employed and no effects of system geometry or the inhomogeneity of the flow-field are observed. Such an independence of the rheological properties from these effects does not occur for low molecular weight liquid crystals and this is, perhaps, the first time this has been reported for polymeric lyotropic liquid crystals; the physical basis for this major difference is discussed briefly. A Semi-empirical constant in eq. (18), N/m² - c rod concentration, rods/m³ - c ^* critical rod concentration at which the isotropic phase becomes unstable, rods/m³ - C interaction potential in the Doi theory defined in eq. (3) - d rod diameter, m - D semi-empirical constant in eq. (19), s^–1 - D _r lumped rotational diffusivity defined in eq. (4), s^–1 - rotational diffusivity of rods in a concentrated (liquid crystalline) system, s^–1 - D _ro rotational diffusivity of a dilute solution of rods, s^–1 - f distribution function defining rod orientation - F tensorial term in the Doi theory defined in eq. (7) (or eq. (19)), s^–1 - G tensorial term in the Doi theory defined in eq. (8) - K _B Boltzmann constant, 1.38 × 10^–23 J/K-molecule - L rod length, m - S scalar order parameter - S tensor order parameter defined in eq. (5) - t time, s - T absolute temperature, K - u unit vector describing the orientation of an individual rod - rate of change ofu due to macroscopic flow, s^–1 - v fluid velocity vector, m/s - v velocity gradient tensor defined in eq. (9), s^–1 - V mean field (aligning) potential defined in eq. (2) - x coordinate direction, m - Kronecker delta (= 0 if = 1 if = ) - _r ratio of viscosity of suspension to that of the solvent at the same shear stress - _s solvent viscosity, Pa · s - ^* viscosity at the critical concentrationc ^*, Pa · s - v ₁, v₂ numerical factors in eqs. (3) and (4), respectively - deviatoric stress tensor, N/m² - volume fraction of rods - ⁰ constant in eq. (16) - ^* volume fraction of rods at the critical concentrationc ^* - average over the distribution functionf(u, t) (= d ²u f(u, t)) - gradient operator - d ²u integral over the surface of the sphere (|u| = 1)     

Some nonlinear problems in servomechanisms

Summary A mathematical model for a hydraulic servomechanism is constructed. It is shown that the model, in general, reduces to a nonlinear third-order equation of the formxx+(₁+xx+–₂ x=p(t). Under certain conditions imposed on the constants involved, it is proved that above equation possesses a periodic solution.     

Optical co-ordinates in general relativity 2 — the problem of distance

Summary Let denote the congruence of null geodesics associated with a given optical observer inV ₄. We prove that determines a unique collection of vector fieldsM₍₎ ( =1, 2, 3) and ₍₀₎ overV ₄, satisfying a weak version of Killing's conditions.This allows a natural interpretation of these fields as the infinitesimal generators of spatial rotations and temporal translation relative to the given observer. We prove also that the definition of the fieldsM₍₎ and ₍₀₎ is mathematically equivalent to the choice of a distinguished affine parameter f along the curves of, playing the role of a retarded distance from the observer.The relation between f and other possible definitions of distance is discussed.

---

Sommario Sia la congruenza di geodetiche nulle associata ad un osservatore ottico assegnato nello spazio-tempoV ₄. Dimostriamo che determina un'unica collezione di campi vettorialiM₍₎ ( =1, 2, 3) e ₍₀₎ inV ₄ che soddisfano una versione in forma debole delle equazioni di Killing. Ciò suggerisce una naturale interpretazione di questi campi come generatori infinitesimi di rotazioni spaziali e traslazioni temporali relative all'osservatore assegnato. Dimostriamo anche che la definizione dei campiM₍₎, ₍₀₎ è matematicamente equivalente alla scelta di un parametro affine privilegiato f lungo le curve di, che gioca il ruolo di distanza ritardata dall'osservatore. Successivamente si esaminano i legami tra f ed altre possibili definizioni di distanza in grande.

---

---

Work performed in the sphere of activity of: Gruppo Nazionale per la Fisica Matematica del CNR.     

Stationary travelling waves on a film flowing down an inclined plane

At small flow rates, the study of long-wavelength perturbations reduces to the solution of an approximate nonlinear equation that describes the change in the film thickness [1–3]. Steady waves can be obtained analytically only for values of the wave numbers close to the wave number _n that is neutral in accordance with the linear theory [1, 2]. Periodic solutions were constructed numerically for the finite interval of wave numbers 0.5_{n n} in [4]. In the present paper, these solutions are found in almost the complete range of wave numbers 0 _n that are unstable in the linear theory. In particular, soliton solutions of this equation are obtained. The results were partly published in [5].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 142–146, July–August, 1980.     

On surface waves in generalized thermoelasticity

Knowles' representation theorem for harmonically time-dependent free surface waves on a homogeneous, isotropic elastic half-space is extended to include harmonically time-dependent free processes for thermoelastic surface waves in generalized thermoelasticity of Lord and Shulman and of Green and Lindsay.r , , r , , .This work was done when author was unemployed.     

Response of a viscous annular liquid layer in zero-gravity to different axial excitations of the rigid boundary plates

A cylindrical annular liquid layer between two plates and around a rigid center-core consisting of incompressible and viscous liquid is subjected to different axial excitations, such as one-sided, counter-directional and double-sided unequal excitations. The response of the free liquid surface, the velocity- and pressure-distribution has been determined.

---

Zusammenfassung Eine zylindrische Flüssigkeitsschicht bestehend aus inkompressibler und viskoser Flüssigkeit wurde verschiedenen harmonischen Anregungsformen ausgesetzt. Dabei wurden die Fälle einseitiger, doppelseitiger entgegengesetzter und ungleicher doppelseitiger Anregung mit Phase behandelt. Die Vergrößerungsfunktionen für die freie Flüssigkeitsoberfläche, für die Geschwindigkeits- und Druckverteilung wurden bestimmt.

---

List of symbols a radius of liquid layer - b radius of inner cylindrical core - (a–b) thickness of layer - e _r , e , k unit vectors in the radial, angular and axial direction resp. - h length of layer - I _m , K _m modified Bessel functions of first and second kind and order m - diameter ratio - p pressure - q 2na/h - q* na/h - r, , z cylindrical coordinates - complex frequency - S sa ²/ - t time - u, w velocity components in the radial- and axial direction - ₀ excitation amplitude - abbreviation - surface tension parameter - surface tension - dynamic viscosity - kinematic viscosity - density of liquid - free liquid surface elevation - dimensionless time - _rz shear stress - reduced forcing frequency - forcing frequency - stream function - _mn natural frequency of non-viscous liquid     

Asymptotic Nusselt numbers for Newtonian flow through a parallel-plate channel with recycle

General expressions for evaluating the asymptotic Nusselt number for a Newtonian flow through a parallel-plate channel with recycle at the ends have been derived. Numerical results with the ratio of thicknesses as a parameter for various recycle ratios are obtained. A regression analysis shows that the results can be expressed by log Nur^0.83=0.3589 (log)² -0.2925 (log) + 0.3348 forR 3, 0.1 0.9; logNu=0.5982(log)² +0.3755 × 10^–2 (log) +0.8342 forR 10^–2, 0.1 0.9.

---

Asymptotische Nusselt-Zahlen für die Newtonsche Strömung durch einen Kanal aus parallelen Platten mit Rückführung

Zusammenfassung In dieser Untersuchung wurden allgemeine Ausdrücke hergeleitet um die asymptotische Nusselt-Zahl für eine Newtonsche Strömung durch einen Kanal aus parallelen Platten mit Rückführung an den Enden berechnen zu können. Es wurden numerische Ergebnisse mit den Dicken-Verhältnissen, als Parameter für verschiedene Rückführungs-verhältnisse, erhalten. Eine Regressionsanalyse zeigt, daß die Ergebnisse wie folgt ausgedrückt werden können: log Nur^0,83=0,3589 (log)² -0,2925 (log) + 0,3348 fürR 3, 0,1 0,9; logNu=0,5982(log)² +0,3755 × 10^–2 (log) + 0,8342 fürR 10^–2, 0,1 0,9.

---

Nomenclature A₁ shooting value,d(0)/d - A₂ shooting value,d(1)/d - B channel width - Gz Graetz number, U_bW²/L - h _m logarithmic average convective heat transfer coefficient - h _x average local convective heat transfer coefficient - k thermal conductivity - L channel length - Nu average local Nusselt number, 2 h_xW/k - Nu _m logarithmic average Nusselt number, 2h_mW/k - R recycle ratio, reverse volume flow rate divided by input volume flow rate - T temperature of fluid - T _m bulk temperature, Eq. (8) - T ₀ temperature of feed stream - T _s wall temperature - U velocity distribution - U _b reference velocity,V/BW - V input volume flow rate - v dimensionless velocity distribution, U/U_b - W channel thickness - x longitudinal coordinate - y transversal coordinate - Z₁-z₆ functions defined in Eq. (A1) - thermal diffusivity - least squares error, Eq. (A7) - weight, Eqs. (A8), (A9) - dimensionless coordinate,y/W - dimensionless coordinate,x/GzL - function, Eq. (7)     

Boundary layer effect on thermal NO decomposition behind incident shock waves

The thermal decomposition of nitric oxide (diluted in Argon) has been measured behind incident shock waves by means of IR diode laser absorption spectroscopy. In two independent runs the diode laser was tuned to the=0 =1² _3/2 R(18.5)-rotational vibrational transition and the=1 =2² _3/2 R(20.5)-rotational vibrational transition of nitric oxide, respectively. These two transitions originating from the vibrational ground state (=0) and the first excited vibrational state (=1) were selected in order to probe the homogeneity along the absorption path. The measured NO decomposition could satisfactorily be described by a chemical reaction mechanism after taking into account boundary layer corrections according to the theory of Mirels. The study forms a further proof of Mirels' theory including his prediction of the laminar-turbulent transition. It also shows, that the inhomogeneities from the boundary layer do not affect the IR linear absorption markedly.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Die Temperaturabhängigkeit der Oberflächenspannung von reinen Kältemitteln vom Tripelpunkt bis zum kritischen Punkt

Zusammenfassung Die Oberflächenspannung von sechs reinen Substanzen — SF₆, CCl₃F, CCl₂F₂, CClF₃, CBrF₃ und CHClF₂ — wurde mit Hilfe einer modifizierten Kapillarmethode gemessen. Die zur Berechnung der Oberflächenspannung erforderlichen Sättigungsdichten und wurden teils aus vorhandenen Zustandsgleichungen, teils aus ebenfalls gemessenen Brechungsindizes bestimmt. Die Temperaturabhängigkeit der Oberflächenspannung läßt sich durch einen erweiterten Ansatz nach van der Waals =_O (T_c-T)(1+...) darstellen, wobei bei einfachen Stoffen ein eingliedriger, bei polaren und assoziierenden Stoffen ein zweigliedriger Ansatz notwendig und ausreichend ist. Für den kritischen Exponenten der Oberflächenspannung wurde ein von der molekularen Substanz weitgehend unabhängiger Wert von =1.284±0.005 gefunden.

---

Temperature dependence of surface tension of pure refrigerants from triple point up to the critical point

The surface tension of six fluids (SF₆, CCl₃F, CCl₂F₂, CClF₃, CBrF₃, CHClF₂) have been measured by means of a modified capillary rise method. The liquid vapor densities, which are needed to calculate the surface tension, have partly been determined by means of refractive indices simultaneously measured in the same apparatus. The temperature dependence of the surface tension is described by an extended van der Waals power law =_O(T_c-T)(1+...). For simple fluids one term and for polar and associating fluids two terms are necessary and sufficient. The critical exponent is found to be 1.284 ± 0.005 and nearly independent of the molecular structure.

---

Formelzeichen a² Laplace-Koeffizient - a Parameter - B_O, B_on Koeffizient der Koexistenzkurve - g Erdbeschleunigung - H Höhe, kapillare Steighöhe - LL Lorentz-Lorenz-Funktion oder Refraktionskonstante - M molare Masse - M Zahl der Meßwerte - N Zahl der unbekannten Parameter - n Brechungsindex - p Druck - R,r Radius - s Entropie - SD Standardabweichung - T, t Temperatur - u innere Energie Griechische Formelzeichen Exponent des Laplace-Koeffizienten - Exponent der Koexistenzkurve - 2. Exponent der Oberflächenspannung - Wellenlänge des Lichts - Exponent der Oberflächenspannung - _D Dipolmoment - , Dichte der Flüssigkeit bzw. des Dampfes - Oberflächenspannung - reduzierte Temperatur (1-T/T_c) - ² gewichtete Varianz Indizes c kritischer Zustand - D Differenz - m Mittelwert - T Isotherme - t Zustand am Tripelpunkt - S Zustand am Schmelzpunkt - bezogen auf Oberfläche     

Viscoelastic swirling flow with free surface in cylindrical chambers

The flow of a viscoelastic liquid driven by the steadily rotating bottom cover of a cylindrical cup is investigated. The flow field and the shape of the free surface are determined at the lowest significant orders of the regular domain perturbation in terms of the angular velocity of the bottom cap. The meridional field superposed on a primary azimuthal field shows a structure of multiple cells. The velocity field and the shape of the free surface are strongly effected by the cylinder aspect ratio and the elasticity of the liquid. The use of this flow configuration as a free surface rheometer to determine the first two Rivlin-Ericksen constants is shown to be promising.Nomenclature R, ,Z Coordinates in the physical domain D - , , Coordinates in the rest stateD ₀ - r, ,z Dimensionless coordinates in the rest stateD ₀ - Angular velocity - Zero shear viscosity - Surface tension coefficient - Density - Dimensionless surface tension parameter - ₁, ₂ The first two Rivlin-Ericksen constants - Stream function - Dimensionless second order meridional stream function - ^* Dimensionless second normal stress function - 2 Dimensionless sum of the first and second normal stress functions - N ₁,N ₂ The first and second normal stress functions - n Unit normal vector - D Stretching tensor - A _n nth order Rivlin-Ericksen tensor - S Extra-stress - u Velocity field - U Dimensionless second order meridional velocity field - V Dimensionless first order azimuthal velocity field - p Pressure - Modified pressure field - P Dimensionless second order pressure field - J Mean curvature - a Cylinder radius - d Liquid depth at rest - D Dimensionless liquid depth at rest - h Free surface height - H Dimensionless free surface height at the second order     

Flow in porous media III: Deformable media

Stokes flow in a deformable medium is considered in terms of an isotropic, linearly elastic solid matrix. The analysis is restricted to steady forms of the momentum equations and small deformation of the solid phase. Darcy's law can be used to determine the motion of the fluid phase; however, the determination of the Darcy's law permeability tensor represents part of the closure problem in which the position of the fluid-solid interface must be determined.Roman Letters A interfacial area of the- interface contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A _e area of entrances and exits for the-phase contained within the macroscopic system, m² - A ^* interfacial area of the- interface contained within a unit cell, m² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - E Young's modulus for the-phase, N/m² - e i unit base vectors (i = 1, 2, 3) - g gravity vector, m²/s - H height of elastic, porous bed, m - k unit base vector (=e ₃) - characteristic length scale for the-phase, m - L characteristic length scale for volume-averaged quantities, m - n unit normal vector pointing from the-phase toward the-phase (n = -n ) - p pressure in the-phase, N/m² - P p – g·r, N/m² - r ₀ radius of the averaging volume, m - r position vector, m - t time, s - T total stress tensor in the-phase, N/m² - T ⁰ hydrostatic stress tensor for the-phase, N/m² - u displacement vector for the-phase, m - V averaging volume, m₃ - V volume of the-phase contained within the averaging volume, m³ - v velocity vector for the-phase, m/s Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - shear coefficient of viscosity for the-phase, Nt/m² - first Lamé coefficient for the-phase, N/m² - second Lamé coefficient for the-phase, N/m² - bulk coefficient of viscosity for the-phase, Nt/m² - T –T ⁰ , a deviatoric stress tensor for the-phase, N/m²     

Berechnung des Temperaturfeldes um eine im Erdboden verlegte Punkt- und Linienquelle bei Randbedingung III. Art an der Erdoberfläche

Zusammenfassung Für den Fall, daß sich in einem halbunendlichen Körper in der Tiefe L eine Punkt- bzw. Linienquelle befindet und daß an der Oberfläche des Körpers ein örtlich und zeitlich konstanter Wärmeübergangskoeffizient herrscht, wird das stationäre Temperaturfeld analytisch berechnet. Beim Vergleich mit einer Näherungslösung (Hilfsschicht) zeigt sich, daß nicht so sehr die Biot-Zahl Bi= · L/ als vielmehr der größte Winkel zwischen Wandnormale und Wärmestromdichte in der Hilfsschicht ein Maß für die Genauigkeit der Näherungslösung ist.

---

Calculation of the temperature field around a buried point- and linesource, respectively, when the boundary condition is Newton's law

The steady state temperature field in a semiinfinite body caused by a buried point- and linesource, respectively, has been analytically calculated. The comparison with a simple approach (additional-layer) shows that the greatest angle between the normal of the wall and the heat flux density in the additional-layer, describes the quality of the approach better than the Biot-number Bi=L/ does.

---

Formelzeichen A Fläche - Bi Biot-Zahl - C Eigenwertfunktion - E₁ Exponentialintegral - exp Exponentialfunktion - i komplexe Einheit - J₀ Besselfunktion nullter Ordnung und 1. Grades - L Verlegungstiefe der Punkt- bzw. Linienquelle - Q Quellstärke - r Radius - Re Realteil eines Ausdruckes - T Temperatur - t Integrationsvariable - x, y, z Ortskoordinaten - Wärmeübergangskoeffizienten an der Erdoberfläche - Laplace-Operator - Wärmeleitfähigkeit des Erdbodens - dimensionslose Temperatur - Integrationsvariable - dimensionsloser Radius - komplexe Ortskoordination Indizes 0 Erdoberfläche, senkrecht über der Quelle - 1 Lösung für das 1. Randwertproblem - 3 Lösung für das 3. Randwertproblem - 13 Zusatzfunktion - w Erdoberfläche - Umgebungstemperatur - Näherungslösung     

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.