Evaporation of single liquid drops in an immiscible liquid at elevated pressures: equipment and preliminary results

This paper describes a piece of equipment which has been specially designed to permit the study of the evaporation of single liquid drops in an immiscible liquid at elevated pressures up to the order of 1 MPa. It is equipped with a piezoelectric discharger for yielding nucleation in each drop and with a dilatometer for detecting the change in volume of the drop in the course of evaporation succeeding the nucleation. Some preliminary results, obtained with the equipment, for n-pentane drops evaporating in a water medium at pressures up to 0.5 MPa (5 atm) are also presented.List of symbols C _D drag coefficient - D ₀ initial, equivalent spherical diameter of drop (= (6 V ₀/)^1/3) - H vertical position of two-phase bubble measured from the nozzle tip - H _t column height required for drop to vaporize completely - Nu Nusselt number related to instantaneous two-phase bubble diameter and thermal conductivity of the continuous phase - p ^* pressure at the nozzle tip - Pe Peclet number related to instantaneous diameter and rise velocity of two-phase bubble and to thermal diffusivity of the continuous phase - Re Reynolds number related to instantaneous diameter and rise velocity of two-phase bubble and to viscosity of the continuous phase - t time - t _v time required for drop to vaporize completely - T temperature in the bulk of the medium - T excess, assumed in theoretical model, of T above the temperature of two-phase bubble - T ₁ ^* excess of T above the saturation temperature of the dispersed-phase fluid corresponding to p ^* - T ₂ ^* excess of T above the temperature at which the sum of the saturated vapor pressures of the dispersed- and the continuous-phase fluids is equal to p ^* - V volume of two-phase bubble - V ₀ initial value of V     

Theoretical derivation of tangential velocity profiles in a flat vortex chamber-influence of turbulence and wall friction

Summary As part of a study on the hydrodynamics of a cyclone separator, a theoretical investigation of the flow pattern in a flat box cyclone (vortex chamber) has been carried out. Expressions have been derived for the tangential velocity profile as influenced by internal friction (eddy viscosity) and wall friction. The most important parameter controlling the tangential velocity profile is = –u ₀ R/(v+ ), where u ₀ is the radial velocity at the outer radius R of the cyclone, the kinematic liquid viscosity and is the kinematic eddy viscosity. For values of greater than about 10 the tangential velocity profile is nearly hyperbolic, for smaller than 1 the tangential velocity even decreases towards the centre. It is shown how and also the wall friction coefficient may be obtained from experimental velocity profiles with the aid of suitable graphs. Because of the close relation between eddy viscosity and eddy diffusion, measurements of velocity profiles in flat box cyclones will also provide information on the eddy motion of particles in a cyclone, a motion reducing its separation efficiency.List of symbols A cross-sectional area of cyclone inlet - h height of cyclone - p static pressure in cyclone - p static pressure difference in cyclone between two points on different radius - r radius in cyclone - r ₁ radius of cyclone outlet - R radius of cyclone circumference - u radial velocity in cyclone - u ₀ radial velocity at circumference of flat box cyclone - v tangential velocity - v ₀ tangential velocity at circumference of flat box cyclone - w axial velocity - z axial co-ordinate in cyclone - friction coefficient in flat box cyclone (for definition see § 5) - ₁ value of friction coefficient for ₁<< ₂ - ₂ value of friction coefficient for ₂<<1 - = - ₁ value of for ₁<< ₂ - ₂ value of for ₂<<1 - thickness of laminar boundary layer - =/h - turbulent kinematic viscosity - ratio of z to h - k ratio of height of cyclone to radius R of cyclone - parameter describing velocity profile in cyclone =–u ₀ R/(+) - kinematic viscosity of fluid - density of fluid - ratio of r to R - ₁ value of at outlet of cyclone - ₂ value of at inner radius of cyclone inlet - _w shear stress at cyclone wall - angular momentum in cyclone/angular momentum in cyclone inlet - ₁ value of at = ₁ - ₂ value of at = ₂     

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     

Motion of large gas bubbles ascending in a liquid

