Qualitative behavior of dissipative wave equations on bounded domains

The qualitative behavior of solutions of the mixed problem u_tt = u-a(x)u_t in IR x , u=0 on IR x , is studied in the case when a>0 and IRⁿ is bounded. Roughly speaking, if aa_min>0, then solutions decay at least as fast as exp t( –1/2a_min), with the possible exception of a finite dimensional set of smooth solutions whose existence is associated with a phenomenon of overdamping. If a_max is sufficiently small, depending on , then no overdamping occurs.Partially supported by NSF grant NSF GP 34260.This work was partially supported by the National Science Foundation under Grant No. GP 34260     

Experimental research on the laws of propagation of an underexpanded jet in a subsonic counterplan

The fields are obtained of total pressures and velocities along the axis of a jet flowing out into a subsonic counterflow from an isolated conical nozzle with semiaperture angle _a = 10°. The nozzle was located at the vertex of a cone with half-angle = 75° at the vertex and ratio of the radii at the midsection and the nozzle equal to 10. It is established that for an underexpanded jets in sub- and supersonic counterflows the main dimensionless number generalizing the data on the distribution of the gas-dynamic parameters in the jet at various pressure-ratio numbers n, and Mach numbers of the counterflow M and of the jet M_a, is the ratio of the momentum flux density J_a of the jet to the velocity head q: J=J_a/q. Generalizing dependences are obtained for the distribution of the total pressure and the velocity along the axis of the jet, and also for the jet range.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 186–188, May–June, 1986.     

Eine Methode zur Bestimmung des Transportkoeffizienten bei nichtlinearen Transportvorgängen

Zusammenfassung Es wird gezeigt, daß man nichtlineare Transportkoeffizienten aus Potentialprofilmessungen auch ohne die Lösung der nichtlinearen Transportgleichung bestimmen kann. Die Methode wird auf zwei Ziegeltone angewendet, deren Feuchtigkeitsdiffusionskoeffizienten recht komplizierte Funktionen des Feuchtegehaltes sind. Die Ergebnisse werden mit den Ergebnissen aus einer numerischen Lösung für einen empirischen Ansatz nach Fujita verglichen.

---

A method for evaluation of conductive properties in nonlinear transport phenomena

A prediction procedure for the conductive properties based upon the measured data on potential field has been developed. The method is particularly related to the nonlinear phenomena when the analytical solution is not available. The approach proposed is applied to the determination of the moisture diffusivity for two kinds of brick clay that posses a pronounced dependance of properties on moisture content. The results are used for Fujita's numerical solution which is compared with measured moisture field. The comparison has been performed and an agreement was found.

---

Bezeichnungen u [kg_w/kg_s] Feuchtigkeit bezogen auf trockene Substanz - w [kg_w/kg_s] nach Gl. (10) - V allgemeines Potential - _s, _so [kg_s/m³] Dichte des trockenen Skeletts - D, D₀ [m²/s] isotherme Feuchteleitkoeffizient; Diffusionskoeffizient - x [m] Koordinate längs Zylinder - , ₀ [s] Zeit - J_w [kg_w /m²s] Feuchtestrom - , Koeffizienten der allgemeinen Leitungsgleichung (4) - A, B, Konstanten der Gl. (6) - , Konstanten der Gl. (7) auch der Gl. (12) - d [m] Zylinderdurchmesser - PE Petrovaradin - SK Sremski Karlovci     

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

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

A three-parameter model describing the behaviour of a viscoelastic liquid in a tangential annular flow

A three-parameter model describing the shear rate-shear stress relation of viscoelastic liquids and in which each parameter has a physical significance, is applied to a tangential annular flow in order to calculate the velocity profile and the shear rate distribution. Experiments were carried out with a 5000 wppm aqueous solution of polyacrylamide and different types of rheometers. In a shear-rate range of seven decades (5 10^–3 s^–1 < < 1.2 10⁵ s^–1) a good agreement is obtained between apparent viscosities calculated with our model and those measured with three different types of rheometers, i.e. Couette rheometers, a cone-and-plate rheogoniometer and a capillary tube rheometer. a physical quantity defined by:a = {1 – ( / ₀)}/ ₀ (Pa^–1) - C constant of integration (1) - r distancer from the center (m) - r ₁,r ₂ radius of the inner and outer cylinder (m) - v _r local tangential velocity at a distancer from the center (v _r = _r r) (m s^–1) - v ₂ local tangential velocity at a distancer ₂ from the center (m s^–1) - shear rate (s^–1) - local shear rate (s^–1) - ₁ wall shear rate at the inner cylinder (s^–1) - dynamic viscosity (Pa s) - _a apparent viscosity (_a = / ) (Pa s) - _a1 apparent viscosity at the inner cylinder (Pa s) - ₀ zero-shear viscosity (Pa s) - infinite-shear viscosity (Pa s) - shear stress (Pa) - _r local shear stress at a distancer from the center (Pa) - ₀ yield stress (Pa) - ₁, ₂ wall shear-stress at the inner and outer cylinder (Pa) - _r local angular velocity (s^–1) - ₂ angular velocity of the outer cylinder (s^–1)     

