Propagation of a spherical wave in nonlinearly compressible and elastoplastic materials

The problem of spherical wave propagation in soil under the action of an intense uniformly decreasing load ₀(t) applied to the boundary of a cavity with radius r₀ is considered. Soil with a high stress level is modeled either by ideally nonlinearly compressible or elastoplastic material, taking account of linear irreversible unloading for the material. In contrast to [1–7], in order to describe material movement use is made of strain theory [8] with determining functions = (), _i=_i(_i), where , _i, , _i are the first and second invariants of strain and stress tensors. During material loading these functions are presented in the form of polynomials ()=(_i+₂¦¦), _ii)=(_i-_2i)_i, in which constant coefficients _i, _i=1, 2) are determined by experiment, taking account of the triaxial stressed state of soil. Solution of the problem is constructed by an analytically reversible method, with prescribed shape for the shock-wave (SW) surface in the form of a second-degree polynomial relating to time t and a numerical method of characteristics for a prescribed arbitrarily decreasing load _i(t). On the basis of the analytical equations obtained, calculations are carried out for material parameters (including loading profile) in a computer and stresses and mass velocity of plastic and elastoplastic materials are compared.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 95–100, July–August, 1986.The authors express their sincere thanks to Kh. A. Rakhmatulin for discussing the results of this work.     

The theory of a vibrating-rod densimeter

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

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

Zusammenfassung Zur Berechnung turbulenter Strömungen wird das k--Modell im Ansatz für die turbulente Scheinzähigkeit erweitert, so daß es den Querkrümmungs- und Dichteeinfluß auf den turbulenten Transportaustausch erfaßt. Die dabei zu bestimmenden Konstanten werden derart festgelegt, daß die bestmögliche Übereinstimmung zwischen Berechnung und Messung erzielt wird. Die numerische Integration der Grenzschichtgleichungen erfolgt unter Verwendung einer Transformation mit dem Differenzenverfahren vom Hermiteschen Typ. Das erweiterte Modell wird auf rotationssymmetrische Freistrahlen veränderlicher Dichte angewendet und zeigt Übereinstimmung zwischen Rechnung und Experiment.

---

On the influence of transvers-curvature and density in inhomogeneous turbulent free jets

The prediction of turbulent flows based on the k- model is extended to include the influence of transverse-curvature and density on the turbulent transport mechanisms. The empirical constants involved are adjusted such that the best agreement between predictions and experimental results is obtained. Using a transformation the boundary layer equations are solved numerically by means of a finite difference method of Hermitian type. The extended model is applied to predict the axisymmetric jet with variable density. The results of the calculations are in agreement with measurements.

---

Bezeichnungen Wirbelabsorptionskoeffizient - c_i Massenkonzentration der Komponente i - c_D, c_L, c, c₁, c₂ Konstanten des Turbulenzmodells - d Düsendurchmesser - E bezogene Dissipationsrate - f bezogene Stromfunktion - f Korrekturfunktion für die turbulente Scheinzähigkeit - j turbulenter Diffusionsstrom - k Turbulenzenergie - k_i Schrittweite in -Richtung - K dimensionslose Turbulenzenergie - L turbulentes Längenmaß - M_i Molmasse der Komponente i - p Druck - allgemeine Gaskonstante - r Querkoordinate - r_0,5 Halbwertsbreite der Geschwindigkeit - r_0,5c Halbwertsbreite der Konzentration - T Temperatur - u Geschwindigkeitskomponente in x-Richtung - v Geschwindigkeitskomponente in r-Richtung - x Längskoordinate - y allgemeine Funktion - Y_i diskreter Wert der Funktion y - Relaxationsfaktor für Iteration - turbulente Dissipationsrate - transformierte r-Koordinate - kinematische Zähigkeit - Exponent - transformierte x-Koordinate - Dichte - _k, Konstanten des Turbulenzmodells - Schubspannung - allgemeine Variable - Stromfunktion - Turbulente Transportgröße Indizes 0 Strahlanfang - m auf der Achse - r mit Berücksichtigung der Krümmung - t turbulent - mit Berücksichtigung der Dichte - im Unendlichen - Schwankungswert oder Ableitung einer Funktion - – Mittelwert Herrn Professor Dr.-Ing. R. Günther zum 70. Geburtstag gewidmet     

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     

