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

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

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 dynamic and steady shear properties of synovial fluid and of the components making up synovial fluid

Steady-shear and dynamic properties of a pooled sample of cattle synovial fluid have been measured using techniques developed for low viscosity fluids. The rheological properties of synovial fluid were found to exhibit typical viscoelastic behaviour and can be described by the Carreau type A rheological model. Typical model parameters for the fluid are given; these may be useful for the analysis of the complex flow problems of joint lubrication.The two major constituents, hyaluronic acid and proteins, have been successfully separated from the pooled sample of synovial fluid. The rheological properties of the hyaluronic acid and the recombined hyaluronic acid-protein solutions of both equal and half the concentration of the constituents found in the original synovial fluid have been measured. These properties, when compared to those of the original synovial fluid, show an undeniable contribution of proteins to the flow behaviour of synovial fluid in joints. The effect of protein was found to be more prominent in hyaluronic acid of half the normal concentration found in synovial fluid, thus providing a possible explanation for the differences in flow behaviour observed between synovial fluid from certain diseased joints compared to normal joint fluid.Nomenclature A Ratio of angular amplitude of torsion head to oscillation input signal - G Storage modulus - G Loss modulus - I Moment of inertia of upper platen — torsion head assembly - K Restoring constant of torsion bar - N ₁ First normal-stress difference - R Platen radius - S ⁽ⁱ⁾ Geometric factor in the dynamic property analysis - t ₁ Characteristic time parameter of the Carreau model - X, Y Carreau model parameters - Z () Reimann Zeta function of - Carreau model parameter - Shear rate - Apparent steady-shear viscosity - ^* Complex dynamic viscosity - Dynamic viscosity - Imaginary part of the complex dynamic viscosity - ₀ Zero-shear viscosity - ₀ Cone angle - Carreau model characteristic time - Density of fluid - Shear stress - Phase difference between torsion head and oscillation input signals - ₀ Zero-shear rate first normal-stress coefficient - Oscillatory frequency     

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

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

Heat transfer between gas-liquid spray stream flowing perpendicularly to the row of the cylinders

Experimental investigation and analysis of heat transfer process between a gas-liquid spray flow and the row of smooth cylinders placed in the surface perpendicular to the flow has been performed. Among others, there was taken into account in the analysis the phenomenon of droplets bouncing and omitting the cylinder as well as the phenomenon of the evaporation process from the liquid film surface.In the experiments test cylinders were used, which were placed between two other cylinders standing in the row.From the experiments and the analysis the conclusion can be drawn that the heat transfer coefficients values for a row of the cylinders are higher than for a single cylinder placed in the gasliquid spray flow.

---

Wärmeübergang an eine senkrecht anf eine Zylinderreihe auftreffende Gas-Flüssigkeits-Sprüh-Strömung

Zusammenfassung Es wurden Messungen und theoretische Analysen des Wärmeübergangs zwischen einer Gas-FlüssigkeitsSprüh-Strömung und den glatten Oberflächen einer Zylinderreihe durchgeführt, die senkrecht zum Sprühstrahl angeordnet waren. Dabei wurde in der Analyse unter anderem das Phänomen betrachtet, daß die Tropfen die Zylinderwand treffen und verfehlen können und daß sich ein Verdampfungsprozeß aus dem flüssigen Film an der Zylinderoberfläche einstellt.Gemessen wurde an einem zwischen zwei Randzylindern befindlichen Zylinder.Die Experimente und die Analyse gestatten die Schlußfolgerung, daß der Wärmeübergangskoeffizient für eine Zylinderreihe höher ist als für einen einzelnen Zylinder in der Sprühströmung.

---

