Reservoir transport equations by compositional approach

The objective of this paper is to present an overview of the fundamental equations governing transport phenomena in compressible reservoirs. A general mathematical model is presented for important thermo-mechanical processes operative in a reservoir. Such a formulation includes equations governing multiphase fluid (gas-water-hydrocarbon) flow, energy transport, and reservoir skeleton deformation. The model allows phase changes due to gas solubility. Furthermore, Terzaghi's concept of effective stress and stress-strain relations are incorporated into the general model. The functional relations among various model parameters which cause the nonlinearity of the system of equations are explained within the context of reservoir engineering principles. Simplified equations and appropriate boundary conditions have also been presented for various cases. It has been demonstrated that various well-known equations such as Jacob, Terzaghi, Buckley-Leverett, Richards, solute transport, black-oil, and Biot equations are simplifications of the compositional model.Notation List B reservoir thickness - B formation volume factor of phase - Cⁱ mass of component i dissolved per total volume of solution - C ⁱ mass fraction of component i in phase - C heat capacity of phase at constant volume - C_p heat capacity of phase at constant pressure - D ⁱ hydrodynamic dispersion coefficient of component i in phase - D_MTf thermal liquid diffusivity for fluid f - F = F(x, y, z, t) defines the boundary surface - f_p fractional flow of phase - g gravitational acceleration - H_p enthalpy per unit mass of phase - J_p volumetric flux of phase - k_rf relative permeability to fluid f - k₀ absolute permeability of the medium - M_p ⁱ mass of component i in phase - n porosity - N rate of accretion - P_f pressure in fluid f - p_ca capillary pressure between phases and =p-p - Rⁱ rate of mass transfer of component i from phase to phase - Rⁱ _source source rate of component i within phase - S saturation of phase - s gas solubility - T temperature - t time - U displacement vector - u velocity in the x-direction - v velocity in the y-direction - V volume of phase - V_s velocity of soil solids - W_i body force in coordinate direction i - x horizontal coordinate - z vertical coordinate Greek Letters _p volumetric coefficient of compressibility - _T volumetric coefficient of thermal expansion - _ij Kronecker delta - volumetric strain - _m thermal conductivity of the whole matrix - internal energy per unit mass of phase - gf suction head - density of phase - _ij tensor of total stresses - _ij tensor of effective stresses - volumetric content of phase - f viscosity of fluid f     

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.     

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/     

Reflection, transmission and excitation of SH-surface waves by a discontinuity in mass-loading on a semi-infinite elastic medium

The scattering of an SH-wave by a discontinuity in mass-loading on a semi-infinite elastic medium is investigated theoretically. The incident wave is either a plane body wave or a plane SH-surface wave. The problem is reduced to a Wiener-Hopf problem for the scattered wave. In this problem the amplitude spectral density of the particle displacement occurs as unknown function. Special attention is given to the numerical values of the surface wave contributions to the scattered field.Nomenclature x ₁, x ₂, x ₃ Cartesian coordinates - , polar coordinates in x ₁, x ₃-plane - volume mass density - surface mass density of mass-loading - , Lamé constants - U scalar wave function, defined by (2.1) - c _S speed of propagation of uniform shear waves in bulk medium (c _S=(/)^1/2) - angular frequency - t time - k _S wave number of uniform shear waves (k _S=/c _S) - reduced specific acoustic impedance of mass-loading (=k _S /) - k _m wave number of SH-surface wave (k _m=k _S(1+ ²)^1/2) - _1,2,3 partial differentiation with respect to x _1,2,3 - ⁱ angle between x ₃-axis and direction of propagation of incident body wave - ⁱ wave number in horizontal direction of incident body wave ( ⁱ=k _S sin( ⁱ)) - ⁱ wave number in vertical direction of incident body wave ( ⁱ=k _S cos( ⁱ)) - C _1,2 complex amplitude of surface wave excited by a body wave - R reflection factor of surface wave, when surface wave is incident - T transmission factor of surface wave, when surface wave is incident - S particle displacement vector The research presented in this paper has been carried out with partial financial support from the Delfts Hogeschoolfonds.     

Viscous properties of soil

