The virtual origin of riblets in an open channel flow

In this paper, a method using the mean velocity profiles for the buffer layer was developed for the estimation of the virtual origin over a riblets surface in an open channel flow. First, the standardized profiles of the mixing length were estimated from the velocity measurement in the inner layer, and the location of the edge of the viscous layer was obtained. Then, the virtual origins were estimated by the best match between the measured velocity profile and the equations of the velocity profile derived from the mixing length profiles. It was made clear that the virtual origin and the thickness of the viscous layer are the function of the roughness Reynolds number. The drag variation coincided well with other results.Nomenclature f _r skin friction coefficient - f _ro skin friction coefficient in smooth channel at the same flow quantity and the same energy slope - g gravity acceleration - H water depth from virtual origin to water surface - H ⁺ u*H/ - H false water depth from top of riblets to water surface - H ⁺ u*H/ - I _e streamwise energy slope - I _b bed slope - k riblet height - k ⁺ u*k/ - l mixing length - l _s standardized mixing length - Q flow quantity - Re Reynolds number volume flow/unit width/v - s riblet spacing - u mean velocity - u* friction velocity = - u* false friction velocity = - y distance from virtual origin - y distance from top of riblet - y ₀ distance from top of riblet to virtual origin - y _v distance from top of riblet to edge of viscous layer - y ⁺ u*y/ - y ⁺ u*y/ - y ₀ ⁺ u*y ₀/ - u ⁺ u*y/ - shifting coefficient for standardization - thickness of viscous layer=y ₀+y - ⁺ u*/ - ⁺ u*/ - eddy viscosity - ridge angle - v kinematic viscosity - density - shear stress     

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

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

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     

A rotating hot-wire technique for spatial sampling of disturbed and manipulated duct flows

The documentation and control of flow disturbances downstream of various open inlet contractions was the primary focus with which to evaluate a spatial sampling technique. An X-wire probe was rotated about the center of a cylindrical test section at a radius equal to one-half that of the test section. This provided quasi-instantaneous multi-point measurements of the streamwise and azimuthal components of the velocity to investigate the temporal and spatial characteristics of the flowfield downstream of various contractions. The extent to which a particular contraction is effective in controlling ingested flow disturbances was investigated by artificially introducing disturbances upstream of the contractions. Spatial as well as temporal mappings of various quantities are presented for the streamwise and azimuthal components of the velocity. It was found that the control of upstream disturbances is highly dependent on the inlet contraction; for example, reduction of blade passing frequency noise in the ground testing of jet engines should be achieved with the proper choice of inlet configurations.List of symbols K _uv correlation coefficient= - P percentage of time that an azimuthal fluctuating velocity derivative dv/d is found - U streamwise velocity component U=U (, t) - V azimuthal or tangential velocity component due to flow and probe rotation V=V (, t) - mean value of streamwise velocity component - U _m resultant velocity from and - mean value of azimuthal velocity component induced by rotation - u fluctuating streamwise component of velocity u=u(, t) - v fluctuating azimuthal component of velocity v = v (, t) - u phase-averaged fluctuating streamwise component of velocity u=u(0) - v phase-averaged fluctuating azimuthal component of velocity v=v() - û average of phase-averaged fluctuating streamwise component of velocity (u()) over cases I-1, II-1 and III-1 û = û() - average of phase-averaged fluctuating azimuthal component of velocity (v()) over cases I-1, II-1 and III-1 - u fluctuating streamwise component of velocity corrected for non-uniformity of probe rotation and/or phase-related vibration u = u(0, t) - v fluctuating azimuthal component of velocity corrected for non-uniformity or probe rotation and/or phase-related vibration v=v (, t) - u ² rms value of corrected fluctuating streamwise component of velocity - rms value of corrected fluctuating azimuthal component of velocity - phase or azimuthal position of X-probe     

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     

Analysis of suddenly started laminar flow in the entrance region of a circular tube

