The dynamic and steady shear properties of synovial fluid and of the components making up synovial fluid

Steady-shear and dynamic properties of a pooled sample of cattle synovial fluid have been measured using techniques developed for low viscosity fluids. The rheological properties of synovial fluid were found to exhibit typical viscoelastic behaviour and can be described by the Carreau type A rheological model. Typical model parameters for the fluid are given; these may be useful for the analysis of the complex flow problems of joint lubrication.The two major constituents, hyaluronic acid and proteins, have been successfully separated from the pooled sample of synovial fluid. The rheological properties of the hyaluronic acid and the recombined hyaluronic acid-protein solutions of both equal and half the concentration of the constituents found in the original synovial fluid have been measured. These properties, when compared to those of the original synovial fluid, show an undeniable contribution of proteins to the flow behaviour of synovial fluid in joints. The effect of protein was found to be more prominent in hyaluronic acid of half the normal concentration found in synovial fluid, thus providing a possible explanation for the differences in flow behaviour observed between synovial fluid from certain diseased joints compared to normal joint fluid.Nomenclature A Ratio of angular amplitude of torsion head to oscillation input signal - G Storage modulus - G Loss modulus - I Moment of inertia of upper platen — torsion head assembly - K Restoring constant of torsion bar - N ₁ First normal-stress difference - R Platen radius - S ⁽ⁱ⁾ Geometric factor in the dynamic property analysis - t ₁ Characteristic time parameter of the Carreau model - X, Y Carreau model parameters - Z () Reimann Zeta function of - Carreau model parameter - Shear rate - Apparent steady-shear viscosity - ^* Complex dynamic viscosity - Dynamic viscosity - Imaginary part of the complex dynamic viscosity - ₀ Zero-shear viscosity - ₀ Cone angle - Carreau model characteristic time - Density of fluid - Shear stress - Phase difference between torsion head and oscillation input signals - ₀ Zero-shear rate first normal-stress coefficient - Oscillatory frequency     

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/     

A rotating hot-wire technique for spatial sampling of disturbed and manipulated duct flows

The documentation and control of flow disturbances downstream of various open inlet contractions was the primary focus with which to evaluate a spatial sampling technique. An X-wire probe was rotated about the center of a cylindrical test section at a radius equal to one-half that of the test section. This provided quasi-instantaneous multi-point measurements of the streamwise and azimuthal components of the velocity to investigate the temporal and spatial characteristics of the flowfield downstream of various contractions. The extent to which a particular contraction is effective in controlling ingested flow disturbances was investigated by artificially introducing disturbances upstream of the contractions. Spatial as well as temporal mappings of various quantities are presented for the streamwise and azimuthal components of the velocity. It was found that the control of upstream disturbances is highly dependent on the inlet contraction; for example, reduction of blade passing frequency noise in the ground testing of jet engines should be achieved with the proper choice of inlet configurations.List of symbols K _uv correlation coefficient= - P percentage of time that an azimuthal fluctuating velocity derivative dv/d is found - U streamwise velocity component U=U (, t) - V azimuthal or tangential velocity component due to flow and probe rotation V=V (, t) - mean value of streamwise velocity component - U _m resultant velocity from and - mean value of azimuthal velocity component induced by rotation - u fluctuating streamwise component of velocity u=u(, t) - v fluctuating azimuthal component of velocity v = v (, t) - u phase-averaged fluctuating streamwise component of velocity u=u(0) - v phase-averaged fluctuating azimuthal component of velocity v=v() - û average of phase-averaged fluctuating streamwise component of velocity (u()) over cases I-1, II-1 and III-1 û = û() - average of phase-averaged fluctuating azimuthal component of velocity (v()) over cases I-1, II-1 and III-1 - u fluctuating streamwise component of velocity corrected for non-uniformity of probe rotation and/or phase-related vibration u = u(0, t) - v fluctuating azimuthal component of velocity corrected for non-uniformity or probe rotation and/or phase-related vibration v=v (, t) - u ² rms value of corrected fluctuating streamwise component of velocity - rms value of corrected fluctuating azimuthal component of velocity - phase or azimuthal position of X-probe     

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     

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     

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     

Critical phenomena and waves in gas-liquid systems

