Thermal entrance region heat transfer for MHD laminar flow in parallel-plate channels with unequal wall temperatures

The problem of thermal entry heat transfer for Hartmann flow in parallel-plate channels with uniform but unequal wall temperatures considering viscous dissipation, Joule heating and axial conduction effects is approached by the eigenfunction expansion method. The series expansion coefficients for the nonorthogonal eigenfunctions are obtained by using a method for nonorthogonal series described by Kantorovich and Krylov [21]. Numerical results are obtained for the case with entrance condition parameter _o=1 and open circuit condition K=1. The parametric values of Ha=0, 2, 6, 10 and Br=0, –1 are considered for Hartmann and Brinkman numbers, respectively.

---

Zusammenfassung Das Problem der Wärmeübertragung im thermischen Einlauf einer Hartmannströmung im ebenen Spalt mit einheitlichen, aber ungleichen Wandtemperaturen wurde unter Berücksichtigung viskoser Dissipation, Joulescher Heizung und axialer Wärmeleitung mit Hilfe einer Entwicklung nach Eigenfunktionen behandelt. Die Koeffizienten der Entwicklung für nichtorthogonale Eigenfunctionen wurde nach einer Methode für nichtorthogonale Reihen nach Kantorovicz und Krylow [21] berechnet. Numerische Ergebnisse werden für den Eintrittsparameter _o=1 und die Bedingung für den offenen Stromkreis K=1 erhalten. Die Parameterwerte Ha=0, 2, 6, 10 und Br=0, –1 werden für die jeweiligen Werte der Hartmann- und der Brinckman-Zahl betrachtet.

---

Nomenclature a one-half of channel height - ¯B,B₀ magnetic field Induction vector and magnitude of applied magnetic field - Br Brinkman number, _f U_m ²/(k_c) - C_n,D_n coefficients in the series expansion of _e, see eq. 16 - c_p specific heat at constant pressure - ,E₀ electric field intensity vector and component - E_n,O_n even and odd eigenfunctions - Ha Hartmann number, (/_f)^1/2 B_o a - h₁,h₂ local heat transfer coefficients at lower and upper plates - ¯J,J_y electric current density vector and component - K external loading parameter, E_o/(B_o U_m) - k thermal conductivity - Nu₁, Nu₂ local Nusselt numbers, h₁,a/k and h₂a/k, respectively - P fluid pressure - Pe Peclet number, PrRe - Pr Prandtl number, C_{p f}/k - q₁,q₂ rates of heat transfer per unit area,^–k(T/Z)Z=–a'^–k(T/Z) Z=a respectively - Re Reynolds number, U_ma/u_f - T,T₀,T₁,T₂ fluid temperature, uniform entrance temperature, uniform but different lower and upper plate temperatures, respectively - T_b,T_m bulk temperature and (T₁+T₂)/2 - U,U_m,u axial, mean and dimensionless velocities, respectively - ¯V velocity vector - X,Z axial and transverse coordinates - x,z dimensionless coordinates - _n,_n even and odd eigenvalues - ,₀,_b dimensionless fluid, entrance and bulk temperatures, respectively - _c,_e,_f characteristic temperature difference (T₂-T_m), and dimensionless fluid temperatures, defined by eq. (10) - _e,_f magnetic permeability and viscosity of fluid - fluid density - electric conductivity - viscous dissipation function - (1-)/2     

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     

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

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

Asymptotic analysis of equilibrium states for rotating turbulent flows

The equilibrium states of homogeneous turbulence simultaneously subjected to a mean velocity gradient and a rotation are examined by using asymptotic analysis. The present work is concerned with the asymptotic behavior of quantities such as the turbulent kinetic energy and its dissipation rate associated with the fixed point (/kS)=0, whereS is the shear rate. The classical form of the model transport equation for (Hanjalic and Launder, 1972) is used. The present analysis shows that, asymptotically, the turbulent kinetic energy (a) undergoes a power-law decay with time for (P/)<1, (b) is independent of time for (P/)=1, (c) undergoes a power-law growth with time for 1<(P/)<(C ₂–1), and (d) is represented by an exponential law versus time for (P/)=(C ₂–1)/(C ₁–1) and (/kS)>0 whereP is the production rate. For the commonly used second-order models the equilibrium solutions forP/,II, andIII (whereII andIII are respectively the second and third invariants of the anisotropy tensor) depend on the rotation number when (P/kS)=(/kS)=0. The variation of (P/kS) andII versusR given by the second-order model of Yakhot and Orzag are compared with results of Rapid Distortion Theory corrected for decay (Townsend, 1970).     

