Mixed convection flow in narrow vertical ducts

The mixed convection flow in a vertical duct is analysed under the assumption that , the ratio of the duct width to the length over which the wall is heated, is small. It is assumed that a fully developed Poiseuille flow has already been set up in the duct before heat from the wall causes this to be changed by the action of the buoyancy forces, as measured by a buoyancy parameter . An analytical solution is derived for the case when the Reynolds numberRe, based on the duct width, is of 0 (1). This is extended to the case whenRe is 0 (^–1) by numerical integrations of the governing equations for a range of values of representing both aiding and opposing flows. The limiting cases, || 1 andR=Re of 0 (1), andR and both large, with of 0 (R ^1/3) are considered further. Finally, the free convection limit, large with R of 0 (1), is discussed.

---

Mischkonvektion in engen senkrechten Rohren

Zusammenfassung Mischkonvektion in einem senkrechten Rohr wird unter der Voraussetzung untersucht, daß das Verhältnis der Rohrbreite zur Länge, über welche die Wand beheizt wird, klein ist. Es wird angenommen, daß sich bereits eine voll entwickelte Poiseuille-Strömung in dem Rohr eingestellt hat, bevor Antriebskräfte, gemessen mit dem Auftriebsparameter , aufgrund der Wandbeheizung die Strömung verändern. Es wird eine analytische Lösung für den Fall erhalten, daß die mit der Rohrbreite als charakteristische Länge gebildete Reynolds-ZahlRe konstant ist. Dies wird mittels einer numerischen Integration der wichtigsten Gleichungen auf den FallRe =f (^–1) sowohl für Gleich- als auch für Gegenstrom ausgedehnt. Weiterhin werden die beiden Grenzfälle betrachtet, wenn || 1 undR=Re konstant ist, sowieR und beide groß mit proportionalR ^1/3. Schließlich wird der Grenzfall der freien Konvektion, großes mit konstantem R, diskutiert.

---

Nomenclature g acceleration due to gravity - Gr Grashof number - G modified Grashof number - h duct width - l length of the heated section of the duct wall - p pressure - Pr Prandtl number - Q flow rate through the duct - Q ₀ heat transfer on the wally=0 - Q ₁ heat transfer on the wally=1 - Re Reynolds number - R modified Reynolds number - T temperature of the fluid - T ₀ ambient temperature - T applied temperature difference - u, velocity component in thex-direction - v, velocity component in they-direction - x, co-ordinate measuring distance along the duct - y, co-ordinate measuring distance across the duct - buoyancy parameter - ₀ modified buoyancy parameter, ₀=R ^–1/3 - coefficient of thermal expansion - ratio of duct width to heated length, =h/l - (non-dimensional) temperature - _w applied temperature on the wally=0 - kinematic viscosity - density of the fluid - ₀ shear stress on the wally=0 - ₁ shear stress on the wally=1 - stream function     

Transport in ordered and disordered porous media II: Generalized volume averaging

In this paper we develop the averaged form of the Stokes equations in terms of weighting functions. The analysis clearly indicates at what point one must choose a media-specific weighting function in order to achieve spatially smoothed transport equations. The form of the weighting function that produces the cellular average is derived, and some important geometrical theorems are presented.Roman Letters A interfacial area of the- interface associated with the local closure problem, m² - A _e area of entrances and exits for the-phase contained within the averaging system, m² - A _p surface area of a particle, m² - d _p 6V _p/A_p, effective particle diameter, m - g gravity vector, m/s² - I unit tensor - K _m permeability tensor for the weighted average form of Darcy's law, m² - L general characteristic length for volume averaged quantities, m - L _p general characteristic length for volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - l characteristic length (pore scale) for the-phase - l _i i=1, 2, 3 lattice vectors, m - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - m _v special convolution product weighting function associated with the traditional averaging volume - m _g general convolution product weighting function - m _V unit cell convolution product weighting function - m _C special convolution product weighting function for ordered media which produces the cellular average - m _D special convolution product weighting function for disordered media - m _M master convolution product weighting function for ordered and disordered media - n unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p_m superficial weighted average pressure, N/m² - p _m intrinsic weighted average pressure, N/m² - p traditional intrinsic volume averaged pressure, N/m² - p p p _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function, m - r position vector, m - r position vector locating points in the-phase, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - V velocity vector in the-phase, m/s - vm superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - V volume of the-phase contained in the averaging volume, m³ - V _p volume of a particle, m³ - v traditional superficial volume averaged velocity, m/s - v v p _m spatial deviation velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V /V, volume average porosity - _m m * . weighted average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m² - V /V, volume fraction of the-phase     

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