The results of the hydraulic studies of gas-liquid media, wave processes in two-phase media and critical phenomena are described. Some methodological foundations to describe these media and methods to obtain the basic similarity criteria for the hydraulics and gas-dynamics of bubble suspensions are discussed. A detailed consideration is given for the phase transition processes on interfaces and the interface stability. A relation has been revealed between the wave and critical phenomena in two-phase systems.Nomenclature a thermal diffusivity - Ar Archimedes number - B gas constant - C heat capacity - C _p heat capacity at constant pressure - C _v heat capacity at constant volume - c ₀ acoustic velocity in the mixture - c _l acoustic velocity in the liquid - C _f flow resistance coefficient - G mass rate of flow - g gravitational acceleration - L latent heat of evaporation - l initial perturbation width - M Mach number - Nu Nusselt number - P pressure - Pr Prandtl number - R bubble radius - (3P ₀/R ₀ ² _f )^–1 bubble resonance frequency square - T temperature - U medium motion velocity - W heavy phase velocity - W light phase velocity - We Weber number - heat release coefficient - dispersion coefficient - void fraction - adiabatic index - film thickness - dimensionless film thickness - kinematic viscosity coefficient - dynamical viscosity coefficient - dissipation coefficient in the mixture - dispersion parameter - _f liquid phase density - light phase density - heat conductivity - surface tension - frequency, ₀ ² =3P ₀/ _f R ₀ ²     

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     

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)     

Mass transfer effects on flow past a vertical oscillating plate with variable temperature

An exact solution to the flow of a viscous incompressible fluid past an infinite vertical oscillating plate, in the presence of a foreign mass has been derived by the Laplace-transform technique when the plate temperature is linearly varying as time. The velocity profiles are shown on graphs and the numerical values of the skin-friction are listed in a table. It is observed that the skin-friction increases with increasingSc, t orGr but decreases with increasingGm ort, whereSc (the Schmidt number), (frequency),t (time),Gr (the Grashof number) andGm is (the modified Grashof number) andt.

---

Einfluß von Stoffübergang auf die Strömung entlang einer senkrechten oszillierenden Platte veränderlicher Temperaturen

Zusammenfassung Mit Hilfe der Laplace-Transformation wird eine exakte Lösung für die Strömung einer zähen, inkompressiblen Flüssigkeit entlang einer unendlich ausgedehnten, senkrechten, oszillierenden Platte gewonnen, wobei die Einwirkung eines Feststoffs Berücksichtigung findet und die Plattentemperatur linear mit der Zeit veränderlich sein soll. Die Geschwindigkeitsprofile sind in Diagrammen dargestellt und die numerischen Werte der Reibungsschubspannung in einer Tabelle. Letztere wächst mit der Schmidt-ZahlSc, der Grashof-ZahlGr und dem Produkt aus Frequenz und Zeitt; sie nimmt ab, wennGm (die modifizierte Grashof-Zahl) undt zunehmen.

---

Nomenclature C species concentration in the fluid near the plate - C ^– species concentration in the fluid away from the plate - C _W species concentration at the plate - C _p specific heat at constant pressure - D chemical molecular diffusivity - Gm modified Grashof number - Gr Grashof number - g acceleration due to gravity - K thermal conductivity - P Prandtl number - Sc Schmidt number - T Temperature of the fluid near the plate - T _W temperature of the plate - T temperature of the fluid far away from the plate - t time - u velocity of the fluid in the upward direction - U ₀ amplitude of oscillation - x coordinate axis along the plate in the vertically upward direction - y coordinate axis normal to the plate Greek symbols skin-friction - viscosity - coefficient of volume expansion - ^* coefficient of species expansion - density - frequency     

Zum Impuls-, Wärme- und Stoffaustausch in turbulenten Strömungen reagierender Binärgemische

Zusammenfassung Die in Teil I vorgestellten Reynolds 'schen Gleichungen und Transportgleichungen werden für Strömungen mit Grenzschichtcharakter angegeben. Weiter werden Integralbedingungen mitgeteilt. Nach einer Diskussion über die Schließung des Gleichungssystems werden Lösungsverfahren besprochen. Dabei wird speziell auf Integralverfahren eingegangen.

---

About the transfer of momentum, heat and mass in turbulent flows of binary mixturesPart II: Thin shear flow layers

The Reynolds equations and transport equations given in part I are presented for thin shear flow layers. Integral relations are given. After a discussion of the closure problem methods of solution are described. Specially integral methods are discussed.

---