Existence theorem for a minimum problem with free discontinuity set

We study the variational problem Where is an open set in ⁿ ,n2gL ^q () L (), 1q<+, O<, <+ andH _n–1 is the (n–1)-dimensional Hausdorff Measure.     

Flow in porous media II: The governing equations for immiscible,two-phase flow

The Stokes flow of two immiscible fluids through a rigid porous medium is analyzed using the method of volume averaging. The volume-averaged momentum equations, in terms of averaged quantities and spatial deviations, are identical in form to that obtained for single phase flow; however, the solution of the closure problem gives rise to additional terms not found in the traditional treatment of two-phase flow. Qualitative arguments suggest that the nontraditional terms may be important when / is of order one, and order of magnitude analysis indicates that they may be significant in terms of the motion of a fluid at very low volume fractions. The theory contains features that could give rise to hysteresis effects, but in the present form it is restricted to static contact line phenomena.Roman Letters (, = , , and ) A interfacial area of the- interface contained within the macroscopic system, m² - A _e area of entrances and exits for the -phase contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A ^* 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² - g gravity vector, m²/s - H mean curvature of the- interface, m^–1 - H area average of the mean curvature, m^–1 - H – H , deviation of the mean curvature, m^–1 - I unit tensor - K Darcy's law permeability tensor, m² - K permeability tensor for the-phase, m² - K viscous drag tensor for the-phase equation of motion - K viscous drag tensor for the-phase equation of motion - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the-phase, m - n unit normal vector pointing from the-phase toward the-phase (n = –n ) - p _c p – P , capillary pressure, N/m² - p pressure in the-phase, N/m² - p intrinsic phase average pressure for the-phase, N/m² - p – p , spatial deviation of the pressure in the-phase, N/m² - r ₀ radius of the averaging volume, m - t time, s - v velocity vector for the-phase, m/s - v phase average velocity vector for the-phase, m/s - v intrinsic phase average velocity vector for the-phase, m/s - v – v , spatial deviation of the velocity vector for the-phase, m/s - V averaging volume, m³ - V volume of the-phase contained within the averaging volume, m³ Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - viscosity of the-phase, Nt/m² - surface tension of the- interface, N/m - viscous stress tensor for the-phase, N/m² - / kinematic viscosity, m²/s     

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

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

Hydrodynamical changes due to polymer migration in very dilute solutions

The slip hypothesis, based on thermodynamical arguments, has been extended to obtain the flow characteristics of polymer solutions flowing in a nonhomogeneous flow field. An asymptotic analysis, valid for both channel and falling film flows, is presented that predicts the flow enhancement due to polymer migration. Concentration-viscosity coupling is shown to be a critical factor in the hydrodynamic analysis. The analysis, which essentially provides an upper bound on flow enhancement, explicitly accounts for the influence of wall shear stress, initial polymer concentration etc. A comparison with the pertinent experimental data shows reasonable agreement. c concentration - c ₀ concentration in shear-free region - c _i initial concentration - d rate of deformation tensor - g acceleration due to gravity - g ₁ function defined in eq. [13] or [15] - g ₂ function defined in eq. [18] or [20] - H half-channel thickness or film thickness - K gas law constant - L length of the channel or film - q flow rate per unit width - q ^* normalized flow rate - T temperature - v velocity - V mean velocity - y transverse distance - y _c location of solvent layer - _w /µ _s - _w /c ₀ KT - /t convected derivative - dimensionless cenentration,c/c ₀ - _c dimensionless interface concentration - _w dimensionless wall concentration - relaxation time - µ _eff effective viscosity - µ _s solvent viscosity - dimensionless transverse distance,y/H - _c dimensionless interface location - density - stress tensor - _w wall shear stress - c _i KT/ _w - ns no slip NCL-Communication No. 3155     