An equation is derived for the ascent velocity of large gas bubbles in a liquid. This velocity is assumed to be governed by the propagation of a wavelike perturbation caused by the bubble in the liquid.Notation w bubble (or drop) velocity - specific gravity - dynamic viscosity - kinematic viscosity - r bubble (or drop) radius - surface tension - coefficient of friction - g gravitational acceleration - D bubble (or drop) diameter - p pressure - c propagation velocity of the wavelike perturbation - wavelength     

Parameter determination for thermodynamical models of the pseudoelastic behaviour of shape memory alloy by resistivity measurements and infrared thermography

Two thermodynamical models of pseudoelastic behaviour of shape memory alloys have been formulated. The first corresponds to the ideal reversible case. The second takes into account the hysteresis loop characteristic of this shape memory alloys.Two totally independent techniques are used during a loading-unloading tensile test to determine the whole set of model parameters, namely resistivity and infrared thermography measurements. In the ideal case, there is no difficulty in identifying parameters.Infrared thermography measurements are well adapted for observing the phase transformation thermal effects.Notations 1 austenite 2 martensite - ₍₎ Macroscopic infinitesimal strain tensor of phase - ₍₂₎ ^f Traceless strain tensor associated with the formation of martensite phase - Macroscopic infiniesimal strain tensor - Macroscopic infinitesimal strain tensor deviator - _f Trace - Equivalent strain - ^pe Macroscopic pseudoelastic strain tensor - x Distortion due to parent (austenite =1)product (martensite =2) phase transformation (traceless symmetric second order tensor) - M Total mass of a system - M() Total mass of phase - V Total volume of a system - V() Total volume of phase - z=M(2)/M Weight fraction of martensite - 1-z=M(1)/M Weight fraction of austenite - u ₀ ^* () Specific internal energy of phase (=1,2) - s ₀ ^* () Specific internal entropy of phase - Specific configurational energy - Specific configurational entropy - ₀ ^f (T) Driving force for temperature-induced martensitic transformation at stress free state ( ₀ ^f T) = T ^*–Ts ^*) - Kirchhoff stress tensor - Kirchhoff stress tensor deviator - Equivalent stress - Cauchy stress tensor - Mass density - K Bulk moduli (K ₀=K) - L Elastic moduli tensor (order 4) - E Young modulus - Energetic shear (₀ = ) - Poisson coefficient - M _s ^o (M _F ^o ) Martensite start (finish) temperature at stress free state - A _s ^o (A _F ^o ) Austenite start (finish) temperature at stress free state - C _v Specific heat at constant volume - k Conductivity - Pseudoelastic strain obtained in tensile test after complete phase transformation (AM) (unidimensional test) - ₀ Thermal expansion tensor - r Resistivity - 1MPa 10⁶ N/m ² - ⁽⁾ Specific free energy of phase - _n Specific free energy at non equilibrium (R model) - _n ^eq Specific free energy at equilibrium (R model) - _n ^v Volumic part of ^eq - Specific free energy at non equilibrium (R _L model) - _conf Specific coherency energy (R _L model) - _c Specific free energy at constrained equilibria (R _L model) - _it (T) Coherency term (R _L model)     

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

Effects of MHD free convection and mass transfer on the flow past an oscillating infinite coaxial vertical circular cylinder

The effects of MHD free convection and mass transfer are taken into account on the flow past oscillating infinite coaxial vertical circular cylinder. The analytical expressions for velocity, temperature and concentration of the fluid are obtained by using perturbation technique.

---

Einwirkungen von freier MHD-Konvektion und Stoffübertragung auf eine Strömung nach einem schwingenden unendlichen koaxialen vertikalen Zylinder

Zusammenfassung Die Einwirkungen der freien MHD-Konvektion und Stoffübertragung auf eine Strömung nach einem schwingenden, unendlichen, koaxialen, vertikalen Zylinder wurden untersucht. Die analytischen Ausdrücke der Geschwindigkeit, Temperatur und Fluidkonzentration sind durch die Perturbationstechnik erhalten worden.

---

