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

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

Nomenclature

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

Pulsatile flow in a tube with wall injection and suction

This paper studies similarity solutions for pulsatile flow in a tube with wall injection and suction. The Navier-Stokes equations are reduced to a system of three ordinary differential equations. Two of the equations represent the effects of suction and injection on the steady flow while the third represents the effects of suction and injection on pulsatile flow. Since the equations for steady flow have been studied previously, the analysis centers on the third equation. This equation is solved numerically and by the method of matched asymptotic expansions. The exact numerical solutions compare well with the asymptotic solutions.The effects of suction and injection on pulsatile flow are the following: a) Small values of suction can cause a resonance-like effect for low frequency pulsatile flow. b) The annular effect still occurs but for large injection or suction the frequency at which this effect becomes dominant depends on the cross-flow Reynolds number. c) The maximum shear stress at the wall is decreased by injection, but may be increased or decreased by suction.Nomenclature a radius of the tube - a ₀ ² i ² - A₀, B₀, C₀, D₀, E₀ constant coefficients appearing in the expression for pressure - b a non-dimensionalized length - b ₀ ² i ^{2 2} - b _k complex coefficients of a power series - B - C ₁, C ₂, D complex constants - d - D _1,2 - f() F(a ^1/2)/aV - f ₀,f ₁,... functions of order one used in asymptotic expansions of f() - F(r) rv _r - g() - G(r) a steady component of velocity in axial direction - h() 4/C₀ a ² H(a ^1/2) - h ₀,h ₁,h ₂,...;l ₀,l ₁,l ₂,... functions of order one used in asymptotic expansions for h() in outer regions - H(r) complex valued function giving unsteady component of velocity - H ₀, H ₁, H ₂, ... K ₀, K ₁, K ₂, ...; L ₀, L ₁, L ₂, ... functions of order one used in asymptotic expansions for h() in inner regions - i - J ₀, J ₁, Y ₀, Y ₁ Bessel functions of first and second kind - k - K Rk/2b ² - O order symbol - p pressure - p ₁(z, t) arbitrary function related to pressure - r radial coordinate - r ₀ (1+16 ^{4 4})^1/4 - R Va/, the crossflow Reynolds number - t time - u() G(r)/V - v _r radial velocity - v _z axial velocity - V constant velocity at which fluid is injected or extracted - z axial coordinate - ² a ²/4 - 4.196 - small parameter; =–2/R (Sect. 4); =–R/2 (Sect. 5); =2/R(Sect. 6) - r ²/a ² - * 0.262 - Arctan (4 ^{2 2}) - , inner variables - kinematic viscosity - b - * zero of g() - density - (r, t) arbitrary function related to axial velocity - frequency     

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 natural convection from a vertical plate with a prescribed surface heat flux in porous media

This paper presents a theoretical and numerical investigation of the natural convection boundary-layer along a vertical surface, which is embedded in a porous medium, when the surface heat flux varies as (1 +x ²)), where is a constant andx is the distance along the surface. It is shown that for > -1/2 the solution develops from a similarity solution which is valid for small values ofx to one which is valid for large values ofx. However, when -1/2 no similarity solutions exist for large values ofx and it is found that there are two cases to consider, namely < -1/2 and = -1/2. The wall temperature and the velocity at large distances along the plate are determined for a range of values of .Notation g Gravitational acceleration - k Thermal conductivity of the saturated porous medium - K Permeability of the porous medium - l Typical streamwise length - q _w Uniform heat flux on the wall - Ra Rayleigh number, =gK(q _w /k)l/(v) - T Temperature - Too Temperature far from the plate - u, v Components of seepage velocity in the x and y directions - x, y Cartesian coordinates - Thermal diffusivity of the fluid saturated porous medium - The coefficient of thermal expansion - An undetermined constant - Porosity of the porous medium - Similarity variable, =y(1+x ) ^/3/x ^1/3 - A preassigned constant - Kinematic viscosity - Nondimensional temperature, =(T – T )Ra^1/3 k/qw - Similarity variable, = =y(log_e x)^1/3/x ^2/3 - Similarity variable, =y/x ^2/3 - Stream function     

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²     