On surface waves in generalized thermoelasticity

Knowles' representation theorem for harmonically time-dependent free surface waves on a homogeneous, isotropic elastic half-space is extended to include harmonically time-dependent free processes for thermoelastic surface waves in generalized thermoelasticity of Lord and Shulman and of Green and Lindsay.r , , r , , .This work was done when author was unemployed.     

On the properties of compressible gas flow in a porous media

The flow of an adiabatic gas through a porous media is treated analytically for steady one- and two-dimensional flows. The effect on a compressible Darcy flow by inertia and Forchheimer terms is studied. Finally, wave solutions are found which exhibit a cut-off frequency and a phase shift between pressure and velocity of the gas, with the velocity lagging behind the pressure.Nomenclature A area of tube for one-dimensional flow - B drag coefficient associated with Forchheimer term - c speed of sound - M Mach number - p ^* gas pressure - p dimensionless gas pressure - s coordinate along the axis of tube - t ^* time variable - t dimensionless time variable - V^* gas velocity in the porous media - V dimensionless gas velocity Greek Letters ratio of specific heat capacities - phase angle between gas pressure and velocity for linear waves - parameter indicating the importance of the inertia term - viscosity - _p natural frequency of the porous media - ^* gas density - dimensionless gas density - parameter indicating the importance of the Forchheimer term - porosity of porous media - velocity potential - stream function     

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     

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

Transport in ordered and disordered porous media I: The cellular average and the use of weighting functions

In this work we consider transport in ordered and disordered porous media using singlephase flow in rigid porous mediaas an example. We defineorder anddisorder in terms of geometrical integrals that arise naturally in the method of volume averaging, and we show that dependent variables for ordered media must generally be defined in terms of thecellular average. The cellular average can be constructed by means of a weighting function, thus transport processes in both ordered and disordered media can be treated with a single theory based on weighted averages. Part I provides some basic ideas associated with ordered and disordered media, weighted averages, and the theory of distributions. In Part II a generalized averaging procedure is presented and in Part III the closure problem is developed and the theory is compared with experiment. Parts IV and V provide some geometrical results for computer generated porous media.Roman Letters A interfacial area of the- interface contained within the macroscopic region, m² - A_e area of entrances and exits for the-phase contained within the macroscopic system, m² - g gravity vector, m/s² - I unit tensor - K traditional Darcy's law permeability tensor, m² - L general characteristic length for volume averaged quantities, m - characteristic length (pore scale) for the-phase - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - N unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p0 reference pressure in the-phase, N/m² - p traditional intrinsic volume averaged pressure, N/m² - r0 radius of a spherical averaging volume, m - r position vector, m - r position vector locating points in the-phase, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - v velocity vector in the-phase, m/s - v traditional superficial volume averaged velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V/V, volume average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m²     

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     

Perturbation analysis for periodic heat transfer in radiating fins

