Thermal and diffusive relaxation of an evaporating drop with internal heat liberation

The problem of nonsteady-state evaporation or growth of a radiating drop with uniformly distributed internal heat sources is considered. The Reynolds R=ua/v 1 and Peclet P_D= ua/D 1 numbers are assumed to be small (a is the radius of the drop, u the velocity of its relative motion, andv, D, the coefficients of viscosity, diffusion and thermal diffusivity of the vapor-gas medium). This enables the convective transfer of vapor and heat to be neglected, and the concentration and temperature fields to be regarded as spherically symmetric [1]. In view of the fact that the density of saturated vapor is less than the density of liquid the convective flow caused by the change in radius of the drop is not taken into account [2]. It has already been shown [3,4], that for _r (, _r are the coefficients of molecular and radiative thermal conductivity) there exists a bounded region ryo (1/) /_r ( is the absorption coefficient for radiation in the gas), in which the effect of radiation on the temperature relaxation of the vapor-gas medium is negligible. If the conditiona (1/) /_r is satisfied, then the temperature at the outer boundary of this region will be practically the same as the temperature at infinity T=T. This means that terms in the energy equation connected with energy transferred by radiation can be neglected. It is assumed that the free path of molecules in the gas is less than the radius of the drop, and so concentration and temperature discontinuities close to the surface of the drop can be neglected [2].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, Vol. 11, No. 1, pp. 78–87, January–February, 1970.The authors are grateful to V. G. Levich for discussing the results of the paper.     

Free-forced convection from a heated longitudinal horizontal cylinder embedded in a porous medium

The flow and heat transfer from a heated semi-infinite horizontal circular cylinder which is moving with a constant speed into a porous medium is considered. It is assumed that the Grashof and Reynolds numbers are large so that the governing equations are the three dimensional boundary-layer equations. A numerical procedure for solving these equations is described and the asymptotic solutions which are valid both near and distant from the leading edge of the cylinder are presented. The range of validity of these asymptotic solutions is discussed and the results are compared in detail with the full numerical solution. The problem is of practical importance, for example in the drilling of pipes into a geothermal reservoir.

---

Freie erzwungene Konvektion von einem beheizten schlanken horizontalen Zylinder, eingebettet in ein poröses Medium

Zusammenfassung Es wird die Strömung und der Wärmeübergang an einem beheizten, halbunendlichen horizontalen Kreiszylinder betrachtet, der mit konstanter Geschwindigkeit sich in ein poröses Medium bewegt. Dabei wird angenommen, daß die Grashof- und Reynolds-Zahlen groß sind, so daß die Bestimmungsgleichungen von den dreidimensionalen Grenzschichtgleichungen gebildet werden. Es wird ein numerisches Verfahren zur Lösung dieser Gleichungen beschrieben und eine asymptotische Lösung präsentiert, die sowohl in der Nähe als auch in großem Abstand von dem vorderen Ende des Zylinders gültig ist. Der Gültigkeitsbereich dieser asymptotischen Lösungen wird diskutiert und die Ergebnisse werden im Detail mit vollständigen numerischen Lösungen verglichen. Das Problem ist z.B. beim Eindringen von Rohrleitungen in geothermische Reservoire von praktischer Wichtigkeit.

---