Suddenly started laminar flow in the entrance region of a circular tube, with constant inlet velocity, is investigated analytically by using integral momentum approach. A closed form solution to the integral momentum equation is obtained by the method of characteristics to determine boundary layer thickness, entrance length, velocity profile, and pressure gradient.Nomenclature M(, , ) a function - N(, , ) a function - p pressure - p* p/1/2U ², dimensionless pressure - Q(, , ) a function - R radius of the tube - r radial distance - Re 2RU/, Reynolds number - t time - U inlet velocity, constant for all time, uniform over the cross section - u velocity in the boundary layer - u* u/U, dimensionless velocity - u ₁ velocity in the inviscid core - x axial distance - y distance perpendicular to the axis of the tube - y* y/R, dimensionless distance perpendicular to the axis - boundary layer thickness - * displacement thickness - /R, dimensionless boundary layer thickness - momentum thickness - absolute viscosity of the fluid - /, kinematic viscosity of the fluid - x/(R Re), dimensionless axial distance - density of the fluid - tU/(R Re), dimensionless time - _w wall shear stress     

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.     

Shear softening and thixotropic properties of wheat flour doughs in dynamic testing at high shear strain

Shear softening and thixotropic properties of wheat flour doughs are demonstrated in dynamic testing with a constant stress rheometer. This behaviour appears beyond the strictly linear domain (strain amplitude ₀ 0.2%),G,G and |^*| decreasing with ₀, the strain response to a sine stress wave yet retaining a sinusoidal shape. It is also shown thatG recovers progressively in function of rest time. In this domain, as well as in the strictly linear domain, the Cox-Merz rule did not apply but() and | ^*())| may be superimposed by using a shift factor, its value decreasing in the former domain when ₀ increases. Beyond a strain amplitude of about 10–20%, the strain response is progressively distorted and the shear softening effects become irreversible following rest.     

On shock layer method for the hypersonic flow problems

Chernyi’s series method[1] is not proper for the case that(γ-l)/(γ+l)<<2/(γ+1)×M²sin²β (γ=c_p/c_v-adiabatic index number, M-Much number, β-shock incidence). In this paper, we only suppose that in the neighbour of the shock, there exists a shock layer in which the density of the gas is very big, but we do not remove the case that (γ-1)/(γ+1)<<2/(γ+1)M²sin²β.     

Effects of temperature and pressure on the rheological behavior of montmorillonite sols

Summary The effects of high temperatures up to 180 C and high pressures up to 560 kg/cm² on the rheological properties of pure montmorillonite suspensions with well-defined base-exchange cations have been investigated. The suspensions behave very much asBingham plastics according to the equation=+, in which is the shear stress,D the shear rate, the plastic viscosity, and the yield stress which is largely a measure of residual flocculation.The observed effects depend strongly on the interparticle forces that govern the colloidal stability and the rheological behavior of the suspensions. One can distinguish between two categories of suspensions:P-type sols in which the clay particles are associated through Coulombic attraction between positive edges and negative faces and are located in a primary potential energy minimum, andS-type sols in which the particles are associated edge-to-edge and are located in a weaker secondary potential minimum obtained by the summation ofvan der Waals attraction and double layer repulsion.Both theBingham yield stress and the plastic viscosity of theP-type sols decrease with increasing temperature. The temperature dependence of follows theArrhenius equation. TheP-type suspensions are either weakly or not at all thixotropic at room temperature and are definitely non-thixotropic at higher temperatures. Pressure slightly increases both the plastic viscosity and the yield stress.TheS-type sols, on the other hand, display an increase in yield stress and degree of thixotropy with increasing temperature and generally a decrease in the plastic viscosity. This behavior is modified in the case of Ca-montmorillonite suspensions, in which both and pass through a maximum at 150C, followed by a decline. The maximum can be explained by disaggregation of face-to-face aggregated clay packets.Pressure causes a decrease both in and the degree of thixotropy in theS-type suspensions, while it causes a slight increase in the plastic viscosity. This behavior is a consequence of the destruction of the hydration shell caused by high pressure.

---