Flow birefringence of polymer melts: Application to the investigation of time dependent rheological properties

Summary Transient stresses including normal stresses, which are developed in a polymer melt by a suddenly imposed constant rate of shear, are investigated by mechanical measurement and, indirectly, with the aid of the flow birefringence technique. For the latter purpose use is made of the so-called stress-optical law, which is carefully checked.It appears that the essentially linear model of the rubberlike liquid, as proposed byLodge, is capable of describing the behaviour of polymer melts rather well, if the applied total shear does not exceed unity. In order to describe also steady state values of the stresses successfully, one should extend measurements to extremely low shear rates.These statements are verified with the aid of a method which was originally designed bySchwarzl andStruik for the practical calculation of interrelations between linear viscoelastic functions. In the present paper dynamic shear moduli are used as reference functions.

---

Zusammenfassung Mit der Zeit anwachsende Spannungen, darunter auch Normalspannungen, wie sie sich nach dem plötzlichen Anlegen einer konstanten Schergeschwindigkeit in einer Polymerschmelze entwickeln, werden mit Hilfe mechanischer Messungen und indirekt mit Hilfe der Strömungsdoppelbrechung untersucht. Für den letzteren Zweck wird das sogenannte spannungsoptische Gesetz herangezogen, dessen Gültigkeit sorgfältig überprüft wird.Es ergibt sich, daß das im Wesen lineare Modell der gummiartigen Flüssigkeit, wie es vonLodge vorgeschlagen wurde, sich recht gut zur Beschreibung des Verhaltens von Polymerschmelzen eignet, solange der im ganzen angelegte Schub den Wert Eins nicht überschreitet. Um auch stationäre Werte der Spannungen in die Beschreibung erfolgreich einzubeziehen, sollte man die Messungen bis zu extrem niedrigen Schergeschwindigkeiten ausdehnen.Die gemachten Feststellungen werden mit Hilfe einer Methode verifiziert, die vonSchwarzl undStruik ursprünglich für die praktische Berechnung von Beziehungen zwischen Zustandsfunktionen entwickelt wurde, die dem linear viskoelastischen Verhalten entsprechen. In der vorliegenden Veröffentlichung dienen die dynamischen Schubmoduln als Bezugsfunktionen.

---

a _T shift factor - B _ij Finger deformation tensor - C stress-optical coefficient, (m²/N) - f (p _jl ) undetermined scalar function - G shear modulus, (N/m²) - G(t) time dependent shear modulus, (N/m²) - G() shear storage modulus, (N/m²) - G() shear loss modulus, (N/m²) - G _r reduced shear storage modulus, (N/m²) - G _r reduced shear loss modulus, (N/m²) - H() shear relaxation time spectrum, (N/m²) - k Boltzmann constant, (Nm/°K) - n _ik refractive index tensor - p undetermined hydrostatic pressure, (N/m²) - p _ij ,p _ik stress tensor, (N/m²) - p ₂₁ shear stress, (N/m²) - p ₁₁ –p ₂₂ first normal stress difference, (N/m²) - p ₂₂ –p ₃₃ second normal stress difference, (N/m²) - q shear rate, (s^–1) - t, t time, (s) - T absolute temperature, (°K) - T ₀ reference temperature, (°K) - x the ratiot/ - x position vector of a material point after deformation, (m) - x position vector of a material point before deformation, (m) - ₀, ₁ constants in eq. [37] - ₀, ₁ constants in eq. [37] - shear deformation - (t, t) time dependent shear deformation - _ij unity tensor - n flow birefringence in the 1–2 plane - (q) non-Newtonian shear viscosity, (N s/m²) - ^* () complex dynamic viscosity, (N s/m²) - | ^* ()| absolute value of complex dynamic viscosity, (N s/m²) - () real part of complex dynamic viscosity, (N s/m²) - () imaginary part of complex dynamic viscosity, (N s/m²) - (t — t) memory function, (N/m² · s) - v number of effective chains per unit of volume, (m^–3) - temperature dependent density, (kg/m³) - ₀ density at reference temperatureT ₀, (kg/m³) - relaxation time, (s) - integration variable, (s) - (x) approximate intensity function - ₁ (x) error function - extinction angle - _m orientation angle of the stress ellipsoid - circular frequency, (s^–1) - 1 direction of flow - 2 direction of the velocity gradient - 3 indifferent direction - t time dependence The present investigation has been carried out under the auspices of the Netherlands Organization for the Advancement of Pure Research (Z. W. O.).North Atlantic Treaty Organization Science Post Doctoral Fellow.Research Fellow, Delft University of Technology.With 11 figures and 2 tables     