Dynamic problems connected with the wave propagation in soils not saturated with water and with wave interaction with obstacles and structural elements at the present time are solved on the basis of models in which plastic but not viscous soil properties are taken into account [1–5]. An analysis of experimental data and their comparison with the calculated results [4, 5] confirms that it is permissible to apply the model of an elasticplastic medium to soils in problems concerning the interaction of waves and structures. At the same time plane-wave damping in soils takes place more intensively than would follow from calculations carried out on the basis of models of an elastic-plastic medium. For example, if in a section of a poured sandy soil, taken as the initial section, the maximum stress in the wave is _m=ll kgf/cm² and its duration is 6=8 msec, then at a distance of 25 cm the calculations give _m=9.5 kgf/cm², while the experiment gives _m= 5 kgf/cm². If in the initial section _m= 20 kgf/cm² and =6 msec, then at a distance of 35 cm the calculation gives _m= l7 kgf/cm², while the experiment gives _m= 9 kgf/cm². In the calculations it was assumed that unloading takes place with a constant strain. This deviation of the calculated results from the experiment can be explained, in the first place, by the dependence of the () on the strain rate , which is not taken into account in the model of an elastic-plastic medium. The viscous properties cause additional energy losses and a more intensive damping of the waves. Experimentally the dependence of the () curves on the strain rate has been investigated for many soils [5–8]. The dynamic load on the test sample was produced by a body falling from a height or being accelerated by some method. Below we present test results of viscous soil properties when the test sample is compressed by an air shock wave. Compression curves and approximate numerical values of the coefficient of viscosity are obtained.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, Vol. 9, No. 4, pp. 68–71, July–August, 1968.The author thanks A. I. Shishikin for his participation in the experiments.     

Unsteady convective diffusion in a rotating parallel plate channel

The paper presents an exact analysis of the dispersion of a passive contaminant in a viscous fluid flowing in a parallel plate channel driven by a uniform pressure gradient. The channel rotates about an axis perpendicular to its walls with a uniform angular velocity resulting in a secondary flow. Using a generalized dispersion model which is valid for all time, we evaluate the longitudinal dispersion coefficientsK _i (i=1, 2, ...) as functions of time. It is shown thatK ₁=0 andK ₃,K ₄, ... decay rapidly in comparison withK ₂. ButK ₂ decreases with increasing (the dimensionless rotation parameter) for values of upto approximately =2.2. ThereafterK ₂ increases with further increase in and its value gets saturated for large values of (say, 500) and does not change any further with increase in . A physical explanation of this anomalous behaviour ofK ₂ is given.

---

Instationäre konvektive Diffusion in einem rotierenden Parallelplattenkanal

Zusammenfassung In dieser Untersuchung wird eine exakte Analyse der Ausbreitung eines passiven Kontaminierungsstoffes in einer zähen Flüssigkeit gegeben, die, befördert durch einen gleichförmigen Druckgradienten, in einem Parallelplattenkanal strömt. Der Kanal rotiert mit gleichförmiger Winkelgeschwindigkeit um eine zu seinen Wänden senkrechte Achse, wodurch sich eine Sekundärströmung ausbildet. Unter Verwendung eines generalisierten, für alle Zeiten gültigen Dispersionsmodells werden die longitudinalen DispersionskoeffizientenK _i (i=1, 2, ...) als Funktionen der Zeit ermittelt. Es wird gezeigt, daßK ₁=0 gilt und dieK ₃,K ₄, ... gegenüberK ₂ schnell abnehmen.K ₂ nimmt ab, wenn , der dimensionslose Rotationsparameter, bis etwa zum Wert 2,2 ansteigt. Danach wächstK ₂ mit bis auf einem Endwert an, der etwa ab =500 erreicht wird. Dieses anomale Verhalten vonK ₂ findet eine physikalische Erklärung.

---

List of symbols C solute concentration - D molecular diffusivity - K _i longitudinal dispersion coefficients - 2L depth of the channel - P ₀ dimensionless pressure gradient along main flow - Pe Péclet number - q velocity vector - Q _x,Q _y mass flux along the main flow and the secondary flow directions - dimensionless average velocity along the main flow direction - (x, y, z) Cartesian co-ordinates Greek symbols dimensionless rotation parameter - the inclination of side walls withx-axis - kinematic viscosity - fluid density - dimensionless time - angular velocity of the channel - dimensionless distance along the main flow direction - dimensionless distance along the vertical direction - dimensionless solute concentration - integral of the dispersion coefficientK ₂() over a time interval     

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     