Pressure pulsations during interaction between a subsonic jet and the subsonic section of a supersonic jeton a flat obstacle

Pressure pulsations were measured during in-leakage of a subsonic jet and the subsonic section of a heated supersonic jet on a flat obstacle. Data have been obtained on the total and spectrum levels of the pressure pulsations at different spacings X of the obstacle from the nozzle exit. It is shown that when the obstacle is disposed at the section of the jet where the local velocity is subsonic, the pulsation levels outside the dependence on the conditions at the nozzle exit (Mach number Maxa 0 a 3.0; stagnation temperature T₀=280–1200K) vary in direct proportion to the local velocity head q. The ratio between the total level and q is (/g)=0.2–0.3. It is established that for a subsonic velocity ahead of the obstacle, all the spectra obtained for different values of M _a , T₀, d _a and X in the coordinates Sh=f(d/V) and (₁*/q)(V/d) will lie on a single generalized spectrum. Here ₁* is the pulsation level in a 1-Hz band, and d and V are, respectively, the jet diameter and velocity directly in front of the obstacle.Translated from Izvestiya Akademiya Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 172–174, September–October, 1975.     

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     

Determination of the steady-state shear viscosity from measurements of the apparent viscosity for some common types of viscometers

Considering a number of model fluids, the relation between the (measurable) apparent viscosity _a and the (true) shear viscosity is studied for some commonly used viscometers, like capillary, slit, plate-plate and concentric cylinders (including the influence of the bottom of the cylinder), as well as for one laboratory type of viscometer. As long as is a purely monotonic function, a shift factor < 1 allows one to deduce from _a . Though in general variable, it frequently suffices for practical purposes to use a constant shift factor (the constant being characteristic of the type of viscometer used). This does not apply to dilute solutions or any fluids with two plateau values for . For plastic fluids, it is shown that Casson or Bingham behavior can — if valid at all — only describe the high shear stress limit of _a .     

A Strongly Degenerate Quasilinear Equation: the Parabolic Case

We prove the existence and uniqueness of entropy solutions of the Neumann problem for the quasilinear parabolic equation u_t=÷ a(u, Du), where a(z,)=f(z,), and f is a convex function of with linear growth as ||||, satisfying other additional assumptions. In particular, this class includes the case where f(z,)=(z)(), >0, and is a convex function with linear growth as ||||.     

On laminar flow through a uniformly porous pipe

Numerous investigations ([1] and [4–9]) have been made of laminar flow in a uniformly porous circular pipe with constant suction or injection applied at the wall. The object of this paper is to give a complete analysis of the numerical and theoretical solutions of this problem. It is shown that two solutions exist for all values of injection as well as the dual solutions for suction which had been noted by previous investigators. Analytical solutions are derived for large suction and injection; for large suction a viscous layer occurs at the wall while for large injection one solution has a viscous layer at the centre of the channel and the other has no viscous layer anywhere. Approximate analytic solutions are also given for small values of suction and injection.

Nomenclature

General r distance measured radially - z distance measured along axis of pipe - u velocity component in direction of z increasing - v velocity component in direction of r increasing - p pressure - density - coefficient of kinematic viscosity - a radius of pipe - V velocity of suction at the wall - r ²/a ² - R wall or suction Reynolds number, Va/ - f() similarity function defined in (6) - u ₀() eigensolution - U(0) a velocity at z=0 - K an arbitrary constant - B _K Bernoulli numbers Particular Section 5 perturbation parameter, –2/R - ² a constant, –K - x / - g(x) f()/ Section 6 perturbation parameter, –R/2 - ² a constant, –K - g() f() - g _c ()=g() near centre of pipe - * point where g()=0 Section 7 2/R - ² K - t (1–)/ - w(t, ) [1–f(t)]/ - ₀, ₁ constants - g() f()– ₀ - ₀/ - ₀ a constant - * point where f()=0     

