On the closure problem for Darcy's law

In a previous derivation of Darcy's law, the closure problem was presented in terms of an integro-differential equation for a second-order tensor. In this paper, we show that the closure problem can be transformed to a set of Stokes-like equations and we compare solutions of these equations with experimental data. The computational advantages of the transformed closure problem are considerable.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 _e area of entrances and exits for the-phase contained within the averaging volume, m² - B second-order tensor used to respresent the velocity deviation - b vector used to represent the pressure deviation, m^–1 - C second-order tensor related to the permeability tensor, m^–2 - D second-order tensor used to represent the velocity deviation, m² - d vector used to represent the pressure deviation, m - g gravity vector, m/s² - I unit tensor - K C ^–1,–D, Darcy's law permeability tensor, m² - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the-phase, m - l _i i=1, 2, 3, lattice vectors, m - 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 - p pressure in the-phase, N/m ² - p intrinsic phase average pressure, N/m² - p –p , spatial deviation of the pressure in the-phase, N/m² - r position vector locating points in the-phase, m - r ₀ radius of the averaging volume, m - t time, s - v velocity vector in the-phase, m/s - v intrinsic phase average velocity in the-phase, m/s - v phase average or Darcy velocity in the \-phase, m/s - v –v , spatial deviation of the velocity in the-phase m/s - V averaging volume, m³ - V volume of the-phase contained in 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²     

Existence and Uniqueness of Time-Periodic Physically Reasonable Navier-Stokes Flow Past a Body

Let be a three-dimensional exterior domain of class C^2,, 0<<1. Assume that a Navier-Stokes liquid is moving in under the action of a body force F that is time-periodic of period T, and that the velocity of the liquid is zero at spatial infinity. In this paper we show that, if F satisfies suitable conditions, and its norm, in appropriate function spaces, is sufficiently small, there is at least one time-periodic strong solution. Furthermore, the velocity field v of such a solution decays to zero for large |x| as |x|^–1 and its spatial gradient decays as |x|^–2, both uniformly in time. In addition, the pressure p decays like |x|^–2 and its gradient like |x|^–3, for almost all t[0,T]. In the special case where F is time-independent, these solutions are also time-independent and coincide with Finns physically reasonable solutions [4]. Moreover, we show that our time-periodic solutions are unique in a very large class, namely, the class of time-periodic weak solutions satisfying the energy inequality and with corresponding pressure fields verifying mild summability conditions in ×[0,T].     

On Perturbative Expansions to the Stochastic Flow Problem

When analyzing stochastic steady flow, the hydraulic conductivity naturally appears logarithmically. Often the log conductivity is represented as the sum of an average plus a stochastic fluctuation. To make the problem tractable, the log conductivity fluctuation, f, about the mean log conductivity, lnK _G, is assumed to have finite variance, _f ². Historically, perturbation schemes have involved the assumption that _f ²<1. Here it is shown that _f may not be the most judicious choice of perturbation parameters for steady flow. Instead, we posit that the variance of the gradient of the conductivity fluctuation, _f ², is more appropriate hoice. By solving the problem withthis parameter and studying the solution, this conjecture can be refined and an even more appropriate perturbation parameter, , defined. Since the processes f and f can often be considered independent, further assumptions on f are necessary. In particular, when the two point correlation function for the conductivity is assumed to be exponential or Gaussian, it is possible to estimate the magnitude of f in terms of _f and various length scales. The ratio of the integral scale in the main direction of flow ( _x ) to the total domain length (L^*), _x ²=_x/L^*, plays an important role in the convergence of the perturbation scheme. For _x smaller than a critical value _c, _x < _c, the scheme's perturbation parameter is =_f/_x for one- dimensional flow, and =_f/_x ² for two-dimensional flow with mean flow in the x direction. For _x > _c, the parameter =_f/_x ³ may be thought as the perturbation parameter for two-dimensional flow. The shape of the log conductivity fluctuation two point correlation function, and boundary conditions influence the convergence of the perturbation scheme.     

Existence of L 1-Connections Between Equilibria of a Semilinear Parabolic Equation