Experimental investigation of the creeping motion of a drop in a vertical tube

New experimental data regarding the motion of a drop along the axis of a vertical tube, filled with another highly viscous liquid, are obtained. The experiments are realised with sufficiently large drops for an internal circulation to develop and also for different pairs of fluids; the preponderant role of the gravity on the drop shape and consequently on its terminal velocity is pointed out. Moreover, by means of a visualization technique, details on the flow both inside and outside the drop are given.List of symbols g gravity acceleration - r distance from the drop center - R equivalent radius of the drop, i.e. the radius of the sphere having the same volume as the drop - R _EQ radius of the equatorial section of the drop - R _T tube radius - L _AX half length of the drop - U ₀ terminal velocity of the drop - P _s Poiseuille number= U _{0 e} /4 g R ² - Fr Fronde number = U ₀ ² _e /2 g R - Re Reynolds number = 2 U ₀ R _e / _e - E _o Eötvös number = 4g R ²/ - deformation parameter = _e U ₀/ - apparent density of the suspended liquid= | _i – _e | - _i viscosity of the suspended liquid - _e viscosity of the suspending liquid - drop-to-tube radius ratio = R/R _T - _EQ equatorial drop-to-tube radius ratio = R EQ/R _T - interfacial tension     

Transport in ordered and disordered porous media IV: Computer generated porous media for three-dimensional systems

In the method of volume averaging, the difference between ordered and disordered porous media appears at two distinct points in the analysis, i.e. in the process of spatial smoothing and in the closure problem. In theclosure problem, the use of spatially periodic boundary conditions isconsistent with ordered porous media and the fields under consideration when the length-scale constraint,r ₀L is satisfied. For disordered porous media, spatially periodic boundary conditions are an approximation in need of further study.In theprocess of spatial smoothing, average quantities must be removed from area and volume integrals in order to extractlocal transport equations fromnonlocal equations. This leads to a series of geometrical integrals that need to be evaluated. In Part II we indicated that these integrals were constants for ordered porous media provided that the weighting function used in the averaging process contained thecellular average. We also indicated that these integrals were constrained by certain order of magnitude estimates for disordered porous media. In this paper we verify these characteristics of the geometrical integrals, and we examine their values for pseudo-periodic and uniformly random systems through the use of computer generated porous media.

Nomenclature

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 _i i=1, 2, 3 gaussian probability distribution used to locate the position of particles - I unit tensor - L general characteristic length for volume averaged quantities, m - L characteristic length for , m - L characteristic length for , m - characteristic length for the -phase particles, m - ⁰ reference characteristic length for the-phase particles, m - characteristic length for the-phase, m - _i i=1, 2, 3 lattice vectors, m - m convolution product weighting function - m _v special convolution product weighting function associated with the traditional volume average - n _i i=1, 2, 3 integers used to locate the position of particles - n unit normal vector pointing from the-phase toward the-phase - n _e outwardly directed unit normal vector at the entrances and exits of the-phase - r _p position vector locating the centroid of a particle, m - r gaussian probability distribution used to determine the size of a particle, m - r ₀ characteristic length of an averaging region, m - r position vector, m - r _m support of the weighting functionm, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume,, m³ - x positional vector locating the centroid of an averaging volume, m - x ₀ reference position vector associated with the centroid of an averaging volume, m - y position vector locating points 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 - /L, small parameter in the method of spatial homogenization - standard deviation ofa _i - _r standard deviation ofr - _r intrinsic phase average of      

The rheology of polymeric liquid crystals