Nomenclature a distance between axes of cylinders, m - c _l specific heat capacity of liquid, J/kg K - c _g specific heat capacity of gas, J/kg K - D cylinder diameter, m - g _l mass velocity of liquid, kg/m²s - ¯k average volume ratio of liquid entering film to amount of liquid directed at the cylinder in gas-liquid spray flow, dimensionless - k() local volume ratio of liquid entering film to amount of liquid directed at the cylinder in gas-liquid spray flow, dimensionless - L specific latent heat of vaporisation, J/kg - m mass fraction of water in gas-liquid spray flow, dimensionless - M constant in Eq. (9) - p pressure, Pa - p _g statical pressure of gas, Pa - p _w pressure of gas on the cylinder surface, Pa - p external pressure on the liquid film surface, Pa - r cylindrical coordinate, m - R radius of cylinder, m - T temperature, K, °C - T _l, t_l liquid temperature in the gas-liquid spray, K, °C - T _w,t_w temperature of cylinder surface, K, °C - T temperature of gas-liquid film interface, K - U liquid film velocity, m/s - w gas velocity on cylinder surface, m/s - w _g gas velocity in free stream, m/s - W _l liquid vapour mass ratio in free stream, dimensionless - W liquid vapour mass ratio at the edge of a liquid film, dimensionless - x coordinate, m - y coordinate, m - z complex variable, dimensionless - average heat transfer coefficient, W/m²K - local heat transfer coefficient, W/m² K - average heat transfer coefficient between cylinder surface and gas, W/m² K - _g, local heat transfer coefficient between cylinder surface and gas, W/m² K - mass transfer coefficient, kg/m²s - liquid film thickness, m - _lg dynamic diffusion coefficient of liquid vapour in gas, kg/m s - pressure distribution function on a cylinder surface - function defined by Eq. (3) - _l liquid dynamic viscosity, kg/m s - _g gas dynamic viscosity, kg/m s - cylindrical coordinate, rad, deg - _l thermal conductivity of liquid, W/m K - _g thermal conductivity of gas, W/m K - mass transfer driving force, dimensionless - _l density of liquid, kg/m³ - _g density of gas, kg/m³ - _w shear stress on the cylinder surface, N/m² - _w shear stress exerted by gas at the liquid film surface, N/m² - air relative humidity, dimensionless - T -T _w - _w =T _wT_l Dimensionless parameters I= enhancement factor of heat transfer - m ^*=M _l/M_g molar mass of liquid to the molar mass of gas ratio - Nu _g= D/ _g gas Nusselt number - Pr _g=c _{g g}/_g gas Prandtl number - Pr _l=c_lll liquid Prandtl number - ¯r=(r–R)/ dimensionless coordinate - Re _g=w_gD _g/_g gas Reynolds number - Re _g,max=w_g,max D _g/_g gas Reynolds number calculated for the maximal gas velocity between the cylinders - Sc=m ^* _g/_l–g Schmidt number =/R dimensionless film thickness     

A non-dualistic unified field theory of gravitation,electromagnetism and spin

The wisdom of classicalunified field theories in the conceptual framework of Weyl, Eddington, Einstein and Schrödinger has often been doubted and in particular there does not appear to be any empirical reason why the Einstein-Maxwell (E-M) theory needs to be geometrized. The crux of the matter is, however not whether the E-M theory is aesthetically satisfactory but whether it answers all the modern questions within the classical context. In particular, the E-M theory does not provide a classical platform from which the Dirac equation can be derived in the way Schrödinger's equation is derived from classical mechanics via the energy equation and the Correspondence Principle. The present paper presents a non-dualistic unified field theory (UFT) in the said conceptual framework as propounded by M. A. Tonnelat. By allowing the metric formds ²=g dx ^v x ^v and the non-degenerate two-formF=(1/2> _l) dx ^vdx ^vto enter symmetrically into the theory we obtain a UFT which contains Einstein's General Relativity and the Born-Infeld electrodynamics as special cases. Above all, it is shown that the Dirac equation describing the electron in an external gravito-electromagnetic field can be derived from the non-dualistic Einstein equation by a simple factorization if the Correspondence Principle is assumed.     