Transport in ordered and disordered porous media III: Closure and comparison between theory and experiment

In this paper we examine the closure problem associated with the volume averaged form of the Stokes equations presented in Part II. For both ordered and disordered porous media, we make use of a spatially periodic model of a porous medium. Under these circumstances the closure problem, in terms of theclosure variables, is independent of the weighting functions used in the spatial smoothing process. Comparison between theory and experiment suggests that the geometrical characteristics of the unit cell dominate the calculated value of the Darcy's law permeability tensor, whereas the periodic conditions required for thelocal form of the closure problem play only a minor role.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² - A interfacial area of the- interface associated with the local closure problem, m² - A _p surface area of a particle, m² - b vector used to represent the pressure deviation, m^–1 - B ⁰ B+I, a second order tensor that maps v _m ontov - B second-order tensor used to represent the velocity deviation - d _p 6V _p/A_p, effective particle diameter, m - d a vector related to the pressure, m - D a second-order tensor related to the velocity, m² - g gravity vector, m/s² - I unit tensor - K traditional Darcy's law permeability tensor calculated on the basis of a spatially periodic model, m² - 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 characteristic length for the volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - characteristic length (pore scale) for the-phase - _i i=1, 2, 3 lattice vectors, m - weighting function - m(-y) , convolution product weighting function - 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 - 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 _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function - r position vector, m - r position vector locating points in the-phase, m. - V averaging volume, m³ - B 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 _m superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - v traditional superficial volume averaged velocity, m/s - v – v _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²     

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

Free convection on an arbitrarily inclined plate with uniform surface heat flux

This paper proposed a proper inclination parameter and transformation variables for the analysis of free convection from an inclined plate with uniform surface heat flux to fluids of any Prandtl number. Very accurate numerical results and a simple correlation equation are obtained for arbitrary inclination from the horizontal to the vertical and for 0.001 Pr. Maximum deviation between the correlated and calculated data is less than 1.2%.

---

Freie Konvektion an einer beliebig geneigten Platte mit erheblicher Wärmestromdichte an der Oberfläche

Zusammenfassung Für die Berechnung von freier Konvektion von Fluiden mit beliebiger Prandtl-Zahl an einer geneigten Platte mit einheitlicher Wärmestromdichte an der Oberfläche werden ein zweckmäßiger Neigungsparameter und Transformationsvariablen eingeführt. Sehr genaue numerische Ergebnisse und eine einfache Korrelationsgleichung wurden für beliebige Neigungen zwischen der Horizontalen und der Vertikalen und für 0.001Pr erhalten. Die größte Abweichung zwischen Korrelations- und berechneten Daten liegt bei weniger als 1.2%.

---

Nomenclature f reduced stream function - g gravitational acceleration - h local heat transfer coefficient - k thermal conductivity - Nu local Nusselt number - p static pressure difference - Pr Prandtl number - q _w wall heat flux - Ra* modified local Rayleigh number,g(q _w x/k)x ³/ - T fluid temperature - T temperature of ambient fluid - u velocity component inx-direction - v velocity component iny-direction - x coordinate parallel to the plate - y coordinate normal to the plate Greek symbols thermal diffusivity - thermal expansion coefficient - (Ra* |sin|)^1/5/( Ra* cos)^1/6 - ( Ra* cos)^1/6+( Ra*|sin|)^1/5 - (y/x) - dimensionless temperature, (T–T )/(q _w x/k) - kinematic viscosity - [1+( Ra* cos)^1/6/( Ra*|sin|)^1/5]^–1 - density of fluid - Pr/(1+Pr) - _w wall shear stress - angle of inclination measured from the horizontal - stream function - dimensionless static pressure difference, p x ²/ ⁴     

Stability characteristics of condensate films