The molecular theory of Doi has been used as a framework to characterize the rheological behavior of polymeric liquid crystals at the low deformation rates for which it was derived, and an appropriate extension for high deformation rates is presented. The essential physics behind the Doi formulation has, however, been retained in its entirety. The resulting four-parameter equation enables prediction of the shearing behavior at low and high deformation rates, of the stress in extensional flows, of the isotropic-anisotropic phase transition and of the molecular orientation. Extensional data over nearly three decades of elongation rate (10^–2–10¹) and shearing data over six decades of shear rate (10^–2–10⁴) have been correlated using this analysis. Experimental data are presented for both homogeneous and inhomogeneous shearing stress fields. For the latter, a 20-fold range of capillary tube diameters has been employed and no effects of system geometry or the inhomogeneity of the flow-field are observed. Such an independence of the rheological properties from these effects does not occur for low molecular weight liquid crystals and this is, perhaps, the first time this has been reported for polymeric lyotropic liquid crystals; the physical basis for this major difference is discussed briefly. A Semi-empirical constant in eq. (18), N/m² - c rod concentration, rods/m³ - c ^* critical rod concentration at which the isotropic phase becomes unstable, rods/m³ - C interaction potential in the Doi theory defined in eq. (3) - d rod diameter, m - D semi-empirical constant in eq. (19), s^–1 - D _r lumped rotational diffusivity defined in eq. (4), s^–1 - rotational diffusivity of rods in a concentrated (liquid crystalline) system, s^–1 - D _ro rotational diffusivity of a dilute solution of rods, s^–1 - f distribution function defining rod orientation - F tensorial term in the Doi theory defined in eq. (7) (or eq. (19)), s^–1 - G tensorial term in the Doi theory defined in eq. (8) - K _B Boltzmann constant, 1.38 × 10^–23 J/K-molecule - L rod length, m - S scalar order parameter - S tensor order parameter defined in eq. (5) - t time, s - T absolute temperature, K - u unit vector describing the orientation of an individual rod - rate of change ofu due to macroscopic flow, s^–1 - v fluid velocity vector, m/s - v velocity gradient tensor defined in eq. (9), s^–1 - V mean field (aligning) potential defined in eq. (2) - x coordinate direction, m - Kronecker delta (= 0 if = 1 if = ) - _r ratio of viscosity of suspension to that of the solvent at the same shear stress - _s solvent viscosity, Pa · s - ^* viscosity at the critical concentrationc ^*, Pa · s - v ₁, v₂ numerical factors in eqs. (3) and (4), respectively - deviatoric stress tensor, N/m² - volume fraction of rods - ⁰ constant in eq. (16) - ^* volume fraction of rods at the critical concentrationc ^* - average over the distribution functionf(u, t) (= d ²u f(u, t)) - gradient operator - d ²u integral over the surface of the sphere (|u| = 1)     

Numerical computation of high Rayleigh number natural convection and prediction of hot radiator induced room air motion

The two-dimensional stationary turbulent buoyant flow and heat transfer in a cavity at high Rayleigh numbers was computed numerically. The k– turbulence model was used. The time-averaged equations for momentum, energy and continuity, which are coupled to the turbulence equations, were solved using a finite difference formulation. In order to validate the computer code, a comparison exercise was carried out. The test results are in good agreement with the internationally accepted benchmark solution. Grid-refinement shows the necessity of a very fine grid at high Rayleigh numbers with especially small grid-distances in the near-wall region. The computed boundary layer velocity profiles are in excellent agreement with available experimental data. The local heat transfer in the turbulent part of the boundary layers is predicted 20% too high. Computations were carried out for the natural convective flow in a room induced by a hot radiator and a cold window. Various radiator configurations and types of thermal boundary conditions were applied including thermal radiation interaction between surfaces.Nomenclature a thermal diffusivity (m²/s) - C constant in _t expression - D cavity dimensions (m) - g acceleration of gravity (m/s²) - G _k production/destruction of k by buoyancy (kg/ms³) - h enthalpy (J/kg) - IX index of grid point - k turbulent kinetic energy (m²/s²) - m dimensionless stratification parameter - Nu overall Nusselt number - Nu _y local Nusselt number - NX total number of grid points - p pressure (N/m²) - P _k production of k by shear stress (kg/ms³) - Q heat flux through wall (W/m) - Ra overall Rayleigh number - Ra _y local Rayleigh number - Re _t turbulent Reynolds number - S source term in -equation (kg/ms⁴) - S source term for - T _c, T _h temperatures of cold and hot walls (K) - T _s (y) stratification temperature on vertical mid-line (K) - T ₀ mean cavity temperature (K) - u, v horizontal and vertical velocity components (m/s) - u ₀ Brunt-Vaisälä velocity scale (m/s) - x, y horizontal and vertical coordinates (m) - non-linearity parameter for grid - coefficient of thermal expansion (l/K) - jet angle (°) - diffusivity for - S dissipation rate for turbulent kinetic energy (m²/s³) - variable to be solved - thermal conductivity (W/mK) - , _t kinematic and eddy viscosities (m²/s) - stream function (kg/ms) - density (kg/m³) - _k, , _t constants in k– model     