Flow through a helically coiled annulus

In this paper the flow is studied of an incompressible viscous fluid through a helically coiled annulus, the torsion of its centre line taken into account. It has been shown that the torsion affects the secondary flow and contributes to the azimuthal component of velocity around the centre line. The symmetry of the secondary flow streamlines in the absence of torsion, is destroyed in its presence. Some stream lines penetrate from the upper half to the lower half, and if is further increased, a complete circulation around the centre line is obtained at low values of for all Reynolds numbers for which the analysis of this paper is valid, being the ratio of the torsion of the centre line to its curvature.Nomenclature A =constant - a outer radius of the annulus - b unit binormal vector to C - C helical centre line of the pipe - D rL - g 1000 - K Dean number=Re² - L 1+r sin - M (L ²+ ² r ²)^1/2 - n unit normal vector to C - P, P pressure and nondimensional pressure - p ₀, p pressures of O(1) and O() - Re Reynolds number=aW ₀/ - (r, , s), (r, , s) coordinates and nondimensional coordinates - nonorthogonal unit vectors along the coordinate directions - r ₀ radius of the projection of C - t unit tangent vector to C - V _r, V , V _s velocity components along the nonorthogonal directions - V_r, V, V _s nondimensional velocity components along - W ₀ average velocity in a straight annulus Greek symbols , curvature and nondimensional curvature of C - U, V, W lowest order terms for small in the velocity components along the orthogonal directions t - _r, , _s first approximations to V _r , V, V _s for small - =/=/ - kinematic viscosity - density of the fluid - , torsion and nondimensional torsion of C - , stream function and nondimensional stream function - nondimensional streamfunction for U, V - a inner radius of the annulus After this paper was accepted for publication, a paper entitled On the low-Reynolds number flow in a helical pipe, by C.Y. Wang, has appeared in J. Fluid. Mech., Vol 108, 1981, pp. 185–194. The results in Wangs paper are particular cases of this paper for =0, and are also contained in [9].     

Determination of permeability tensors for two-phase flow in homogeneous porous media: Theory

In this paper we continue previous studies of the closure problem for two-phase flow in homogeneous porous media, and we show how the closure problem can be transformed to a pair of Stokes-like boundary-value problems in terms of pressures that have units of length and velocities that have units of length squared. These are essentially geometrical boundary value problems that are used to calculate the four permeability tensors that appear in the volume averaged Stokes' equations. To determine the geometry associated with the closure problem, one needs to solve the physical problem; however, the closure problem can be solved using the same algorithm used to solve the physical problem, thus the entire procedure can be accomplished with a single numerical code.Nomenclature a a vector that maps V onto , m^-1. - A a tensor that maps V onto . - A area of the - interface contained within the macroscopic region, m². - A area of the -phase entrances and exits contained within the macroscopic region, m². - A area of the - interface contained within the averaging volume, m². - A area of the -phase entrances and exits contained within the averaging volume, m². - Bo Bond number (= (=(–)g²/). - Ca capillary number (= v/). - g gravitational acceleration, m/s². - H mean curvature, m^-1. - I unit tensor. - permeability tensor for the -phase, m². - viscous drag tensor that maps V onto V. - ^* dominant permeability tensor that maps onto v , m². - ^* coupling permeability tensor that maps onto v , m². - characteristic length scale for the -phase, m. - l characteristic length scale representing both and , m. - L characteristic length scale for volume averaged quantities, m. - n unit normal vector directed from the -phase toward the -phase. - n unit normal vector representing both n and n . - n unit normal vector representing both n and n . - P pressure in the -phase, N/m². - p superficial average pressure in the -phase, N/m². - p intrinsic average pressure in the -phase, N/m². - p –p , spatial deviation pressure for the -phase, N/m². - r ₀ radius of the averaging volume, m. - r position vector, m. - t time, s. - v fluid velocity in the -phase, m/s. - v superficial average velocity in the -phase, m/s. - v intrinsic average velocity in the -phase, m/s. - v –v , spatial deviation velocity in the -phase, m/s. - V volume of the -phase contained within the averaging volmue, m³. - averaging volume, m³. Greek Symbols V /, volume fraction of the -phase. - viscosity of the -phase, Ns/m². - density of the -phase, kg/m³. - surface tension, N/m. - (v +v ^T ), viscous stress tensor for the -phase, N/m².     