Formelzeichen c Massenkonzentration der Komponente - c_t charakteristische Konzentrationsschwankung - c_o Bezugskonzentration - c spezifische Wärme bei konstantem Druck - c_f Reibungsbeiwert - c_D Dissipationsintegral - c_E Entrainment-Funktion - c Schubspannungsintegral - D binsrer Diffusionskoeffizient - H Formparameter - H₁₂ Formparameter - H₃₂ Formparameter - j Kassendiffusionsstrom - L Bezugslänge - p Druck - p_t charakteristische Druckschwankung - p_o Bezugsdruck - Pr Prandtl-Zahl - q Wärmestrom - q²/2 kinetische Energie der Schwankungsbewegung - Re_L mit L gebildete Reynolds-Zahl - Re mit gebildete Reynolds-Zahl - Re₂ mit ₂ gebildete Reynolds-Zahl - Sc Schmidt-Zahl - T absolute Temperatur - T_t charakteristische TemperaturSchwankung - T_o Bezugstemperatur - u,v,w Geschwindigkeitskomponenten - u_t charakteristische Geschwindigkeitsschwankung - u_o Bezugsgeschwindigkeit - U=/ü dimensionslose. x-Komponente der Geschwindigkeit - x,y,z Komponenten des Ortsvektors Griechische Symbole Grenzschichtdicke - ₁ Verdrängungsdicke - ₂ Impulsverlustdicke - ₃ Energieverlustdicke - _T Enthalpieverlustdicke - _c Konzentrationsverlustdicke - =d/dx Parameter für die Grenzschichtabsch:atzung - turbulente Impulsaustauschgröße - _D turbulente Stoffaustauschgröße - _q turbulente Energieaustauschgröße - Dissipationsfunktion - Wärmeleitfähigkeit - dynamische Viskosität - v=/ kinematische Viskosität - Dichte - Produktionsdichte - Schubspannung Indizes mol molekularer Anteil - tur turbulenter Anteil - res resultierender Anteil - Außenrand der Grenzschicht - w Wand     

Vaporization of single liquid drops in an immiscible liquid part II: Heat transfer characteristics

Heat transfer characteristics during the vaporization process of a pentane or furan drop in an aqueous glycerol of high viscosity has been studied. With the progress of vaporization, the overall heat transfer coefficient related to the liquid-liquid interfacial area of a two-phase bubble increases monotonically, and influences of initial drop diameter and temperature difference reduce. Some convection or circulation seems to occur in the unvaporized-liquid phase.

---

Verdampfung einzelner Flüssigkeitstropfen in einer nicht mischbaren Flüssigkeit. Teil II: Der Wärmeübergang

Zusammenfassung In dieser Arbeit wird der Wärmeübergang während der Verdampfung von Pentan- und Furan-Tropfen in einer wässerigen Glyzerinlösung hoher Viskosität untersucht. Mit fortschreitender Verdampfung steigt der Wärmeübergangskoeffizient, bezogen auf die Grenzfläche flüssig-flüssig der zweiphasigen Blase monoton an, wobei Einflüsse des anfänglichen Tropfendurchmessers und der Temperaturdifferenz abnehmen. In der nichtverdampften Flüssigkeitsphase scheint Konvektion oder Zirkulation aufzutreten.

---

Nomenclature A total surface area of two-phase bubble - A_L liquid-liquid interfacial area of two-phase bubble - D equivalent spherical diamter of two-phase bubble - D_i initial drop diameter - h average overall heat transfer coefficient related to A - h_c average outside heat transfer coefficient related to A - q local outside heat transfer coefficient - h_L average overall heat transfer coefficient related to A_L - h_Lc average outside heat transfer coefficient related to A_L - k_c thermal conductivity of continuous-phase liquid - k_dl thermal conductivity of dispersed-phase liquid - k_v correction factor of velocity [cf. Eq.(2)] - Nu_c =h_c D/k - Nu_c =h_c D/k_c - Pe_c =UD/_c - Pr_c =_c/_c - Q cumulative heat transferred into two-phase bubble - q local heat flux - r radial distance in spherical co-ordinates - R radius of two-phase bubble - T temperature - T_L interface temperature between continuousphase and dispersed-phase component in liquid phase - T bulk temperature - T temperature difference - T nominal temperature difference - U velocity of rise of two-phase bubble - u velocity gradient in r direction [cf. Eq.(9)] - u_r velocity component in r direction - u velocity component in direction - V volume of two-phase bubble - V_dl volume of dispersed-phase component in liquid phase - X defined in Eq.(7) - Y defined in Eq.(8) - Z defined in Eq.(12) - _c thermal diffusivity of continuous-phase liquid - half opening angle of vapor phase in two-phase bubble - average thickness of dispersed-phase component in liquid phase [cf. Eq.(22)] - angle in spherical co-ordinates - vaporization ratio - time     

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     

Nonlinear analysis of the forced response of a beam with three mode interaction