A perturbation analysis is presented for periodic heat transfer in radiating fins of uniform thickness. The base temperature is assumed to oscillate around a mean value. The perturbation expansion is carried out in terms of dimensionless amplitude of the base temperature oscillation. The zero-order problem which is nonlinear, and corresponds to the steady state fin behaviour, is solved by quasilinearization. A method of complex combination is used to reduce both the first and the second order problems to two, coupled linear boundary value problems which are subsequently solved by a noniterative numerical scheme. The second-order term is composed of an oscillatory component with twice the frequency of base temperature oscillation and a time-independent term which causes a net change in the steady state values of temperature and heat transfer rate. Within the range of parameters used, the net effect is to decrease the mean temperature and increase the mean heat transfer rate. This is in constrast to the linear case of convecting fins where the mean values are unaffected by base temperature oscillations. Detailed numerical results are presented illustrating the effects of fin parameter N and dimensionless frequency B on temperature distribution, heat transfer rate, and time-average fin efficiency. The time-average fin efficiency is found to reduce significantly at low N and high B.

---

Störungsanalyse für periodische Wärmeübertragung an Strahlungsrippen

Zusammenfassung Eine Störungsanalyse wird für periodische Wärmeübertragung in Strahlungsrippen gleicher Dicke vorgelegt. Die Fußtemperatur wird als um einen Mittelwert schwingend angenommen. Die Störungsentwicklung wird in Termen einer dimensionslosen Amplitude e dieser Schwingung angesetzt. Das Problem nullter Ordnung, das nichtlinear ist und dem stationären Verhalten der Rippe entspricht, wird durch Quasilinearisierung gelöst. Eine Methode der komplexen Kombination wird angewandt, um die Probleme erster und zweiter Ordnung auf zwei gekoppelte Grenzwertprobleme zu reduzieren, die nacheinander nach einem nichtiterativen Schema gelöst werden. Der Term zweiter Ordnung besteht aus einer Schwingungskomponente mit der doppelten Frequenz der Schwingung der Fußtemperatur und einem zeitunabhängigen Term, der eine Nettoänderung der stationären Werte der Temperatur und der Wärmeübertragung verursacht. Im verwendeten Bereich der Parameter tritt eine Abnahme der mittleren Temperatur und eine Zunahme der mittleren Wärmeübertragung auf. Das steht im Gegensatz zum linearen Fall der Konvektionsrippe, bei dem die Mittelwerte durch Schwingungen der Fußtemperatur nicht beeinflußt werden. Detaillierte numerische Ergebnisse zeigen die Einflüsse des Rippenparameters N und der dimensionslosen Frequenz B auf Temperatur Verteilung, Wärmeübertragung und zeitliches Mittel des Rippengütegrades. Dieses zeitliche Mittel nimmt merklich ab bei kleinem N und hohem B.

---

Nomenclature b fin thickness - B dimensionless frequency, L²/ - E emissivity - f₀, f₁ functions of X - g₀, g₁, g₂ functions of X - h₀, h₁, h₂ functions of X - k thermal conductivity - L fin Length - N fin parameter, 2EL²T_bm/bk - q heat transfer rate - Q dimensionless heat transfer rate, qL/kbT_bm - t time - T temperature - T_b fin base temperature - T_S effective sink temperature - T_bm mean fin base temperature - x axial distance - X dimensionless axial distance, x/L - dimensionless amplitude of base temperature (s. Eq.2) - thermal diffusivity - instantaneous fin efficiency - time-average fin efficiency - _ss steady state fin efficiency - dimensionless temperature, T/T_bm - ₀ zero-order approximation - ₁ first-order approximation - ₂ second-order approximation - _2s steady component of ₂ - , ₁, ₂ constants - complex function of X - ₁ real part of - ₂ imaginary part of - complex function of X - ₁ real part of Y - ₂ imaginary part of - dimensionless time, t/L² - frequency of base temperature oscillation     

Evaluating Herschel-Bulkley fluids with the back extrusion (annular pumping) technique