Nomenclature C _p Specific heat at constant temperature - C the species concentration near the circular cylinder - C _w the species concentration of the circular cylinder - C the species concentration of the fluid at infinite - * dimensionless species concentration - D chemical molecular diffusivity - g acceleration due to gravity - Gr Grashof number - Gm modified Grashof number - K thermal conductivity - Pr Prandtl number - r _a ,r _b radius of inner and outer cylinder - a, b dimensionless inner and outer radius - r,r coordinate and dimensionless coordinate normal to the circular cylinder - Sc Schmidt number - t time - t dimensionless time - T temperature of the fluid near the circular cylinder - T _w temperature of the circular cylinder - T temperature of the fluid at infinite - u velocity of the fluid - u dimensionless velocity of the fluid - U ₀ reference velocity - z,z coordinate and dimensionless coordinate along the circular cylinder - coefficient of volume expansion - * coefficient of thermal expansion with concentration - dimensionless temperature - H ₀ magnetic field intensity - coefficient of viscosity - _e permeability (magnetic) - kinematic viscosity - electric conductivity - density - M Hartmann number - dimensionless skin-friction - frequency - dimensionless frequency     

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 ₀ ²     

The contribution of internal degrees of freedom to the non-Newtonian rheology of model polymer fluids

We report non-equilibrium molecular dynamics simulations of rigid and non-rigid dumbbell fluids to determine the contribution of internal degrees of freedom to strain-rate-dependent shear viscosity. The model adopted for non-rigid molecules is a modification of the finitely extensible nonlinear elastic (FENE) dumbbell commonly used in kinetic theories of polymer solutions. We consider model polymer melts — that is, fluids composed of rigid dumbbells and of FENE dumbbells. We report the steady-state stress tensor and the transient stress response to an applied Couerte strain field for several strain rates. We find that the rheological properties of the rigid and FENE dumbbells are qualitatively and quantitatively similar. (The only exception to this is the zero strain rate shear viscosity.) Except at high strain rates, the average conformation of the FENE dumbbells in a Couette strain field is found to be very similar to that of FENE dumbbells in the absence of strain. The theological properties of the two dumbbell fluids are compared to those of a corresponding fluid of spheres which is shown to be the most non-Newtonian of the three fluids considered.Symbol Definition b dimensionless time constant relating vibration to other forms of motion - F force on center of mass of dumbbell - F _i force on bead i of dumbbell - F force between center of masses of dumbbells and - F _ij force between beads i and j - h vector connecting bead to center of mass of dumbbell - H dimensionless spring constant for dumbbells, in units of / ² - I moment of inertia of dumbbell - J general current induced by applied field - k _B Boltzmann's constant - L angular momentum - m mass of bead, (= m/2) - M mass of dumbbell, g - N number of dumbbells in simulation cell - P translational momentum of center of mass of dumbbell - P pressure tensor - P _xy xy component of pressure tensor - Q separation of beads in dumbbell - Q _eq equilibrium extension of FENE dumbbell and fixed extension of rigid dumbbell - Q ₀ maximum extension of dumbbell - r _ij vector connecting beads i and j - r position vector of center of mass dumbbell - R vector connecting centers of mass of two dumbbells - t time - t ^* dimensionless time, in units of m/ - T ^* dimensionless temperature, in units of /k - u potential energy - u velocity vector of flow field - u _x x component of velocity vector - V volume of simulation cell - X general applied field - strain rate, s^–1 - ^* dimensionless shear rate, in units of /m ² - general transport property - Lennard-Jones potential well depth - friction factor for Gaussian thermostat - shear viscosity, g/cms - ^* dimensionless shear viscosity, in units of m/ ² - ^* dimensionless number density, in units of ^–3 - Lennard-Jones separation of minimum energy - relaxation time of a fluid - angular velocity of dumbbell - orientation angle of dumbbell      

Eine analytische Näherungslösung für die eingefrorene laminare Grenzschichtströmung eines binären Gemisches