Investigation of supersonic viscous flow past a sphere in the presence of subsonic and sonic injection

Supersonic flow past a sphere with a given rate of gas injection along the generator is investigated numerically on the range Re=10²–10⁴. Calculations have been made on the interval 0 90°, where is the angle between the axis of symmetry and the normal to the surface. It is shown that for high subsonic and sonic injection rates it is possible to observe qualitatively new features in the flow structure and in the distribution of the local supersonic flow characteristics around the perimeter of the sphere not previously noted in [9]. In the case of sonic injection the changes in flow structure occur only in the supersonic zone. In the neighborhood of the transition from a subsonic to sonic injection velocity the heat flux has a local maximum, which in absolute value does not exceed the heat flux in the absence of injection. It is shown that there may be qualitative differences in the pressure distribution over the surface of the body with increase in the injection parameter depending on the distribution and value of the injected gas flow rate and, moreover, the number Re.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 83–89, January–February, 1988.     

Two-phase flow in heterogeneous porous media II: Numerical experiments for flow perpendicular to a stratified system

In this work, we make use of numerical experiments to explore our original theoretical analysis of two-phase flow in heterogeneous porous media (Quintard and Whitaker, 1988). The calculations were carried out with a two-region model of a stratified system, and the parameters were chosen be consistent with practical problems associated with groundwater flows and petroleum reservoir recovery processes. The comparison between theory (the large-scaled averaged equations) and experiment (numerical solution of the local volume averaged equations) has allowed us to identify conditions for which the quasi-static theory is acceptable and conditions for which a dynamic theory must be used. Byquasi-static we mean the following: (1) The local capillary pressure,everywhere in the averaging volume, can be set equal to the large-scale capillary pressure evaluated at the centroid of the averaging volume and (2) the large-scale capillary pressure is given by the difference between the large-scale pressures in the two immiscible phases, and is therefore independent of gravitational effects, flow effects and transient effects. Bydynamic, we simply mean a significant departure from the quasi-static condition, thus dynamic effects can be associated with gravitational effects, flow effects and transient effects. To be more precise about the quasi-static condition we need to refer to the relation between the local capillary pressure and the large-scale capillary pressure derived in Part I (Quintard and Whitaker, 1990). Herep _c ¦_y represents the local capillary pressure evaluated at a positiony relative to the centroid of the large-scale averaging volume, and {p _c }¦_x represents the large-scale capillary pressure evaluated at the centroid.In addition to{p _c } ^c being evaluated at the centroid, all averaged terms on the right-hand side of Equation (1) are evaluated at the centroid. We can now write the equations describing the quasi-static condition as , , This means that the fluids within an averaging volume are distributed according to the capillary pressure-saturation relationwith the capillary pressure held constant. It also means that the large-scale capillary pressure is devoid of any dynamic effects. Both of these conditions represent approximations (see Section 6 in Part I) and one of our main objectives in this paper is to learn something about the efficacy of these approximations. As a secondary objective we want to explore the influence of dynamic effects in terms of our original theory. In that development only the first four terms on the right hand side of Equation (1) appeared in the representation for the local capillary pressure. However, those terms will provide an indication of the influence of dynamic effects on the large-scale capillary pressure and the large-scale permeability tensor, and that information provides valuable guidance for future studies based on the theory presented in Part I.Roman Letters A scalar that maps {}*/t onto - A scalar that maps {}*/t onto - A interfacial area between the -region and the -region contained within, m² - A interfacial area between the -region and the -region contained within, m² - A interfacial area between the -region and the -region contained within, m² - a vector that maps ({}*/t) onto , m - a vector that maps ({}*/t) onto , m - b vector that maps ({p}– g) onto , m - b vector that maps ({p}– g) onto , m - B second order tensor that maps ({p}– g) onto , m² - B second order tensor that maps ({p}– g) onto , m² - c vector that maps ({}*/t) onto , m - c vector that maps ({}*/t) onto , m - C second order tensor that maps ({}*/t) onto , m² - C second order tensor that maps ({}*/t) onto . m² - D third order tensor that maps ( ) onto , m - D third order tensor that maps ( ) onto , m - D second order tensor that maps ( ) onto , m² - D second order tensor that maps ( ) onto , m² - E third order tensor that maps () onto , m - E third order tensor that maps () onto , m - E second order tensor that maps () onto - E second order tensor that maps () onto - p _c =(), capillary pressure relationship in the-region - p _c =(), capillary pressure relationship in the-region - g gravitational vector, m/s² - largest of either or - - - i unit base vector in thex-direction - I unit tensor - K local volume-averaged-phase permeability, m² - K local volume-averaged-phase permeability in the-region, m² - K local volume-averaged-phase permeability in the-region, m² - {K } large-scale intrinsic phase average permeability for the-phase, m² - K –{K }, large-scale spatial deviation for the-phase permeability, m² - K –{K }, large-scale spatial deviation for the-phase permeability in the-region, m² - K –{K }, large-scale spatial deviation for the-phase permeability in the-region, m² - K ^* large-scale permeability for the-phase, m² - L characteristic length associated with local volume-averaged quantities, m - characteristic length associated with large-scale averaged quantities, m - I _i i = 1, 2, 3, lattice vectors for a unit cell, m - l characteristic length associated with the-region, m - ; characteristic length associated with the-region, m - l _H characteristic length associated with a local heterogeneity, m - - n unit normal vector pointing from the-region toward the-region (n =–n ) - n unit normal vector pointing from the-region toward the-region (n =–n ) - p pressure in the-phase, N/m² - p local volume-averaged intrinsic phase average pressure in the-phase, N/m² - {p } large-scale intrinsic phase average pressure in the capillary region of the-phase, N/m² - p local volume-averaged intrinsic phase average pressure for the-phase in the-region, N/m² - p local volume-averaged intrinsic phase average pressure for the-phase in the-region, N/m² - p –{p }, large scale spatial deviation for the-phase pressure, N/m² - p –{p }, large scale spatial deviation for the-phase pressure in the-region, N/m² - p –{p }, large scale spatial deviation for the-phase pressure in the-region, N/m² - P _c p –{p }, capillary pressure, N/m² - {p_c}^c large-scale capillary pressure, N/m² - r ₀ radius of the local averaging volume, m - R ₀ radius of the large-scale averaging volume, m - r position vector, m - , m - S /, local volume-averaged saturation for the-phase - S ^* {}*{}*, large-scale average saturation for the-phaset time, s - t time, s - u , m - U , m² - v -phase velocity vector, m/s - v local volume-averaged phase average velocity for the-phase in the-region, m/s - v local volume-averaged phase average velocity for the-phase in the-region, m/s - {v } large-scale intrinsic phase average velocity for the-phase in the capillary region of the-phase, m/s - {v } large-scale phase average velocity for the-phase in the capillary region of the-phase, m/s - v –{v }, large-scale spatial deviation for the-phase velocity, m/s - v –{v }, large-scale spatial deviation for the-phase velocity in the-region, m/s - v –{v }, large-scale spatial deviation for the-phase velocity in the-region, m/s - V local averaging volume, m³ - V volume of the-phase in, m³ - V large-scale averaging volume, m³ - V capillary region for the-phase within, m³ - V capillary region for the-phase within, m³ - V _c intersection of m³ - V volume of the-region within, m³ - V volume of the-region within, m³ - V () capillary region for the-phase within the-region, m³ - V () capillary region for the-phase within the-region, m³ - V () , region in which the-phase is trapped at the irreducible saturation, m³ - y position vector relative to the centroid of the large-scale averaging volume, m Greek Letters local volume-averaged porosity - local volume-averaged volume fraction for the-phase - local volume-averaged volume fraction for the-phase in the-region - local volume-averaged volume fraction for the-phase in the-region - local volume-averaged volume fraction for the-phase in the-region (This is directly related to the irreducible saturation.) - {} large-scale intrinsic phase average volume fraction for the-phase - {} large-scale phase average volume fraction for the-phase - {}* large-scale spatial average volume fraction for the-phase - –{}, large-scale spatial deviation for the-phase volume fraction - –{}, large-scale spatial deviation for the-phase volume fraction in the-region - –{}, large-scale spatial deviation for the-phase volume fraction in the-region - a generic local volume-averaged quantity associated with the-phase - mass density of the-phase, kg/m³ - mass density of the-phase, kg/m³ - viscosity of the-phase, N s/m² - viscosity of the-phase, N s/m² - interfacial tension of the - phase system, N/m - , N/m - , volume fraction of the-phase capillary (active) region - , volume fraction of the-phase capillary (active) region - , volume fraction of the-region ( + =1) - , volume fraction of the-region ( + =1) - {p }– g, N/m³ - {p }– g, N/m³     