Theoretical investigation of the effect of an internal heat source on turbulent liquid metal heat transfer in a channel with uniform heat flux

Summary The effect of an internal heat source on the heat transfer characteristics for turbulent liquid metal flow between parallel plates is studied analytically. The analysis is carried out for the conditions of uniform internal heat generation, uniform wall heat flux, and fully established temperature and velocity profiles. Consideration is given both to the uniform or slug flow approximation and the power law approximation for the turbulent velocity profile. Allowance is made for turbulent eddying within the liquid metal through the use of an idealized eddy diffusivity function. It is found that the Nusselt number is unaffected by the heat source strength when the velocity profile is assumed to be uniform over the channel cross section. In the case of a 1/7-power velocity expression, the Nusselt numbers are lower than those in the absence of internal heat generation, and decrease with diminishing eddy conduction. Nusselt numbers, in the absence of an internal heat source, are compared with existing calculations, and indications are that the present results are adequate for preliminary design purposes.Nomenclature A hydrodynamic parameter - a half height of channel - a ₁ a constant, 1+0.01 Pr Re ^0.9 - a ₂ a constant, 0.01 Pr Re ^0.9 - C _p specific heat at constant pressure - D _h hydraulic diameter of channel, 4a - h heat transfer coefficient, q _w/(t _w–t _b) - I ₁ integral defined by (17) - I ₂ integral defined by (18) - k diffusivity parameter, (1+0.01 Pr Re ^0.9)^1/2 - m exponent in power velocity expression - Nu Nusselt number, hD _h/ - Nu ₀ Nusselt number in absence of internal heat generation - Pr Prandtl number, / - Q heat generation rate per volume - q _w wall heat flux - Re Reynolds number for channel, 2/ - s ratio of heat generation rate to wall heat flux, Qa/q _w - T dimensionless temperature, (t _w–t)/(t _w–t _b) - t fluid temperature, t _w wall temperature, t _b fluid bulk temperature - u fluid velocity in x direction, , fluid mean velocity - x longitudinal coordinate measured from channel entrance - x ⁺ dimensionless longitudinal coordinate, 2(x/a)/Pr Re - y transverse coordinate measured from channel centerline - z transverse coordinate measured from channel wall, a–y - molecular diffusivity of heat, /C _p - dummy variable of integration - dummy variable of integration - _H eddy diffusivity of heat - _M eddy diffusivity of momentum - dummy variable of integration - fluid thermal conductivity - _T dimensionless diffusivity, Pr ( _H/) - fluid kinematic viscosity - dummy variable of integration - fluid density - dummy variable of integration - ratio of eddy diffusivity for heat transfer to that for momentum transfer, _H/ _M - average value of - dimensionless velocity distribution, u/     

Magnetohydrodynamic natural convection heat transfer from radiate vertical porous surfaces

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

Asymptotic Nusselt numbers for Newtonian flow through a parallel-plate channel with recycle