Turbulent heat and mass transfer for Pr ≫ 1 and law of turbulent diffusion damping at a solid wall

The question of determining the law of damping for the turbulent diffusion coefficient at a smooth wall according to data on mass and heat transfer for Pr 1 is discussed. It is proved that the hypothesis that this law is determined by the first member of the Taylor series expansion of , namely, / = yⁿ ₊ is valid in the Pr range from 10³ to 10⁵ only under the assumption that the subsequent terms in the expansion have smaller coefficients. A statistical analysis of electrochemical and other experiments devoted to this problem shows that apparently n = 3, but singularities in the experimental results do not permit making a final conclusion. Requirements on a conclusive experiment are formulated on the basis of the analysis made.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 172–175, March–April, 1977.     

Chaos and integrability in a nonlinear wave equation

We consider the parametrized family of equations _tt ,u- _xx u-au+u ₂ ² u=O,x(0,L), with Dirichlet boundary conditions. This equation has finite-dimensional invariant manifolds of solutions. Studying the reduced equation to a four-dimensional manifold, we prove the existence of transversal homoclinic orbits to periodic solutions and of invariant sets with chaotic dynamics, provided that =2, 3, 4,.... For =1 we prove the existence of infinitely many first integrals pairwise in involution.     

Ü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     

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     

A noninvasive method for the measurement of flow-induced surface displacement of a compliant surface

A noninvasive optical method is described which allows the measurement of the vertical component of the instantaneous displacement of a surface at one or more points. The method has been used to study the motion of a passive compliant layer responding to the random forcing of a fully developed turbulent boundary layer. However, in principle, the measurement technique described here can be used equally well with any surface capable of scattering light and to which optical access can be gained. The technique relies on the use of electro-optic position-sensitive detectors; this type of transducer produces changes in current which are linearly proportional to the displacement of a spot of light imaged onto the active area of the detector. The system can resolve displacements as small as 2 m for a point 1.8 mm in diameter; the final output signal of the system is found to be linear for displacements up to 200 m, and the overall frequency response is from DC to greater than 1 kHz. As an example of the use of the system, results detailing measurements obtained at both one and two points simultaneously are presented.List of symbols C _t elastic transverse wave speed = (G/)^1/2 - d ⁺ spot diameter normalized by viscous length scale - G frequency average of G() - G() shear storage modulus - G() shear loss modulus - l. viscous length scale = v/u _* - N total number of sampled data values - r separation vector for 2-point measurements = (, ) - rms root-mean-square value - R momentum thickness Reynolds number = U _t8/v - t time - u (y) mean streamwise component of velocity in boundary layer - u _* friction velocity = (t _w/)^1/2 - U free-stream velocity - x, y, z longitudinal, normal and spanwise directions - y _o undisturbed surface position - vertical component of compliant surface displacement - ₉₉ boundary layer thickness for which u(y) = 0.99 U _t8 - _l viscous sublayer thickness 5 l _* - frequency average of G()/ - boundary layer momentum thicknes = - fluid dynamic viscosity - v fluid kinematic viscosity = / - , longitudinal, spanwise components of separation vector r - fluid density - time delay - _w wall shear stress     

Propagation of a spherical wave in nonlinearly compressible and elastoplastic materials

The problem of spherical wave propagation in soil under the action of an intense uniformly decreasing load ₀(t) applied to the boundary of a cavity with radius r₀ is considered. Soil with a high stress level is modeled either by ideally nonlinearly compressible or elastoplastic material, taking account of linear irreversible unloading for the material. In contrast to [1–7], in order to describe material movement use is made of strain theory [8] with determining functions = (), _i=_i(_i), where , _i, , _i are the first and second invariants of strain and stress tensors. During material loading these functions are presented in the form of polynomials ()=(_i+₂¦¦), _ii)=(_i-_2i)_i, in which constant coefficients _i, _i=1, 2) are determined by experiment, taking account of the triaxial stressed state of soil. Solution of the problem is constructed by an analytically reversible method, with prescribed shape for the shock-wave (SW) surface in the form of a second-degree polynomial relating to time t and a numerical method of characteristics for a prescribed arbitrarily decreasing load _i(t). On the basis of the analytical equations obtained, calculations are carried out for material parameters (including loading profile) in a computer and stresses and mass velocity of plastic and elastoplastic materials are compared.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 95–100, July–August, 1986.The authors express their sincere thanks to Kh. A. Rakhmatulin for discussing the results of this work.     