On the effects of uniform high suction on the rotationally symmetric flow of a conducting liquid near a stationary disc in the presence of a transverse magnetic field

In the present paper an attempt has been made to find out effects of uniform high suction in the presence of a transverse magnetic field, on the motion near a stationary plate when the fluid at a large distance above it rotates with a constant angular velocity. Series solutions for velocity components, displacement thickness and momentum thickness are obtained in the descending powers of the suction parameter a. The solutions obtained are valid for small values of the non-dimensional magnetic parameter m (= 4 _e ² H ₀ ² /) and large values of a (a2).Nomenclature a suction parameter - E electric field - E _r , E , E _z radial, azimuthal and axial components of electric field - F, G, H reduced radial, azimuthal and axial velocity components - H magnetic field - H _r , H , H _z radial, azimuthal and axial components of magnetic field - H ₀ uniform magnetic field - H* displacement thickness and momentum thickness ratio, */ - h induced magnetic field - h _r , h , h _z radial, azimuthal and axial components of induced magnetic field - J current density - m nondimensional magnetic parameter - p pressure - P reduced pressure - R Reynolds number - U ₀ representative velocity - V velocity - V _r , V , V _z radial, azimuthal and axial velocity components - w ₀ uniform suction through the disc. - density - electrical conductivity - kinematic viscosity - _e magnetic permeability - a parameter, (/)^1/2 z - a parameter, a - * displacement thickness - momentum thickness - angular velocity     