General expressions for evaluating the asymptotic Nusselt number for a Newtonian flow through a parallel-plate channel with recycle at the ends have been derived. Numerical results with the ratio of thicknesses as a parameter for various recycle ratios are obtained. A regression analysis shows that the results can be expressed by log Nur^0.83=0.3589 (log)² -0.2925 (log) + 0.3348 forR 3, 0.1 0.9; logNu=0.5982(log)² +0.3755 × 10^–2 (log) +0.8342 forR 10^–2, 0.1 0.9.

---

Asymptotische Nusselt-Zahlen für die Newtonsche Strömung durch einen Kanal aus parallelen Platten mit Rückführung

Zusammenfassung In dieser Untersuchung wurden allgemeine Ausdrücke hergeleitet um die asymptotische Nusselt-Zahl für eine Newtonsche Strömung durch einen Kanal aus parallelen Platten mit Rückführung an den Enden berechnen zu können. Es wurden numerische Ergebnisse mit den Dicken-Verhältnissen, als Parameter für verschiedene Rückführungs-verhältnisse, erhalten. Eine Regressionsanalyse zeigt, daß die Ergebnisse wie folgt ausgedrückt werden können: log Nur^0,83=0,3589 (log)² -0,2925 (log) + 0,3348 fürR 3, 0,1 0,9; logNu=0,5982(log)² +0,3755 × 10^–2 (log) + 0,8342 fürR 10^–2, 0,1 0,9.

---

Nomenclature A₁ shooting value,d(0)/d - A₂ shooting value,d(1)/d - B channel width - Gz Graetz number, U_bW²/L - h _m logarithmic average convective heat transfer coefficient - h _x average local convective heat transfer coefficient - k thermal conductivity - L channel length - Nu average local Nusselt number, 2 h_xW/k - Nu _m logarithmic average Nusselt number, 2h_mW/k - R recycle ratio, reverse volume flow rate divided by input volume flow rate - T temperature of fluid - T _m bulk temperature, Eq. (8) - T ₀ temperature of feed stream - T _s wall temperature - U velocity distribution - U _b reference velocity,V/BW - V input volume flow rate - v dimensionless velocity distribution, U/U_b - W channel thickness - x longitudinal coordinate - y transversal coordinate - Z₁-z₆ functions defined in Eq. (A1) - thermal diffusivity - least squares error, Eq. (A7) - weight, Eqs. (A8), (A9) - dimensionless coordinate,y/W - dimensionless coordinate,x/GzL - function, Eq. (7)     

Natural convection heat transfer through an enclosed horizontal layer of supercritical carbon dioxide

In natural convection heat transfer through a thin horizontal layer of carbon dioxide, maxima in the equivalent thermal conductivities are obtained in the vicinity of the respective pseudocritical temperatures at pressures of 75.8, 89.6 and 103.4 bar. The maxima are the more pronounced, the closer the critical point is approached.Comparison of experimental results with Nusselt equations shows good agreement except for the immediate vicinity of the pseudocritical temperature.In visual observations a distinct change in flow structure appears in the immediate vicinity of the pseudocritical temperature. A steady state polygon pattern and a boiling-like action could not be observed in this geometry.

---

Zusammenfassung Beim Wärmetransport durch freie Konvektion in einer dünnen waagerechten Schicht von Kohlendioxid ergaben sich Maxima der scheinbaren Wärmeleitfähigkeit in der Nähe der pseudokritischen Temperaturen bei Drükken von 75,8, 89,6 und 103,4 bar. Die Maxima sind um so ausgeprägter, je mehr man sich dem kritischen Punkt nähert.Ein Vergleich der Versuchsergebnisse mit Nusseltbeziehungen ergibt gute Übereinstimmung außer in unmittelbarer Umgebung der pseudokritischen Temperatur. Direkte Beobachtungen der Konvektionsmuster zeigen in unmittelbarer Umgebung der pseudokritischen Temperatur eine deutliche Strukturänderung. Ein stationäres Zellmuster und siedeähnliche Vorgänge konnten in dieser Anordnung nicht beobachtet werden.

---