A mathematical model was developed to describe the behavior of Herschel-Bulkley fluids in a back extrusion (annular pumping) device. A technique was also developed to determine the rheological properties (yield stress, flow behavior index, and consistency coefficient) of these fluids. Mathematical terms were expressed in four dimensionless terms, and graphical aids and tables were prepared to facilitate the handling of the expressions.Nomenclature a radius of the plunger, m - dv/dr shear rate, s^–1 - F force applied to the plunger, N - F _b buoyancy force, N - F _cb force corrected for buoyancy, N - F _T recorded force just before the plunger is stopped, N - F _Te recorded force after the plunger is stopped, N - g acceleration due to gravity, m/s² - H(t) momentary height between plunger and container bottom, m - K a/R, dimensionless - L length of annular region, m - L(t) depth of plunger penetration, m - n flow behavior index, dimensionless - p static pressure, Pa - P _L pressure in excess of hydrostatic pressure at the plunger base, Pa - p ₀ pressure at entrance to annulus, Pa - P pressure drop per unit of length, Pa/m - Q total volumetric flow rate through the annulus, m³/s - r radial coordinate, measured from common axis of cylinder forming annulus, m - R radius of outer cylinder of annulus, m - s reciprocal of n, dimensionless - t time, s - T dimensionless shear stress, defined in Eq. (3) - T ₀ dimensionless yield stress, defined in Eq. (4) - T _w dimensionless shear stress at the plunger wall - _p velocity of plunger, m/s - velocity, m/s - mass density of fluid, kg/m³ - Newtonian viscosity, Pa s - P p ₀ –p _L , Pa - consistency coefficient, Pa sⁿ - value of where shear stress is zero - _–, ₊ limits of the plug flow region (Fig. 1) - r/R - shear stress, Pa - _y yield stress, Pa - _w shear stress at the plunger wall, Pa - dimensionless flow rate defined in Eq. (24) - dimensionless velocity defined by Eq. (5) - _–, ₊ dimensionless velocity outside the plug flow region - _max dimensionless maximum velocity in the plug flow region - _p dimensionless velocity at the plunger wall     

Normal forms for random diffeomorphisms

Given a dynamical system (,, ,) and a random diffeomorphism (): ^{d d} with fixed point at x=0. The normal form problem is to construct a smooth near-identity nonlinear random coordinate transformation h() to make the random diffeomorphism ()=h()^–1() h() as simple as possible, preferably linear. The linearization D(, 0)=:A() generates a matrix cocycle for which the multiplicative ergodic theorem holds, providing us with stochastic analogues of eigenvalues (Lyapunov exponents) and eigenspaces. Now the development runs pretty much parallel to the deterministic one, the difference being that the appearance of turns all problems into infinite-dimensional ones. In particular, the range of the homological operator is in general not closed, making the conceptof-normal form necessary. The stochastic versions of resonance and averaging are developed. The case of simple Lyapunov spectrum is treated in detail.     

Interpolation formulas for maxwell viscosity of certain metals as a function of shear-strain intensity and temperature

The purpose of this study is the construction of interpolation formulas for the dependence of Maxwell viscosity, a quantity which is the reciprocal of shear-strain relaxation time , on shear-strain intensity and temperature for several metals: iron, aluminum, copper, and lead. This function was interpolated in various temperature and deformation velocity ranges in accordance with available experimental data for iron (0 10⁷ sec^–1, 200 ° T 1500 °); aluminum (0 10⁷ sec^–1, 300 ° T 900 °); copper (0 10⁵ sec^–1, 300 ° T 1300 °); lead (0 10⁶ sec^–1, 90 ° T 400 °); temperatures in °K.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 114–118, July–August, 1974.     

New concepts in relative permeabilities at high capillary numbers for surfactant flooding

