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²     

Hyperbolic phenomena in a strongly degenerate parabolic equation

We consider the equation u _t =((u) (u _x )) _x , where >0 and where is a strictly increasing function with lim _s = <. We solve the associated Cauchy problem for an increasing initial function, and discuss to what extent the solution behaves qualitatively like solutions of the first-order conservation law u _t = ((u)) _x . Equations of this type arise, for example, in the theory of phase transitions where the corresponding free-energy functional has a linear growth rate with respect to the gradient.     

Instability of a plane axisymmetric flow with a shear discontinuity in the velocity profile

It is proposed to investigate the stability of a plane axisymmetric flow with an angular velocity profile (r) such that the angular velocity is constant when r < r_O – L and r > r_O + L but varies monotonically from ₁ to ₂ near the point r_O, the thickness of the transition zone being small L r_O, whereas the change in velocity is not small ¦₂ – ₁¦ ₂, ₁. Obviously, as L O short-wave disturbances with respect to the azimuthal coordinate (k=m/r_O 1/r_O) will be unstable with a growth rate-close to the Kelvin—Helmholtz growth rate. In the case L=O (i.e., for a profile with a shear-discontinuity) we find the instability growth rate _O and show that where the thickness of the discontinuity L is finite (but small) the growth rate does not differ from _O up to terms proportional to kL 1 and 1/m 1. Using this example it is possible to investigate the effect of rotation on the flow stability. It is important to note that stabilization (or destabilization) of the flow in question by rotation occurs only for three-dimensional or axisymmetric perturbations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 111–114, January–February, 1985.     

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      

Analysis of slip and temperature jump coefficients in a binary gas mixture

The possibility of simplifying the formulas obtained by the Maxwell-Loyalka method for the velocity _u, temperature _T and diffusion _d slip coefficients and the temperature jump coefficient in a binary gas mixture with frozen internal degrees of freedom of the molecules is considered. Special attention is paid to gases not having sharply different physicochemical properties. The formulas are written in a form convenient for use without linearization in the thermal diffusion coefficient. They are systematically analyzed for mixtures of inert gases, N₂, O₂, CO₂, and H₂ at temperatures extending from room temperature to 2500°K. It is shown that for the molecular weight ratios m^* = m₂/m₁ considered the expressions for _u and can be radically simplified. With an error acceptable for practical purposes (up to 10%) it is possible to employ expressions of the same structural form as for a single-component gas: for _u if 1 m^* 6, and for if 1 m^* 3. When 1 m^* 2 the expression for _T can be simplified with a maximum error of 5%. Within the limits of accuracy of the method the expression for _t can be linearized in the thermal diffusion coefficient. An approximate expression convenient for practical calculations is proposed for _d Finally, the , _u, and _T for a single-component polyatomic gas with easy excitation of the internal degrees of freedom of the molecules are similarly analyzed; it is shown that these expressions can be considerably simplified.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 152–159, November–December, 1990.     

Deformation to breaking of deep water gravity waves

The results of laboratory observations of the deformation of deep water gravity waves leading to wave breaking are reported. The specially developed visualization technique which was used is described. A preliminary analysis of the results has led to similar conclusions than recently developed theories. As a main fact, the observed wave breaking appears as the result of, first, a modulational instability which causes the local wave steepness to approach a maximum and, second, a rapidly growing instability leading directly to the breaking.List of symbols L total wave length - H total wave height - crest elevation above still water level - trough depression below still water level - wave steepness =H/L - crest steepness =/L - trough steepness =/L - F ₁ forward horizontal length from zero-upcross point (A) to wave crest - F ₂ backward horizontal length from wave crest to zero-downcross point (B) - crest front steepness =/F ₁ - crest rear steepness =/F ₂ - vertical asymmetry factor=F ₂/F ₁ (describing the wave asymmetry with respect to a vertical axis through the wave crest) - µ horizontal asymmetry factor=/H (describing the wave asymmetry with respect to a horizontal axis: SWL) - T ₀ wavemaker period - L ₀ theoretical wave length of a small amplitude sinusoïdal wave generated at T _inf0 ^sup–1 frequency - ₀ average wave height     

An application of nonsymmetric Lax-Milgram lemma to nonassociated plasticity

