The rheology of Hydroxy-propylguar (HPG) solutions and its influence on the flow through a porous medium and turbulent tube flow,respectively (Part 1)

The rheology of aqueous HPG solutions in the range 100 wppm to 5000 wppm is investigated. The flow through a porous medium and turbulent tube flow, respectively, of these solutions is studied as well. Especially with respect to the higher concentrations, the data correlate nicely only after the effect of shear is extracted, i.e., after the variable viscosity is taken into account. This is accomplished by working with an apparent viscosity _c , defined such that, the Hagen Poiseuille law (with _c ) holds in laminar tube flow.     

Sand Erosion in Axial Flow Conditions

A mathematical model of sand erosion in axial flow conditions is presented. The basic mass balance equations and sand erosion constitutive equation were given in Vardoulakis et al. (1996). As opposed to reference Vardoulakis et al. (1996), we consider here the extreme case where convection is null and hydrodynamic dispersion dominates. In addition, Brinkman's extension of Darcy's law is adopted to account for a smooth transition between channel flow and Darcian flow. The set of governing PDE's is presented in dimensionless form and is solved numerically. In concordance with the basic constitutive equation for erosion kinetics, the analysis shows that erosion progresses in time as a front of high transport concentration. This result is justified by the highly non-linear character of the erosion source term which dominates in the diffusion-like governing equation.     

Invasion Percolation on Correlated and Elongated Lattices: Implications for the Interpretation of Residual Saturations in Rock Cores

The invasion percolation model is used to investigate the effect of correlated heterogeneity on capillary dominated displacements in porous media. The breakthrough and residual saturations are shown to be strongly influenced by the correlations. Correlated heterogeneity leads to lower residual saturations than those observed in random systems and the scatter commonly observed in laboratory core measurements of the residual saturations can be attributed to the presence of such heterogeneity at the pore scale. Invasion percolation computations on elongated lattices, those with a geometry of L ^d–1 × nL where n denotes the aspect ratio, show that residual saturations for systems with correlated heterogeneity exhibit a strong dependence on aspect ratio. This effect is not considered by experimentalists who advocate the use of long (high aspect ratio) cores in order to overcome end-effects in experiments on shorter cores. A new scaling law is proposed for the residual saturations in elongated systems with correlated heterogeneity, and is confirmed by numerical simulations.     

On the properties of compressible gas flow in a porous media

The flow of an adiabatic gas through a porous media is treated analytically for steady one- and two-dimensional flows. The effect on a compressible Darcy flow by inertia and Forchheimer terms is studied. Finally, wave solutions are found which exhibit a cut-off frequency and a phase shift between pressure and velocity of the gas, with the velocity lagging behind the pressure.Nomenclature A area of tube for one-dimensional flow - B drag coefficient associated with Forchheimer term - c speed of sound - M Mach number - p ^* gas pressure - p dimensionless gas pressure - s coordinate along the axis of tube - t ^* time variable - t dimensionless time variable - V^* gas velocity in the porous media - V dimensionless gas velocity Greek Letters ratio of specific heat capacities - phase angle between gas pressure and velocity for linear waves - parameter indicating the importance of the inertia term - viscosity - _p natural frequency of the porous media - ^* gas density - dimensionless gas density - parameter indicating the importance of the Forchheimer term - porosity of porous media - velocity potential - stream function     

Elongational flow opto-rheometry for polymer melts

Polymer melt elongation is one of the most important procedures in polymer processing. To understand its molecular mechanisms, we constructed an elongational flow opto-rheometer (EFOR) in which a high precision birefringence apparatus of reflection-double path type was installed into a Meissner's new elongational rheometer of a gas cushion type (commercialized as RME from Rheometric Scientific) just by mounting a small reflecting mirror at the center of the RME's sample supporting table. The EFOR enabled us to achieve simultaneous measurements of tensile stress (t) and birefringence n(t) as a function of time t under a given constant strain rate within the range of 0.001 to 1.0s^–1. (t) can be monitored upto the maximum Hencky strain (t) of 7 as attained, in principle, with RME, while the measurable range of the phase difference in the birefringence was 0 to 250 (0 to 79 100 nm for He-Ne laser light) within the accuracy of ±0.1 (±31.6 nm) up to (t) 4. The performance was tested on an anionically polymerized polystyrene (PS) and a low density polyethylene (LDPE). For both polymers (t) first followed the linear viscoelasticity rule in that the elongational viscosity, , is three times the steady shear viscosity, 3 _o(t), at low shear rate , but the _E (t) tended to deviate upward after a certain Hencky strain was attained. The birefringence n(t) was a function of both Hencky strain and strain rate in such a way that the stress-optical law holds with the stress-optical coefficient C(t) = n(t)/(t) being equal to the ones reported from shear flow experiments. Interestingly, however, for PS elongated at low strain rates the C(t) vs (t) relation exhibited a strong nonlinearity as soon as (t) reached steady state. This implies that the tensile stress reaches the steady state but the birefringence continues to increase in the low strain-rate elongation. For the PS melt elongated at high strain rates, on the other hand, C(t) was nearly a constant in the entire range observed. For LDPE with long-chain branchings, (t) exhibited tendency of strain-induced hardening after certain critical strain, but C(t) was nearly a constant in the entire range of (t) observed.     