We give a necessary and sufficient condition for the existence of L ¹-connections between equilibria of a semilinear parabolic equation. By an L ¹-connection from an equilibrium ^– to an equilibrium ⁺ we mean a function u(, t) which is a classical solution on the interval (–, T) for some T and blows up at t = T but continues to exist in the space L ¹ for t [T, ) and satisfies u(, t) ^± (in a suitable sense) as t ±. The main tool in our analysis is the zero number.     

Testing viscometric predictions of the internal viscosity model,using dilute viscous theta solutions

A stress-symmetrized internal viscosity (I.V.) model for flexible polymer chains, proposed by Bazua and Williams, is scrutinized for its theoretical predictions of complex viscosity ^* () = – i and non-Newtonian viscosity (), where is frequency and is shear stress. Parameters varied are the number of submolecules,N (i.e., molecular weightM = NM _s ); the hydrodynamic interaction,h ^*; and/f, where andf are the I.V. and friction coefficients of the submolecule. Detailed examination is made of the eigenvalues _p (N, h ^*) and how they can be estimated by various approximations, and property predictions are made for these approximations.Comparisons are made with data from our preceding companion paper, representing intrinsic properties [], [], [] in very viscous theta solutions, so that theoretical foundations of the model are fulfilled. It is found that [ ()] data can be predicted well, but that [ ()] data cannot be matched at high. The latter deficiency is attributed in part to unrealistic predictions of coil deformation in shear.     

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

Experimental study of mutual interference between two spheres placed on a plane boundary

This paper describes an experimental study of the mutual interference between two spheres placed on a plane boundary. The experiment was carried out in an N. P. L. type wind-tunnel having a working section of 500×500×2000 mm³ in size at a Reynolds number of 4.74×10⁴. The surface-pressure distributions of two spheres were measured for the various relative positions of two spheres and the drag, side-force, and lift coefficients were determined from surface-pressure distributions. The separation of the flow and the formation of vortices were observed by the method of visualization. The distributions of velocities, and turbulent intensities of the flow past two spheres were measured. The experimental results for two spheres were compared with those of a single sphere.List of symbols C _D drag coefficient - C _L lift coefficient - C _p surface-pressure coefficient of sphere=(P-P )(qU ² ) - C _s coefficients of side force - D diameter of sphere [mm] - P static pressure [Pa] - P static pressure in free stream [Pa] - Re Reynolds number= DU/v - S spacing between the centers of two adjoining spheres in plane view [mm] - U time-mean velocity in X-direction [m/s] - [m/s] free stream velocity [m/s] - u, v, w X, Y and Z-components of velocity fluctuation [m/s] - X, Y, Z coordinate axes with origin at the bottom center of test sphere, X, Y, Z axis being taken in the streamwise, lateral and vertical directions respectively [mm] (Fig. 1) - latitude angle [°] - longitude angle [°] - angle between the line connected with the centers of two spheres and wind direction [°] (Fig. 2) - kinematic viscosity of air [m²/s] - density of air [N/m³] This paper was presented at the 10th Symposium on Turbulence, University of Missouri-Rolla, Sept. 22–24, 1986     

Symplectic approach to the theory of quantized fields. II

We prove that the set D of vector fields on the configuration space B of a field whose 1-parameter groups locally associated are groups of fibre-preserving transformations of B that leave invariant that field in the sense of variational theory, is a Lie algebra with respect to ordinary addition, multiplication by real numbers and Lie brackets. We see that this Lie algebra structure can be carried over to the corresponding set of Noether invariants, which then becomes a Lie algebra in a natural way.Further, we define the n-form of Poincaré-Cartan of a field, and we use it to generalize the Lie algebras D and in a reasonable way. The algebras D and are subalgebras of the new Lie algebras D and introduced. A main result in this connection is the following: the differential d of the n-form of Poincaré-Cartan is –(d+f), where (, d+f) are the field equations on the vertical bundle B.The symplectic manifold of solutions associated with a field is introduced in a formal way and the former Lie algebras D, , D, are interpreted on this manifold. In imitation of the case of analytical dynamics, the main results in this direction are: a) Every vector field of the Lie algebra D defines, in a canonical way, a vector field on the manifold of solutions such that its polar 1-form with respect to the symplectic metric ₂ is the differential of its corresponding Noether invariant, and b) the Lie bracket [, ] of two Noether invariants , is the Noether invariant given by ₂(D, D), where D, D are the vector fields on the manifold of solutions defined, in the sense a), by two infinitesimal generators of , , respectively. This will allow us to regard the Lie algebra as the analogous object in field theory to the Poisson algebra of analytic dynamics.We apply the general formalism to the relativistic theory of non-linear scalar fields, and we compare our results with the formalism developed by I. Segal for this case.     