Zusammenfassung Der Einfluß von hohen Temperaturen (bis 180 C) und hohen Drucken (bis zu 560 kp/cm²) auf die rheologischen Eigenschaften reiner Suspensionen von Montmorillonit mit definierten Austauschionen wurde untersucht. Die Suspensionen folgen dem Binghamschen Flie\gesetz gemÄ\ der Gleichung=+, worin die Schubspannung,D das GeschwindigkeitsgefÄlle, dieBingham- ViskositÄt und die Flie\spannung darstellen. Letztere zeigt in erster Linie das Ausma\ der Koagulation an.Die beobachteten Effekte hÄngen stark von den KrÄften zwischen den Teilchen ab, welche die StabilitÄt der Kolloide und das rheologische Verhalten der Suspensionen bestimmen. Man kann zwei Kategorien von Suspensionen unterscheiden: Sole vom sog.P-Typ, in denen die Tonteilchen durch Coulombsche AnziehungskrÄfte zwischen positiven Kanten und negativen FlÄchen assoziiert sind und sich in einem primÄren Potentialminimum befinden, und Sole vomS-Typ, in denen die Teilchen Kante-zu-Kante assoziiert sind und sich in einem flacheren sekundÄren Potentialminimum befinden, welches durch das Zusammenwirken vonvan der Waalsschen AnziehungskrÄften und von Absto\ungskrÄften infolge der Wechselwirkung der elektrischen Doppelschichten entsteht.Sowohl die Flie\grenze als auch dieBingham-ViskositÄt der Sole vomP-Typ nehmen mit wachsender Temperatur ab. Die TemperaturabhÄngigkeit von folgt derArrhenius-Gleichung. Die Suspensionen vomP-Typ sind bei Zimmertemperatur entweder gar nicht oder nur schwach thixotrop, wÄhrend sie bei hoher Temperatur in keinem Fall thixotrop sind. Mit wachsendem Druck erhöht sich sowohl dieBingham-ViskositÄt als auch die Flie\grenze ein wenig.Bei den Solen vomS-Typ andererseits steigt die Flie\spannung und der Grad der Thixotropie mit steigender Temperatur, wÄhrend dieBingham-ViskositÄt im allgemeinen abnimmt. Bei Ca-Montmorillonit-Suspensionen ist das Verhalten etwas anders: Sowohl als auch erreichen bei 150 C ein Maximum und fallen dann wieder ab. Das Maximum kann durch Desaggregation flÄchenhaft aggregierter Tonpartikel erklÄrt werden. Bei steigendem Druck fÄllt sowohl als auch der Grad der Thixotropie in den Suspensionen vomS-Typ ab, wÄhrend dieBingham-YiskositÄt leicht ansteigt. Dieses Verhalten ist eine Folge der dann eintretenden Zerstörung der Solvathüllen.

     

The dynamic performance of the Weissenberg rheogoniometer III. Large amplitude oscillatory shearing — harmonic analysis

The harmonic content of the nonlinear dynamic behaviour of 1% polyacrylamide in 50% glycerol/water was studied using a standard Model R 18 Weissenberg Rheogoniometer. The Fourier analysis of the Oscillation Input and Torsion Head motions was performed using a Digital Transfer Function Analyser.In the absence of fluid inertia effects and when the amplitude of the (fundamental) Oscillation Input motion _I is much greater than the amplitudes of the Fourier components of the Torsion Head motion _Tn empirical nonlinear dynamic rheological propertiesG _n (, ⁰),G _n (, ⁰) and/or _n (, ⁰), _n (, ⁰) may be evaluated without a-priori-knowledge of a rheological constitutive equation. A detailed derivation of the basic equations involved is presented.Cone and plate data for the third harmonic storage modulus (dynamic rigidity)G ₃ (, ⁰), loss modulusG ₃ (, ⁰) and loss angle ₃ (, ⁰) are presented for the frequency range 3.14 × 10^–2 1.25 × 10² rad/s at two strain amplitudes, _CP ⁰ = 2.27 and 4.03. Composite cone and plate and parallel plates data for both the third and fifth harmonic dynamic viscosities ₃ (, ⁰), _S (, ⁰) and dynamic rigiditiesG ₃ (, ⁰),G ₅ (, ⁰) are presented for strain amplitudes in the ranges 1.10 _CP ⁰ 4.03 and 1.80 _PP ⁰ 36 for a single frequency, = 3.14 × 10^–1 rad/s. Good agreement was obtained between the results from both geometries and the absence of significant fluid inertia effects was confirmed by the superposition of the data for different gap widths.     

An “ideal equivalent gas method” for the study of shock waves in supersonic real gas flows