Über die Stabilität viskometrischer Strömungen

Zusammenfassung Die Stabilität der ebenen Couette- und der ebenen Poiseuille-Strömung nicht-newtonscher Fluide wird für kleine Störungen in der viskometrischen Ebene untersucht. Der Einfluß der Relaxationszeit der Störungen wird vernachlässigt. Es wird gezeigt, daß die ebene Couette-Strömung unabhängig von der ReZahl instabil wird, fallsd(N)/d > 4 _>d gilt. Hier bedeuten die Schergeschwindigkeit,N den ersten Normalspannungskoeffizienten, die Viskosität und _d die differentielle Viskosität ( _d =d/d). Das gleiche Kriterium gilt mit den Daten an der Kanalwand auch für die Poiseuille-Strömung. In diesem Fall oszillieren die Eigenfunktionen in einer sehr dünnen, wandnahen Schicht und klingen im Flüssigkeitsinnern sehr rasch ab.

---

Summary The stability of plane Couette and plane Poiseuille flow of a non-Newtonian fluid is investigated for small perturbations in the viscometric plane. The influence of the relaxation time of the perturbations is neglected. It is shown that plane Couette flow will become unstable independently of Reynolds number ifd(N)/d > 4 _d holds. Here are the rate of shear velocity,N the first normal stress coefficient, the viscosity and _d the differential viscosity ( _d =d/d). The same criterion holds also for plane Poiseuille flow with the data taken at the wall. In this case the eigenfunctions are oscillating in a very thin layer near the wall and decaying very rapidly in the inner region of the flow field.

Mit 11 Abbildungen     

Interior Differentiability of Weak Solutions to the Equations of Stationary Motion of a Class of Non-Newtonian Fluids

In this paper, we consider weak solutions to the equations of stationary motion of a class of non-Newtonian fluids the constitutive law of which includes the power law model as special case. We prove the existence of second order derivatives of weak solutions to these equations.     

Correlations between porous medium flow data and turbulent tube flow results for aqueous hydroxypropylguar solutions (Part 2)

Turbulent tube flow and the flow through a porous medium of aqueous hydroxypropylguar (HPG) solutions in concentrations from 100 wppm to 5000 wppm is investigated. Taking the rheological flow curves into account reveals that the effectiveness in turbulent tube flow and the efficiency for the flow through a porous medium both start at the same onset wall shear stress of 1.3 Pa. The similarity of the curves = ( _w ) and = ( _w ), respectively, leads to a simple linear relation / =k, where the constantk or proportionality depends uponc. This offers the possibility to deduce (for turbulent tube flow) from (for flow through a porous medium). In conjunction with rheological data, will reveal whether, and if yes to what extent, drag reduction will take place (even at high concentrations).The relation of our treatment to the model-based Deborah number concept is shown and a scale-up formula for the onset in turbulent tube flow is deduced as well.     

Simulation of viscoelastic stagnation flow

A finite element simulation for the steady flow in a planar stagnation die was used to compute the velocity, pressure and stress fields. It is predicted that a region surrounds the stagnation point where the flow approximates a planar extension. This region is circular for the Newtonian liquid and becomes an ellipse for the Maxwell fluid. An isotropic point in the stress field is found for the Newtonian case as well as for the Maxwell fluid. Lubrication of the die wall, modeled as a finite slip, increases the size of extensional flow region by as much as 100%, while causing a migration of the isotropic stress point towards the die wall.A slight increase in the apparent planar extensional viscosity occurs before the numerical scheme fails at deformation rates well below the extensional singularity in the Maxwell model. Slip at the walls does not significantly alter the convergence behavior, which appears to be limited by effects in the entry region. In this region, at a Weissenberg number of 1.16, spatial oscillations in the pressure, deformation rate and stress develop. The stress normal to the main flow direction in the entry region is compressive for all values of the slip coefficient. a Length of hyperbolic region in the stagnation die - A total area for flow in thexy-plane - C die constant;XY = C/h ² at the die wall - da differential area element - ds differential contour length along the boundary - ;A boundary of the flow - h half-width at die entry - I unit tensor - l _i inlet length in the stagnation die - l _o outlet length in the stagnation die - n unit outward normal to the boundary or to a streamline - n _i unit outward normal to a finite element at the boundary - P pressure; normalized byV/h - P _i nodal value of pressure in an element - Q volumetric flow rate through the die - T total stress tensor; normalized byV/h - t unit tangent to the boundary or to a streamline - t _i unit tangent to a finite element at the boundary - V maximum speed at the centerline of inlet or outlet - v velocity vector with components (u, ); normalized byV - v _i nodal velocities in a finite element - v _s velocity of the solid surface; zero in this study - Ws Weissenberg number;V/h - W width of the die normal to thexy-plane - X position vector with components (X, Y) - slip coefficient; normalized byh/ - rate of deformation tensor; normalized byV/h - elongation rate; dimensionless - relaxation time of main fluid - shear viscosity of main fluid - _ex extensional viscosity (_xx — _yy) - finite element shape function for velocity - finite element shape function for pressure and stress - extra stress; normalized byV/h - _i nodal value of extra stress in an element - partial differentiation - gradient operator     