Material functions generated by the complex viscosity

Based on the complex viscosity model various steady-state and transient material functions have been completed. The model is investigated in terms of a corotational frame reference. Also, BKZ-type integral constitutive equations have been studied. Some relations between material functions have been derived. C ^–1 Finger tensor - F[], (F ^–1[]) Fourier (inverse) transform - rate of deformation tensor in corotating frame - h(I, II) Wagner's damping function - J (x) Bessel function - m parameter inh (I, II) - m(s) memory function - m _k, n_k integers (powers in complex viscosity model) - P principal value of the integral - parameter in the complex viscosity model - rate of deformation tensor - shear rates - [], [] incomplete gamma function - (a) gamma function - steady-shear viscosity - ^* complex viscosity - , real and imaginary parts of ^* - ₀ zero shear viscosity - ⁺, ₁ ⁺ stress growth functions - ^–, ₁ ^- stress relaxation functions - (s) relaxation modulus - ₁(s) primary normal-stress coefficient - ø(a, b; z) degenerate hypergeometric function - ₁, ₂ time constants (parameters of ^*) - frequency - extra stress tensor     

Interaction of an oscillating vortex with a turbulent boundary layer

An oscillating vortex embedded within a turbulent boundary layer was generated experimentally by forcing a periodic lateral translation of a half-delta wing vortex generator. The objective of the experiment was to investigate the possibility that a natural oscillation, or meander, might be responsible for flattened vortex cores observed in previous work, which could also have contaminated previous turbulence measurements. The effect of this forced oscillation was characterized by comparison of measurements of the mean velocities and Reynolds stresses at two streamwise stations, for cases with and without forcing. The Reynolds stresses, especially w, were affected significantly by the forced oscillation, mainly through contributions from the individual production terms, provided the vortex was not too diffuse.List of Symbols a amplitude of forced vortex motion - f frequency of forced vortex generator motion - l vortex generator root chord - L flow length scale - R _Y , R _Z vortex core radial dimensions in vertical and spanwise directions, respectively - R_r vortex circulation Reynolds number R = / - u, v, w instantaneous velocity components in X, Y, Z directions - U, V, W mean velocities; shorthand notation for u, , w - X, Y, Z right-hand Cartesian streamwise, vertical, and spanwise coordinate directions - boundary-layer thickness - overall circulation - air kinematic viscosity - _x streamwise vorticity, _X = W/Y–V/d+t6Z - ( )₀ reference value (measured at X = 10 cm) - ( )c refers to vortex center - ( ) _max maximum value for a particular crossflow plane - ( ) (overbar) time average - ( ) (prime) fluctuating component, e.g., u=U+u     

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.     

Flow over a body with development of a steady separation zone as re → ∞

This paper discusses formulation of the total problem of flow of an incompressible liquid over a body, with formation of a closed stationary separation zone as Re . The scheme used is based on the method of matched asymptotic expansions [1]. Following [1], it is postulated that the separated zone is developed (i.e., it is not infinitely fragmented and does not vanish as Re ), and the flow inside it has a definite degree of regularity with respect to Re. With these hypotheses we can use the Prandtl-Batchelor theorem [2], which states that, in the limit as Re , a region of circulating flow becomes vortex flow of an inviscid liquid with constant vorticity . Therefore, a basis for constructing matched asymptotic expansions is the vortex-potential problem (the problem of determining a stream function , satisfying the equation = 0 in the region of translational motion and the equation = in a certain region, unknowna priori, of circulating motion). In the general case the solution of the vortex-potential problem depends on two parameters: the total pressure p_o and the vorticity in the separated zone. These parameters appear in the condition for matching the solutions of the first and second boundary-layer approximations (at the boundary of the separated zone for the end Re values) with the corresponding solutions for the inviscid flow. It is shown in the present paper that the conditions for matching the cyclic boundary layer with the external translational flow are the same additional relations which allow us to close the total problem. Thus, in using the method of matched asymptotic expansions to solve the problem of flow over a body with closed stationary separation zones one must simultaneously consider no less than two approximations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 28–37, March–April, 1978.The authors thank G. Yu. Stepanov for discussion of the paper and valuable comments.     