Some remarks on differential equations of quadratic type

In this paper we study differential equations of the formx(t) + x(t)=f(x(t)), x(0)=x ₀ C HereC is a closed, bounded convex subset of a Banach spaceX,f(C) C, and it is often assumed thatf(x) is a quadratic map. We study the differential equation by using the general theory of nonexpansive maps and nonexpansive, non-linear semigroups, and we obtain sharp results in a number of cases of interest. We give a formula for the Lipschitz constant off: C C, and we derive a precise explicit formula for the Lipschitz constant whenf is quadratic,C is the unit simplex inR ⁿ, and thel ₁ norm is used. We give a new proof of a theorem about nonexpansive semigroups; and we show that if the Lipschitz constant off: CC is less than or equal to one, then lim_tf(x(t))–x(t)=0 and, if {x(t):t 0} is precompact, then lim_tx(t) exists. Iff¦C=L¦C, whereL is a bounded linear operator, we apply the nonlinear theory to prove that (under mild further conditions on C) lim_t f(x(t))–x(t)=0 and that lim_t x(t) exists if {x(t):t 0} is precompact. However, forn 3 we give examples of quadratic mapsf of the unit simplex ofR ⁿ into itself such that lim_t x(t) fails to exist for mostx ₀ C andx(t) may be periodic. Our theorems answer several questions recently raised by J. Herod in connection with so-called model Boltzmann equations.     

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.     

Necessary and sufficient conditions for isotropic rank-one convex functions in dimension 2

We provide some new necessary and sufficient conditions for regular isotropic rank-one convex functions on M ₂ ⁺={2×2 matrices such that det M0}. It is well known that isotropic functions W (M) can be written as W (M)=G(₁, ₂) where _i are the singular values of M. One of these conditions allows us to understand better the gap between the rank-one convexity and the quasiconvexity.     

Mixed convection on vertical plate for fluids of any prandtl number