Zusammenfassung Für die eingefrorene laminare Grenzschichtströmung eines teilweise dissoziierten binären Gemisches entlang einer stark gekühlten ebenen Platte wird eine analytische Näherungslösung angegeben. Danach läßt sich die Wandkonzentration als universelle Funktion der Damköhler-Zahl der Oberflächenreaktion angeben. Für das analytisch darstellbare Konzentrationsprofil stellt die Damköhler-Zahl den Formparameter dar. Die Wärmestromdichte an der Wand bestehend aus einem Wärmeleitungs- und einem Diffusionsanteil wird angegeben und diskutiert. Das Verhältnis beider Anteile läßt sich bei gegebenen Randbedingungen als Funktion der Damköhler-Zahl ausdrücken.

---

An analytical approximation for the frozen laminar boundary layer flow of a binary mixture

An analytical approximation is derived for the frozen laminar boundary layer flow of a partially dissociated binary mixture along a strongly cooled flat plate. The concentration at the wall is shown to be a universal function of the Damkohler-number for the wall reaction. The Damkohlernumber also serves as a parameter of shape for the concentration profile which is presented in analytical form. The heat transfer at the wall depending on a conduction and a diffusion flux is derived and discussed. The ratio of these fluxes is expressed as a function of the Damkohler-number if the boundary conditions are known.

---

Formelzeichen A Atom - A₂ Molekül - C Konstante in Gl. (20) - c₁=1/(2C) Konstante in Gl. (35) - c_p spezifische Wärme bei konstantem Druck - D binärer Diffusionskoeffizient - Ec=u ² /(2h_f) Eckert-Zahl - h spezifische Enthalpie - h_t=h+u²/2 totale spezifische Enthalpie - h _A ⁰ spezifische Dissoziationsenthalpie - K_w Reaktionsgeschwindigkeitskonstante der heterogenen Wandreaktion - 1= /( ) Champman-Rubesin-Parameter - Le=Pr/Sc Lewis-Zahl - M Molmasse - p statischer Druck - Pr= c_pf/ Prandtl-Zahl - q_w Wärmestromdichte an der Wand - q_cw, q_dw Wärmeleitungsbzw. Diffusionsanteil der Wärmestromdichte an der Wand - universelle Gaskonstante - R=/(2M_a) individuelle Gaskonstante der molekularen Komponente - Re_x= u x/ Reynolds-Zahl - Sc=/( D) Schmidt-Zahl - T absolute Temperatur - T_d=h _A ⁰ /R charakteristische Dissoziationstemperatur - u, v x- und y-Komponenten der Geschwindigkeit - U=u/u normierte x-Komponente der Geschwindigkeit - x, y Koordinaten parallel und senkrecht zur Platte Griechische Symbole - =A/ Dissoziationsgrad - Grenzschichtdicke - ₂ Impulsverlustdicke - Damköhler-Zahl der Oberflächenreaktion - =T/T normierte Temperatur - =y/ normierter Wandabstand - Wärmeleitfähigkeit - dynamische Viskosität - , ^* Ähnlichkeitskoordinaten - Dichte - Schubspannung Indizes A auf ein Atom bezogen - M auf ein Molekül bezogen - f auf den eingefrorenen Zustand bezogen - w auf die Wand bezogen - auf den Außenrand der Grenzschicht bezogen     

Thermal free convection from a solid sphere

Summary Thermal free convection from a sphere has been studied by melting solid benzene spheres in excess liquid benzene (Pr=8,3; 10⁸<Gr<10⁹). Overall heat transfer as well as local heat transfer were investigated. For the effect of cold liquid produced by the melting a correction has been applied. Results are compared with those obtained by other workers who used alternative experimental methods.Nomenclature coefficient of heat transfer - d characteristic length, here diameter of sphere - thermal conductivity - g acceleration of free fall - cubic expansion coefficient - T temperature difference between wall and fluid at infinity - kinematic viscosity - density - c specific heat capacity - a thermal diffusivity (=/c) - D diffusion coefficient - Nu dimensionless Nusselt number (=d/) - Nu* the analogous number for mass transfer (=kd/D) - mean value of Nusselt number - Gr dimensionless Grashof number (=gd ³T/ ²) - Gr* the analogous number for mass transfer (=gd ³x/ ²) - Pr dimensionless Prandtl number (=/a) - Sc dimensionless Schmidt number (=/D)     