Einige Beziehungen zwischen Transportkoeffizienten und thermodynamischen Zustandsgrö\en

Zusammenfassung Eine früher mitgeteilte Beziehung [1] zwischen der ViskositÄt und dem isenthalpen Joule-Thomson-Koeffizienten _h wird für kleine Drücke theoretisch begründet und an sieben Stoffen nachgeprüft. Die WärmeleitfÄhigkeit wird als Funktion von c_v für drei Stoffe dargestellt.

---

Some relations between transport coefficients and thermodynamical properties

A formerly given relation [1] between viscosity and isenthalpic Joule-Thomson-coefficient _h is proved theoretically for small pressures and checked with seven substances. The heat conductivity is presented as a function ofc_v for three substances.

---

Bezeichnungen B ^* dimensionsloser zweiter Virialkoeffizient - B ₁ ^* Ableitung vonB ^* nachT ^*.B ₁ ^*=T ^* dB ^*/dT ^* - c _v isochore spezifische WärmekapazitÄt - C _p ^o isobare molare WärmekapazitÄt fürp 0 - h spezifische Enthalpie - k Boltzmann-Konstante.k=R/N _A - M molare Masse - N _A Avogadro-Konstante - p Druck - R molare Gaskonstante - R _i spezifische Gaskonstante des Stoffesi - it Celsius-Temperatur - T Kelvin-Temperatur - T ^* dimensionslose Temperatur.T ^*=kT/ - _h isenthalper Joule-Thomson-Koeffizient. _h=(T/p)_h - , Konstanten der Potentialfunktion - ViskositÄt - WärmeleitfÄhigkeit - (^2,2)^* dimensionsloses Sto\integral     

Concentration dependence of the critical viscoelastic properties of gelatin at the gel point

Gelatin gel properties have been studied through the evolution of the storage [G()] and the loss [G()] moduli during gelation or melting near the gel point at several concentrations. The linear viscoelastic properties at the percolation threshold follow a power-law G()G() and correspond to the behavior described by a rheological constitutive equation known as the Gel Equation. The critical point is characterized by the relation: tan = G/G = cst = tan ( · /2) and it may be precisely located using the variations of tan versus the gelation or melting parameter (time or temperature) at several frequencies. The effect of concentration and of time-temperature gel history on its variations has been studied. On gelation, critical temperatures at each concentration were extrapolated to infinite gel times. On melting, critical temperatures were determined by heating step by step after a controlled period of aging. Phase diagrams [T = f(C)] were obtained for gelation and melting and the corresponding enthalpies were calculated using the Ferry-Eldridge relation. A detailed study of the variations of A with concentration and with gel history was carried out. The values of which were generally in the 0.60–0.72 range but could be as low as 0.20–0.30 in some experimental conditions, were compared with published and theoretical values.     

Normal forms for random diffeomorphisms

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

Supersonic air flow past a sphere with account for vibrational relaxation

A numerical solution is obtained for the problem of air flow past a sphere under conditions when nonequilibrium excitation of the vibrational degrees of freedom of the molecular components takes place in the shock layer. The problem is solved using the method of [1]. In calculating the relaxation rates account was taken of two processes: 1) transition of the molecular translational energy into vibrational energy during collision; 2) exchange of vibrational energy between the air components. Expressions for the relaxation rates were computed in [2]. The solution indicates that in the state far from equilibrium a relaxation layer is formed near the sphere surface. A comparison is made of the calculated values of the shock standoff with the experimental data of [3].Notation uV_max, vV_max velocity components normal and tangential to the sphere surface - V_max maximal velocity - P V _max ² pressure - density - TT temperature - e_viRT vibrational energy of the i-th component per mole (i=–O₂, N₂) - =rb–1 shock wave shape - a _f the frozen speed of sound - HRT/m gas total enthalpy     

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

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

