The effect of surface mass transfer on buoyancy induced flow in a variable porosity medium adjacent to a vertical heated plate

The effect of surface mass transfer on buoyancy induced flow in a variable porosity medium adjacent to a heated vertical plate is studied for high Rayleigh numbers. Similarity solutions are obtained within the frame work of boundary layer theory for a power law variation in surface temperature,T _Wx and surface injectionv _Wx(–1/2). The analysis incorporates the expression connecting porosity and permeability and also the expression connecting porosity and effective thermal diffusivity. The influence of thermal dispersion on the flow and heat transfer characteristics are also analysed in detail. The results of the present analysis document the fact that variable porosity enhances heat transfer rate and the magnitude of velocity near the wall. The governing equations are solved using an implicit finite difference scheme for both the Darcy flow model and Forchheimer flow model, the latter analysis being confined to an isothermal surface and an impermeable vertical plate. The influence of the intertial terms in the Forchheimer model is to decrease the heat transfer and flow rates and the influence of thermal dispersion is to increase the heat transfer rate.

---

Der Effekt des Oberflächenstoffaustausches bei auftriebsinduzierter Strömung in einem variablen porösen Medium, das an eine vertikale, beheizte Platte angrenzt

Zusammenfassung Es wird der Effekt des Oberflächenstoffaustausches in auftriebsinduzierter Strömung in einem variablen porösen Medium, das an eine vertikale, beheizte Platte angrenzt, für große Reynoldszahlen untersucht. Ähnliche Lösungen werden im Rahmen der Grenzschicht-Theorie, durch Variation des Potenzansatzes der Oberflächentemperatur,T _Wx , und der Oberflächengeschwindigkeit,v _Wx(–1/2), erreicht. Die Analyse vereinigt sowohl den Ausdruck, der Porösität und Permeabilität verbindet, als auch den Ausdruck, der Porösität und Wärmeleitfähigkeit miteinander verbindet. Der Einfluß der Temperaturverteilung auf Strömung und Wärmeübergangskennzahlen wird ebenfalls im Detail analysiert. Als Ergebnis der vorliegenden Untersuchung ergibt sich die Tatsache, daß variable Porösität Wärmeübertragungsrate und Betrag der Geschwindigkeit in Wandnähe steigert. Die bestimmenden Gleichungen, sowohl für das Darcysche Strömungsmodell als auch für das Forchheimersche Strömungsmodell, werden mit Hilfe eines implizierten Differenzenschemas gelöst. Die Berechnung wird für die beiden Fälle, isotherme Oberfläche und undurchlässige vertikale Platte, angewandt. Der Einfluß der Terme für die Trägheitskräfte im Forchheimerschen Modell senkt Wärmeübergangs- und Durchgangsrate, wogegen die Wärmeübergangsrate durch den Einfluß der Temperaturverteilung erhöht wird.

---

Nomenclature a constant defined by Eq. (12) - A constant defined by Eq. (12) - B constant defined by Eq. (3) - b _s/_f ratio of thermal conductivities - C constant defined by Eq. (1) - C _P specific heat of the convective fluid - d particle diameter - f dimensionless function defined by Eq. (14) - f _w lateral mass flux parameter - g acceleration due to gravity - k ₀ mean permeability of the mediumk ₀= ₀ ³ d ²/150 (1– ₀)² k ₀=1.75d/(1– ₀) 150 (Inertia parameter) - L length of the source or sink - m mass transfer - n constant defined in Eq. (12) - k (y) permeability of the porous medium - k (y) interial coefficient in the Ergun expression - Gr modified Grashof numberGr=(g k ₀ k ₀ (T _w–))/ ² - R _a Rayleigh number (g k ₀ x T _w–)/ - R _ad modified Rayleigh number (g k ₀ d|T _w–|)/ - N _u Nusselt number - s x/d - Q overall heat transfer rate - T temperature - T _w surface temperature - T ambient fluid temperature - u velocity in vertical direction - v velocity in horizontal direction - x vertical coordinate - y horizontal coordinate Greek symbols ₀ mean thermal diffusivity _f/ C_p - coefficient of thermal expansion - constant defined in Eq. (4) - ratio of particle to bed diameter - _e effective thermal conductivity - f thermal conductivity of fluid - _s thermal conductivity of solid - dimensionless similarity variable in Eq. (13) - value of at the edge of the boundary layer - constant defined in Eq. (1) - _e effective molecular thermal diffusivity - (y) porosity of the medium - ₀ mean porosity of the medium - viscosity of the fluid - ₀ density of the convective fluid - stream function - w condition at the wall - condition at infinity     