Momentum and energy balances for dispersed two-phase flow

Summary The equations of motion and the mechanical energy balances for two-phase flow systems are derived by integration over a volume containing a large number of elements of the dispersed phase.List of symbols A, A boundary of volumes V, V - dA, dA surface element of A, A - A _s boundary of particles in V - dA _s surface element of A _s - F force per unit volume of the system - g –gz=gravity vector - g acceleration by gravity - I unit tensor - p pressure - Q dissipation in the continuous phase - Q _s dissipation in the dispersed phase - R compression work in the continuous phase - R _s compression work in the dispersed phase - t time - u velocity of continuous phase - u _s velocity of dispersed phase - u magnitude of u - u _s magnitude of u _s - V volume in the two-phase system - V part of V occupied by the continuous phase - W work done by F - z vertical coordinate - local volume fraction of the dispersed phase - pI–=stress tensor of the continuous phase - _s turbulent particle stress tensor - density of the continuous phase - _s density of the dispersed phase - shearing-stress tensor of the continuous phase - _s turbulent particle shearing-stress tensor - nabla operator - u, u _s velocity gradient tensor - substantial derivative (Shell Internationale Research Maatschappij N.V.)(Bataafse Internationale Petroleum Maatschappij N.V.)     

The theory of a vibrating-rod viscometer

The paper presents a complete theory for a viscometer based upon the principle of a circular-section rod, immersed in a fluid, performing transverse oscillations perpendicular to its axis. The theory is established as a result of a detailed analysis of the fluid flow around the rod and is subject to a number of criteria which subsequently constrain the design of an instrument. Using water as an example it is shown that a practical instrument can be designed so as to enable viscosity measurement with an accuracy of ±0.1%, although it is noted that many earlier instruments failed to satisfy one or more of the newly-established constraints.Nomenclature A, D constants in equation (46) - A _m , B _m , C _m , D _m constants in equations (50) and (51) - A _j , B _j constants in equation (14) - a _j ⁺ , a _j ^– wavenumbers given by equation (15) - C _f drag coefficient defined in equation (53) - c speed of sound - D _b drag force of fluid b - D ₀ coefficient of internal damping - E extensional modulus - f(z) initial deformation of rod - f(), F _m () functions of defined in equation (41) - F force in the rod - force per unit length near t=0 - F dimensionless force per unit length near t=0 - g _m amplitude of transient force - G modulus of rigidity - h, h* functions defined by equations (71) and (72) - H functions defined by equation (69) and (70) - I second moment of area - I _0,1, J _0,1, K _0,1 modified Bessel functions - k, k functions defined in equations (2) - L half-length of oscillator - Ma Mach number - m _b added mass per unit length of fluid b - m _s mass per unit length of solid - n _j eigenvalue defined in equations (15) and (16) - R radius of rod - R _c radius of container - 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 - y ₀ initial lateral displacement - y ₁, y ₂ successive maximum lateral displacement - z axial coordinate - dimensionless tension - dimensionless mass of fluid - dimensionless drag of fluid - amplification factor - logarithmic decrement in a fluid - _a , _b logarithmic decrement in fluids a and b - ₀ logarithmic decrement in vacuo - _j logarithmic decrement in mode j in a fluid - spatial resolution of amplitude - _v voltage resolution - r, , , _s, , increments in R, , , _s , , - dimensionless amplitude of oscillation - dimensionless axial coordinate - angular coordinate - _f thermal conductivity of fluid - viscosity of fluid - viscosity of fluid calculated on assumption that * - _a , _b viscosity of fluids a and b - _m constants in equation (10) - dimensionless displacement - _j j the component of - density of fluid - _a , _b density of fluids a and b - _s density of tube or rod material - dimensionless radial coordinate - * dimensionless radius of container - dimensionless times - spatial component of defined in equation (11) - _j , _tm jth, mth component of - dimensionless streamfunction - ⁰, ¹ components of in series expansion in powers of - streamfunction - dimensionless frequency (based on ) - angular frequency - ₀ angular frequency in absence of fluid and internal damping - _j angular frequency in mode j in a fluid - _a , _b frequencies in fluids a and b     

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.     

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     