Messung des binären Diffusionskoeffizienten in einem Entmischungssystem mit Hilfe der Photonen-Korrelationsspektroskopie

Zusammenfassung Die dynamische Lichtstreuung in Form der Photonen-Korrelationsspektroskopie wird in einem Entmischungssystem beispielhaft zur Messung des binären Diffusionskoeffizienten eingesetzt. Mit einem Versuchsaufbau wird nahe dem kritischen Entmischungspunkt mit der Homodyn-Methode und weit davon entfernt mit der Heterodyn-Methode gearbeitet. Ein Verfahren ermöglicht die Korrektur der Störeinflüsse des Heterodyn-Anteils bei der Homodyn-Messung. Fürn-Hexan/Nitrobenzol wird der Diffusionskoeffizient für vier unterschiedliche Konzentrationswerte als Funktion der Temperatur ausgemessen. Bei der kritischen Konzentration zeigt er bei gleicher Temperatur den kleinsten Wert und läßt sich nahe dem kritischen Entmischungspunkt durch einen einfachen Potenzansatz mit der reduzierten Temperaturdifferenz =T (– T_c)/T_c ausdrücken. Die statistische Genauigkeit ist besser als 1 %. Die Übereinstimmung mit Literaturwerten ist gut.

---

Measuring of the binary diffusion coefficient in a separation system with photon-correlation-spectroscopy

The dynamic light scattering in form of photon-correlation spectroscopy is examplary used in a separation system for measuring the binary diffusion coefficient. In a test setup the homodyntechnique is used near the critical separation point and in distance the heterodyn-technique is used. A special method allows the correction of the disturbing influences of the heterodyn-part using the homodyn-measuring. Forn-hexane/nitrobenzene the diffusion coefficient is measured for four different concentration values as a function of temperature. At the critical concentration with constant temperature the coefficient shows the minimum value and it is expressed near the critical separation point with an elementary exponential equation with the reduced temperature difference =T (– T_c)/T_c. The statistical precision is better than 1%. The conformity with the literature is well.

---

Abbreviations

Formelzeichen A Konstante - B Konstante - b Konstante - b ₁ Konstante - b₂ Konstante - C Konstante - c Konzentration - c _c kritische Konzentration - D Konstante - D ₁₂ binärer Diffusionskoeffizient - E ₀ elektrisches Feld des Laserlichts - E _s elektrisches Feld des Streulichts - G () Korrelationsfunktion - I ₀ Intensität des Referenzlichts - _s mittlere Intensität des Streulichts - k ₀ Wellenvektor des Laserlichts - k _s Wellenvektor des Streulichts - n Brechungsindex der zu untersuchenden Flüssigkeit - p Druck - q Streuvektor - R Ortsvektor - r Ortsvektor - T Temperatur, Zeit - t Zeit - x Molenbruch Griechische Buchstaben reduzierte Temperaturdifferenz; Dielektrizität - Frequenz des Laserlichts - statistische Schwankungen der Dielektrizitätskonstante - ₀ Wellenlänge des Laserlichts - Streuwinkel - Zeit - _c Zeitkonstante - kritischer Exponent     

Comparison of a new internal viscosity model with other constrained-connector theories of dilute polymer solution rheology

Complex viscosity ^* = -i predictions of the Dasbach-Manke-Williams (DMW) internal viscosity (IV) model for dilute polymer solutions, which employs a mathematically rigorous formulation of the IV forces, are examined in the limit of infinite IV over the full range of frequency number of submolecules N, and hydrodynamic interaction h ^*. Although the DMW model employs linear entropic spring forces, infinite IV makes the submolecules rigid by suppressing spring deformations, thereby emulating the dynamics of a freely jointed chain of rigid links. The DMW () and () predictions are in close agreement with results for true freely jointed chain models obtained by Hassager (1974) and Fixman and Kovac (1974 a, b) with far more complicated formalisms. The infinite-frequency dynamic viscosity predicted by the DMW infinite-IV model is also found to be in remarkable agreement with the calculations of Doi et al. (1975). In contrast to the other freely jointed chain models cited above, however, the DMW model yields a simple closed-form solution for complex viscosity expressed in terms of Rouse-Zimm relaxation times.     

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