Transport in ordered and disordered porous media V: Geometrical results for two-dimensional systems

In this paper we continue the geometrical studies of computer generated two-phase systems that were presented in Part IV. In order to reduce the computational time associated with the previous three-dimensional studies, the calculations presented in this work are restricted to two dimensions. This allows us to explore more thoroughly the influence of the size of the averaging volume and to learn something about the use of anon-representative region in the determination of averaged quantities.

Nomenclature

Roman Letters A interfacial area of the– interface associated with the local closure problem, m² - a _i i=1, 2, gaussian probability distribution used to locate the position of particles - l unit tensor - characteristic length for the-phase particles, m - ⁰ reference characteristic length for the-phase particles, m - characteristic length for the-phase, m - _i i=1,2,3 lattice vectors, m - m convolution product weighting function - m _V special convolution product weighting function associated with a unit cell - n _i i=1, 2 integers used to locate the position of particles - n unit normal vector pointing from the-phase toward the-phase - r _p position vector locating the centroid of a particle, m - r gaussian probability distribution used to determine the size of a particle, m - r ₀ characteristic length of an averaging region, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume,V, m³ - x position of the centroid of an averaging area, m - x ₀ reference position of the centroid of an averaging area, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters V /V, volume average porosity - _{a i} standard deviation ofa _i - _r standard deviation ofr - intrinsic phase average of      

On the cross-effects of coupled macroscopic transport equations in porous media

In this paper, the derivation of macroscopic transport equations for this cases of simultaneous heat and water, chemical and water or electrical and water fluxes in porous media is presented. Based on themicro-macro passage using the method of homogenization of periodic structures, it is shown that the resulting macroscopic equations reveal zero-valued cross-coupling effects for the case of heat and water transport as well as chemical and water transport. In the case of electrical and water transport, a nonsymmetrical coupling was found.Notations b mobility - c concentration of a chemical - D rate of deformation tensor - D molecular diffusion coefficient - D _ij ^eff macroscopic (or effective) diffusion tensor - electric field - E ₀ initial electric field - k _ij molecular tensor - j, j ^*, current densities - K _ij macroscopic permeability tensor - l characteristic length of the ERV or the periodic cell - L characteristic macroscopic length - L _ijkl coupled flows coefficients - n _i unit outward vector normal to - p pressure - q _t ,q _t ⁺ , heat fluxes - q _c ,q _c ⁺ , chemical fluxes - s specific entropy or the entropy density - S entropy per unit volume - t time variable - t _ij local tensor - T absolute temperature - v _i velocity - V ₀ initial electric potential - V electric potential - x macroscopic (or slow) space variable - y microscopic (or fast) space variable - _i local vectorial field - _i local vectorial field - electric charge density on the solid surface - , bulk and shear viscosities of the fluid - _ij local tensor - _ij local tensor - _i local vector - _ij molecular conductivity tensor - _ij ^eff effective conductivity tensor - homogenization parameter - fluid density - ₀ ion-conductivity of fluid - _ij dielectric tensor - _i ¹ , _i ² , _i ³ local vectors - ⁴ local scalar - _S solid volume in the periodic cell - _L volume of pores in the periodic cell - boundary between _S and _L - ^s rate of entropy production per unit volume - total volume of the periodic cell - _l volume of pores in the cell On leave from the Politechnika Gdanska; ul. Majakowskiego 11/12, 80-952, Gdask, Poland.     

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/     

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     

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     

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     

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     

Reflection, transmission and excitation of SH-surface waves by a discontinuity in mass-loading on a semi-infinite elastic medium