Theoretische Untersuchung der Strömung und des Wärmetransports bei erzwungener turbulenter Konvektion im Rechteckkanal unter Berücksichtigung von Sekundäreffekten

Zusammenfassung Im Bereich des abgeschlossenen Einlaufs eines Rechteckkanals wird bei erzwungener turbulenter Konvektion das Strömungsfeld und das Temperaturfeld berechnet. Die Analyse behandelt die Einflüsse der Richtungsabhängigkeit der turbulenten Austauschgrößen und der infolge der Anisotropie der Turbulenz existierenden Sekundärströmung auf die axiale Komponente des Strömungsfeldes und auf das Temperaturfeld. Die thermischen Effekte der Kanalwand sowie Wärmetransport durch Leitung und Temperaturstrahlung werden berücksichtigt. Die Analyse zeigt, daß nur eine einzige Drehung der Sekundärströmung innerhalb eines TrapezSymmetrieelementes des Rechteckkanals existiert. Die Sekundäreffekte bewirken, daß bei einseitiger oder zweiseitig gegenüberliegender Beheizung die maximale Wandtemperatur von der Kanal ecke zur Mitte der beheizten Wand verschoben wird.

---

An analytical study of momentum and heat transfer on turbulent forced convection in rectangular channels with allowance for secondary effects

A theoretical study was performed to investigate turbulent forced-convective momentum and heat transport in a rectangular channel under fully developed flow and heat transfer conditions. Main emphasis has been devoted to analyse the effects of the anisotropic turbulent transport properties and the turbulence-induced secondary flow on the main flow field and the temperature distribution. The effects of peripheral wall conduction as well as radiation within the channel are included in the analysis. The analysis reveals that only a single secondary current occurs in the trapezoidal symmetry element of the rectangular duct. Furthermore, when the heating extends over one or two oppositely located sides only, the location of the maximum wall temperature is shifted from the corner to the center of the wall.

---