Hydrodynamical changes due to polymer migration in very dilute solutions

The slip hypothesis, based on thermodynamical arguments, has been extended to obtain the flow characteristics of polymer solutions flowing in a nonhomogeneous flow field. An asymptotic analysis, valid for both channel and falling film flows, is presented that predicts the flow enhancement due to polymer migration. Concentration-viscosity coupling is shown to be a critical factor in the hydrodynamic analysis. The analysis, which essentially provides an upper bound on flow enhancement, explicitly accounts for the influence of wall shear stress, initial polymer concentration etc. A comparison with the pertinent experimental data shows reasonable agreement. c concentration - c ₀ concentration in shear-free region - c _i initial concentration - d rate of deformation tensor - g acceleration due to gravity - g ₁ function defined in eq. [13] or [15] - g ₂ function defined in eq. [18] or [20] - H half-channel thickness or film thickness - K gas law constant - L length of the channel or film - q flow rate per unit width - q ^* normalized flow rate - T temperature - v velocity - V mean velocity - y transverse distance - y _c location of solvent layer - _w /µ _s - _w /c ₀ KT - /t convected derivative - dimensionless cenentration,c/c ₀ - _c dimensionless interface concentration - _w dimensionless wall concentration - relaxation time - µ _eff effective viscosity - µ _s solvent viscosity - dimensionless transverse distance,y/H - _c dimensionless interface location - density - stress tensor - _w wall shear stress - c _i KT/ _w - ns no slip NCL-Communication No. 3155     

Rheological behavior and birefringence investigations on drag-reducing surfactant solutions of tallow-(tris-hydroxyethyl)-ammonium acetate/sodiumsalicylate mixtures

Rheological and flow birefringent properties of a drag-reducing mixture of tallow-(tris-hydroxiethyl)-ammonium acetate (ETHOQUAD T/13-50) and sodiumsalicylate (NaSal) have been studied as a function of the concentration and of the salt/surfactant molar ratio x. The optimum molar ratio x for drag reduction is around 2.5. It is shown that shear-induced supramicellar structures (SIS) which are believed to be responsible for friction reduction in turbulent pipe flow develop in the presence of NaSal. It was observed that SIS are also formed even if the concentration c exceeds c ^*, i.e., the concentration where the volumes of rotation of the individual rodlike micelles start to overlap. The validity of the stress optical law is discussed. A switch from a reptation-controlled stress relaxation to a kinetically controlled mechanism takes place at x 2.5 for this system.     

Newtonian Mechanics on Kähler Manifold

In this paper we discuss Newtonian Mechanics on Kahler Manifold, and also givefoe complex mathematical aspects of Newton’s law, the law of kinetic energy, the lawof kinetic quantity,the equation of motion and the "general equation of dynamics",and so on.     

Deformed particle image pattern matching in Particle Image Velocimetry

The deformation of particle image patterns due to velocity gradients causes errors of velocity measurements and false velocity detections in PIV (Particle Image Velocimetry). A novel technique to overcome those limitations inherent in the conventional PIV by correcting the particle image pattern according to the local velocity gradients in two dimensional flows, i.e. u/x, u/y, v/x and v/y, is proposed and successfully applied to a water flow downstream of a backward facing step.     

A minimal composite theory for stability of non-parallel compressible boundary-layer flow