The linear stability theory is used to study stability characteristics of laminar condensate film flow down an arbitrarily inclined wall. A critical Reynolds number exists above which disturbances will be amplified. The magnitude of the critical Reynolds number is in all practical situations so small that a laminar gravity-induced condensate film can be expected to be unstable. Several stabilizing effects are acting on the film flow; at an inclined wall these effects are due to surface tension, gravity and condensation mass transfer.

---

Zusammenfassung Mit Hilfe der linearen Stabilitätstheorie werden die Stabilitätseigenschaften laminarer Kondensatfilme an einer geneigten Wand untersucht. Es zeigt sich, daß Kondensatfilme in jedem praktischen Fall ein unstabiles Verhalten aufweisen. Der stabilisierende Einfluß von Oberflächenspannung, Schwerkraft und Stoffübertragung durch Kondensation bewkkt jedoch, daß Störungen in bestimmten Wellenlängenbereichen gedämpft werden.

---

Nomenclature c=c^*/u₀ complex wave velocity, celerity, dimensionless - c^*=c _r ^* + i c _i ^* complex wave velocity, celerity, dimensional - c_p specific heat at constant pressure - g gravitational acceleration - h_fg latent heat - k thermal conductivity of liquid - p^* pressure - p=p^*/u₀ ² dimensionless pressure - Pe=Pr Re^* Peclet number - Pr Prandtl number - Re^*=u₀ / Reynolds number (defined with surface velocity) - S temperature perturbation amplitude - t^* time - t=t^* u₀/ dimensionless time - T temperature - T_s saturation temperature - T_w wall temperature - T=T_s-T_w temperature drop across liquid film - u^*, v^* velocity components - u=u^*/u₀ dimensionless velocity components - v=v^*/u₀ dimensionless velocity components - u₀ surface velocity of undisturbed film flow - v _g ^* vapor velocity - x^*, y^* coordinates - x=x^*/ dimensionless coordinates - y=y^*/ dimensionless coordinates Greek Symbols =^* wave number, dimensionless - ^*=2 /^* wave number dimensional - ^* wave length, dimensional - =^*/ wave length, dimensionless - local thickness of undisturbed condensate film - kinematic viscosity, liquid - density, liquid - _g density vapor - surface tension - = (1 +) film thickness of disturbed film, Fig. 1 - stream function perturbation amplitude - angle of inclination Base flow quantities are denoted by^–, disturbance quantities are denoted by.     

Oscillating sedimentation of spheres in viscoplastic fluids

In applications of chemical engineering often the sedimentation is used to separate disperse particles from liquid phases. Some real liquids, e.g., polymer fluids, paints, and skin creams show viscoplastic flow behavior, i.e., they have a yield stress. In such fluids it is possible that suspended particles do not move under action of gravity although the density of the particles is greater than the fluid density. A possibility to sediment stuck spherical particles is shown. The fluid is set in sinusoidal vibration so that the particles undergo forced oscillations. This effect is investigated for single spheres. A model is given and several theoretical results are discussed. A criterion is presented that allows one to predict the combinations of the vibration parameters (amplitude and frequency) which are needed to sediment the spheres. The theoretical investigations are confirmed by experiments. The motion of several glass and steel spheres in an oscillating tube filled with aqueous carbopol solutions are detected. The comparison between theory and experiment shows good agreement.Nomenclature C Stokes drag coefficient - D strain rate tensor - E unit tensor - F, F _w external resp. drag force - f body force vector - G Green deformation tensor - G dimensionless (shear) modulus - g acceleration of gravity - K abbreviation (Eq. (15)) - p pressure - R, spherical coordinates - R ₀, R _a sphere resp. body radius - T extra stress tensor - t time - S stress tensor - U ₀, W ₀ displacement amplitude - u displacement vector - velocity vector - V viscoelastic number - V _k sphere volume - V steady sink velocity - Y yield stress parameter - Y _g limiting value - Z Stokes number - density ratio - second invariant of D - radii ratio - y differential viscosity - , Lamé constants - ^* = + i complex (shear) modulus - f, _k fluid resp. sphere density - second invariant of T - _f yield stress - phase angle - frequency     

Nonlocal dispersion in media with continuously evolving scales of heterogeneity