Viscoelastic swirling flow with free surface in cylindrical chambers

The flow of a viscoelastic liquid driven by the steadily rotating bottom cover of a cylindrical cup is investigated. The flow field and the shape of the free surface are determined at the lowest significant orders of the regular domain perturbation in terms of the angular velocity of the bottom cap. The meridional field superposed on a primary azimuthal field shows a structure of multiple cells. The velocity field and the shape of the free surface are strongly effected by the cylinder aspect ratio and the elasticity of the liquid. The use of this flow configuration as a free surface rheometer to determine the first two Rivlin-Ericksen constants is shown to be promising.Nomenclature R, ,Z Coordinates in the physical domain D - , , Coordinates in the rest stateD ₀ - r, ,z Dimensionless coordinates in the rest stateD ₀ - Angular velocity - Zero shear viscosity - Surface tension coefficient - Density - Dimensionless surface tension parameter - ₁, ₂ The first two Rivlin-Ericksen constants - Stream function - Dimensionless second order meridional stream function - ^* Dimensionless second normal stress function - 2 Dimensionless sum of the first and second normal stress functions - N ₁,N ₂ The first and second normal stress functions - n Unit normal vector - D Stretching tensor - A _n nth order Rivlin-Ericksen tensor - S Extra-stress - u Velocity field - U Dimensionless second order meridional velocity field - V Dimensionless first order azimuthal velocity field - p Pressure - Modified pressure field - P Dimensionless second order pressure field - J Mean curvature - a Cylinder radius - d Liquid depth at rest - D Dimensionless liquid depth at rest - h Free surface height - H Dimensionless free surface height at the second order     

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 in partitioned enclosed spaces of solar collector

The two-dimensional, steady, laminar natural convection phenomena in partitioned enclosure of solar collector has been studied numerically. Heat conduction along the partition is considered. An iterative finite-difference scheme is employed to solve the governing equations in the flow field. The effects of Rayleigh number, thermal conductivity ratio, partition angle, tilt angle, and aspect ratio on both the local and average heat transfer coefficients of the solar collector have been discussed. The range of Rayleigh number tested was up to 5 × 10⁴, the thermal conductivity ratio of 4.5 and 30, partition angle from 10 deg to 170 deg, tilt angle from 10 deg to 90 deg, and aspect ratio varied between 0.2 and 10. The results indicate that the convective heat transfer is strongly affected with the aspect ratio of the enclosures.

---

Freie Konvektion in unterteilten Kammern von Solarkollektoren

Zusammenfassung Die zweidimensionale, stetige, laminare freie Konvektion in unterteilten Kammern von Solarkollektoren wurde numerisch untersucht. Die Wärmeübertragung entlang dieser Kammern wurde betrachtet. Ein iteratives Finite-Differenzen-Schema wurde angewandt um die Gleichungen, welche das Strömungsfeld beschreiben, zu lösen. Der Einfluß der Rayleigh-Zahl, der thermische Leitfähigkeit, des Kammerwinkels, des Neigungswinkels und der Längenverhältnisse auf die örtlichen und durchschnittlichen Wärmeübertragungskoeffizienten des Solarkollektors wurde diskutiert. Der Bereich der Rayleigh-Zahl erstreckte sich bis zu 5 × 10⁴, das Verhältnis der thermischen Leitfähigkeit betrug 4.5 und 30, der Kammerwinkel lag zwischen 10° und 170°, der Neigungswinkel zwischen 10° und 90° und das Längenverhältnis variierte zwischen 0.2 und 10. Die Ergebnisse beinhalten, daß die konvektive Wärmeübertragung sehr stark durch das Längenverhältnis der Kammern beeinflußt wird.

---