The scattering of an SH-wave by a discontinuity in mass-loading on a semi-infinite elastic medium is investigated theoretically. The incident wave is either a plane body wave or a plane SH-surface wave. The problem is reduced to a Wiener-Hopf problem for the scattered wave. In this problem the amplitude spectral density of the particle displacement occurs as unknown function. Special attention is given to the numerical values of the surface wave contributions to the scattered field.Nomenclature x ₁, x ₂, x ₃ Cartesian coordinates - , polar coordinates in x ₁, x ₃-plane - volume mass density - surface mass density of mass-loading - , Lamé constants - U scalar wave function, defined by (2.1) - c _S speed of propagation of uniform shear waves in bulk medium (c _S=(/)^1/2) - angular frequency - t time - k _S wave number of uniform shear waves (k _S=/c _S) - reduced specific acoustic impedance of mass-loading (=k _S /) - k _m wave number of SH-surface wave (k _m=k _S(1+ ²)^1/2) - _1,2,3 partial differentiation with respect to x _1,2,3 - ⁱ angle between x ₃-axis and direction of propagation of incident body wave - ⁱ wave number in horizontal direction of incident body wave ( ⁱ=k _S sin( ⁱ)) - ⁱ wave number in vertical direction of incident body wave ( ⁱ=k _S cos( ⁱ)) - C _1,2 complex amplitude of surface wave excited by a body wave - R reflection factor of surface wave, when surface wave is incident - T transmission factor of surface wave, when surface wave is incident - S particle displacement vector The research presented in this paper has been carried out with partial financial support from the Delfts Hogeschoolfonds.     

Mixed convection on vertical plate for fluids of any prandtl number