General nonlocal diffusive and dispersive transport theories are derived from molecular hydrodynamics and associated theories of statistical mechanical correlation functions, using the memory function formalism and the projection operator method. Expansion approximations of a spatially and temporally nonlocal convective-dispersive equation are introduced to derive linearized inverse solutions for transport coefficients. The development is focused on deriving relations between the frequency-and wave-vector-dependent dispersion tensor and measurable quantities. The resulting theory is applicable to porous media of fractal character.Nomenclature C _v (t) particle velocity correlation function - C _v ,(t) particle fluctuation velocity correlation function - C _j (x,t) current correlation function - D(x,t) dispersion tensor - D(x,t) fluctuation dispersion tensor - f ₀(x,p) equilibrium phase probability distribution function - f(x, p;t) nonequilibrium phase probability distribution function - G(x,t) conditional probability per unit volume of finding a particle at (x,t) given it was located elsewhere initially - (k,t) Fourier transform ofG(x,t) - G(x,t) fluctuation conditional probability per unit volume of finding a particle at (x,t) given it was located elsewhere initially - k wave vector - K(t) memory function - L Liouville operator - m mass - p(t) particle momentum coordinate - P = (0)( , (0)) projection operator - Q =I-P projection operator - s real Laplace space variable - S(k, ) time-Fourier transform of(k,t) - t time - v(t) particle velocity vector - v(t) particle fluctuation velocity vector - V phase space velocity - time-Fourier variable - ^(itn)(k) frequency moment of(k,t) - x(t) particle displacement coordinate - x(t) particle displacement fluctuation coordinate - friction coefficient - (t) normalized correlation function General Functions () Dirac delta function - () Gamma function Averages ₀ Equilibrium phase-space average - Nonequilibrium phase-space average - (,) L ² inner product with respect tof ₀     

Dynamic behaviour of the constant temperature anemometer due to thermal inertia of the wire

The time dependent differential equation for the local wire temperature of a constant temperature anemometer is solved by a perturbation method in case of a harmonically changing heat transfer coefficient. The time dependent power supply to the wire follows from the condition of constant mean temperature imposed by the anemometer circuit. The influence of thin supporting wires, or copper-plated wire ends, is evaluated also. Numerical results are given for a number of cases that are of practical interest.Nomenclature c specific heat - D diameter of the wire - D _u diameter of the copper-plated ends of the wire - f D - g I ² r ₀ - I electric current - L length of the wire - P 1/4D ² c - q 1/4D ² - r resistance of the wire per unit length at temperature T' - r ₀ resistance of the wire per unit length at temperature T - T T' – T - T' local temperature of the wire - T ambient temperature - T _w constant mean temperature imposed by the anemometer circuit - T _u local difference between the temperature of the supporting wire and the ambient temperature - t time - x axial coordinate with the origin in the middle of the wire - heat transfer coefficient - temperature coefficient of the resistance - small parameter - time constant = cD ²/4D - _u time constant of the copper-plated ends cD _u ² /4D _u - thermal conductivity of wire material - _u thermal conductivity of the copper-plated wire ends - density - circular frequency     

Mean and turbulent velocity measurements of supersonic mixing layers