Concentration-dependent behaviour of the shear viscosity of coal-fuel oil suspensions

A function correlating the relative viscosity of a suspension of solid particles in liquids to their concentration is derived here theoretically using only general thermodynamic ideas, with out any consideration of microscopic hydrodynamic models. This function ( _r = exp (1/2B ^* C ²)) has a great advantage over the many different functions proposed in literature, for it depends on a single parameter,B ^*, and is therefore concise. To test the validity of this function, a least-squares regression analysis was undertaken of available data on the viscosity and concentration of suspensions of coal particles in fuel oil, which promise to be a useful alternative to fuel oil in the near future. The proposed function was found to accurately describe the concentration-dependent behaviour of the relative viscosity of these suspensions. Furthermore, an attempt was made to obtain information about the factors affecting the value ofB ^*, however the results were only qualitative because of, among other things, the inaccuracy of the viscosity measurements in such highly viscous fluids. shear viscosity of the suspension - ₀ shear viscosity of the Newtonian suspending medium - _r = /₀ relative viscosity - solid volume concentration - c solid weight concentration - _m maximum attainable volume concentration of solids - solid volume concentration at which the relative viscosity of the suspension becomes infinite - c _m maximum attainable solid weight concentration - _s density of the solid phase - _l density of the liquid phase - _m density of the suspension - k _n coefficients of theø-power series expansion of _r - { _j } sets of parameters specifying the thermodynamic state of the solid phase of a suspension - T absolute temperature (K) - f (c, T, _j) formal expression for the relative variation of the viscosity with concentration = [1 / (/c)] _T,j - d median size of the granulometric distribution - _B plastic or Bingham viscosity - K consistency factor - n flow index - g ([c _m –c],T, _j ) function including an asymptotic divergence asc tends toc _m , formally describing the concentration dependent behaviour of the shear viscosity of a suspension - A (T, _j) regression analysis parameters - B (T, _j) regression analysis parameters - B ^* (T, _j ) regression analysis parameters     

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     

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

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

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     

Rheological properties of reactive polymer melts.

A new method for describing the rheological properties of reactive polymer melts, which was presented in an earlier paper, is developed in more detail. In particular, a detailed derivation of the equation of a first-order rheometrical flow surface is given and a procedure for determining parameters and functions occurring in this equation is proposed. The experimental verification of the presented approach was carried out using our data for polyamide-6.Notation E Dimensionless reduced viscosity, eq. (34) - E ₀ Newtonian asymptote of the function (36) - E power-law asymptote of the function (36) - E _{= 1} the value ofE at = 1 - k degradation reaction rate constant, s^–1 - k ₁ rate constant of function (t), eq. (26), s^–1 - k ₂ rate constant of function (t), eq. (29), s^–1 - K(t) residence-time-dependent consistency factor, eq. (22) - M _w weight-average molecular weight - M _x x-th moment of the molecular weight distribution - R gas constant - S _x M _x /M _w - t residence time in molten state, s - t _j thej-th value oft, s - T temperature, K - % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xd9vqpe0x% c9q8qqaqFn0dXdir-xcvk9pIe9q8qqaq-xir-f0-yqaqVeLsFr0-vr% 0-vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaaieGaceWFZo% Gbaiaaaaa!3B4E!\[\dot \gamma \] shear rate, s^–1 - _i thei-th value of , s^–1 - _r =1 the value of at = 1, s^–1 - ^* reduced shear rate, eq. (44), s^–1 - dimensionless reduced shear rate, eq. (35) - viscosity, Pa · s - shear-rate and residence-time dependent viscosity, Pa · s - zero-shear-rate degradation curve - degradation curve at - _{t0 (t)} zero-residence-time flow curve - Newtonian asymptote of the RFS - instantaneous flow curve - power-law asymptote of the RFS - _0,0 zero-shear-rate and zero-residence-time viscosity, Pa · s - _E=1 value of viscosity atE=1, Pa · s - ^* reduced viscosity, eq. (43), Pa · s - zero-residence-time rheological time constant, s - density, kg/m³ - (t),(t) residence time functions     