Nomenclature a slope of the partition with respect to the horizontal - A H/L=cell aspect ratio - A _w t/L=wall aspect ratio - g acceleration due to gravity - h local heat transfer coefficient - average heat transfer coefficient - H cell length - k thermal conductivity of fluid within the cell - k _w thermal conductivity of the cell wall - L plate spacing - Nu _f h L/k=local cell Nusselt number - L/k=average cell Nusselt number - overall average Nusselt number - qL/k _w t(T _h–T _c)=average wall Nusselt number - Pr /=Prandtl number - q heat transfer in the cell wall from the hot to cold plate per unit depth - Ra g L ³ T/=Rayleigh number - R _k k _w/k=ratio of wall thermal conductivity to that of the fluid - t thickness of cell wall - T _c cold plate temperature - T _f temperature in cell - T _h hot plate temperature - T _w temperature in cell wall - u, U dimension and dimensionless velocities inx-direction - v, V dimension and dimensionless velocities iny-direction - x distance measured from the bottom of the cell (Fig. 1) - X x/L=normalized distance ofx - y distance measured from hot plate (Fig. 1) - Y y/L=normalized distance ofy - x ₁ distance measured in wall (Fig. 1) - X ₁ x/L=normalized distance ofx ₁ Greek symbols thermal diffusivity of fluid - coefficient of volumetric expansion of fluid - partition angle with respect to the hot plate - _f T _f–T _c/T _h–T _c=dimensionless temperature in cell - _w T _w–T _c/T _h–T _c=dimensionless temperature in cell wall - kinematic viscosity of fluid - enclosure tilt angle from horizontal - dimensional vorticity - L ²/=dimensionless vorticity - dimensionless streamline     

Impacts of local dispersion and first-order decay on solute transport in randomly heterogeneous porous media

Stochastic subsurface transport theories either disregard local dispersion or take it to be constant. We offer an alternative Eulerian-Lagrangian formalism to account for both local dispersion and first-order mass removal (due to radioactive decay or biodegradation). It rests on a decomposition of the velocityv into a field-scale componentv , which is defined on the scale of measurement support, and a zero mean sub-field-scale componentv _s , which fluctuates randomly on scales smaller than. Without loss of generality, we work formally with unconditional statistics ofv _s and conditional statistics ofv . We then require that, within this (or other selected) working framework,v _s andv be mutually uncorrelated. This holds whenever the correlation scale ofv is large in comparison to that ofv _s . The formalism leads to an integro-differential equation for the conditional mean total concentration c which includes two dispersion terms, one field-scale and one sub-field-scale. It also leads to explicit expressions for conditional second moments of concentration cc. We solve the former, and evaluate the latter, for mildly fluctuatingv by means of an analytical-numerical method developed earlier by Zhang and Neuman. We present results in two-dimensional flow fields of unconditional (prior) mean uniformv . These show that the relative effect of local dispersion on first and second moments of concentration dies out locally as the corresponding dispersion tensor tends to zero. The effect also diminishes with time and source size. Our results thus do not support claims in the literature that local dispersion must always be accounted for, no matter how small it is. First-order decay reduces dispersion. This effect increases with time. However, these concentration moments c and cc of total concentrationc, which are associated with the scale below, cannot be used to estimate the field-scale concentrationc directly. To do so, a spatial average over the field measurement scale is needed. Nevertheless, our numerical results show that differences between the ensemble moments ofc and those ofc are negligible, especially for nonpoint sources, because the ensemble moments ofc are already smooth enough.     

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.     

Die Temperaturabhängigkeit der Oberflächenspannung von reinen Kältemitteln vom Tripelpunkt bis zum kritischen Punkt

Zusammenfassung Die Oberflächenspannung von sechs reinen Substanzen — SF₆, CCl₃F, CCl₂F₂, CClF₃, CBrF₃ und CHClF₂ — wurde mit Hilfe einer modifizierten Kapillarmethode gemessen. Die zur Berechnung der Oberflächenspannung erforderlichen Sättigungsdichten und wurden teils aus vorhandenen Zustandsgleichungen, teils aus ebenfalls gemessenen Brechungsindizes bestimmt. Die Temperaturabhängigkeit der Oberflächenspannung läßt sich durch einen erweiterten Ansatz nach van der Waals =_O (T_c-T)(1+...) darstellen, wobei bei einfachen Stoffen ein eingliedriger, bei polaren und assoziierenden Stoffen ein zweigliedriger Ansatz notwendig und ausreichend ist. Für den kritischen Exponenten der Oberflächenspannung wurde ein von der molekularen Substanz weitgehend unabhängiger Wert von =1.284±0.005 gefunden.