Nomenclature a radius of cylinder - Gr Grashof number (=g(T_w-Ta/²) - g acceleration due to gravity - permeability in the porous medium - Nu local Nusselt number - total heat flux from cylinder - q _w heat flux from cylinder - r radial co ordinate - Ra Rayleigh number (=g (T_w - T_t8) a/ ) - Re Reynolds number (=U _t8 a/) - T temperature - u, v, w speeds inx, , r directions - x axial co ordinate - equivalent thermal diffusivity - thermal expansion coefficient - ratioGr/Re - similarity variable - dimensionless temperature (=(T- T)/(T _w- T) - kinematic viscosity - azimuthal co ordinate - w cylinder surface - free stream     

Anwendung des Übertragungsverfahrens bei stark abklingenden Lösungsfunktionen

Übersicht Bei stark abklingenden Funktionen wird die Übertragungsmatrix U() aufgespalten in die Anteilc U ₁() e^– und U ₂() e. Der zweite Term spielt am Rand = 0 keinc Rolle. Die unbekannten Anfangswerte sind über die Matrix U ₁(0) an die bekannten gebunden und eindeutig bestimmbar.

---

Summary For strongly decaying solution functions the transfer matrix U() is splitted into the parts U ₁() e^– and U ₂() e. The second term does not influence at the boundary = 0. The unknown initial values are related by the matrix U ₁(0) to the known values and they can be uniquely determined.

     

Calculation of shear influences on the phase separation of polymer blends exhibiting upper critical solution temperatures

Calculations were performed on the basis of a generalized Gibbs energy of mixing G , which is the sum of the Gibbs energy of mixing of the stagnant system and E _s, the energy stored in the system during stationary flow. With increasing shear rate , the demixing temperatures shift to lower values (shear-induced mixing; diminution of the heterogeneous area), then to higher values (shear-induced demixing), and finally to lower values again before the effects fade out. The details of the rather complex phase diagrams resulting for a given shear rate are primarily determined by a band in the T/ plane ( = mole fraction) within which (² E _s/²) _T <0 (i.e., E _S acts towards phase separation). There are two ranges of within which closed miscibility gaps can exist: The more common outer islands are partly or totally situated outside the equilibrium gap (and within the above mentioned band). As is raised they break away from the mainland at the upper end of the first region of shear-induced mixing and shift to T>UCST where they submerge. Bound to a suitable choice of parameters, a second kind of closed miscibility gaps, the inner islands, which always remain within the equilibrium solubility gap (and outside the band of negative curvature of E _S) is additionally observed. This time the islands break away from the mainland at the lower end of the first region of shear-induced mixing where they also submerge. The present findings are compared with the results of previous calculations for LCSTs.     

Mass transfer with chemical reaction in a laminar falling film

Zusammenfassung Es wird eine analytische Lösung für die Absorption in einem laminaren Rieselfilm mit homogener und heterogener chemischer Reaktion 1. Ordnung vorgestellt, wobei der Stofftransportwiderstand auf der Gasseite liegt. Die Lösung ist eine Funktion von drei dimensionslosen ParameternBi, und, welche die BiotZahl und einen homogenen bzw. heterogenen Reaktionsparameter darstellen. Es wird gezeigt, daß für feste Werte vonBi und die Absorptionsrate (bezogen auf die Breite 1 des Rieselfilms) über eine gewisse Länge (dimensionslos) des Rieselfilms unabhängig von ist, wenn, < 0,6 ist. Die laufende Länge wird von der Stelle aus gemessen, an der die Absorption beginnt. Für b 0,6 nimmt der FlußQ mit zu, erreicht aber einen Sättigungswert bei=10, wonachQ nurmehr sehr langsam anwächst. Jedoch für ein gegebenes und ohne Übergangswiderstand im Film (Bi ) nimmtQ mit für alle 0 zu.

---

Mass transfer with chemical reaction in a laminar falling film

An analytical solution is presented for gas absorption in a laminar falling film with first-order homogeneous and heterogeneous chemical reaction and external gas-phase mass transfer resistance. The solution depends on three dimensionless parametersBi, and, wich represent the Biot number, homogeneous and heterogeneous reaction parameters, respectively. It is shown that for fixed values ofBi and, the rate of gas absorption (per unit breadth) over a certain length; (dimensionless) along the falling film measured from the point where surface absorption begins is independent of if < 0.6. For 0.6, this fluxQ increases with but reaches a saturation value at=10 beyond whichQ increases very slowly. But for given and zero gas film resistance (Bi ),Q increases with for all 0.

     

Nonisothermal non-Newtonian chemical reactions in a tubular reactor under conditions of variable viscosity

The nonisothermal non-Newtonian chemical reactions in a tubular reactor are investigated. The non-Newtonian fluid is assumed to be characterized by the Ostwald-de Waele power-law model, which represents the majority of laminar flow of food products and many polymer melts and solutions. The temperature effect on the viscosity is considered and is found to be very significant. The effects of other important dimensionless parameters on the reactor performances are examined.Nomenclature c mass fraction of reactant - c ₀ inlect mass fraction of reactant - C _p heat capacity, J/kg K - C dimensionless concentration of reactant, c/c ₀ - C _b dimensionless bulk concentration of reactant - D molecular diffusivity, m²/s - E activation energy, J/kg - H heat of reaction, J/m³ - k ₁ frequency factor, s^–1 - k _t heat conductivity, J/m K kg - K fluid consistency, kg s ^n–1/m - K ₁ dimensionless frequency factor, k ₁ r ₀ ² c ^m–1 exp(–₁)/D - K ₀ constant in Eq. (6) - m order of chemical reaction - n rheological parameter - p pressure, kg/m s² - r radial coordinate, m - r ₀ radius of reactor, m - R dimensionless radial coordinate, r/r ₀ - R _g gas constant, J/kg K - T temperature, K - T ₀ inlet temperature, K - u velocity, m/s - u _b bulk velocity, m/s - U dimensionless velocity, u/u _b - x axial coordinate, m - X dimensionless axial coordinate, xD/r ₀ ² u _b Greek symbols dimensionless parameter, - dimensionless parameter, ₀T₀ - ₁ dimensionless activation energy, E/R _g T ₀ - ₂ dimensionless heat generation, - dimensionless temperature, (T–T ₀)/T ₀ - _b dimensionless bulk temperature - liquid density, kg/m³     

Die Torsionsspannungen in Wellen mit tiefen Längsnuten

Übersicht MitF(x, y) als Spannungsfunktion einer Welle ohne Nut und(, y) als Potentialfunktion des Quelle-Senke-Systems erhält man Spannungsfunktionen(, y) =F(x, y) –(, y) für Wellen mit tiefen Längsnuten. Es wird gezeigt, daß sich damit die Schubspannungen in den Läufern von Schraubenverdichtern ermitteln lassen.

---

Shearing stresses in shafts with deep longitudinal grooves

Summary The stress functions(, y) of shafts with deep longitudinal grooves may be represented by(, y) =F(x, y) –(, y) whereF(x, y) is the stress function of a cylindrical shaft without grooves and(, y) denotes the potential function of the source-sink system. It is shown that the shearing stresses in rotors of screw-compressors may be obtained in this way.

     

Flow of rarefied gas past a sweptback cylinder

A numerical investigation has been made of the hypersonic flow of a rarefied monatomic gas past the windward part of the side surface of an infinite circular cylinder. The calculation was made by direct statistical Monte Carlo modeling for freestream Mach number M_t8=20, ratio of the surface temperature of the body to the stagnation temperature equal to t_tw =T _tw/T _t0 = 0.03, sweep angle 75°, and Reynolds number Re_t0 30.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 146–154, January–February, 1992.     

On the Ellipticity of the Middle Surface of a Shell and its Application to the Asymptotic Analysis of ‘Membrane’ Shells

We consider a surface S = (), where ² is a bounded, connected, open set with a smooth boundary and : ³ is a smooth map; let () denote the components of the two-dimensional linearized strain tensor of S and let 0 with length 0 > 0. We assume the the norm ,|| ()||0, in the space V0() = { H¹() × H¹() × L²(); = 0 on 0 } is equivalent to the usual product norm on this space. We then establish that this assumption implies that the surface S is uniformly elliptic and that we necessarily have 0 = .     

Heteroclinic orbits in singular systems: A unifying approach

We consider singularly perturbed systems , such that=f(, _o, 0). _o ^m , has a heteroclinic orbitu(t). We construct a bifurcation functionG(, ) such that the singular system has a heteroclinic orbit if and only ifG(, )=0 has a solution=(). We also apply this result to recover some theorems that have been proved using different approaches.     

The Conjugacy of Stochastic and Random Differential Equations and the Existence of Global Attractors

We consider stochastic differential equations in d-dimensional Euclidean space driven by an m-dimensional Wiener process, determined by the drift vector field f₀ and the diffusion vector fields f₁,...,f_m, and investigate the existence of global random attractors for the associated flows . For this purpose is decomposed into a stationary diffeomorphism given by the stochastic differential equation on the space of smooth flows on R^d driven by m independent stationary Ornstein Uhlenbeck processes z¹,...,z^m and the vector fields f₁,...,f_m, and a flow generated by the nonautonomous ordinary differential equation given by the vector field (_t/x)^–1[f₀(_t)+ _i=1 ¹ f_i(_t)z _t ⁱ ]. In this setting, attractors of are canonically related with attractors of . For , the problem of existence of attractors is then considered as a perturbation problem. Conditions on the vector fields are derived under which a Lyapunov function for the deterministic differential equation determined by the vector field f₀ is still a Lyapunov function for , yielding an attractor this way. The criterion is finally tested in various prominent examples.     

Diffuse approximation method for solving natural convection in porous media

The diffuse approximation is presented and applied to natural convection problems in porous media. A comparison with the control volume-based finite-element method shows that, overall, the diffuse approximation appears to be fairly attractive.Nomenclature H height of the cavities - I functional - K permeability - p(M _i ,M) line vector of monomials - p ^T p-transpose - M current point - Nu Nusselt number - Ri inner radius - Ro outer radius - Ra Rayleigh number - x, y cartesian coordinates - u, v velocity components - T temperature - M vector of estimated derivatives - _t thermal diffusivity - coefficient of thermal expansion - practical aperture of the weighting function - scalar field - (M, M _i ) weighting function - streamfunction - kinematic viscosity     

An optical method for the measurement of unsteady film thickness

This paper describes a measurement technique to quantify temporal variations in the thickness of an unsteady liquid film with a resolution which is independent of the thickness. The optical transformation function has been derived for fringes of equal inclination and, for a temporally varying film, allows the unsteady component of film thickness to be measured in terms of frequency modulated signal analysis of light intensity variations. As it does not require calibration, the method is suited to in-situ measurements of complex and rapidly varying films as encountered in engineering two-phase flows. It requires the inversion of the frequency time-series of the light intensity observed by a photodetector which represents the absolute values of the time-derivative of the thickness variation.The technique has been used to measure the thickness of the film formed as a result of impingement of a pulsating two-phase jet onto a heated flat plate with surface temperatures of 150 °C and 240 °C and located 143 nozzle exit diameters downstream of the nozzle. The angle between the jet axis and the surface normal was 20 degrees and the injection frequency was 16.7 Hz corresponding to a flow rate of 7.2 mm³ per injection. The results along the line of incidence showed that the ensemble-averaged space-time structure of the film was qualitatively independent of the plate temperature with three peaks, two of which occurred at large radial distances and disappeared in less than 10 ms. The third peak was close to the impingement region and persisted for more than 50 ms due to the small velocities of the incoming two-phase jet as the nozzle needle closed and the low momentum wall jet which was unable to transport the droplets radially outwards. At the higher surface temperature, the rate of evaporation and the amplitude variation of the unsteady component of the overall film thickness increased, and the film covered a smaller area.List of symbols A amplitude of electric field - E electric field produced at a point - f _i (t) instantaneous frequency - f _PMT (t) frequency observed by photomultiplier tube - h thickness resolution - h _m minimum thickness that can be measured - h (t, x) unsteady film thickness - h _s (x) steady film thickness - h _t (t,x) overall film thickness - H nozzle-to-plate distance - i - I light intensity at a point - interference term - k wavenumber of the monochromatic point source - n refractive index of the air - refractive index of thin film material - r radial distance from the geometrical impingement point - rms root-mean-square - t time - t ₀ zero-crossing time defined in Fig. 2 - t _r rise-time defined in Fig. 2 - T _w wall temperature - x position vector in 3-D space Greek symbols angle of impingement - angle of reflection - integration constant defined in Eq. (15) - _c integration constant at r = 0 - phase angle - P optical path difference - P minimum optical path difference - azimuthal coordinate defined in Fig. 3 a - initial phase angle defined as = h _s (x) - wavelength of the illuminating source - initial phase of an electric field produced by a source - _f (t) ensemble-averaged rms series for f _PMT (t) - angle of incidence - angle of refraction - angular frequency of the electric field produced by a source - optical transformation function     

A class of nonlinear,completely integrable abstract wave equations

IfL is a positive self-adjoint operator on a Hubert spaceH, with compact inverse, the second-order evolution equation int,u+Lu+u _H ² u=0 has an infinite number of first integrals, pairwise in involution. It follows from this that no nontrivial solution tends weakly to 0 inH ast. Under an additional separation assumption on the eigenvalues ofL, all trajectories (u,u) are relatively compact inD(L ^1/2)×H. Finally, if all the eigenvalues are simple, the set of initial values of quasi-periodic solutions is dense in the ball B=(u ₀,u ₀ )D(L ^1/2)×H; L^1/2 u _{0 H} ² +u ² < for sufficiently small.     

A class of hypo-elastic non-elastic materials and their thermodynamics

Thermodynamics is developed for a class of thermo-hypo-elastic materials. It is shown that materials of this class obey the laws of thermodynamics, but are not elastic.

Table of Symbols

Latin Letters A _ijkl tensor-valued function of t _ij appearing in hypo-elastic constitutive relation - B _ijkl another tensor-valued function. See equation (4.2) - B the square of - d _ij rate of deformation tensor - d _ij deviator of rate of deformation - f, k functions of pressure, p - g, h functions of the invariant - p pressure - q _i heat flux vector - s _ij stress deviator - _ij co-rotational derivative of stress deviator - t time - t ₁ t ₂ specific values of time - t _ij stress tensor - t _ij ⁰ a specific value of stress - T Temperature - T ₀ a specific value of temperature - u _i velocity - V(t) a material volume as a function of time, t - V ₀ a material volume at a reference configuration - W work (W = work done in a deformation—section 5) Sript Letters Specific internal energy - Specific Helmholtz free energy - G Specific Gibbs function Greek Letters an invariant of the stress deviator—see eq. (2.4) - _ij kroneker delta - (W = work done in a deformation—section 5) - specific entropy - hypo-elastic potential - hypo-elastic potential - mass density - ₀ mass density in a reference configuration - specific volume = 1/ - a function of p - _ijkl a constant tensor—see eq. (2.5) - –G/ - _ij rate of rotation tensor This work is dedicated to Jerald L. Ericksen, without whose influence it would not have been possible     

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

Modelling of heat transfer by conduction in a transition region between a porous medium and an external fluid

We study the modelling of purely conductive heat transfer between a porous medium and an external fluid within the framework of the volume averaging method. When the temperature field for such a system is classically determined by coupling the macroscopic heat conduction equation in the porous medium domain to the heat conduction equation in the external fluid domain, it is shown that the phase average temperature cannot be predicted without a generally negligible error due to the fact that the boundary conditions at the interface between the two media are specified at the macroscopic level.Afterwards, it is presented an alternative modelling by means of a single equation involving an effective thermal conductivity which is a function of point inside the interfacial region.The theoretical results are illustrated by means of some numerical simulations for a model porous medium. In particular, temperature fields at the microscopic level are presented.Roman Letters _sf interfacial area of thes-f interface contained within the macroscopic system m² - A _sf interfacial area of thes-f interface contained within the averaging volume m² - C _p mass fraction weighted heat capacity, kcal/kg/K - g vector that maps to _s , m - h vector that maps to _f , m - K _eff effective thermal conductivity tensor, kcal/m s K - l _s,l _f microscopic characteristic length m - L macroscopic characteristic length, m - n _fs outwardly directed unit normal vector for thef-phase at thef-s interface - n outwardly directed unit normal vector at the dividing surface. - R ₀ REV characteristic length, m - T _i macroscopic temperature at the interface, K - error on the external fluid temperature due to the macroscopic boundary condition, K - T ^* macroscopic temperature field obtained by solving the macroscopic Equation (3), K - V averaging volume, m³ - V _s,V _f volume of the considered phase within the averaging volume, m³. - _mp volume of the porous medium domain, m³ - _ex volume of the external fluid domain, m³ - _s , _f volume of the considered phase within the volume of the macroscopic system, m³ - dividing surface, m² - x, z spatial coordinates Greek Letters _s, _f volume fraction - ratio of the effective thermal conductivity to the external fluid thermal conductivity - ^* macroscopic thermal conductivity (single equation model) kcal/m s K - _s, _f microscopic thermal conductivities, kcal/m s K - spatial average density, kg/m³ - microscopic temperature, K - ^* microscopic temperature corresponding toT ^*, K - spatial deviation temperature K - error in the temperature due to the macroscopic boundary conditions, K - ^* _i macroscopic temperature at the interface given by the single equation model, K - spatial average - ^s , ^f intrinsic phase average.     

Some remarks on differential equations of quadratic type

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

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³     

Stabilitätsrechnung der Rohrströmung einer inkompressiblen viskoelastischen Flüssigkeit zweiter Ordnung

Zusammenfassung Für die Durchsatzströmung im Rohr wird mit Hilfe der klassischen hydrodynamischen Stabilitätstheorie gezeigt, daß die inkompressible Flüssigkeit zweiter Ordnungs = –pI + 2(d + 2t ₁ d ² –t ₀ d) stabil ist gegenüber kleinen rotationssymmetrischen Störungen.

---

Summary For Poiseuille pipe flow it is shown by means of the classical theory of hydrodynamic stability, that the incompressible second-order fluids = –pI + 2(d + 2t ₁ d ² –t ₀ d) is stable with respect to small disturbances of rotational symmetry.

---

Nomenklatur a _n Koeffizienten der Reihenentwicklung - c = /k komplexe Wellengeschwindigkeit - d Deformationsgeschwindigkeitstensor - D, D dimensionsloser Deformationsgeschwindigkeitstensor (Grund- und Störtensor) - e _i kovariante Basis - g Vektor der Erdbeschleunigung - I Einheitstensor - k Wellenzahl - M, O, S, Q, T Funktion vonk, Re, ₀ - p, P, p Gesamt-, Grund-, Stördruck - r, (r, , z) dimensionsloser Ortsvektor (Zylinderkoordinaten) - R Rohrradius - Re =U _M R/ Reynoldszahl - s(s ^*=s*pI) Spannungstensor (Isotroper Anteil des ) - t ₀,t ₁ Stoffzeiten, Parameter der Flüssigkeit zweiter Ordnung - t Zeit - u, U, u Vektor der Gesamt-, Grund-, Störgeschwindigkeit - U _M Maximale Grundgeschwindigkeit - v, V, v Vektor der dimensionslosen Gesamt-, Grund-, Störgeschwindigkeit - w Rotationsgeschwindigkeitstensor - W, W Rotationsgeschwindigkeitstensor, dimensionslos (Grund-, Störtensor) - x (x _r ,x ,x _z ) Ortsvektor (Zylinderkoordinaten) - Viskosität - ₀, ₁ dimensionslose Stoffzeiten - dimensionsloser Druck - Dichte - dimensionslose Zeit - Stromfunktion, dimensionslos - komplexe Frequenz, dimensionslos - = e ⁱ /x _i Nablaoperator (e ⁱ kontravariante Basis) - _* Nablaoperator, dimensionslos - R, I Real-, Imaginärteil Mit 4 Abbildungen