Summary A method developed by the author for the systematic study of the thermodynamic and dynamic properties of the gas behind a shock wave is reported.The method is applicable to supersonic flow regimes for which the excitation, dissociation and ionization effects invalidate the usually adopted hypothesis of ideal gas.An Ideal Equivalent Gas, having the ratio of the specific heats _s dependent on Mach number and altitude of flight is postulated.On the basis of the mass, momentum and energy conservation equations, valid through the shock wave, the relations defining the thermodynamic and dynamic state of the gas behind the shock wave are derived. These relations establish an extension of the classic relations valid for the ideal gas and reduce to them identically for _s=.The dependence of the ratio of specific heats _s of the Ideal Equivalent Gas on Mach number and altitude has been established, over a wide range, on the basis of the real gas solutions derived by Huber.

---

Sommario Nella presente nota viene esposto un metodo sviluppato dall'autore per lo studio sistematico dello stato termodinamico e dinamico del gas a valle di un'onda d'urto in regime supersonico, allorchè cioè gli effetti dell'eccitazione dei gradi di libertà vibrazionali delle molecole e della loro dissociazione e successiva ionizzazione invalidano l'ipotesi di gas ideale generalmente adottata.Viene definito un gas ideale equivalente avente rapporto dei calori specifici _s funzione del numero di Mach e della quota di volo ed in base alle equazioni di conservazione della massa, della quantità di moto e dell'energia, valide attraverso all'onda d'urto, vengono derivate delle relazioni definenti lo stato termodinamico e dinamico del gas a valle dell'onda d'urto. Tali relazioni costituiscono una estensione delle classiche relazioni dell'urto valide per il gas ideale alle quali si riducono per _s=.La dipendenza del rapporto dei calori specifici _s del gas ideale equivalente, dal numero di Mach e dalla quota è stata stabilita sulla base delle soluzioni ottenute da Huber per il gas reale.

---

     

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

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

Über den Einfluß von Querkrümmung und Dichte bei inhomogenen turbulenten Freistrahlen

Zusammenfassung Zur Berechnung turbulenter Strömungen wird das k--Modell im Ansatz für die turbulente Scheinzähigkeit erweitert, so daß es den Querkrümmungs- und Dichteeinfluß auf den turbulenten Transportaustausch erfaßt. Die dabei zu bestimmenden Konstanten werden derart festgelegt, daß die bestmögliche Übereinstimmung zwischen Berechnung und Messung erzielt wird. Die numerische Integration der Grenzschichtgleichungen erfolgt unter Verwendung einer Transformation mit dem Differenzenverfahren vom Hermiteschen Typ. Das erweiterte Modell wird auf rotationssymmetrische Freistrahlen veränderlicher Dichte angewendet und zeigt Übereinstimmung zwischen Rechnung und Experiment.

---

On the influence of transvers-curvature and density in inhomogeneous turbulent free jets

The prediction of turbulent flows based on the k- model is extended to include the influence of transverse-curvature and density on the turbulent transport mechanisms. The empirical constants involved are adjusted such that the best agreement between predictions and experimental results is obtained. Using a transformation the boundary layer equations are solved numerically by means of a finite difference method of Hermitian type. The extended model is applied to predict the axisymmetric jet with variable density. The results of the calculations are in agreement with measurements.

---

Bezeichnungen Wirbelabsorptionskoeffizient - c_i Massenkonzentration der Komponente i - c_D, c_L, c, c₁, c₂ Konstanten des Turbulenzmodells - d Düsendurchmesser - E bezogene Dissipationsrate - f bezogene Stromfunktion - f Korrekturfunktion für die turbulente Scheinzähigkeit - j turbulenter Diffusionsstrom - k Turbulenzenergie - k_i Schrittweite in -Richtung - K dimensionslose Turbulenzenergie - L turbulentes Längenmaß - M_i Molmasse der Komponente i - p Druck - allgemeine Gaskonstante - r Querkoordinate - r_0,5 Halbwertsbreite der Geschwindigkeit - r_0,5c Halbwertsbreite der Konzentration - T Temperatur - u Geschwindigkeitskomponente in x-Richtung - v Geschwindigkeitskomponente in r-Richtung - x Längskoordinate - y allgemeine Funktion - Y_i diskreter Wert der Funktion y - Relaxationsfaktor für Iteration - turbulente Dissipationsrate - transformierte r-Koordinate - kinematische Zähigkeit - Exponent - transformierte x-Koordinate - Dichte - _k, Konstanten des Turbulenzmodells - Schubspannung - allgemeine Variable - Stromfunktion - Turbulente Transportgröße Indizes 0 Strahlanfang - m auf der Achse - r mit Berücksichtigung der Krümmung - t turbulent - mit Berücksichtigung der Dichte - im Unendlichen - Schwankungswert oder Ableitung einer Funktion - – Mittelwert Herrn Professor Dr.-Ing. R. Günther zum 70. Geburtstag gewidmet     