The dynamic performance of the Weissenberg rheogoniometer III. Large amplitude oscillatory shearing — harmonic analysis

The harmonic content of the nonlinear dynamic behaviour of 1% polyacrylamide in 50% glycerol/water was studied using a standard Model R 18 Weissenberg Rheogoniometer. The Fourier analysis of the Oscillation Input and Torsion Head motions was performed using a Digital Transfer Function Analyser.In the absence of fluid inertia effects and when the amplitude of the (fundamental) Oscillation Input motion _I is much greater than the amplitudes of the Fourier components of the Torsion Head motion _Tn empirical nonlinear dynamic rheological propertiesG _n (, ⁰),G _n (, ⁰) and/or _n (, ⁰), _n (, ⁰) may be evaluated without a-priori-knowledge of a rheological constitutive equation. A detailed derivation of the basic equations involved is presented.Cone and plate data for the third harmonic storage modulus (dynamic rigidity)G ₃ (, ⁰), loss modulusG ₃ (, ⁰) and loss angle ₃ (, ⁰) are presented for the frequency range 3.14 × 10^–2 1.25 × 10² rad/s at two strain amplitudes, _CP ⁰ = 2.27 and 4.03. Composite cone and plate and parallel plates data for both the third and fifth harmonic dynamic viscosities ₃ (, ⁰), _S (, ⁰) and dynamic rigiditiesG ₃ (, ⁰),G ₅ (, ⁰) are presented for strain amplitudes in the ranges 1.10 _CP ⁰ 4.03 and 1.80 _PP ⁰ 36 for a single frequency, = 3.14 × 10^–1 rad/s. Good agreement was obtained between the results from both geometries and the absence of significant fluid inertia effects was confirmed by the superposition of the data for different gap widths.     

Flow of second-order fluids over an enclosed rotating disc