Injection moulding of thermoplastics: Freezing of variable-viscosity fluids.

The injection moulding of thermoplastics involves, during mould filling, flows of hot polymer melts into mould networks, the walls of which are so cold that frozen layers form on them. An analytical study of such flows is presented here for the case when the Graetz and Nahme numbers are large and the Pearson number is small. Thus the flows are developing and temperature differences due to heat generation by viscous dissipation are sufficiently large to cause significant variations in viscosity (but the difference between the entry temperature of the polymer to a specific part of the mould network and the melting temperature of the polymer is not). Br Brinkman number - Gz Graetz number - h half-height of channel or disc - h ^* half-height of polymer melt region in channel or disc - L length of channel or pipe - m viscosity shear-rate exponent - Na Nahme number - p pressure - P pressure drop - Pe Péclet number - Pn Pearson number - Q volumetric flowrate - r radial coordinate in pipe or disc - R radius of pipe - Re Reynolds number - R _i inner radius of disc - R _o outer radius of disc - R ^* radius of polymer melt region in pipe - T temperature - T _ad adiabatic temperature rise - T _e entry polymer melt temperature - T _m melting temperature of polymer - T _max maximum temperature - T ₀ reference temperature - T _w wall temperature - flow-average temperature rise - u _r radial velocity in pipe or disc - u _x axial velocity in channel - u _y transverse velocity in channel or disc - u _z axial velocity in pipe - w width of channel - x axial coordinate in channel or modified radial coordinate in disc - y transverse coordinate in channel or disc - z axial coordinate in pipe - thermal conductivity of molten polymer - thermal conductivity of frozen polymer - scaled dimensionless axial coordinate in channel or pipe or radial coordinate in disc - ₀ undetermined integration constant - heat capacity of molten polymer - viscosity temperature exponent - dimensionless transverse coordinate in channel or disc - ^* dimensionless half-height of polymer melt region in channel or disc - H ^* scaled dimensionless half-height of polymer melt region in channel or disc or radius of polymer melt region in pipe - dimensionless temperature - ^* dimensionless wall temperature - scaled dimensionless temperature - numerical constant - µ viscosity of molten polymer - µ ₀ consistency of molten polymer - dimensionless pressure gradient - scaled dimensionless pressure gradient - density of molten polymer - dimensionless radial coordinate in pipe or disc - _i dimensionless inner radius of disc - ^* dimensionless radius of polymer melt region in pipe - dimensionless streamfunction - scaled dimensionless streamfunction - dummy variable - streamfunction - similarity variable - similarity variable     

Squeezing flows of Newtonian liquid films an analysis including fluid inertia

A theoretical study is made of the flow behavior of thin Newtonian liquid films being squeezed between two flat plates. Solutions to the problem are obtained by using a numerical method, which is found to be stable for all Reynolds numbers, aspect ratios, and grid sizes tested. Particular emphasis is placed on including in the analysis the inertial terms in the Navier-Stokes equations.Comparison of results from the numerical calculation with those from Ishizawa's perturbation solution is made. For the conditions considered here, it is found that the perturbation series is divergent, and that in general one must use a numerical technique to solve this problem.Nomenclature a half of the distance, or gap, between the two plates - a ₀ the value of a at time t=0 - adot da/dt - ä d² a/dt ² - d³ a/dt ³ - a ⁱ components of a contravariant acceleration vector - f unknown function of z ₀ and t defined in (6) - f ⁱ function defined in (9) f ¹=r ₀ g(z ₀, t) f ²= ₀ f ³=f(z ₀, t) - F force applied to the plates - g unknown function of z ₀ and t defined in (6) - g g/z ₀ - h grid dimension in the z ₀ direction (see Fig. 5) - Christoffel symbol - i, j, k, l indices - k grid dimension in the t direction (see Fig. 5) - r radial coordinate direction defined in Fig. 1 - r ₀ radial convected coordinate - R radius of the circular plates - t time - v _r fluid velocity in the r direction - v _z fluid velocity in the z direction - v fluid velocity in the direction - x ⁱ cylindrical coordinate x ¹=r x²= x³=z - z vertical coordinate direction defined in Fig. 1 - z ₀ vertical convected coordinate - tangential coordinate direction - ₀ tangential convected coordinate - viscosity - kinematic viscosity, / - ⁱ convected coordinate ¹=r₀ ²=₀ ³=z₀ - density     