Based on a general assumption for plastic potential and yield surface, some properties of the nonassociated plasticity are studied, and the existence and uniqueness of the distribution of incremental stress and displacement for work-hardening materials are proved by using nonsymmetric Lax-Milgram lemma, when the work-hardening parameter A>F/Q/–F/, Q/.     

A Strongly Degenerate Quasilinear Equation: the Parabolic Case

We prove the existence and uniqueness of entropy solutions of the Neumann problem for the quasilinear parabolic equation u_t=÷ a(u, Du), where a(z,)=f(z,), and f is a convex function of with linear growth as ||||, satisfying other additional assumptions. In particular, this class includes the case where f(z,)=(z)(), >0, and is a convex function with linear growth as ||||.     

Über die Übertragungseinheiten und ein „Vierfelderdiagramm” zur Berechnung von Wärmeübertragern

Zusammenfassung Es wird dargelegt, wie man nach Einführung des bekannten Begriffes Übertragungseinheit und mit Hilfe eines Vierfelderdiagramms zu einer relativ einfachen und durchsichtigen Berechnung von Gleich- und Gegenstrom-Wärmeübertragern gelangt.

---

On transfer units and a four-quadrant-diagram for the calculation of heat exchangers

It is shown that the introduction of the known concept of a transfer unit and the use of a four-quadrant diagram leads to a relatively simple and clear computation method for co-current and counter-current heat exchangers.

     

Nonlinear analysis of the forced response of a beam with three mode interaction

An analysis is presented for the primary resonance of a clamped-hinged beam, which occurs when the frequency of excitation is near one of the natural frequencies,_n . Three mode interaction (₂ 3₁ and _{3 1} + 2₂) is considered and its influence on the response is studied. The case of two mode interaction (₂ 3₁) is also considered to compare it with the case of three mode interaction. The straight beam experiencing mid-plane stretching is governed by a nonlinear partial differential equation. By using Galerkin's method the governing equation is reduced to a system of nonautonomous ordinary differential equations. The method of multiple scales is applied to solve the system. Steady-state responses and their stability are examined. Results of numerical investigations show that there exists no significant difference between both modal interactions' influences on the responses.     

Calculation of heat transfer accompanying hypersonic three-dimensional flow of air over slender blunt-nosed cones