Relatively Compact Orbits and Compact Attractors for a Class of Nonlinear Evolution Equations

In the present paper we consider the nonlinear evolution equation u+AuG(u), where A:D(A)XX is m-accretive with (I+A)^–1 compact for some >0, and is continuous, and we prove that the orbit is relatively compact if and only if u is uniformly continuous, and both u and G^u are bounded on . In the same spirit, we derive conditions for orbits of bounded sets to have compact attractors. Some consequences and an example from age-structured population dynamics illustrate the effectiveness of the abstract result.     

Conversion of oscillatory shear data from highly filled polymers by direct numerical integration

A procedure was developed to enable the direct numerical integration of the Fourier integral transform equation relating G(t) to G() by considering integration limits that vary as a function of time and which define a range of discrete sub-intervals within the complete frequency domain data set. The method provides results that are in very close agreement to results determined from a relaxation spectrum. However, at low values of time the solution to the variable limit integral transform is sensitive to the absence of a contribution beyond the upper experimental limit of the frequency domain data. G(t) results determined from the conversion of shifted master G() experimental data using the variable limit integral transformation and the relaxation spectrum compared favourably with actual shifted master G(t) experimental data. The former curves were characterised by the same form and trend as the experimental results, confirming that the underlying viscoelastic behaviour is well represented. While the variable limit Fourier integral transform procedure provides a good approximation to relaxation spectrum results, the latter is clearly the more robust method of converting data from the frequency to the time domain.It was observed that the time-temperature superposition procedure used in the construction of shifted master curves can magnify potential differences between the shifted G(t) values determined from the conversion of G() data and the actual experimental G(t) results, when compared to data that has not been shifted to a master curve.     

The mirror relations and nonlinear viscoelasticity of polymer melts

An attempt is made to incorporate into a quasilinear viscoelastic constitutive equation of the Boltzmann superposition type the two mirror relations of Gleissle, as well as his relation between the steady-state first normal-stress difference and the shear viscosity curve. It is shown that the three relations can hold separately within this constitutive model, but not simultaneously, because they require a different nonlinear strain measure, namelyS ₁₂ () = – a ( – 1) (a = 0 for 1,a = 1 for 1) for the mirroring of the viscosities,S ₁₂ () = – a (–k ²/) (a = 0 for k, a = 1 for k) for the mirroring of the first normal-stress coefficients, and for the third relation. Here denotes the shear strain and erf the error function. Experimental data on melts of a low-density polyethylene, a high-density polyethylene and a polypropylene show that the mirror relations are passable approximations, but that the third relation meets reality surprisingly close if the right value ofk is used.     

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²     

Determination of permeability tensors for two-phase flow in homogeneous porous media: Theory