---

Temperature dependence of surface tension of pure refrigerants from triple point up to the critical point

The surface tension of six fluids (SF₆, CCl₃F, CCl₂F₂, CClF₃, CBrF₃, CHClF₂) have been measured by means of a modified capillary rise method. The liquid vapor densities, which are needed to calculate the surface tension, have partly been determined by means of refractive indices simultaneously measured in the same apparatus. The temperature dependence of the surface tension is described by an extended van der Waals power law =_O(T_c-T)(1+...). For simple fluids one term and for polar and associating fluids two terms are necessary and sufficient. The critical exponent is found to be 1.284 ± 0.005 and nearly independent of the molecular structure.

---

Formelzeichen a² Laplace-Koeffizient - a Parameter - B_O, B_on Koeffizient der Koexistenzkurve - g Erdbeschleunigung - H Höhe, kapillare Steighöhe - LL Lorentz-Lorenz-Funktion oder Refraktionskonstante - M molare Masse - M Zahl der Meßwerte - N Zahl der unbekannten Parameter - n Brechungsindex - p Druck - R,r Radius - s Entropie - SD Standardabweichung - T, t Temperatur - u innere Energie Griechische Formelzeichen Exponent des Laplace-Koeffizienten - Exponent der Koexistenzkurve - 2. Exponent der Oberflächenspannung - Wellenlänge des Lichts - Exponent der Oberflächenspannung - _D Dipolmoment - , Dichte der Flüssigkeit bzw. des Dampfes - Oberflächenspannung - reduzierte Temperatur (1-T/T_c) - ² gewichtete Varianz Indizes c kritischer Zustand - D Differenz - m Mittelwert - T Isotherme - t Zustand am Tripelpunkt - S Zustand am Schmelzpunkt - bezogen auf Oberfläche     

Rotationsviskosimetrie mit kugelförmigen Spindeln

Zusammenfassung Die Rheodynamik der stationären viskometrischen Drehströmung um eine rotierende Kugel wird mit Methoden der Variationsrechnung untersucht. Neben iterativen numerischen Lösungsmethoden, die zu exakten Resultaten führen, wird auch eine approximative Ein-Gradienten-Lösung konstruiert, die durch Quadraturen dargestellt wird. Ausgehend von dieser analytischen Approximation werden einfache Methoden zur Auswertung von Experimentaldaten vorgeschlagen, die mit Hilfe von Eintauch-Rotationsviskosimetern mit kugelförmigen Meßspindeln gewonnen wurden.

---

Summary The rotational viscometric flow around a rotating sphere has been studied by variational methods. The exact numerical, as well as an approximate analytical solutions are given. Employing the analytical approximation, a simple method of evaluating viscometric data from immersional (portable) viscometers with a rotating sphere is proposed.

---