A mixed convection parameter=(Ra) ^1/4/(Re)^1/2, with=Pr/(1+Pr) and=Pr/(1 +Pr)^1/2, is proposed to replace the conventional Richardson number, Gr/Re², for combined forced and free convection flow on an isothermal vertical plate. This parameter can readily be reduced to the controlling parameters for the relative importance of the forced and the free convection,Ra ^1/4/(Re ^1/2 Pr ^1/3) forPr 1, and (RaPr)^1/2/(RePr ^1/2 forPr 1. Furthermore, new coordinates and dependent variables are properly defined in terms of, so that the transformed nonsimilar boundary-layer equations give numerical solutions that are uniformly valid over the entire range of mixed convection intensity from forced convection limit to free convection limit for fluids of any Prandtl number from 0.001 to 10,000. The effects of mixed convection intensity and the Prandtl number on the velocity profiles, the temperature profiles, the wall friction, and the heat transfer rate are illustrated for both cases of buoyancy assisting and opposing flow conditions.

---

Mischkonvektion an einer vertikalen Platte für Fluide beliebiger Prandtl-Zahl

Zusammenfassung Für die kombinierte Zwangs- und freie Konvektion an einer isothermen senkrechten Platte wird ein Mischkonvektions-Parameter=( Ra) ^1/4 (Re)^1/2, mit=Pr/(1 +Pr) und=Pr/(1 +Pr)^1/2 vorgeschlagen, den die gebräuchliche Richardson-Zahl, Gr/Re², ersetzen soll. Dieser Parameter kann ohne weiteres auf die maßgebenden Kennzahlen für den relativen Einfluß der erzwungenen und der freien Konvektion reduziert werden,Ra ^1/4/(Re ^1/2 Pr ^1/3) fürPr 1 und (RaPr)^1/4/(RePr)^1/2 fürPr 1. Weiterhin werden neue Koordinaten und abhängige Variablen als Funktion von definiert, so daß für die transformierten Grenzschichtgleichungen numerische Lösungen erstellt werden können, die über den gesamten Bereich der Mischkonvektion, von der freien Konvektion bis zur Zwangskonvektion, für Fluide jeglicher Prandtl-Zahl von 0.001 bis 10.000 gleichmäßig gültig sind. Der Einfluß der Intensität der Mischkonvektion und der Prandtl-Zahl auf die Geschwindigkeitsprofile, die Temperaturprofile, die Wandreibung und den Wärmeübergangskoeffizienten werden für die beiden Fälle der Strömung in und entgegengesetzt zur Schwerkraftrichtung dargestellt.

---

Nomenclature C _f local friction coefficient - C _p specific heat capacity - f reduced stream function - g gravitational acceleration - Gr local Grashoff number,g T _w –T )x³/v² - Nu local Nusselt number - Pr Prandtl number,v/ - Ra local Rayleigh number,g T _w –T x ³/( v) - Re local Reynolds number,u x/v - Ri Richardson number,Gr/Re ² - T fluid temperature - T _w wall temperature - T free stream temperature - u velocity component in thex direction - u free stream velocity - v velocity component in they direction - x vertical coordinate measuring from the leading edge - y horizontal coordinate Greek symbols thermal diffusivity - thermal expansion coefficient - mixed convection parameter (Ra)1/4/Re)^1/2 - pseudo-similarity variable,(y/x) - ₀ conventional similarity variable,(y/x)Re ^1/2 - dimensionless temperature, (T–T T _W –T - unified mixed-flow parameter, [(Re) ^1/2 + (Ra)^1/4] - dynamic viscosity - kinematic viscosity - stretched streamwise coordinate or mixed convection parameter, [1 + (Re)^1/2/(Ra) ^1/4]^–1=/(1 +) - density - Pr/(1 + Pr) _w wall shear stress - stream function - Pr/(l+Pr)^1/3 This research was supported by a grand from the National Science Council of ROC     

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.     

The dynamics of solid-solid phase transitions 2. Incoherent interfaces

Incoherent phase transitions are more difficult to treat than their coherent counterparts. The interface, which appears as a single surface in the deformed configuration, is represented in its undeformed state by a separate surface in each phase. This leads to a rich but detailed kinematics, one in which defects such as vacancies and dislocations are generated by the moving interface. In this paper we develop a complete theory of incoherent phase transitions in the presence of deformation and mass transport, with phase interface structured by energy and stress. The final results are a complete set of interface conditions for an evolving incoherent interface.Frequently used symbols A_i,C_i generic subsurface of S_t - B_i undeformed phase-i region - C configurational bulk stress, Eshelby tensor - F deformation gradient - G inverse deformation gradient - H relative deformation gradient - J bulk Jacobian of the deformation - ¯K, K_i total (twice the mean) curvature of and S_i - Lin (U, V) linear transformations from U into V - Lin⁺ linear transformations of ³ with positive determinant - Orth⁺ rotations of ³ - Q^a external bulk mass supply of species a - ¯S bulk Cauchy stress tensor - S bulk Piola-Kirchhoff stress tensor - S_i undeformed phase i interface - U_i relative velocity of S_i - Unim⁺ linear transformations of ³ with unit determinant - ¯V, V_i normal velocity of and S_i - intrinsic edge velocity of S and A _i S - W_i volume flow across the phase-i interface - X material point - b external body force - e internal bulk configurational force - f_i external interfacial force (configurational) - ¯g external interfacial force (deformational) - grad, div spatial gradient and divergence - gradient and divergence on - h relative deformation - h^a, diffusive mass flux of species a and list of mass fluxes - ¯m outward unit normal to a spatial control volume - ¯n, n_i unit normal to and S_i - n subspace of ³ orthogonal to n - ¯q^a external interfacial mass supply of species a - s ......... - ¯v, v_i compatible velocity fields of and S_i - ¯w, w_i compatible edge velocity fields for and A_i - x spatial point - y_i deformation or motion of phase i - y^. material velocity - generic subsurfaces of - , _i deformed body and deformed phase-i region - () energy supplied to by mass transport - symmetry group of the lattice - _i, surface jacobians - lattice - () power expended on - spatial control volume - S deformed phase interface - lattice point density - interfacial power density - , A total surface stress - C configurational surface stress for phase 1 (material) - ¯C_i configurational surface stress (spatial) - F_i tangential deformation gradient - G_i inverse tangential deformation gradient - H incoherency tensor - ¯1(x), 1_i(X) inclusions of ¯n(x) and n _i (X) into ³ - K configurational surface stress for phase 2 (material) - ¯L, l_i curvature tensor of and S_i - ¯P(x), P_i(X) projections of ³ onto ¯n(x) and n_i (X) - ¯S, S deformational surface stress (spatial and material) - ¯a, a normal part of total surface stress - c normal part of configurational surface stress for phase 1 (material) - e_i internal interfacial configurational force - ¯v, v_i unit normal to and A _i - (x),_i(X) projections of ³ onto ¯n(x) and n _i (X) - _i normal internal force (material) - bulk free energy - slip velocity - _i=(–1)ⁱ _i ......... - ^a, chemical potential of species a and list of potentials - ^a, bulk molar density of species a and list of molar densities - _i normal internal force (spatial) - surface tension - , _i effective shear - referential-to-spatial transform of field - interfacial energy - grand canonical potential - l unit tensor in ³ - x, vector and tensor product in ³ - (...)^., _t(...) material and spatial time derivative - , Div material gradient and divergence - gradient and divergence on S_i - (...), (...) normal time derivative following and S_i - (...) limit of a bulk field asx ,x_i - [...],...> jump and average of a bulk field across the interface - (...)_ext extension of a surface tensor to ³ - tangential part of a vector (tensor) on and S_i     

Experimental and numerical analysis of oscillatory tube flow of viscoelastic fluids represented at the example of human blood

Summary The 4 constant parameters of an Oldroyd type constitutive equation for normal human blood of 45% haematocrit were determined by means of the steady flow curve and the material functions () and (), measured at 2 Hz in an oscillatory capillary viscometer. It was found, as other authors did before, that the phenomenological behaviour of blood flow can be reproduced qualitatively thereby. The quantitative behaviour, however, cannot be described by thus developed parameters. The parameters of the constitutive equationµ ₀, ₁ and ₂ were therefore generalized to become dependent of shear rate and frequency respectively.In itself this is nothing but a transformation of the material functions ( ),() and (), but these can be used as parameters in a constitutive equation though having lost the property of constancy.In this way the linear region of viscoelasticity and the steady flow curve can be reproduced quantitatively. A computer simulation of oscillatory flow for large amplitudes shows another tendency for the phase shift between pressure and flow than the experiment in the oscillatory capillary rheometer does.The applicability of a constitutive equation modified in this manner for other than oscillatory flow should be further examined especially for pulsatile flow.

---

Zusammenfassung Mit Hilfe der Fließkurve für die stationäre Strömung und Messungen der Materialfunktion () im Oszillations-Kapillarrheometer bei 2 Hz wurden die Konstanten der 4-Konstanten-Oldroyd-Stoffgleichung für gesundes Humanblut von 45% Hämatokrit bestimmt. Es zeigte sich, daß sich — wie auch schon von anderen Autoren mitgeteilt wurde — mit den so ermittelten Modellkonstanten die Phänomenologie des Fließverhaltens des Blutes qualitativ gut beschreiben läßt. Zur quantitativen Beschreibung reicht dieses Modell jedoch nicht aus, wie man an der Wiedergabe der stationären Fließkurve und der Materialfunktionen der linearen Viskoelastizität erkennt. Aus diesem Grund wurden die Parameterµ ₀, ₁ und ₂ in Abhängigkeit der Schergeschwindigkeit bzw. der Frequenz angesetzt.Dies bedeutet zunächst nichts anderes als eine Transformation der Materialfunktionen, jedoch in einer Art, daß letztere als Parameter in einer Stoffgleichung verwendet werden können, was allerdings mit dem Verlust der Konstanz der Parameter verbunden ist. Mit einer derart modifizierten Stoffgleichung lassen sich der Bereich der linearen Viskoelastizität und die stationäre Fließkurve quantitativ beschreiben. Eine Computersimulation der oszillierenden Rohrströmung zeigt für große Amplituden eine andere Tendenz für die Phasenverschiebung zwischen Druck und Volumenstrom, als sie sich bei Messungen am Oszillations-Kapillarrheometer ergibt.Die Anwendbarkeit der modifizierten Stoffgleichung für andere Strömungsformen, wie z.B. pulsierende Rohrströmungen, muß noch geprüft werden.

---

---

Paper, presented at the Annual Conference of the Deutsche Rheologische Gesellschaft in Berlin, May 8–10, 1978.

---

With 8 figures