In this paper we continue previous studies of the closure problem for two-phase flow in homogeneous porous media, and we show how the closure problem can be transformed to a pair of Stokes-like boundary-value problems in terms of pressures that have units of length and velocities that have units of length squared. These are essentially geometrical boundary value problems that are used to calculate the four permeability tensors that appear in the volume averaged Stokes' equations. To determine the geometry associated with the closure problem, one needs to solve the physical problem; however, the closure problem can be solved using the same algorithm used to solve the physical problem, thus the entire procedure can be accomplished with a single numerical code.Nomenclature a a vector that maps V onto , m^-1. - A a tensor that maps V onto . - A area of the - interface contained within the macroscopic region, m². - A area of the -phase entrances and exits contained within the macroscopic region, m². - A area of the - interface contained within the averaging volume, m². - A area of the -phase entrances and exits contained within the averaging volume, m². - Bo Bond number (= (=(–)g²/). - Ca capillary number (= v/). - g gravitational acceleration, m/s². - H mean curvature, m^-1. - I unit tensor. - permeability tensor for the -phase, m². - viscous drag tensor that maps V onto V. - ^* dominant permeability tensor that maps onto v , m². - ^* coupling permeability tensor that maps onto v , m². - characteristic length scale for the -phase, m. - l characteristic length scale representing both and , m. - L characteristic length scale for volume averaged quantities, m. - n unit normal vector directed from the -phase toward the -phase. - n unit normal vector representing both n and n . - n unit normal vector representing both n and n . - P pressure in the -phase, N/m². - p superficial average pressure in the -phase, N/m². - p intrinsic average pressure in the -phase, N/m². - p –p , spatial deviation pressure for the -phase, N/m². - r ₀ radius of the averaging volume, m. - r position vector, m. - t time, s. - v fluid velocity in the -phase, m/s. - v superficial average velocity in the -phase, m/s. - v intrinsic average velocity in the -phase, m/s. - v –v , spatial deviation velocity in the -phase, m/s. - V volume of the -phase contained within the averaging volmue, m³. - averaging volume, m³. Greek Symbols V /, volume fraction of the -phase. - viscosity of the -phase, Ns/m². - density of the -phase, kg/m³. - surface tension, N/m. - (v +v ^T ), viscous stress tensor for the -phase, N/m².     

Mixed convection flow in narrow vertical ducts

The mixed convection flow in a vertical duct is analysed under the assumption that , the ratio of the duct width to the length over which the wall is heated, is small. It is assumed that a fully developed Poiseuille flow has already been set up in the duct before heat from the wall causes this to be changed by the action of the buoyancy forces, as measured by a buoyancy parameter . An analytical solution is derived for the case when the Reynolds numberRe, based on the duct width, is of 0 (1). This is extended to the case whenRe is 0 (^–1) by numerical integrations of the governing equations for a range of values of representing both aiding and opposing flows. The limiting cases, || 1 andR=Re of 0 (1), andR and both large, with of 0 (R ^1/3) are considered further. Finally, the free convection limit, large with R of 0 (1), is discussed.

---

Mischkonvektion in engen senkrechten Rohren

Zusammenfassung Mischkonvektion in einem senkrechten Rohr wird unter der Voraussetzung untersucht, daß das Verhältnis der Rohrbreite zur Länge, über welche die Wand beheizt wird, klein ist. Es wird angenommen, daß sich bereits eine voll entwickelte Poiseuille-Strömung in dem Rohr eingestellt hat, bevor Antriebskräfte, gemessen mit dem Auftriebsparameter , aufgrund der Wandbeheizung die Strömung verändern. Es wird eine analytische Lösung für den Fall erhalten, daß die mit der Rohrbreite als charakteristische Länge gebildete Reynolds-ZahlRe konstant ist. Dies wird mittels einer numerischen Integration der wichtigsten Gleichungen auf den FallRe =f (^–1) sowohl für Gleich- als auch für Gegenstrom ausgedehnt. Weiterhin werden die beiden Grenzfälle betrachtet, wenn || 1 undR=Re konstant ist, sowieR und beide groß mit proportionalR ^1/3. Schließlich wird der Grenzfall der freien Konvektion, großes mit konstantem R, diskutiert.

---

Nomenclature g acceleration due to gravity - Gr Grashof number - G modified Grashof number - h duct width - l length of the heated section of the duct wall - p pressure - Pr Prandtl number - Q flow rate through the duct - Q ₀ heat transfer on the wally=0 - Q ₁ heat transfer on the wally=1 - Re Reynolds number - R modified Reynolds number - T temperature of the fluid - T ₀ ambient temperature - T applied temperature difference - u, velocity component in thex-direction - v, velocity component in they-direction - x, co-ordinate measuring distance along the duct - y, co-ordinate measuring distance across the duct - buoyancy parameter - ₀ modified buoyancy parameter, ₀=R ^–1/3 - coefficient of thermal expansion - ratio of duct width to heated length, =h/l - (non-dimensional) temperature - _w applied temperature on the wally=0 - kinematic viscosity - density of the fluid - ₀ shear stress on the wally=0 - ₁ shear stress on the wally=1 - stream function