A Achsenschnitt durch den Bereich der Strömung - B - b, c anpaßbare empirische Konstanten - C Kalibrierungsoperator - D Schergeschwindigkeit der viskosimetrischen Strömung - D _ij Komponenten des Deformationsgeschwindigkeitstensors - D _I, _I Stoffkonstanten der VF des Ellis-Modells - g metrischer Koeffizient - H() Funktional der Ein-Gradienten-Approximation, Gl. [27] - J[] energetisches Potential - J _a[] Ein-Gradienten-Approximation fürJ - K Konsistenzkoeffizient, Parameter der VF des Potenzmodells - m Parameter des Ellis-Modells - M Drehmoment - n Parameter des Potenzmodells - n, n Differentialindices der VF, Gl. [20c, d] - n*,n** Differentialindices der RC, Gl. [9], [13] - r, , z polare Zylinderkoordinaten - R Spindelhalbmesser - rheometrischer Operator - S Spindeloberfläche - U(D) energetische Funktion nachBird, Gl. [20e] - v _i physikalische Komponenten der Geschwindigkeit - Z() transformierte VF, Gl. [20f] - (n) durch Gl. [35] definierte Funktion - k Verhältnis der Radien von Spindel und Wand - ( durch Gl. [43] definierte Funktion - natürliche (Radial-)Koordinate - Schubspannung der viskosimetrischen Strömung - _ij Komponenten des Spannungstensors - _S() Spannungsprofil an der Spindeloberfläche - _M Maximalspannung an der Spindeloberfläche - mittlere Spannung an der Spindeloberfläche, Gln. [3], [22] - natürliche (Meridional-) Koordinate - Winkelgeschwindigkeit in der Flüssigkeit - Winkelgeschwindigkeit der Spindelrotation - ( rheometrische Charakteristik Mit 4 Abbildungen und 3 Tabellen     

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

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      

Dielectric properties of the composite system polyurethane — sodium chloride

Summary The dielectric properties of the composite system polyurethane-sodium chloride have been measured at frequencies between 10^–4 Hz and 3 · 10⁵ Hz in the temperature range from –150 °C up to +90 dgC. Three dielectric loss mechanisms have been found; they are indicated by ₁, ₂ and. The activation energy of the ₁-transition is 35 kcal/mole, that of the-transition 6.7 kcal/mole. The ₂-loss peak was only observed at frequencies of 10³ Hz and higher, forming one broad peak with the ₁-loss peak at lower frequencies. In the composite materials, the- and ₂-loss peaks measured at fixed frequencies were found at the same temperature. The ₂-loss peak, however, was shifted to a lower temperature, due to the sodium chloride filler. Comparison of experimental data of and tan with theoretical predictions concerning the dielectric properties of composite systems showed only partial agreement. The difference mainly consisted in. the temperature shift in the tan-peak of the ₁-transition.

---

Zusammenfassung Die dielektrischen Eigenschaften des Verbundssystems Kochsalz-Polyurethankautschuk wurden im Frequenzgebiet zwischen 10^–4 Hz und 3.10⁵ Hz und im Temperaturbereich von –150 °C bis +90 °C gemessen. Es wurden drei dielektrische Verlustmaxima gefunden, die mit ₁, ₂ und angedeutet werden. Die Aktivierungsenergie des ₁-Überganges beträgt 35 kcal/Mol, die des-Überganges 6.7 kcal/Mol. Das ₂-Maximum konnte nur bei Frequenzen höher als 10³Hz vom ₁-Maximum gesondert erfaßt werden. Die Lage der ₂- und-Maxima war vom Füllgrad unabhängig. Das ₁-Maximum verschiebt sich mit steigendem Füllgrad zu niedrigeren Temperaturen. Die gemessenen Werte von und tan stimmen nur teilweise mit den Aussagen einer Theorie der dielektrischen Eigenschaften von Mischkörpern überein.

     

On the localness of the spectral energy transfer in turbulence

By utilizing available experimental data for net energy transfer spectra for homogeneous turbulence, contributions P(, ) to the energy transfer at a wavenumber from various other wavenumbers are calculated. This is done by fitting a truncated power-exponential series in and to the experimental data for the net energy transfer T(), and using known properties of P(, ). Although the contributions P(, ) obtained by using this procedure are not unique, the results obtained by using various assumptions do not differ significantly. It seems clear from the results that for a region where the energy entering a wavenumber band dominates that leaving, much of the energy entering the band comes from wavenumbers which are about an order of magnitude smaller. That is, the energy transfer is rather nonlocal. This result is not significantly dependent on Reynolds number (for turbulence Reynolds numbers based on microscale from 3 to 800). For lower wavenumbers, where more energy leaves than enters a wavenumber band, the energy transfer into the band is more local, but much of the energy then leaves at distant wavenumbers.     

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      

On forced vibration of anisotropic cylinders

Summary The dynamic response of a circular cylinder with thick walls of transverse curvilinear isotropy subjected to a uniformly distributed pressure varying periodically with time is analyzed by means of the Laplace transformation, and the exact solution is obtained in closed form. The previously obtained solutions for forced vibrations with isotropy, and free vibrations with transverse curvilinear isotropy are included as special cases of the general results reported here.Nomenclature t time - r, , z cylindrical coordinates - _ii components of normal strain - _ii components of normal stress - u radial displacement - c _ij elastic constant - mass density - c ² c ₁₁/ - ² c ₂₂/c ₁₁ - a, b inner, outer radius of the cylinder - , A, B constants - forced angular frequency - function defined by (9) - p, real, complex variables - constant defined by (14) - real number - , Lamé elastic constants - J (x) Bessel function of first kind - Y (x) Bessel function of second kind - I (x) modified Bessel function of first kind - K (x) modified Bessel function of second kind     

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.     

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

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     

Nonstationary vibration of a rotating shaft with nonlinear spring characteristics during acceleration through a critical speed (A critical speed of a summed-and-differential harmonic oscillation)

Nonstationary vibration of a flexible rotating shaft with nonlinear spring characteristics during acceleration through a critical speed of a summed-and-differential harmonic oscillation was investigated. In numerical simulations, we investigated the influence of the angular acceleration , the initial angular position of the unbalance _n and the initial rotating speed on the maximum amplitude. We also performed experiments with various angular accelerations. The following results were obtained: (1) the maximum amplitude depends not only on but also on _n and : (2) when the initial angular position _n changes. the maximum amplitude varies between two values. The upper and lower bounds of the maximum amplitude do not change monotonously for the angular acceleration: (3) In order to always pass the critical speed with finite amplitude during acceleration. the value of must exceed a certain critical value.Nomenclature O-xyz rectangular coordinate system - , ₁, ₁ inclination angle of rotor and its projections to thexy- andyz-planes - I _r polar moment of inertia of rotor - I diametral moment of inertia of rotor - i _r ratio ofI _r toI - dynamic unbalance of rotor - directional angle of fromx-axis - c damping coefficient - spring constant of shaft - N _nt ,N _nt nonlinear terms in restoring forees in ₁ and ₁ directions - ₄ representative angle - a small quantity - V. V _u .V _N potential energy and its components corresponding to linear and nonlinear terms in the restoring forees - directional angle - _n coefficients of asymmetrical nonlinear terms - _n coefficients of symmetrical nonlinear terms - coefficients of asymmetrical nonlinear terms experessed in polar coordinates - coefficients of symmetrical nonlinear terms expressed in polar coordinates - rotating speed of shaft - t time - _n initial angular position of att=0 - p natural frequency - p ₁.p _t natural frequencies of forward and backward precessions - , ₁, ₁ total phases of harmonic, forward precession and backward precession components in summed-and-differential harmonic oscillation - , ₁, ₁ phases of harmonic, forward precession and backward precession components in summed-and-differential harmonic oscillation - P, R _t ,R _b amplitudes of harmonic, forward precession and backward precession components in summed-and-differential harmonic oscillation - difference between phases ( = _f – _u) - acceleration of rotor - initial rotating speed - t _t ,r _b amplitudes of nonstationary oscillation during acceleration - (r _t )_max, (r _b )_max maximum amplitudes of nonstationary oscillation during acceleration - (r ₁ ¹ )_max, (r _b ¹ )_max maximum value of angular acceleration of non-passable case - ₀ critical value over which the rotor can always pass the critical speed - p ₁,p ₂,p ₃,p ₄ natural frequencies of experimental apparatus     

Optical co-ordinates in general relativity 2 — the problem of distance

Summary Let denote the congruence of null geodesics associated with a given optical observer inV ₄. We prove that determines a unique collection of vector fieldsM₍₎ ( =1, 2, 3) and ₍₀₎ overV ₄, satisfying a weak version of Killing's conditions.This allows a natural interpretation of these fields as the infinitesimal generators of spatial rotations and temporal translation relative to the given observer. We prove also that the definition of the fieldsM₍₎ and ₍₀₎ is mathematically equivalent to the choice of a distinguished affine parameter f along the curves of, playing the role of a retarded distance from the observer.The relation between f and other possible definitions of distance is discussed.

---

Sommario Sia la congruenza di geodetiche nulle associata ad un osservatore ottico assegnato nello spazio-tempoV ₄. Dimostriamo che determina un'unica collezione di campi vettorialiM₍₎ ( =1, 2, 3) e ₍₀₎ inV ₄ che soddisfano una versione in forma debole delle equazioni di Killing. Ciò suggerisce una naturale interpretazione di questi campi come generatori infinitesimi di rotazioni spaziali e traslazioni temporali relative all'osservatore assegnato. Dimostriamo anche che la definizione dei campiM₍₎, ₍₀₎ è matematicamente equivalente alla scelta di un parametro affine privilegiato f lungo le curve di, che gioca il ruolo di distanza ritardata dall'osservatore. Successivamente si esaminano i legami tra f ed altre possibili definizioni di distanza in grande.

---

---

Work performed in the sphere of activity of: Gruppo Nazionale per la Fisica Matematica del CNR.