Darcy-Forchheimer natural,forced and mixed convection heat transfer in non-Newtonian power-law fluid-saturated porous media

The governing equation for Darcy-Forchheimer flow of non-Newtonian inelastic power-law fluid through porous media has been derived from first principles. Using this equation, the problem of Darcy-Forchheimer natural, forced, and mixed convection within the porous media saturated with a power-law fluid has been solved using the approximate integral method. It is observed that a similarity solution exists specifically for only the case of an isothermal vertical flat plate embedded in the porous media. The results based on the approximate method, when compared with existing exact solutions show an agreement of within a maximum error bound of 2.5%.Nomenclature A cross-sectional area - b _i coefficient in the chosen temperature profile - B ₁ coefficient in the profile for the dimensionless boundary layer thickness - C coefficient in the modified Forchheimer term for power-law fluids - C ₁ coefficient in the Oseen approximation which depends essentially on pore geometry - C _i coefficient depending essentially on pore geometry - C _D drag coefficient - C _t coefficient in the expression forK ^* - d particle diameter (for irregular shaped particles, it is characteristic length for average-size particle) - f _p resistance or drag on a single particle - F _R total resistance to flow offered byN particles in the porous media - g acceleration due to gravity - g _x component of the acceleration due to gravity in thex-direction - Grashof number based on permeability for power-law fluids - K intrinsic permeability of the porous media - K ^* modified permeability of the porous media for flow of power-law fluids - l _c characteristic length - m exponent in the gravity field - n power-law index of the inelastic non-Newtonian fluid - N total number of particles - Nux,_D,F local Nusselt number for Darcy forced convection flow - Nux,_D-F,F local Nusselt number for Darcy-Forchheimer forced convection flow - Nux,_D,M local Nusselt number for Darcy mixed convection flow - Nux,_D-F,M local Nusselt number for Darcy-Forchheimer mixed convection flow - Nux,_D,N local Nusselt number for Darcy natural convection flow - Nux,_D-F,N local Nusselt number for Darcy-Forchheimer natural convection flow - pressure - p exponent in the wall temperature variation - _Pe c characteristic Péclet number - _Pe x local Péclet number for forced convection flow - _Pe x modified local Péclet number for mixed convection flow - _Ra c characteristic Rayleigh number - _Ra x local Rayleigh number for Darcy natural convection flow - _Ra x local Rayleigh number for Darcy-Forchheimer natural convection flow - Re convectional Reynolds number for power-law fluids - Reynolds number based on permeability for power-law fluids - T temperature - T _e ambient constant temperature - T w,_ref constant reference wall surface temperature - T _w(X) variable wall surface temperature - T _w temperature difference equal toT w,_ref–T _e - T ₁ term in the Darcy-Forchheimer natural convection regime for Newtonian fluids - T ₂ term in the Darcy-Forchheimer natural convection regime for non-Newtonian fluids (first approximation) - T _N term in the Darcy/Forchheimer natural convection regime for non-Newtonian fluids (second approximation) - u Darcian or superficial velocity - u ₁ dimensionless velocity profile - u _e external forced convection flow velocity - u _s seepage velocity (local average velocity of flow around the particle) - u _w wall slip velocity - U c _M characteristic velocity for mixed convection - U c _N characteristic velocity for natural convection - x, y boundary-layer coordinates - x ₁,y ₁ dimensionless boundary layer coordinates - X coefficient which is a function of flow behaviour indexn for power-law fluids - effective thermal diffusivity of the porous medium - shape factor which takes a value of/4 for spheres - shape factor which takes a value of/6 for spheres - ₀ expansion coefficient of the fluid - _T boundary-layer thickness - T ₁ dimensionless boundary layer thickness - porosity of the medium - similarity variable - dimensionless temperature difference - coefficient which is a function of the geometry of the porous media (it takes a value of 3 for a single sphere in an infinite fluid) - ₀ viscosity of Newtonian fluid - ^* fluid consistency of the inelastic non-Newtonian power-law fluid - constant equal toX(2 ^2–n )/ - density of the fluid - dimensionless wall temperature difference     

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

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