Nomenclature A area of the heating or cooling plate - C constant in the correlation - g acceleration of gravity - h heat transfer coefficient - k thermal conductivity of fluid in the gap - k _e equivalent thermal conductivity - m, n exponents of dimensionless numbers - q heat flux - T _C,PC absolute temperature; critical C, pseudocritical PC - Gr Grashof numberg ( _h– _c) ³/ ² - Nu Nusselt numberh/k - Pr Prandtl number/ - thermal diffusivity - coefficient of volume expansion - width of gap - _c,h temperature of cooling (c)-, heating (h)-plate - _m arithmetic mean temperature ( _c+ _h)/2 - kinematic viscosity - _c,h fluid density at the temperature of the cooling (c)- or heating (h)-plate - heat flow rate through the gap     

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     

Helical flow of molten polymers in a cylindrical annulus

Summary The paper is concerned with an analytical investigation of helical flow of a non-Newtonian fluid through an annulus with a rotating inner cylinder. The shear dependence of viscosity is described by a power law and the temperature dependence by an exponential function.Velocity and temperature profiles, energy input and shear along the stream lines, pressure drop, and torque are presented for the range of input parameters encountered in polymer extrusion.The results of the study can be applied to a mixing element in a screw extruder and for a device to control extrudate temperature and output.Nomenclature a thermal diffusivity [m²/s] - b temperature coefficient [K^–1], see eq. [4] - c heat capacity [J/kg K] - h slot width [m] - I ₁,I ₂,I ₃ invariants of the rate of deformation tensor, see eq. [5] - k thermal conductivity [J/m s K] - l, L = 1/h length of the slot - l _T ,l _K thermal and kinematic entrance length - m power law exponent, see eq. [3] - M torque [m N] - p pressure [N/m²] - P dimensionless pressure gradient, see eq. [24] - P _R,P _RZ dimensionless components of the shear stress tensor, see eq. [25] and eq. [26] - r, R = r/r _wa radial coordinate - r _wa, r_wi outer and inner radius of annulus [m] - t time [s]; dwell time in the annulus - T temperature [K] - v , v_r, V_z velocity components [m/s] - v ₀ angular velocity at inner wall [m/s] - average velocity inz-direction [m/s] - V , V_R, V_Z dimensionless velocity components,v /v₀, v_r/v₀, v_z/v₀ - V _z velocity ratio, helical parameter - Y coordinate inr-direction, see eq. [20] - z, Z = z/h Pe axial coordinate - deformation - rate of deformation tensor [s^–1] - apparent viscosity [N s/m²], see eq. [3] - dimensionless temperature,b (T – T ₀) - azimuth coordinate - ratio of radii,r _wi/r_wa - density [kg/m³] - , kl shear stress tensor [N/m²] - fluidity [m²w/Nw s], see eq. [4] - Gf Griffith number, see eq. [12] - Pe Péclet number, see eq. [13] - Re Reynolds number, - 0 initial state, reference state - equilibrium state - e entrance - wi, wa at surface of inner or outer wall - r, R, z, Z, coordinates - i, j radial and axial position of nodal point in the grid - k, l tensor components Presented at Euromech 37, Napoli 6. 20–23. 1972.With 15 figuresDedicated to Prof. Dr.-Ing. G. Schenkel on his 60th birthday     

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

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

Instantaneous visualization of surface flows

An experimental system is described for visualizing the surface flow of a wing, using an oil smoke tracer technique. The method leads to the determination of the instantaneous velocity direction at the output of surface injectors. A preliminary investigation is made on a flat plate to optimize the conditions of oil smoke injection. Then, the visualization is performed on the upperside of a sweptback wing in the vicinity of the reattachment of the vortex flow. This visualization technique can be applied to other types of wall flows — separated or not — around various bodies.List of symbols b wing span - c _n normal (to leading edge) chord - c _r streamwise (or root) chord - d diameter of the injectors - distance from the apex along the leading edge - relative distance from the apex along the leading edge ( = /C _d) - sweep angle - e injector geometric parameter (e = d/l) - angle of attack - K injection parameter - l length of the injectors - v kinematic viscosity - P _t, P_s total and static pressure of the flow - P _inj injection pressure - P _r reduced pressure (P _r = (P_inj – P_t)/(P_t – P_s)) - Re flow Reynolds number (Re = V ·c _n/v) - Re _i injector Reynolds number (Re = V ·d/v) - s curvilinear distance along c _d - s relative curvilinear distance along c _d(s = s/c _d) - V infinite upstream flow velocity     