Summary This paper is devoted to a study of the flow of a second-order fluid (flowing with a small mass rate of symmetrical radial outflow m, taken negative for a net radial inflow) over a finite rotating disc enclosed within a coaxial cylinderical casing. The effects of the second-order terms are observed to depend upon two dimensionless parameters ₁ and ₂. Maximum values ₁ and ₂ of the dimensionless radial distances at which there is no recirculation, for the cases of net radial outflow (m>0) and net radial inflow (m<0) respectively, decrease with an increase in the second-order effects [represented by T(=₁+₂)]. The velocities at ₁ and ₂ as well as at some other fixed radii have been calculated for different T and the associated phenomena of no-recirculation/recirculation discussed. The change in flow phenomena due to a reversal of the direction of net radial flow has also been studied. The moment on the rotating disc increases with T.Nomenclature , , z coordinates in a cylindrical polar system - z ₀ distance between rotor and stator (gap length) - =/z ₀, dimensionless radial distance - =z/z ₀, dimensionless axial distance - _s = _s/z₀, dimensionless disc radius - V =(u, v, w), velocity vector - dimensionless velocity components - uniform angular velocity of the rotor - , p fluid density and pressure - P =p/(₂ z ₀₂ ² , dimensionless pressure - ₁, ₂, ₃ kinematic coefficients of Newtonian viscosity, elastico-viscosity and cross-viscosity respectively - ₁, _{2 2}/z ₀ ² , resp. ₃/z ₀ ² , dimensionless parameters representing the ratio of second-order and inertial effects - m = , mass rate of symmetrical radial outflow - l a number associated with induced circulatory flow - Rm =m/(z ₀₁), Reynolds number of radial outflow - R _l =l/(z ₀₁), Reynolds number of induced circulatory flow - Rz =z ₀ ² /₁, Reynolds number based on the gap - ₁, ₂ maximum radii at which there is no recirculation for the cases Rm>0 and Rm<0 respectively - ₁(T), ₂(T) ₁ and ₂ for different T - U _1(T) ⁽⁺⁾ = dimensionless radial velocity, Rm>0 - V _1(T) ⁽⁺⁾ = , dimensionless transverse velocity, Rm>0 - U _2(T) ^(–) = , dimensionless radial velocity, Rm=–Rn<0, m=–n - V _2(T) ^(–) = , dimensionless transverse velocity, Rm<0 - C _m moment coefficient     

A visualization study of the interaction of a free vortex with the wake behind an airfoil

The two-dimensional interaction of a single vortex with a thin symmetrical airfoil and its vortex wake has been investigated in a low turbulence wind tunnel having velocity of about 2 m/s in the measuring section. The flow Reynolds number based on the airfoil chord length was 4.5 × 10³. The investigation was carried out using a smoke-wire visualization technique with some support of standard hot-wire measurements. The experiment has proved that under certain conditions the vortex-airfoil-wake interaction leads to the formation of new vortices from the part of the wake positioned closely to the vortex. After the formation, the vortices rotate in the direction opposite to that of the incident vortex.List of symbols c test airfoil chord - C vortex generator airfoil chord - TA test airfoil - TE test airfoil trailing edge - TE _G vortex generator airfoil trailing edge - t dimensionless time-interval measured from the vortex passage by the test airfoil trailing edge: gDt=(T-T- _TE)·U/c - T time-interval measured from the start of VGA rotation - U free stream velocity - U vortex induced velocity fluctuation - VGA vortex generator airfoil - y distance in which the vortex passes the test airfoil - Z vortex circulation coefficient: Z=/(U · c/2) - vortex generator airfoil inclination angle - vortex circulation - vortex strength: =/2     

A distributed electrical analog for transient diffusion-techniques of model fabrication and calibration

Simulation of transient two-dimensional diffusion by means of a distributed electrical analog is discussed. After present techniques of analog model construction and calibration are reviewed, an improved calibration technique is presented and a convenient method of analog fabrication, not previously reported, is described. The proposed new method allows complete access to any point on the analog model during a test. Frequency response and step response measurements indicate that an adequate simulation is provided by this particular type of analog model.

---

Zusammenfassung Die Simulierung eines nichtstationären, zweidimensionalen Diffusionsvorganges mittels eines kontinuierlichen, elektro-thermischen Analogapparates wird besprochen. Eine Übersicht der gegenwärtigen Methoden für die Konstruktion und Kalibrierung von elektro-thermischen Analogmodellen wird gegeben. Ein verbessertes Verfahren für die Kalibrierung und eine handliche Fertigungsmethode von Analogmodellen, die noch nicht in der Literatur beschrieben wurden, werden dargestellt. Das vorgeschlagene neue Verfahren gestattet vollständigen Zugang zu jedem Punkt im Analogmodell während des Experiments. Meßbeobachtungen des periodischen Frequenzverhaltens und des nichtstationären Verhaltens zeigen, daß dieses spezielle Analogmodell eine ausreichende Simulierung des Diffusionsvorganges gestattet.

---

Nomenclature A ₀ Potential amplitude atx=0 (see Fig. 5) - A _x Potential amplitude at locationx (see Fig. 5) - A _L Potential amplitude atx=L (see Fig. 5) - C Capacitance per unit area - j Frequency - L Characteristic length - R Resistance per square - t Time - x Coordinate - X Dimensionless distance,x/L - y Coordinate - Y Dimensionless distance,y/L - Diffusivity - _x Phase angle at locationx (see Fig. 5) - _L Phase angle atx=L (see Fig. 5) - _d Thickness of dielectric sheet - _r Thickness of resistance sheet - Dielectric constant - Resistivity - Dimensionless time,t/L ² - Potential - ₁ Reference potential - ₂ Reference potential - Dimensionless potential, ( – ₁)/( ₂ – ₁) - Angular frequency, 2f - Dimensionless frequency,L ²/ The investigation was performed while the first author was Visiting Associate Professor at Purdue University during 1967/68.     

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/