The effective length method [1, 2] has been used to make systematic calculations of the heat transfer for laminar and turbulent boundary layers on slender blunt-nosed cones at small angles of attack ( + 5° in a separationless hypersonic air stream dissociating in equilibrium (half-angles of the cones 0 20°, angles of attack 0 15°, Mach numbers 5 M 25). The parameters of the gas at the outer edge of the boundary layer were taken equal to the inviscid parameters on the surface of the cones. Analysis of the results leads to simple approximate dependences for the heat transfer coefficients.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 173–177, September–October, 1981.     

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

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

---

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

---

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

Nonlinear rheological behavior of a concentrated spherical silica suspension

Linear and nonlinear viscoelastic properties were examined for a 50 wt% suspension of spherical silica particles (with radius of 40 nm) in a viscous medium, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture. The effective volume fraction of the particles evaluated from zero-shear viscosities of the suspension and medium was 0.53. At a quiescent state the particles had a liquid-like, isotropic spatial distribution in the medium. Dynamic moduli G^* obtained for small oscillatory strain (in the linear viscoelastic regime) exhibited a relaxation process that reflected the equilibrium Brownian motion of those particles. In the stress relaxation experiments, the linear relaxation modulus G(t) was obtained for small step strain (0.2) while the nonlinear relaxation modulus G(t, ) characterizing strong stress damping behavior was obtained for large (>0.2). G(t, ) obeyed the time-strain separability at long time scales, and the damping function h() (–G(t, )/G(t)) was determined. Steady flow measurements revealed shear-thinning of the steady state viscosity () for small shear rates (< ^–1; = linear viscoelastic relaxation time) and shear-thickening for larger (>^–1). Corresponding changes were observed also for the viscosity growth and decay functions on start up and cessation of flow, + (t, ) and _– (t, ). In the shear-thinning regime, the and dependence of ₊(t,) and _–(t,) as well as the dependence of () were well described by a BKZ-type constitutive equation using the G(t) and h() data. On the other hand, this equation completely failed in describing the behavior in the shear-thickening regime. These applicabilities of the BKZ equation were utilized to discuss the shearthinning and shear-thickening mechanisms in relation to shear effects on the structure (spatial distribution) and motion of the suspended particles.Dedicated to the memory of Prof. Dale S. Parson     

Calculation of intense blowing on a blunt body and an airfoil

Various aspects of the problem of intense blowing through the surface of bodies have, been theoretically studied by a number of authors, within the framework of inviscid flow theory. A detailed bibliography on this topic is given, e.g., in [1, 2]. The well-known approaches to solution of this problem have a limited area of application. For example, asymptotic methods can be used for hypersonic flow regimes only at relatively low levels of the blown gas momentum ( = ² = _ov_o ²/ V ² 1). The same limitation applies to the numerical method of straight lines [2]. The forward Eulerian calculation schemes [3, 4] smear the contact discontinuity severely, and cannot handle the case where the blown gas and the gas in the incident flow have different thermodynamic properties (_o ). This paper presents results of a numerical investigation of supersonic flow over two-dimensional and axisymmetric bodies with intense blowing on the forward surface, performed using a time-dependent finite-difference method [5] with an explicit definition of the contact interface between the two cases. The calculations encompass a family of elliptic cylinders with semiaxis ratio 0.5 4, a flat-face cylinder, and a flat plate with rounding near the midsection, with variations in the blowing law, the incident flow Mach number M (3 M 10), the adiabatic indices, and the blowing parameter 0 0.5.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 117–124, January–February, 1977.In conclusion, the authors thank T. S. Novikov and I. D. Sandomirskii, who took part In the present calculations.     

Heat transfer from laminar flow in ducts with elliptic cross-section

Summary The cooling of a hot fluid in laminar Newtonian flow through cooled elliptic tubes has been calculated theoretically. Numerical data have been computed for the two values 1.25 and 4 of the axial ratio of the elliptic cross-section . For =1.25 the influence of non-zero thermal resistance between outmost fluid layer and isothermal surroundings has also been investigated. Special attention has been given to the distribution of heat flux around the perimeter; when increases the flux varies more with the position at the circumference. This positional dependence becomes less pronounced, however, as the (position-independent) thermal resistance of the wall increases.Flattening of the conduit, while maintaining its cross-sectional area constant, improves the cooling. Comparison with rectangular pipes shows that this improvement is not as marked with elliptic as with rectangular pipes.Nomenclature A _k =A _{m, n} coefficients of expansion (6) - a, b half-axes of ellipse, b<a - a _p =a _{r, s} coefficients of representation (V) - D hydraulic diameter, = 4S/P; S = cross-sectional area, P = perimeter - D _e equivalent diameter, according to (13) - n coordinate (outward) normal to the tube wall - T temperature of fluid - T _i temperature of fluid at the inlet - T _s temperature of surroundings - v ₀ mean velocity of fluid - v _z longitudinal velocity of fluid - x, y carthesian coordinates coinciding with axes of ellipse - z coordinate in flow direction - , dimensionless half-axes of ellipse, =a/D and =b/D - _t heat transfer coefficient from fluid at bulk temperature to surroundings; equation (11) - _w heat transfer coefficient at the wall; equation (3) - axial ratio of ellipse, = a/b = / - , , , dimensionless coordinates; =x/D, =y/D, =z/D, =n/D - dimensionless temperature, = (T–T _s)/(T _i–T _s) - ₀ cup-mixing mean value of ; equation (10) - thermal conductivity of fluid - _m,n = _k eigenvalue - c volumetric heat capacity of fluid - _{m, n} = _k = _k eigenfunction; equations (6) and (I) - Nu total Nusselt number, = _t D/ - Nusselt number at large distance from the inlet - Nu _w wall Nusselt number, = _w D/, based on _w - Pé Péclet number, = ₀ Dc/     

Anlaufkorrektur und Elastizität bei der Kapillarströmung von Lösungen

Zusammenfassung Aus der Anlaufkorrektur kann man nach einer Rechnung vonFromm, die auf dem Maxwellschen Modell basiert, eine Relaxationszeit und eine korrigierte Viskosität _c ermitteln. Der Quotient _c/ stellt einen Schermodul dar. Diese Größe wird für Lösungen von Cellulosetrinitrat in Butylacetat, Polyvinylacetat in Dioxan, Polystyrol in Toluol, Polyacrylamid in Wasser, und Viskose, in Abhängigkeit von der Konzentrationc und dem SchergefälleD ermittelt. Es zeigt sich, daß _c/ etwa im Wendepunkt der Fließkurven eine Art Plateau oder ein flaches Maximum zeigt und in diesem Plateaubereich eine lineare Abhängigkeit von der Konzentration. Die absolute Größe von _c/ ist jedoch um Größenordnungen geringer, als sie nach der Formel vonRouse bzw.Bueche für die erste Relaxationszeit eines Verhängungsnetzwerkes zu erwarten wäre. Das wird so gedeutet, daß bei dem hohen Schergefälle, das bei den Messungen herrschte (D etwa 10⁴ sec^–1), ein Teil der Verhängungen zerstört ist, wodurch die Relaxationszeit vergrößert und der Schermodul verkleinert wird.

---

Summary From the end-correction, according to a calculation byFromm based upon theMaxwell-model, a relaxation time and a corrected viscosity _c can be obtained. The quotient _c/ represents a shear modulus. Its value is determined for solutions of cellulosetrinitrate in butylacetate, polyvinylacetate in dioxane, polystyrene in toluene, polyacryloamide in water, and viscose, in dependence of concentrationc and shear rateD. It is found, that _c/ shows a plateau or a flat maximum at the inflection point of the flow curves. In this range, a linear dependence on concentration is found too. The absolute value of _c/, however, is smaller by orders of magnitude than that to be expected for the first relaxation time of an entanglement network according to the formulas byRouse resp.Bueche. This is explained by a partial disruption of entanglements in the high shear rate prevailing at the experiments (D about 10⁴ sec^–1), which effects an increase of the relaxation time and a decrease of the shear modulus.

---

---

Vorgetragen auf der Jahrestagung der Deutschen Rheologen in Bad Ems vom 18.–19. Mai 1967.     

Boundary layer stability on a rotating disk with corotation of the surrounding fluid

The stability of the steady self-similar flow in the boundary layer on a rotating disk of infinite radius with corotation of the surrounding fluid is analyzed by the normal mode method. The spectral problem for infinitesimal three-dimensional disturbances is solved by a collocation method with expansion of the amplitude functions in Chebyshev polynomials. It is established that for all values of the parameter 0, equal to the ratio of the angular velocities of the fluid and the disk, the lower critical Reynolds number is determined byA-type, waves, whose development is governed by the parallel instability mechanism typical of an Ekman layer. TheB-type instability, associated with the presence of an inflection point on the velocity profile, disappears when 4. The neutral surfaces are calculated for Karman flow (=0) and Bödewadt flow (). It is found that in Karman flowA-type waves may grow at values of the Reynolds number several times smaller than the critical Reynolds number for spiral vortices. The results of the analysis are compared with the available experimental data.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.5, pp. 69–77, September–October, 1992.     

Some self-similar solutions of the Leibenzon equation

Self-similar one-dimensional solutions of the Leibenzon equation c²_t= _zz ^k (z 0, k 2) are considered. Approximate solutions are constructed for the two cases in which the initial value = ₁ = const > 0 and on the boundary either a constant value = ₂ < ₁ is maintained or the flow (directed outwards) is given. In the first problem the dependence of the boundary flow on the governing parameters is determined. A characteristic property of the types of motion in question is the existence near the boundary of a region, expanding with time, in which the flow is almost independent of the coordinate.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 145–150, September–October, 1991.     

Turbulent dispersion of droplets for air flow in a pipe

An optical technique was used to study the dispersion of 50 m, 90 m and 150 m droplets downstream of a source located in the center of a vertical pipe through which turbulent air is flowing. A turbulent dispersion coefficient, _P, and a mean-square of the fluctuations in the turbulent velocity, v _p ² , are determined from the change of the measured mean-square displacement of the droplets over the pipe cross section with time. The interesting aspect of the experiments is that they explored conditions where the inertia of the particles is believed to be a much more important effect than that of the crossing of trajectories associated with the inequality of the average velocities of the particles and the fluid.     

Exact solutions and mathematical properties of boundary‐value problems for dynamic‐diffusion boundary layers

The paper studies boundaryvalue problems for dynamicdiffusion boundary layers occurring near a vertical wall at high Schmidt numbers and for dynamic boundary layers whose inner edge is adjacent to the dynamicdiffusion layers. Exact solutions for boundary layers at small and large times are derived. The wellposedness of the boundaryvalue problem for a steady dynamicdiffusion layer is studied.