Formelzeichen a_i halbe Seitenlänge des Kanals (i=2,3), a₂a₃ - a Funktion - D Durchmesser - e Basis des natürlichen Logarithmus - L turbulentes Längenmaß - 1 normiertes turbulentes Längenmaß, L/R_H - In natürlicher Logarithmus - Nu Nusselt-Zahl - Pe Péclet-Zahl - p Funktion - Q zugeführte Wärme - q_i Wärmefluß (i=1, 2, 3) - R Radius - Re Reynolds-Zahl - Sf Stefan-Zahl - T Temperatur - U Umfang - u_i normierte Geschwindigkeit w_i/w_m, (i=1, 2, 3) - W Strahlungsfluß - W Wandstärke des Kanals - w_i Geschwindigkeit (i=1, 2, 3) - w Schubspannungsgeschwindigkeit an der Wand - x_i kartesische Koordinaten (i=1, 2, 3) - y _i ⁺ normierter Wandabstand, (¦a_i – x_i¦)/R_H (i=2, 3) - _i normierte halbe Seitenlänge des Kanals, ₂=(1 + a₂/a₃)/2, ₃=(1 + a₃/a₂)/2 - i Emissionskoeffizient der Strahlung - _ijkl turbulenter Austauschtensor für Impuls - _qij turbulenter Austauschtensor für Wärme - normierte Temperatur, (T – T_m)_F/(Q_WD_H) - Wandleitungsparameter, W_w/R_HF - Wärmeleitfähigkeit - kinematische Viskosität - _i normierte kartesische Koordinaten, x_i/R_H (i=1, 2, 3) - Dichte - Stefan-Boltzmann-Strahlungskonstante - _W Wandschubspannung - Turbulenzparameter, w_mR_H/ - Stromfunktion - _i normierter Wirbelvektor, _ijku_k/_i Indizes und Überschreibungen a, b Bezeichnung von Flächenelementen - CL Kanalmitte - B Buleev - F Fluid - H hydraulisch - K Kanal - l laminar - m mittel - O Ort, an dem die Turbulenzkorrelation berechnet wird - q Wärme - t turbulent - w Wand - Impuls - - turbulenter Mittelwert - turbulente Fluktuation - pro Längeneinheit - · pro Zeiteinheit Teil der vom Fachbereich für Verfahrenstechnik der Technischen Universität Berlin genehmigten Dissertation des Verfassers.     

Thermal effects in the capillary rheometer

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

---

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

---

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

Supersonic nonequilibrium flow past segmental bodies of a mixture simulating the atmosphere of venus

Calculations of the flow of the mixture 0.94 CO₂+0.05 N₂+0.01 Ar past the forward portion of segmentai bodies are presented. The temperature, pressure, and concentration distributions are given as a function of the pressure ahead of the shock wave and the body velocity. Analysis of the concentration distribution makes it possible to formulate a simplified model for the chemical reaction kinetics in the shock layer that reflects the primary flow characteristics. The density distributions are used to verify the validity of the binary similarity law throughout the shock layer region calculated.The flow of a CO₂+N₂+Ar gas mixture of varying composition past a spherical nose was examined in [1]. The basic flow properties in the shock layer were studied, particularly flow dependence on the free-stream CO₂ and N₂ concentration.New revised data on the properties of the Venusian atmosphere have appeared in the literature [2, 3] One is the dominant CO₂ concentration. This finding permits more rigorous formulation of the problem of blunt body motion in the Venus atmosphere, and attention can be concentrated on revising the CO₂ thermodynamic and kinetic properties that must be used in the calculation.The problem of supersonic nonequilibrium flow past a blunt body is solved within the framework of the problem formulation of [4].Notation V body velocity - shock wave standoff - universal gas constant - ratio of frozen specific heats - hRt/m enthalpy per unit mass undisturbed stream P pressure - density - T temperature - m molecular weight - c_p specific heat at constant pressure - (X) concentration of component X (number of particles in unit mass) - R body radius of curvature at the stagnation point - _j rate of j-th chemical reaction shock layer P V ² pressure - density - TT temperature - mm molecular weight Translated from Izv. AN SSSR. Mekhanika Zhidkosti i Gaza, Vol. 5, No. 2, pp. 67–72, March–April, 1970.The author thanks V. P. Stulov for guidance in this study.     

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     

Evaluation of the dynamic-mechanical properties of solids from a cooperative model for stress relaxation

For many solid materials the stress relaxation process obeys the universal relationF = – (d/d lnt)_max = (0.1 ± 0.01) ( ₀ — _i ), regardless of the structure of the material. Here denotes the stress,t the time, ₀ the initial stress of the experiment and _i the internal stress. A cooperative model accounting for the similarity in relaxation behaviour between different materials was developed earlier. Since this model has a spectral character, the concepts of linear viscoelasticity are used here to evaluate the corresponding prediction of the dynamic mechanical properties, i.e. the frequency dependence of the storageE () and lossE () moduli. Useful numerical approximations ofE () andE () are also evaluated. It is noted that the universal relation in stress relaxation had a counterpart in the frequency dependence ofE (). The theoretical prediction of the loss factor for high-density polyethylene is compared with experimental results. The agreement is good.