The behavior of supersonic mixing layers under three conditions has been examined by schlieren photography and laser Doppler velocimetry. In the schlieren photographs, some large-scale, repetitive patterns were observed within the mixing layer; however, these structures do not appear to dominate the mixing layer character under the present flow conditions. It was found that higher levels of secondary freestream turbulence did not increase the peak turbulence intensity observed within the mixing layer, but slightly increased the growth rate. Higher levels of freestream turbulence also reduced the axial distance required for development of the mean velocity. At higher convective Mach numbers, the mixing layer growth rate was found to be smaller than that of an incompressible mixing layer at the same velocity and freestream density ratio. The increase in convective Mach number also caused a decrease in the turbulence intensity ( _u/U).List of symbols a speed of sound - b total mixing layer thickness between U ₁ – 0.1 U and U ₂ + 0.1 U - f normalized third moment of u-velocity, f u³/(U)³ - g normalized triple product of u² , g u²/(U)³ - h normalized triple product of u ², h u ²/(U)³ - l _u axial distance for similarity in the mean velocity - l _u axial distance for similarity in the turbulence intensity - M Mach number - M _c convective Mach number (for ₁ = ₂), M _c (U ₁ – U ₂)/(a ₁ + a ₂) - P static pressure - r freestream velocity ratio, r U ₂/U ₁ - Re unit Reynolds number, Re U/ - s freestream density ratio, s ₂/₁ - T _t total temperature - u instantaneous streamwise velocity - u deviation of u-velocity, uu – U - U local mean streamwise velocity - U ₁ primary freestream velocity - U ₂ secondary freestream velocity - average of freestream velocities, (U ₁ + U ₂)/2 - U freestream velocity difference, U U ₁ – U ₂ - instantaneous transverse velocity - v deviation of -velocity, – V - V local mean transverse velocity - x streamwise coordinate - y transverse coordinate - y ₀ transverse location of the mixing layer centerline - ensemble average - ratio of specific heats - boundary layer thickness (y-location at 99.5% of free-stream velocity) - similarity coordinate, (y – y ₀)/b - compressible boundary layer momentum thickness - viscosity - density - standard deviation - dimensionless velocity, (U – U ₂)/U - 1 primary stream - 2 secondary stream A version of this paper was presented at the 11th Symposium on Turbulence, October 17–19, 1988, University of Missouri-Rolla     

Production of temperature fluctuations in grid turbulence: Wiskind's experiment revisited

Measurements have been made in nearly-isotropic grid turbulence on which is superimposed a linearly-varying transverse temperature distribution. The mean-square temperature fluctuations, , increase indefinitely with streamwise distance, in accordance with theoretical predictions, and consistent with an excess of production over dissipation some 50% greater than values recorded in previous experiments. This high level of production has the effect of reducing the ratio,r, of the time scales of the fluctuating velocity and temperature fields. The results have been used to estimate the coefficient,C, in Monin's return-to-isotropy model for the slow part of the pressure terms in the temperature-flux equations. An empirical expression by Shih and Lumley is consistent with the results of earlier experiments in whichr 1.5, C 3.0, but not with the present data where r 0.5, C 1.6. Monin's model is improved when it incorporates both time scales.List of symbols C coefficient in Monin model, Eq. (5) - M grid mesh length - m exponent in power law for temperature variance, x ^m - n turbulence-energy decay exponent,q ² x ^-n - p production rate of - p pressure - q ² - R microscale Reynolds number - r time-scale ratiot/t - T mean temperature - U mean velocity - mean-square velocity fluctuations (turbulent energy components) - turbulent temperature flux - x, y, z spatial coordinates - temperature gradient dT/dy - thermal diffusivity - dissipation rate ofq ²/2 - dissipation rate of - Taylor microscale (₂=5q₂/) - temperature microscale - _v temperature-flux correlation coefficient, /v - dimensionless distance from the grid,x/M     

Flow in porous media I: A theoretical derivation of Darcy's law