A mixed convection parameter=(Ra) ^1/4/(Re)^1/2, with=Pr/(1+Pr) and=Pr/(1 +Pr)^1/2, is proposed to replace the conventional Richardson number, Gr/Re², for combined forced and free convection flow on an isothermal vertical plate. This parameter can readily be reduced to the controlling parameters for the relative importance of the forced and the free convection,Ra ^1/4/(Re ^1/2 Pr ^1/3) forPr 1, and (RaPr)^1/2/(RePr ^1/2 forPr 1. Furthermore, new coordinates and dependent variables are properly defined in terms of, so that the transformed nonsimilar boundary-layer equations give numerical solutions that are uniformly valid over the entire range of mixed convection intensity from forced convection limit to free convection limit for fluids of any Prandtl number from 0.001 to 10,000. The effects of mixed convection intensity and the Prandtl number on the velocity profiles, the temperature profiles, the wall friction, and the heat transfer rate are illustrated for both cases of buoyancy assisting and opposing flow conditions.

---

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

Zusammenfassung Für die kombinierte Zwangs- und freie Konvektion an einer isothermen senkrechten Platte wird ein Mischkonvektions-Parameter=( Ra) ^1/4 (Re)^1/2, mit=Pr/(1 +Pr) und=Pr/(1 +Pr)^1/2 vorgeschlagen, den die gebräuchliche Richardson-Zahl, Gr/Re², ersetzen soll. Dieser Parameter kann ohne weiteres auf die maßgebenden Kennzahlen für den relativen Einfluß der erzwungenen und der freien Konvektion reduziert werden,Ra ^1/4/(Re ^1/2 Pr ^1/3) fürPr 1 und (RaPr)^1/4/(RePr)^1/2 fürPr 1. Weiterhin werden neue Koordinaten und abhängige Variablen als Funktion von definiert, so daß für die transformierten Grenzschichtgleichungen numerische Lösungen erstellt werden können, die über den gesamten Bereich der Mischkonvektion, von der freien Konvektion bis zur Zwangskonvektion, für Fluide jeglicher Prandtl-Zahl von 0.001 bis 10.000 gleichmäßig gültig sind. Der Einfluß der Intensität der Mischkonvektion und der Prandtl-Zahl auf die Geschwindigkeitsprofile, die Temperaturprofile, die Wandreibung und den Wärmeübergangskoeffizienten werden für die beiden Fälle der Strömung in und entgegengesetzt zur Schwerkraftrichtung dargestellt.

---

Nomenclature C _f local friction coefficient - C _p specific heat capacity - f reduced stream function - g gravitational acceleration - Gr local Grashoff number,g T _w –T )x³/v² - Nu local Nusselt number - Pr Prandtl number,v/ - Ra local Rayleigh number,g T _w –T x ³/( v) - Re local Reynolds number,u x/v - Ri Richardson number,Gr/Re ² - T fluid temperature - T _w wall temperature - T free stream temperature - u velocity component in thex direction - u free stream velocity - v velocity component in they direction - x vertical coordinate measuring from the leading edge - y horizontal coordinate Greek symbols thermal diffusivity - thermal expansion coefficient - mixed convection parameter (Ra)1/4/Re)^1/2 - pseudo-similarity variable,(y/x) - ₀ conventional similarity variable,(y/x)Re ^1/2 - dimensionless temperature, (T–T T _W –T - unified mixed-flow parameter, [(Re) ^1/2 + (Ra)^1/4] - dynamic viscosity - kinematic viscosity - stretched streamwise coordinate or mixed convection parameter, [1 + (Re)^1/2/(Ra) ^1/4]^–1=/(1 +) - density - Pr/(1 + Pr) _w wall shear stress - stream function - Pr/(l+Pr)^1/3 This research was supported by a grand from the National Science Council of ROC     

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     

Temperature distribution in a thin rotating disk

The problem of heat conduction in a thin rotating disk with heat input at a fixed point is considered. The disk is cooled by forced convection from its lateral surfaces. By defining a complex temperature, the temperature throughout the disk is presented as a series of Bessel functions of complex argument. Results are given for a range of rotational speeds.Nomenclature R radial coordinate - angular coordinate - a radius of disk - b thickness of disk - T temperature - T ambient temperature - rotational speed of disk - q heat flux into disk - k thermal conductivity of disk - density of disk - c specific heat of disk - h coefficient of convective heat transfer - r dimensionless radial coordinate, R/a - T* characteristic temperature, q ₀ a/ k - t dimensionless temperature, (T–T )/T* - C ₁, C ₂ dimensionless parameters defined in (3)     

On the closure problem for Darcy's law

In a previous derivation of Darcy's law, the closure problem was presented in terms of an integro-differential equation for a second-order tensor. In this paper, we show that the closure problem can be transformed to a set of Stokes-like equations and we compare solutions of these equations with experimental data. The computational advantages of the transformed closure problem are considerable.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 _e area of entrances and exits for the-phase contained within the averaging volume, m² - B second-order tensor used to respresent the velocity deviation - b vector used to represent the pressure deviation, m^–1 - C second-order tensor related to the permeability tensor, m^–2 - D second-order tensor used to represent the velocity deviation, m² - d vector used to represent the pressure deviation, m - g gravity vector, m/s² - I unit tensor - K C ^–1,–D, Darcy's law permeability tensor, m² - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the-phase, m - l _i i=1, 2, 3, lattice vectors, m - n unit normal vector pointing from the-phase toward the-phase - n _e outwardly directed unit normal vector at the entrances and exits of the-phase - p pressure in the-phase, N/m ² - p intrinsic phase average pressure, N/m² - p –p , spatial deviation of the pressure in the-phase, N/m² - r position vector locating points in the-phase, m - r ₀ radius of the averaging volume, m - t time, s - v velocity vector in the-phase, m/s - v intrinsic phase average velocity in the-phase, m/s - v phase average or Darcy velocity in the \-phase, m/s - v –v , spatial deviation of the velocity in the-phase m/s - V averaging volume, m³ - V volume of the-phase contained in 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²     

Unsteady convective diffusion in a rotating parallel plate channel

The paper presents an exact analysis of the dispersion of a passive contaminant in a viscous fluid flowing in a parallel plate channel driven by a uniform pressure gradient. The channel rotates about an axis perpendicular to its walls with a uniform angular velocity resulting in a secondary flow. Using a generalized dispersion model which is valid for all time, we evaluate the longitudinal dispersion coefficientsK _i (i=1, 2, ...) as functions of time. It is shown thatK ₁=0 andK ₃,K ₄, ... decay rapidly in comparison withK ₂. ButK ₂ decreases with increasing (the dimensionless rotation parameter) for values of upto approximately =2.2. ThereafterK ₂ increases with further increase in and its value gets saturated for large values of (say, 500) and does not change any further with increase in . A physical explanation of this anomalous behaviour ofK ₂ is given.

---

Instationäre konvektive Diffusion in einem rotierenden Parallelplattenkanal

Zusammenfassung In dieser Untersuchung wird eine exakte Analyse der Ausbreitung eines passiven Kontaminierungsstoffes in einer zähen Flüssigkeit gegeben, die, befördert durch einen gleichförmigen Druckgradienten, in einem Parallelplattenkanal strömt. Der Kanal rotiert mit gleichförmiger Winkelgeschwindigkeit um eine zu seinen Wänden senkrechte Achse, wodurch sich eine Sekundärströmung ausbildet. Unter Verwendung eines generalisierten, für alle Zeiten gültigen Dispersionsmodells werden die longitudinalen DispersionskoeffizientenK _i (i=1, 2, ...) als Funktionen der Zeit ermittelt. Es wird gezeigt, daßK ₁=0 gilt und dieK ₃,K ₄, ... gegenüberK ₂ schnell abnehmen.K ₂ nimmt ab, wenn , der dimensionslose Rotationsparameter, bis etwa zum Wert 2,2 ansteigt. Danach wächstK ₂ mit bis auf einem Endwert an, der etwa ab =500 erreicht wird. Dieses anomale Verhalten vonK ₂ findet eine physikalische Erklärung.

---

List of symbols C solute concentration - D molecular diffusivity - K _i longitudinal dispersion coefficients - 2L depth of the channel - P ₀ dimensionless pressure gradient along main flow - Pe Péclet number - q velocity vector - Q _x,Q _y mass flux along the main flow and the secondary flow directions - dimensionless average velocity along the main flow direction - (x, y, z) Cartesian co-ordinates Greek symbols dimensionless rotation parameter - the inclination of side walls withx-axis - kinematic viscosity - fluid density - dimensionless time - angular velocity of the channel - dimensionless distance along the main flow direction - dimensionless distance along the vertical direction - dimensionless solute concentration - integral of the dispersion coefficientK ₂() over a time interval     

Magnetohydrodynamic natural convection heat transfer from radiate vertical porous surfaces

Magnetohydrodynamic natural convection heat transfer from radiate vertical surfaces with fluid suction or injection is considered. The nonsimilarity parameter is found to be the conductive fluid injection or suction along the streamwise coordinate = V{4x/² g(T _w – T )}^1/4. Three dimensionless parameters had been found to describe the problem: the magnetic influence number N = B ² _y /V ², the radiation-conduction parameter R _d = k _R /4aT ³ , and the Gebhart number Ge _x = gx/cp to represent the effect of the viscous dissipation. It is found that increasing the magnetic field strength causes the velocity and the heat transfer rates inside the boundary layer to decrease. Its apparent that increasing the radiation-conduction parameter decreases the velocity and enhances the heat transfer rates. The Gebhart number, i.e, the viscous dissipation had no effect on the present problem.Nomenclature a Stefan-Boltzmann constant - B _y Magnetic field flux density Wb/m² - Cf _x Local skin friction factor - c _p Specific heat capacity - f Dimensionless stream function - Ge _x Gebhart number, gx/cp - g Gravitational acceleration - k Thermal Conductivity - L Length of the plate - N Magnetic influence number, B ² _y /V ² - p Pressure - Pr Prandtl number - q _r Radiative heat flux - q _w (x) Local surface heat flux - Q _w (x) Dimensionless Local surface heat flux - R _d Planck number (Radiation-Conduction parameter), k _R /4aT ³ - T Temperature - T Free stream temperature - T _w Wall temperature - u, v Velocity components in x- and y-directions - V Porous wall suction or injection velocity - V _w Porous wall suction or injection velocity - x, y Axial and normal coordinates - Thermal diffusivity Greek symbols _R Roseland mean absorption coefficient, 4/3R _d - Coefficient of thermal expansion - Nonsimilarity parameter, V{4x/² g(T _w – T )}^1/4 - Peseudo-similarity variable - Dimensionless temperature - _w Ratio of surface temperature to the ambient temperature, T _w /T - Dynamice viscosity - Kinemtic viscosity - Fluid density - Electrical conductivity - _w Local wall shear stress - Dimensional stream function     

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