Shear viscosity behavior of highly concentrated suspensions at low and high shear-rates

A method is proposed for determining the shear viscosity behavior of highly concentrated suspensions at low and high shear-rates through the use of a formulation that is a function of three parameters signifying the effects of particle size distribution. These parameters are the intrinsic viscosity [], a parametern that reflects the level of particle association at the initiation of motion and the maximum packing concentration _m. The formulation reduces to the modified Eilers equation withn = 2 for high shear rates. An analytical method was used for the calculation of maximum packing concentration which was subsequently correlated with the experimental values to account for the surface induced interaction of particles with the fluid. The calculated values of viscosities at low and high shear-rates were found to be in good agreement with various experimental data reported in literature. A brief discussion is also offered on the reliability of the methods of measuring the maximum packing concentration. _r = /₀ relative viscosity of the suspension - volumetric concentration of solids - k _n coefficient which characterizes a specific effect of particle interactions - _m maximum packing concentration - _r,0 relative viscosity at low shear-rates - [] intrinsic viscosity - n, n parameter that reflects the level of particle interactions at low and high shear-rates, respectively - _r, relative viscosity at high shear-rates - (_m)_s, (_m)_i, (_m)_l packing factors for small, intermediate and large diameter classes - v _s, v_i, v_l volume fractions of small, intermediate and large diameter classes, respectively - _si, _sl coefficient to be used in relating a smaller to an intermediate and larger particle group, respectively - _is, _il coefficient to be used in relating an intermediate to a smaller and larger particle group, respectively - _ls, _li coefficient to be used in relating a larger to a smaller and intermediate particle group, respectively - _m0 maximum packing concentration for binary mixtures - _m,e measured maximum packing concentration - _m,c calculated maximum packing concentration     

Investigations of mean and turbulent flow characteristics of a two dimensional plane diffuser

This paper reports the investigation of mean and turbulent flow characteristics of a two-dimensional plane diffuser. Both experimental and theoretical details are considered. The experimental investigation consists of the measurement of mean velocity profiles, wall static pressure and turbulence stresses. Theoretical study involves the prediction of downstream velocity profiles and the distribution of turbulence kinetic energy using a well tested finite difference procedure. Two models, viz., Prandtl's mixing length hypothesis and k- model of turbulence, have been used and compared. The nondimensional static pressure distribution, the longitudinal pressure gradient, the pressure recovery coefficient, percentage recovery of static pressure, the variation of U _max/U _bar along the length of the diffuser and the blockage factor have been valuated from the predicted results and compared with the experimental data. Further, the predicted and the measured value of kinetic energy of turbulence have also been compared. It is seen that for the prediction of mean flow characteristics and to evaluate the performance of the diffuser, a simple turbulence model like Prandtl's mixing length hypothesis is quite adequate.List of symbols C ₁ , C ₂ ,C turbulence model constants - F _x body force - k kinetic energy of turbulence - l _m mixing length - L length of the diffuser - u, v, w rms value of the fluctuating velocity - u, v, w turbulent component of the velocity - mean velocity in the x direction - _A average velocity at inlet - U _bar average velocity in any cross section - U _max maximum velocity in any cross section - V mean velocity in the y direction - W local width of the diffuser at any cross section - x, y coordinates - dissipation rate of turbulence - _m eddy diffusivity - Von Karman constant - mixing length constant - _l laminar viscosity - _eff effective viscosity - v kinematic viscosity - density - _k effective Schmidt number for k - effective Schmidt number for - stream function - non dimensional stream function     