Flooding oil reservoirs with surfactant solutions can increase the amount of oil that can be recovered. Macroscopic modelling of the process requires relative permeabilities to be functions of saturation and capillary number. With only limited experimental data, relative permeabilities have usually been assumed to be linear functions of saturation at high capillary numbers. The experimental data is reviewed, some of which suggest that this assumption is not necessarily correct. The basis for the assumption is therefore reviewed and it is concluded that the linear model corresponds to microscopically segregated flow in the porous medium. Based on new but equally plausible complementary assumptions about the flow pattern, a mixed flow model is derived. These models are then shown to be limiting cases of a droplet model which represents the mixing scale within the porous medium and gives a physical basis for interpolating between the models. The models are based on physical concepts of flow in a porous medium and so the approach described here represents a significant improvement in the understanding of high capillary number flow. This is shown by the fact that fewer parameters are needed to describe experimental data.Notation A total cross-sectional area assigned to capillary bundle - A ⁽ⁱ⁾ physical cross-sectional area of tube i - c ⁽ⁱ⁾ ordered configurational label for droplets in tube i - c configuration label for tube i (order not considered) - D defined by Equation (26) - E(...) expectation value with respect to the trinomial distribution - S _r ⁽⁾ fractional flow of phase - k absolute permeability - k _r relative permeability of phase - k _r ⁰ endpoint relative permeability of phase - L capillary tube length in bundle model - m ⁽ⁱ⁾ number of droplets of phase a occupying tube i - n exponent for phase a in Equation (2) - N number of droplets in bundle model - N _c capillary number - p pressure - p(c') probability of configuration c - Q ⁽ⁱ⁾ total volume flow rate in tube i - S saturation of phase - S flowing saturation of phase - S _r residual saturation of phase - S _r ⁽⁾ saturations when fractional flow of phase is 1 in the case of varying residual saturations for three-phase flow ( ) - t _c residence time for droplet configuration c - v ⁽ⁱ⁾ total fluid velocity in bundle tube i - , phase label - p pressure differential across capillary bundle - ⁽ⁱ⁾ tube conductivity defined by Equation (7) - viscosity of phase - interfacial tension - gradient operator - ... average over tube droplet configurations     

Flow in porous media III: Deformable media

Stokes flow in a deformable medium is considered in terms of an isotropic, linearly elastic solid matrix. The analysis is restricted to steady forms of the momentum equations and small deformation of the solid phase. Darcy's law can be used to determine the motion of the fluid phase; however, the determination of the Darcy's law permeability tensor represents part of the closure problem in which the position of the fluid-solid interface must be determined.Roman Letters A interfacial area of the- interface 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 macroscopic system, m² - A ^* interfacial area of the- interface contained within a unit cell, m² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - E Young's modulus for the-phase, N/m² - e i unit base vectors (i = 1, 2, 3) - g gravity vector, m²/s - H height of elastic, porous bed, m - k unit base vector (=e ₃) - characteristic length scale for the-phase, m - L characteristic length scale for volume-averaged quantities, m - n unit normal vector pointing from the-phase toward the-phase (n = -n ) - p pressure in the-phase, N/m² - P p – g·r, N/m² - r ₀ radius of the averaging volume, m - r position vector, m - t time, s - T total stress tensor in the-phase, N/m² - T ⁰ hydrostatic stress tensor for the-phase, N/m² - u displacement vector for the-phase, m - V averaging volume, m₃ - V volume of the-phase contained within the averaging volume, m³ - v velocity vector for the-phase, m/s Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - shear coefficient of viscosity for the-phase, Nt/m² - first Lamé coefficient for the-phase, N/m² - second Lamé coefficient for the-phase, N/m² - bulk coefficient of viscosity for the-phase, Nt/m² - T –T ⁰ , a deviatoric stress tensor for the-phase, N/m²     

Asymptotic behavior of solutions of nonlinear elliptic equations

We study and obtain formulas for the asymptotic behavior as ¦x¦ of C ² solutions of the semilinear equation u=f(x, u), x (*) where is the complement of some ball in ⁿ and f is continuous and nonlinear in u. If, for large x, f is nearly radially symmetric in x, we give conditions under which each positive solution of (*) is asymptotic, as ¦x¦, to some radially symmetric function. Our results can also be useful when f is only bounded above or below by a function which is radially symmetric in x or when the solution oscillates in sign. Examples when f has power-like growth or exponential growth in the variables x and u usefully illustrate our results.     

Nicht-isotherme Rieselfilmströmung einer hochviskosen Flüssigkeit mit stark temperaturabhÄngiger ViskositÄt

Zusammenfassung Es werden Geschwindigkeitsverteilungen und Filmdickenabnahmen von nichtisothermen NEWTONschen und nicht-NEWTONschen (Potenzansatz) Rieselfilmen mit temperaturanhÄngiger ViskositÄt berechnet, wobei die Temperaturverteilung im Film als linear vorausgesetzt wird. An dicken Rieselfilmen mit Re=10^–4... 10^–2 sind Geschwindigkeitsprofile, Filmdicken und OberflÄchentemperaturen gemessen und daraus die thermische EinlauflÄnge bestimmt worden. Ausgehend von der Penetrationstheorie für eine endlich dicke Platte kann man für diese EinlauflÄnge eine Approximationsformel erhalten, die für Strömungen mit Re < 1000 verwendet werden kann.

---

Non-isothermal filmflow of a highly viscous liquid, the viscosity strongly depending on temperature

Velocity distributions and film thicknesses of nonisothermal NEWTONIAN and non-NEWTONIAN (power-law) falling films are computed assuming that the temperature across the film varies linearly. Experimental studies on thick falling films of Re=10^–4...10^–2 had been carried out to measure velocities, film thickness and surface temperature and to calculate the thermal entrance length. One can get for this entrance length a approximation formula which is valid for flows with RePr <1000 by applying the results for the thermal penetration into a material finite plate.

---

Bezeichnungen B dimensionsloser Temperaturkoeffizient - ¯B [K] Temperaturkoeffizient (ln)/(1/T) - c_p [J/kgK] spezif. WÄrme bei konst. Druck - Fo FOURIER-Zahl - g [m/s²] Erdbeschleunigung - H dimensionslose Filmdicke - h [m] Filmdicke - m [Pas^2–n] ViskositÄtskoeffizient im Potenzansatz von OSTWALD-DE WAELE - Nu NUSSELT-Zahl - n Flüssigkeitsexponent im Potenzansatz von OSTWALD-DE WAELE - Pr PRANDTL-Zahl (Gl.3.5) - q [W/m²] WÄrmestromdichte - Re REYNOLDS-Zahl (Gl.3.4) - T [K] Temperatur - t [s] Zeit - U dimensionslose Geschwindigkeit (X-Komponente) - u [m/s] Geschwindigkeitskomponente in x-Richtung - X dimensionslose Koordinate (X=x/h₀) - x [m] LÄnge, Koordinate - Y dimensionslose Koordinate (Y=y/h₀) - y [m] Höhe, Koordinate - [W/m²K] WÄrmeübergangskoeffizient - Plattenneigungswinkel gegen Horizontale - [s^–1] Schergeschwindigkeit - dimensionslose Temperatur (Gl.3.3) - [m²/s] TemperaturleitfÄhigkeit (Gl.3.3) - [W/mK] WÄrmeleitfÄhigkeit - [Pas] ViskositÄt - [kg/m³] spezif. Dichte - [Pa] Schubspannung Indizes a scheinbar (apparent) - 0 bei x=0, auch: isotherm - P auf die Penetrationszeit bezogen - s an der OberflÄche - T bei linearer Temperaturdifferenz T - w an der Wand - 99 auf =0,99 bezogen - gemittelt, Mittelwert - thermisch ausgebildet, bei x - proportional - ¯t ungefÄhr - kleiner oder gleich ungefÄhr     

The virtual origin of riblets in an open channel flow

In this paper, a method using the mean velocity profiles for the buffer layer was developed for the estimation of the virtual origin over a riblets surface in an open channel flow. First, the standardized profiles of the mixing length were estimated from the velocity measurement in the inner layer, and the location of the edge of the viscous layer was obtained. Then, the virtual origins were estimated by the best match between the measured velocity profile and the equations of the velocity profile derived from the mixing length profiles. It was made clear that the virtual origin and the thickness of the viscous layer are the function of the roughness Reynolds number. The drag variation coincided well with other results.Nomenclature f _r skin friction coefficient - f _ro skin friction coefficient in smooth channel at the same flow quantity and the same energy slope - g gravity acceleration - H water depth from virtual origin to water surface - H ⁺ u*H/ - H false water depth from top of riblets to water surface - H ⁺ u*H/ - I _e streamwise energy slope - I _b bed slope - k riblet height - k ⁺ u*k/ - l mixing length - l _s standardized mixing length - Q flow quantity - Re Reynolds number volume flow/unit width/v - s riblet spacing - u mean velocity - u* friction velocity = - u* false friction velocity = - y distance from virtual origin - y distance from top of riblet - y ₀ distance from top of riblet to virtual origin - y _v distance from top of riblet to edge of viscous layer - y ⁺ u*y/ - y ⁺ u*y/ - y ₀ ⁺ u*y ₀/ - u ⁺ u*y/ - shifting coefficient for standardization - thickness of viscous layer=y ₀+y - ⁺ u*/ - ⁺ u*/ - eddy viscosity - ridge angle - v kinematic viscosity - density - shear stress     

Dielectric properties of heterogeneous mixtures with a polar constituent

Summary After defining the boundaries for the dielectric constant of a heterogeneous mixture, the behaviour of such a mixture is studied as a function of the frequency, when one of its components is polar. Deviations from a semicircle are to be expected for the function _m =f( _m ) even when the dielectric properties of the polar constituent can be described with a semicircular Cole-Cole-arc. The relaxation time of the mixture is shorter than that of the polar constituent.     

Approximate solutions for drag coefficient of bubbles moving in shear-thinning elastic fluids

The drag coefficient for bubbles with mobile or immobile interface rising in shear-thinning elastic fluids described by an Ellis or a Carreau model is discussed. Approximate solutions based on linearization of the equations of motion are presented for the highly elastic region of flow. These solutions are in reasonably good agreement with the theoretical predictions based on variational principles and with published experimental data. C _D Drag coefficient - E ^* Differential operator [E ^{* 2} = ²/² + (sin/ ²)/(1/sin /)] - El Ellis number - F _D Drag force - K Consistency index in the power-law model for non-Newtonian fluid - n Flow behaviour index in the Carreau and power-law models - P Dimensionless pressure [=(p – p ₀)/₀ (U /R)] - p Pressure - R Bubble radius - Re ₀ Reynolds number [= 2R U /₀] - Re Reynolds number defined for the power-law fluid [= (2R) ⁿ U ^2–n /K] - r Spherical coordinate - t Time - U Terminal velocity of a bubble - u Velocity - Wi Weissenberg number - Ellis model parameter - Rate of deformation - Apparent viscosity - ₀ Zero shear rate viscosity - Infinite shear rate viscosity - Spherical coordinate - Parameter in the Carreau model - ^* Dimensionless time [=/(U /R)] - Dimensionless length [=r/R] - Second invariant of rate of deformation tensors - ^* Dimensionless second invariant of rate of deformation tensors [=/(U /R)²] - Second invariant of stress tensors - ^* Dimensionless second invariant of second invariant of stress tensor [= / ₀ ² (U /R)²] - Fluid density - Shear stress - ^* Dimensionless shear stress [=/ ₀ (U /R)] - _1/2 Ellis model parameter - ₁ ^2/* Dimensionless Ellis model parameter [= _1/2/ ₀(U /R)] - Stream function - ^* Dimensionless stream function [=/U R ²]