A new minimal composite theory that extends the approach of Govindarajan and Narasimha [1] is proposed here for 2D non-parallel compressible boundary-layer stability subject to 3D disturbances. The mean profiles are obtained from Hortons analysis, which provides a good approximation for a large range of Prandtl numbers at non-zero pressure gradients. In the lowest order, all effects of order lower than O(R^-2/3) anywhere in the boundary-layer are included, R being the local boundary-layer Reynolds number; the resulting non-parallel formulation yields a set of four ordinary differential equations, as compared to the five coupled equations of classical parallel flow theory of Mack [2]. The largest effect on stability of flow non-parallelism is found to be due to the wall-normal advection of velocity and temperature disturbance quantities by the mean flow. The present theory shows that the bulk viscosity, invariably included in compressible stability theories, is irrelevant at the lowest order. In comparison with the full [O(R^-1)] non-parallel theory, the present theory is marginally better than the parallel flow theory. PACS 03.50.De, 04.20-q, 42.65-k     

The rheology of polymeric liquid crystals

The molecular theory of Doi has been used as a framework to characterize the rheological behavior of polymeric liquid crystals at the low deformation rates for which it was derived, and an appropriate extension for high deformation rates is presented. The essential physics behind the Doi formulation has, however, been retained in its entirety. The resulting four-parameter equation enables prediction of the shearing behavior at low and high deformation rates, of the stress in extensional flows, of the isotropic-anisotropic phase transition and of the molecular orientation. Extensional data over nearly three decades of elongation rate (10^–2–10¹) and shearing data over six decades of shear rate (10^–2–10⁴) have been correlated using this analysis. Experimental data are presented for both homogeneous and inhomogeneous shearing stress fields. For the latter, a 20-fold range of capillary tube diameters has been employed and no effects of system geometry or the inhomogeneity of the flow-field are observed. Such an independence of the rheological properties from these effects does not occur for low molecular weight liquid crystals and this is, perhaps, the first time this has been reported for polymeric lyotropic liquid crystals; the physical basis for this major difference is discussed briefly. A Semi-empirical constant in eq. (18), N/m² - c rod concentration, rods/m³ - c ^* critical rod concentration at which the isotropic phase becomes unstable, rods/m³ - C interaction potential in the Doi theory defined in eq. (3) - d rod diameter, m - D semi-empirical constant in eq. (19), s^–1 - D _r lumped rotational diffusivity defined in eq. (4), s^–1 - rotational diffusivity of rods in a concentrated (liquid crystalline) system, s^–1 - D _ro rotational diffusivity of a dilute solution of rods, s^–1 - f distribution function defining rod orientation - F tensorial term in the Doi theory defined in eq. (7) (or eq. (19)), s^–1 - G tensorial term in the Doi theory defined in eq. (8) - K _B Boltzmann constant, 1.38 × 10^–23 J/K-molecule - L rod length, m - S scalar order parameter - S tensor order parameter defined in eq. (5) - t time, s - T absolute temperature, K - u unit vector describing the orientation of an individual rod - rate of change ofu due to macroscopic flow, s^–1 - v fluid velocity vector, m/s - v velocity gradient tensor defined in eq. (9), s^–1 - V mean field (aligning) potential defined in eq. (2) - x coordinate direction, m - Kronecker delta (= 0 if = 1 if = ) - _r ratio of viscosity of suspension to that of the solvent at the same shear stress - _s solvent viscosity, Pa · s - ^* viscosity at the critical concentrationc ^*, Pa · s - v ₁, v₂ numerical factors in eqs. (3) and (4), respectively - deviatoric stress tensor, N/m² - volume fraction of rods - ⁰ constant in eq. (16) - ^* volume fraction of rods at the critical concentrationc ^* - average over the distribution functionf(u, t) (= d ²u f(u, t)) - gradient operator - d ²u integral over the surface of the sphere (|u| = 1)     

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²     

Description of a Flow of a Gas‐Condensate Mixture in an Axisymmetric Capillary Tube by the Density‐Functional Method

An isothermal flow of a twophase multicomponent mixture through a smalldiameter capillary tube is examined by the densityfunctional method. For low ratios of the characteristic radius of the capillary to its length, a general form of the dominating term in the asymptotic solution is found. An improved version of the law of mixture transfer is obtained. The form of possible corrections to the Darcy law for the filtration rates of the phases is discussed.     

Generalized Forchheimer Equation for Two-Phase Flow Based on Hybrid Mixture Theory

In this paper, we derive a Forchheimer-type equation for two-phase flow through an isotropic porous medium using hybrid mixture theory. Hybrid mixture theory consists of classical mixture theory applied to a multiphase system with volume averaged equations. It applies to media in which the characteristic length of each phase is small relative to the extent of the mixture. The derivation of a Forchheimer equation for single phase flow has been obtained elsewhere. These results are extended to include multiphase swelling materials which have nonnegligible interfacial thermodynamic properties.