Flow through a helically coiled annulus

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

Free convection heat transfer from an isothermal plate with arbitrary inclination

This paper presents a new formulation for the laminar free convection from an arbitrarily inclined isothermal plate to fluids of any Prandtl number between 0.001 and infinity. A novel inclination parameter is proposed such that all cases of the horizontal, inclined and vertical plates can be described by a single set of transformed equations. Moreover, the self-similar equations for the limiting cases of the horizontal and vertical plates are recovered from the transformed equations by setting=0 and=1, respectively. Heated upward-facing plates with positive and negative inclination angles are investigated. A very accurate correlation equation of the local Nusselt number is developed for arbitrary inclination angle and for 0.001 Pr .

---

Wärmeübertragung bei freier Konvektion an einer isothermen Platte mit beliebiger Neigung

Zusammenfasssung Diese Untersuchung stellt eine neue Formulierung der laminaren freien Konvektion von Flüssigkeiten mit einer Prandtl-Zahl zwischen 0,001 und unendlich an einer beliebig schräggestellten isothermen Platte dar. Ein neuer Neigungsparameter wird eingeführt, so daß alle Fälle der horizontalen, geneigten oder vertikalen Platte von einem einzigen Satz transformierter Gleichungen beschrieben werden können. Die unabhängigen Gleichungen für die beiden Fälle der horizontalen and vertikalen Platte wurden für=0 und=1 aus den transformierten Gleichungen wieder abgeleitet. Es wurden erwärmte aufwärtsgerichtete Platten mit positiven und negativen Neigungswinkeln untersucht. Eine sehr genaue Gleichung wurde für die lokale Nusselt-Zahl bei beliebigen Neigungswinkeln und für 0,001 Pr entwickelt.

---

Nomenclature C _p specific heat - f reduced stream function - g gravitational acceleration - Gr local Grashof number,g(T _w –T _w ) x³/v² - h local heat transfer coefficient - k thermal conductivity - n constant exponent - Nu local Nusselt number,hx/k - p pressure - Pr Prandtl number, v/ - Ra local Rayleigh number,g(T _w –T )J x³/v - T fluid temperature - T _w wall temperature - T temperature of ambient fluid - u velocity component in x-direction - v velocity component in y-direction - x coordinate parallel to the plate - y coordinate normal to the plate Greek symbols thermal diffusivity - thermal expansion coefficient - (Ra¦sin¦)^1/4/( Ra cos()^1/5 - pseudo-similarity variable, (y/) - dimensionless temperature, (T–T )/(T _w–T ) - ( Ra cos)^1/5+(Rasin)^1/4 - v kinematic viscosity - 1/[1 +(Ra cos)^1/5/( Ra¦sin)^1/4] - density of fluid - Pr/(1+Pr) - _w wall shear stress - angle of plate inclination measured from the horizontal - stream function - dimensionless dynamic pressure     

Gap Estimates for Double-Well Schrödinger Operators and Quantum Stabilization of Anharmonic Crystals

The spectrum of the Schrödinger operator of a one-dimensional quantum anharmonic oscillator of mass m is studied. This spectrum consists of simple (nondegenerate) eigenvalues E _n , $$n\in {\mathbb N}_0$$ such that n ^– E _n + as n + with a certain > 1. The gap parameter =min _n (E _n – E _n-1) is in the center of the study. It is proven that this parameter is a continuous function of m; its small mass and large mass asymptotics are found. The influence of the dependence of on m on the stability of systems of interacting quantum anharmonic oscillators is briefly discussed.