An analysis is presented for the primary resonance of a clamped-hinged beam, which occurs when the frequency of excitation is near one of the natural frequencies,_n . Three mode interaction (₂ 3₁ and _{3 1} + 2₂) is considered and its influence on the response is studied. The case of two mode interaction (₂ 3₁) is also considered to compare it with the case of three mode interaction. The straight beam experiencing mid-plane stretching is governed by a nonlinear partial differential equation. By using Galerkin's method the governing equation is reduced to a system of nonautonomous ordinary differential equations. The method of multiple scales is applied to solve the system. Steady-state responses and their stability are examined. Results of numerical investigations show that there exists no significant difference between both modal interactions' influences on the responses.     

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²     

Turbulent boundary layer over a smooth wall with widely separated transverse square cavities

Laser-Doppler velocimetry (LDV) measurements and flow visualizations are used to study a turbulent boundary layer over a smooth wall with transverse square cavities at two values of the momentum thickness Reynolds number (R =400 and 1300). The cavities are spaced 20 element widths apart in the streamwise direction. Flow visualizations reveal a significant communication between the cavities and the overlying shear layer, with frequent inflows and ejections of fluid to and from cavities. There is evidence to suggest that quasi-streamwise near-wall vortices are responsible for the ejections of fluid out of the cavities. The wall shear stress, which is measured accurately, increases sharply immediately downstream of the cavity. This increase is followed by a sudden decrease and a slower return to the smooth wall value. Integration of the wall shear stress in the streamwise direction indicates that there is an increase in drag of 3.4% at bothR .Nomenclature C _f skin friction coefficient - C _fsw friction coefficient for a continuous smooth wall - k height of the cavity - k ⁺ ku / - R Reynolds number based on momentum thickness (U ₁ /v) - Rx Reynolds number based on streamwise distance (U ₁ x/) - s streamwise distance between two cavities - t time - t ⁺ tu ₂ / - U ₁ freestream velocity - mean velocity inx direction - u,v,w rms turbulent intensities inx,y andz directions - u local friction velocity - u _sw friction velocity for a continuous smooth wall - w width of the cavity - x streamwise co-ordinate measured from the downstream edge of the cavity - y co-ordinate normal to the wall - z spanwise co-ordinate - y ⁺ yu / - boundary layer thickness - ₀ boundary layer thickness near the upstream edge of the cavity - _i thickness of internal layer - kinematic viscosity of water - ⁺ zu / - momentum thickness     

The swirling radial jet

The similarity solution of the radial turbulent jet with weak swirl is discussed and a new solution of the radial turbulent jet with swirl is proposed without restrictions assumed in the weak swirl solution.Nomenclature e swirl parameter - k experimental constant - l non-negative constant - M, M , N, P integral invariants - q velocity component in -direction - q _max maximum velocity component in -direction - u radial velocity component - u _max maximum radial velocity component - v axial velocity component - w peripheral velocity component - w _max maximum peripheral velocity component - x radial coordinate - y transverse coordinate - angle introduced in (28) - characteristic width of a jet - (x, y) similarity variable (scaled x and y coordinate) - molecular kinematic viscosity - _T eddy kinematic viscosity - tangential coordinate - fluid density - turbulent shear stress in -direction - _xy , _y components of turbulent shear stress tensor - (x, y) stream function     

Estimate of the transient conduction of heat in materials with linear thermal properties based on the solution for constant properties

A method of analysis is described which yields quasianalytical solutions for one and multidimensional unsteady heat conduction problems with linearly dependent thermal properties, such as thermal conductivity and volumetric specific heat. The method accomodates rather general thermal boundary conditions including arbitrary variations in surface temperature or in surface heat flux or a convective exchange with a fluid having even varying temperature. Once the solution for the identical problem but with constant properties has been developed, its practical realization is rather direct, being facilitated by a reduced number of iterations. The four applied examples given in this work show that a wide variety of nonlinear heat conduction problems can be tackled by this procedure without much difficulty. These simple solutions compare favorably with more laborious results reported in the archival heat transfer literature.

---

Berechnung nichtstationärer Wärmeleitvorgänge mit linear temperaturabhängigen Stoffwerten aus der Lösung für konstante Stoffwerte