Stokes flow through a rigid porous medium is analyzed in terms of the method of volume averaging. The traditional averaging procedure leads to an equation of motion and a continuity equation expressed in terms of the volume-averaged pressure and velocity. The equation of motion contains integrals involving spatial deviations of the pressure and velocity, the Brinkman correction, and other lower-order terms. The analysis clearly indicates why the Brinkman correction should not be used to accommodate ano slip condition at an interface between a porous medium and a bounding solid surface.The presence of spatial deviations of the pressure and velocity in the volume-averaged equations of motion gives rise to aclosure problem, and representations for the spatial deviations are derived that lead to Darcy's law. The theoretical development is not restricted to either homogeneous or spatially periodic porous media; however, the problem ofabrupt changes in the structure of a porous medium is not considered.Roman Letters A interfacial area of the - interface contained within the macroscopic system, m² - A _e area of entrances and exits for the -phase contained within the macroscopic system, m² - A interfacial area of the - interface contained within the averaging volume, m² - A ^* interfacial area of the - interface contained within a unit cell, m² - Ae area of entrances and exits for the -phase contained within a unit cell, m² - B second order tensor used to represent the velocity deviation (see Equation (3.30)) - b vector used to represent the pressure deviation (see Equation (3.31)), m^–1 - d distance between two points at which the pressure is measured, m - g gravity vector, m/s² - K Darcy's law permeability tensor, m² - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the -phase (see Figure 2), m - characteristic length scale for the -phase (see Figure 2), m - n unit normal vector pointing from the -phase toward the -phase (n =–n ) - n _e unit normal vector for the entrances and exits of the -phase contained within a unit cell - p pressure in the -phase, N/m² - p intrinsic phase average pressure for the -phase, N/m² - p –p , spatial deviation of the pressure in the -phase, N/m² - r ₀ radius of the averaging volume and radius of a capillary tube, m - v velocity vector for the -phase, m/s - v phase average velocity vector for the -phase, m/s - v intrinsic phase average velocity vector for the -phase, m/s - v –v , spatial deviation of the velocity vector for the -phase, m/s - V averaging volume, m³ - V volume of the -phase contained within the averaging volume, m³ Greek Letters V/V, volume fraction of the -phase - mass density of the -phase, kg/m³ - viscosity of the -phase, Nt/m² - arbitrary function used in the representation of the velocity deviation (see Equations (3.11) and (B1)), m/s - arbitrary function used in the representation of the pressure deviation (see Equations (3.12) and (B2)), s^–1     

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     

A three-parameter model describing the experimental relation between shear stress and shear rate for laminar flow of aqueous polymer solutions

Summary A three-parameter model is introduced to describe the shear rate — shear stress relation for dilute aqueous solutions of polyacrylamide (Separan AP-30) or polyethylenoxide (Polyox WSR-301) in the concentration range 50 wppm – 10,000 wppm. Solutions of both polymers show for a similar rheological behaviour. This behaviour can be described by an equation having three parameters i.e. zero-shear viscosity ₀, infinite-shear viscosity , and yield stress ₀, each depending on the polymer concentration. A good agreement is found between the values calculated with this three-parameter model and the experimental results obtained with a cone-and-plate rheogoniometer and those determined with a capillary-tube rheometer.

---

Zusammenfassung Der Zusammenhang zwischen Schubspannung und Schergeschwindigkeit von strukturviskosen Flüssigkeiten wird durch ein Modell mit drei Parametern beschrieben. Mit verdünnten wäßrigen Polyacrylamid-(Separan AP-30) sowie Polyäthylenoxidlösungen (Polyox WSR-301) wird das Modell experimentell geprüft. Beide Polymerlösungen zeigen im untersuchten Schergeschwindigkeitsbereich von ein ähnliches rheologisches Verhalten. Dieses Verhalten kann mit drei konzentrationsabhängigen Größen, nämlich einer Null-Viskosität ₀, einer Grenz-Viskosität und einer Fließgrenze ₀ beschrieben werden. Die Ergebnisse von Experimenten mit einem Kegel-Platte-Rheogoniometer sowie einem Kapillarviskosimeter sind in guter Übereinstimmung mit den Werten, die mit dem Drei-Parameter-Modell berechnet worden sind.

---

a Pa^–1 physical quantity defined by:a = {1 – ( / ₀)}/ ₀ - c l concentration (wppm) - D m capillary diameter - L m length of capillary tube - P Pa pressure drop - R m radius of capillary tube - u m s^–1 average velocity - v _r m s^–1 local axial velocity at a distancer from the axis of the tube - shear rate (–dv _r /dr) - local shear rate in capillary flow - s^–1 wall shear rate in capillary flow - Pa s dynamic viscosity - _a Pa s apparent viscosity defined by eq. [2] - ( _a ) Pa s apparent viscosity in capillary tube at a distanceR from the axis - ₀ Pa s zero-shear viscosity defined by eq. [4] - Pa s infinite-shear viscosity defined by eq. [5] - l ratior/R - kg m^– density - Pa shear stress - ₀ Pa yield stress - _r Pa local shear stress in capillary flow - _R Pa wall shear stress in capillary flow _R = (PR/2L) - _v m³ s^–1 volume rate of flow With 8 figures and 1 table