The interaction of two opposing,asymmetric curved wall jets

An experimental investigation was made of a two dimensional flow formed by the interaction of two asymmetric turbulent curved wall jets past a circular cylinder. Measurements were made of velocity and turbulence intensity profiles of the two curved wall jets before the interaction, and those of the merged jet after the interaction. The location of the interaction region of the two opposing curved wall jets and the flow direction of the merged jet were found to depend primarily on the ratio of initial momentum fluxes. The velocity and turbulence intensity profiles of the merged jet were similar to those of the plane turbulent jet. However, the growth rate of the merged jet was approximately 1.5 times larger than that of the plane jet. The influence of the momentum flux ratio on the growth rate appeared to be insignificant.List of symbols C _f friction coefficient - h slot height - J _p, J _c initial momentum flux of a power jet and of a control jet, respectively - P, Pa wall static and atmospheric pressure, respectively - Re Reynolds number based on slot height - Re _m local Reynolds number U _m y _m /v - U local mean velocity - U _c velocity along the center line of the merged jet - U _m local maximum velocity of the curved wall jet - u r.m.s. value of velocity fluctuations - u _u friction velocity - U ⁺ U/u_t - x distance along the cylinder surface - x distance along the center line of the merged jet - y _1/2, y _1/2 position of y and y where U = U _m /2 and U = U _c /2, respectively - y ⁺ yu _t/V - deflection angle of the merged jet (Fig. 4) - interaction angle (Fig. 4) - merged jet angle (Fig. 4) - angle measured from the center line of the cylinder (Fig. 4) - interception angle (Fig. 8) - , normalized coordinates, y/y _1/2 and y/y _1/2, respectively     

The very slow flow of a Powell-Eyring fluid around a sphere

Summary The very slow flow of a Powell-Eyring type non-Newtonian fluid around a sphere is investigated by a variational technique. The result, a correction factor that is applied to the Stokes' equation, is given as a plot and as an equation which is empirically fit to the plot. Also, a comparison of the very slow flows of a simplified viscoelastic Oldroyd fluid and the Powell-Eyring fluid is made which indicates that in a certain restricted region of the very slow flows, both models give essentially the same results. The Oldroyd and Powell-Eyring model parameters are interrelated by forcing both models to fit the same tube flow viscosity data.Nomenclature B dimensionless quantity, v /R - C dimensionless second invariant - c ₁ constant determined by variational method - D dimensionless variational integral - D _2j , D _j+k position-independent variables used in specification of trial functions - E _2j , E _j+k position-independent variables used in specification of trial functions - f friction factor - f _corr friction factor correction - F _drag drag force on sphere - g, g ₀, g ₁ general trial function; first and second terms in the general trial function - G, H terms in the expression for C - j index - J variational integral - k index - K term in the expression for C - p, q integers - r integer, radial coordinate - R radius of sphere - Re Reynolds number - Re ₀ Reynolds number at point of zero shear rate - Re Reynolds number at infinite distance from sphere - Re _NN Reynolds number based on variable part of viscosity - u, v dimensionless position coordinates - V volume considered - v _i ith velocity component - v _r , v , v _z velocity components in the r, , and z-directions - v approach velocity of the fluid - x/ parameter in Powell-Eyring model - x _i i-position coordinate - parameter in Powell-Eyring model - rate of deformation - , _c , _N , ₀ coefficient of viscosity; cross viscosity; parameter in Powell-Eyring model; viscosity in limit of zero shear rate - spherical coordinate - , _ij rate of deformation tensor; ij-component of rate of deformation tensor - ₁, ₂ parameters in Oldroyd model - Newtonian viscosity - ₁, ₂ parameters in Oldroyd model - dimensionless radial coordinate, r/R - second invariant - fluid density - spherical coordinate - stream function