Zusammenfassung Es werden quasi-analytische Lösungen für ein- und mehrdimensionale nichtstationäre Wärmeleitprobleme mit linear temperaturabhängigen Stoffwerten, wie Wärmeleitfähigkeit und volumetrische Wärmekapazität, mitgeteilt. Die Methode gilt für recht allgemeine Randbedingungen wie beliebige Veränderungen der Oberflächentemperatur, der Wärmestromdichte oder auch konvektiven Wärmeaustausch mit veränderlicher Fluidtemperatur. Ist die Lösung für das identische Problem mit konstanten Stoffwerten bekannt, kann die Methode direkt mit einer begrenzten Zahl von Iterationen angewandt werden. Die vier hier mitgeteilten Beispiele zeigen, daß eine große Zahl nichtlinearer Wärmeleitprobleme auf diese Weise ohne Schwierigkeit angepackt werden können. Die einfachen Lösungen stimmen befriedigend mit komplizierteren Ergebnissen aus der Literatur überein.

---

Nomenclature a side of square bar - B _i0 reference Biot number,hR/k₀ - B _i0 ^T transformed Biot number, equation (16) - c geometric parameter, equation (8) - h convective coefficient - k thermal conductivity - k ₀ value ofk atT ₀ - K dimensionless thermal conductivity,k/k ₀ - K _i value ofK at _i - K _i+1 value ofK at _i+1 - m _k slope of theK- line, equation (3) - m _s slope of theS- line, equation (4) - R characteristic length - s volumetric specific heat - s ₀ value of s at T₀ - S dimensionless volumetric specific heat, s/s₀ - S _i value ofS at _i - S _i+1 value of S at _i+1 - t time - T temperature - T ₀ reference temperature - x, y cartesian coordinates - X, Y dimensionless cartesian coordinates,x/a andy/a - thermal diffusivity - _k transformed time, equation (11) - _s transformed time, equation (37) - _k dimensionless time for variable conductivity, equation (8) - _s dimensionless time for variable specific heat, equation (34) - dimensionless temperature,T/T ₀ - dimensionless coordinate,r/R - ₀ value of at T₀ - _i lower value of the interval (_i, _i+1) - _i+1 upper value of the interval (_i, _i+1     

Deformation to breaking of deep water gravity waves

The results of laboratory observations of the deformation of deep water gravity waves leading to wave breaking are reported. The specially developed visualization technique which was used is described. A preliminary analysis of the results has led to similar conclusions than recently developed theories. As a main fact, the observed wave breaking appears as the result of, first, a modulational instability which causes the local wave steepness to approach a maximum and, second, a rapidly growing instability leading directly to the breaking.List of symbols L total wave length - H total wave height - crest elevation above still water level - trough depression below still water level - wave steepness =H/L - crest steepness =/L - trough steepness =/L - F ₁ forward horizontal length from zero-upcross point (A) to wave crest - F ₂ backward horizontal length from wave crest to zero-downcross point (B) - crest front steepness =/F ₁ - crest rear steepness =/F ₂ - vertical asymmetry factor=F ₂/F ₁ (describing the wave asymmetry with respect to a vertical axis through the wave crest) - µ horizontal asymmetry factor=/H (describing the wave asymmetry with respect to a horizontal axis: SWL) - T ₀ wavemaker period - L ₀ theoretical wave length of a small amplitude sinusoïdal wave generated at T _inf0 ^sup–1 frequency - ₀ average wave height     

Some remarks on differential equations of quadratic type

In this paper we study differential equations of the formx(t) + x(t)=f(x(t)), x(0)=x ₀ C HereC is a closed, bounded convex subset of a Banach spaceX,f(C) C, and it is often assumed thatf(x) is a quadratic map. We study the differential equation by using the general theory of nonexpansive maps and nonexpansive, non-linear semigroups, and we obtain sharp results in a number of cases of interest. We give a formula for the Lipschitz constant off: C C, and we derive a precise explicit formula for the Lipschitz constant whenf is quadratic,C is the unit simplex inR ⁿ, and thel ₁ norm is used. We give a new proof of a theorem about nonexpansive semigroups; and we show that if the Lipschitz constant off: CC is less than or equal to one, then lim_tf(x(t))–x(t)=0 and, if {x(t):t 0} is precompact, then lim_tx(t) exists. Iff¦C=L¦C, whereL is a bounded linear operator, we apply the nonlinear theory to prove that (under mild further conditions on C) lim_t f(x(t))–x(t)=0 and that lim_t x(t) exists if {x(t):t 0} is precompact. However, forn 3 we give examples of quadratic mapsf of the unit simplex ofR ⁿ into itself such that lim_t x(t) fails to exist for mostx ₀ C andx(t) may be periodic. Our theorems answer several questions recently raised by J. Herod in connection with so-called model Boltzmann equations.