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.     

On the localness of the spectral energy transfer in turbulence

By utilizing available experimental data for net energy transfer spectra for homogeneous turbulence, contributions P(, ) to the energy transfer at a wavenumber from various other wavenumbers are calculated. This is done by fitting a truncated power-exponential series in and to the experimental data for the net energy transfer T(), and using known properties of P(, ). Although the contributions P(, ) obtained by using this procedure are not unique, the results obtained by using various assumptions do not differ significantly. It seems clear from the results that for a region where the energy entering a wavenumber band dominates that leaving, much of the energy entering the band comes from wavenumbers which are about an order of magnitude smaller. That is, the energy transfer is rather nonlocal. This result is not significantly dependent on Reynolds number (for turbulence Reynolds numbers based on microscale from 3 to 800). For lower wavenumbers, where more energy leaves than enters a wavenumber band, the energy transfer into the band is more local, but much of the energy then leaves at distant wavenumbers.     

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)     

A Capillary Microstructure of the Wetting Front

This article reports the experimental results of a study of the wetting-front microscale structure formed only by capillary forces in homogeneous and random etched glass capillary models. In the homogeneous model, water propagates through the capillary system, evenly filling the capillaries across the direction of flow. Air is trapped by the pinch-off mechanism inside the pore bodies in the form of individual bubbles. The experiments specified three consecutive steps of the pinch-off mechanism, film flow, snap-off, and interface movement. In the random model, both the bypass and pinch-off, forming bypass/cut-off mechanism, create residual air structure. Bypass traps air inside large capillary-pore aggregates which are bounded by small-diameter capillaries in where pinch-off traps air in the adjacent pores. An analysis of the residual air distribution versus depth below the surface in the homogeneous and random micromodels made it possible to identify three successive zones, namely a transition zone, a transmission zone, and a wetting-and-front zone. In the transition zone, the residual air content increases with depth from zero to the constant value in the transmission zone where it remains practically constant. The capillary processes within the wetting-and-front combined zone govern air replacement with wetting and formation of the transmission zone.     

Nonlinear viscoelastic properties and change in entanglement structure of linear polymer

Simultaneous measurements of stress relaxation and differential dynamic modulus were made at 268 K over a time scale of 10 to 10⁴⁵ s for nearly monodisperse polybutadiene (M _w =2.2x10⁵, 1,2-structure 70%, M _e =3600) and also one having coarse cross-linking (M _c =29000). Static shear strain ranged from 0.1 to 2.0. In a long-time region (t> _k ), the relaxation modulus G (; t) could be expressed by the product G (0; t) h (y). The observed h() agreed well with the Doi-Edwards theory without use of IA approximation. Both the cured and uncured samples showed initial drop of the differential storage modulus G (), ; t) followed by gradual recovery, but did not attain the value before shearing G (, ; t) for the uncured sample showed smaller values than that for the cured one in the whole measured time scale at the higher strain, confirming the two origins of nonlinear viscoelasticity of well entangled polymer; induced chain anisotropy and induced decrement in entanglement density. G (, ; t) curves for the cured sample agreed well with the BKZ predictions. But the curves for the uncured sample agreed well with the BKZ prediction only at the time scale of t< _k . BKZ prediction showed significant upward deviations at t> _k . Such the differences are discussed in terms of the two origins.Dedicated to Prof. John D. Ferry on the occasion of his 85th birthday.     

A New Optimal Growth Constant for the Hele–Shaw Instability

We study the immiscible displacement of the oil from a homogeneous porous medium by using a less viscous fluid (water). We use the Hele–Shaw model, then a sharp interface exists between the fluids. The fingering phenomenon appears, first studied by Saffman and Taylor (1958). Gorell and Homsy (1983) consider an intermediate region (I. R.) between water and oil, containing a polymer mixture. The unknown viscosity in I. R. is a parameter which can improve the stability of the I. R.–oil interface. A numerical optimal viscosity profile in I. R. is given. Carasso and Paa (1998) obtain an explicit formula for an optimal viscosity profile in I. R. An upper estimation of the growth constant is given. In this paper, a very slow viscosity profile in I. R. is defined and an optimal formula for the growth constant is obtained, less than the previous estimation of Carasso and Paa. Moreover, this formula is similar with the Saffman–Taylor result, only the water viscosity is replaced by the limit value of viscosity in I. R. on the interface with the oil. We explain the apparent contradiction between the previous results of Gorell and Homsy (1983) and Paa and Polisevski (1992).     

Solution of regular beam equations in arbitrary emission conditions on a curvilinear surface

The regular beam equations are solved analytically for the case of emission from an arbitrary surface in conditions of total space charge (-mode) and in a given external magnetic field H (§2) for temperature-limited emission (T-mode), in an external magnetic field H (§3); and for emission with nonzero initial velocity (§4). The emitter is taken as the coordinate surface x¹=0 in an orthogonal system x¹ (i = =1,2,3), while the current density J and field on it are given functions j(x², x³), (x², x³. The solution is written as series in (x¹) with coefficients dependent on x², x³, determined from recurrence relations. For emission in the -mode and H 0, =1/3; for temperature-limited emission, =1/2; with nonzero initial velocity, =1. The results are extended to the case of a beam in the presence of a moving background of uniform density (5).     

Ši’lnikov Chaos in the Generalized Lorenz Canonical Form of Dynamical Systems

This paper studies the generalized Lorenz canonical form of dynamical systems introduced by elikovský and Chen [International Journal of Bifurcation and Chaos 12(8), 2002, 1789]. It proves the existence of a heteroclinic orbit of the canonical form and the convergence of the corresponding series expansion. The ilnikov criterion along with some technical conditions guarantee that the canonical form has Smale horseshoes and horseshoe chaos. As a consequence, it also proves that both the classical Lorenz system and the Chen system have ilnikov chaos. When the system is changed into another ordinary differential equation through a nonsingular one-parameter linear transformation, the exact range of existence of ilnikov chaos with respect to the parameter can be specified. Numerical simulation verifies the theoretical results and analysis.     

ON THE UNIFICATION OF THE HAMILTON PRINCIPLES IN NONHOLONOMIC SYSTEM AND IN HOLONOMIC SYSTEM

ＯＮＴＨＥＵＮＩＦＩＣＡＴＩＯＮＯＦＴＨＥＨＡＭＩＬＴＯＮＰＲＩＮＣＩＰＬＥＳＩＮＮＯＮＨＯＬＯＮＯＭＩＣＳＹＳＴＥＭＡＮＤＩＮＨＯＬＯＮＯＭＩＣＳＹＳＴＥＭ（梁立孚）（韦扬）ＯＮＴＨＥＵＮＩＦＩＣＡＴＩＯＮＯＦＴＨＥＨＡＭＩＬＴＯＮＰＲＩＮＣＩＰＬ...     

Transcendentally small transversality in the rapidly forced pendulum

The rapidly forced pendulum equation with forcing sin((t/), where =<₀^p,p = 5, for ₀, sufficiently small, is considered. We prove that stable and unstable manifolds split and that the splitting distanced(t) in the ( ,t) plane satisfiesd(t) = sin(t/) sech(/2) +O( ₀ exp(–/2)) (2.3a) and the angle of transversal intersection,, in thet = 0 section satisfies 2 tan/2 = 2S _s = (/2) sech(/2) +O(( ₀ /) exp(–/2)) (2.3b) It follows that the Melnikov term correctly predicts the exponentially small splitting and angle of transversality. Our method improves a previous result of Holmes, Marsden, and Scheuerle. Our proof is elementary and self-contained, includes a stable manifold theorem, and emphasizes the phase space geometry.     

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.     

Flow in porous media III: Deformable media

Stokes flow in a deformable medium is considered in terms of an isotropic, linearly elastic solid matrix. The analysis is restricted to steady forms of the momentum equations and small deformation of the solid phase. Darcy's law can be used to determine the motion of the fluid phase; however, the determination of the Darcy's law permeability tensor represents part of the closure problem in which the position of the fluid-solid interface must be determined.Roman Letters A interfacial area of the- interface 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 macroscopic system, m² - A ^* interfacial area of the- interface contained within a unit cell, m² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - E Young's modulus for the-phase, N/m² - e i unit base vectors (i = 1, 2, 3) - g gravity vector, m²/s - H height of elastic, porous bed, m - k unit base vector (=e ₃) - characteristic length scale for the-phase, m - L characteristic length scale for volume-averaged quantities, m - n unit normal vector pointing from the-phase toward the-phase (n = -n ) - p pressure in the-phase, N/m² - P p – g·r, N/m² - r ₀ radius of the averaging volume, m - r position vector, m - t time, s - T total stress tensor in the-phase, N/m² - T ⁰ hydrostatic stress tensor for the-phase, N/m² - u displacement vector for the-phase, m - V averaging volume, m₃ - V volume of the-phase contained within the averaging volume, m³ - v velocity vector for the-phase, m/s Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - shear coefficient of viscosity for the-phase, Nt/m² - first Lamé coefficient for the-phase, N/m² - second Lamé coefficient for the-phase, N/m² - bulk coefficient of viscosity for the-phase, Nt/m² - T –T ⁰ , a deviatoric stress tensor for the-phase, N/m²     

The Gibbs-Thompson relation within the gradient theory of phase transitions

This paper discusses the asymptotic behavior as 0+ of the chemical potentials associated with solutions of variational problems within the Van der Waals-Cahn-Hilliard theory of phase transitions in a fluid with free energy, per unit volume, given by ²¦¦²+ W(), where is the density. The main result is that is asymptotically equal to E/d+o(), with E the interfacial energy, per unit surface area, of the interface between phases, the (constant) sum of principal curvatures of the interface, and d the density jump across the interface. This result is in agreement with a formula conjectured by M. Gurtin and corresponds to the Gibbs-Thompson relation for surface tension, proved by G. Caginalp within the context of the phase field model of free boundaries arising from phase transitions.     

The behaviour of a mass-spring system provided with a discontinuous dynamic vibration absorber

Summary The influence of a specific discontinuous dynamic vibration absorber on the motion of a vibrating system is investigated. Attention is paid to the effectiveness in the case of free vibrations and of vibrations due to sinusoidal excitations, where structural damping is also taken into account. For certain configurations numerical results are given.Notation M mass of the vibrating system - m mass of the ball - c spring constant - k structural damping coefficient - g acceleration of gravity - F ₀ amplitude of external force - frequency of external force - 2 distance over which the ball may move - x displacement of centre of tube - y displacement of the ball x and y are measured positive downward from the equilibrium position of the system when no gravitational force would act     

Theoretical derivation of tangential velocity profiles in a flat vortex chamber-influence of turbulence and wall friction

Summary As part of a study on the hydrodynamics of a cyclone separator, a theoretical investigation of the flow pattern in a flat box cyclone (vortex chamber) has been carried out. Expressions have been derived for the tangential velocity profile as influenced by internal friction (eddy viscosity) and wall friction. The most important parameter controlling the tangential velocity profile is = –u ₀ R/(v+ ), where u ₀ is the radial velocity at the outer radius R of the cyclone, the kinematic liquid viscosity and is the kinematic eddy viscosity. For values of greater than about 10 the tangential velocity profile is nearly hyperbolic, for smaller than 1 the tangential velocity even decreases towards the centre. It is shown how and also the wall friction coefficient may be obtained from experimental velocity profiles with the aid of suitable graphs. Because of the close relation between eddy viscosity and eddy diffusion, measurements of velocity profiles in flat box cyclones will also provide information on the eddy motion of particles in a cyclone, a motion reducing its separation efficiency.List of symbols A cross-sectional area of cyclone inlet - h height of cyclone - p static pressure in cyclone - p static pressure difference in cyclone between two points on different radius - r radius in cyclone - r ₁ radius of cyclone outlet - R radius of cyclone circumference - u radial velocity in cyclone - u ₀ radial velocity at circumference of flat box cyclone - v tangential velocity - v ₀ tangential velocity at circumference of flat box cyclone - w axial velocity - z axial co-ordinate in cyclone - friction coefficient in flat box cyclone (for definition see § 5) - ₁ value of friction coefficient for ₁<< ₂ - ₂ value of friction coefficient for ₂<<1 - = - ₁ value of for ₁<< ₂ - ₂ value of for ₂<<1 - thickness of laminar boundary layer - =/h - turbulent kinematic viscosity - ratio of z to h - k ratio of height of cyclone to radius R of cyclone - parameter describing velocity profile in cyclone =–u ₀ R/(+) - kinematic viscosity of fluid - density of fluid - ratio of r to R - ₁ value of at outlet of cyclone - ₂ value of at inner radius of cyclone inlet - _w shear stress at cyclone wall - angular momentum in cyclone/angular momentum in cyclone inlet - ₁ value of at = ₁ - ₂ value of at = ₂     

Über singuläre Lastfälle in einer linearen Schalentheorie und ihre finite Behandlung

Übersicht Ausgehend von bekannten Fundamentallösungen für Platten bzw. Scheiben wird die Erstellung singulärer Ansatzfunktionen gezeigt, wie sie für finite Näherungsverfahren benötigt werden, die von den Funktionalen der totalen Energie bzw. der komplementären Energie ausgehen. Das Vorgehen wird eingehend an Kreiszylinderschalen erläutert.

---

Summary Starting from known fundamental solutions of plates the construction of singular basic estimate functions is shown. These are necessary in finite approximation methods basing on the functionals of total energy or complementary energy. The proceeding is explained in detail in a cylindrical shell analysis.

     

Der Spannungs- und Verschiebungszustand drehsymmetrischer Membrane mit beliebig gekrümmter Meridiankurve

Zusammenfassung Die Einführung von Zylinderkoordinaten (x, r, ) in die Gleichgewichtsbedingungen der Schnittkräfte bzw. in die Beziehungen zwischen Verzerrung und Verschiebungen am differentialen Schalenabschnitt ermöglicht die Berechnung des Spannungs- und Verschiebungszustandes von drehsymmetrischen Membranen mit beliebig gekrümmter Meridiankurve auf die Integration einer einfachen, linearen partiellen Differentialgleichung zweiter Ordnung für eine charakteristische FunktionF bzw. zurückzuführen. Eine geschlossene Lösung und damit eine Darstellung der Schnittkräfte und Verschiebungen durch explizite Formeln ist bei harmonischer Belastung cosn für zwei Funktionsgruppen=x ² und=x ^–3 möglich. Im Sonderfall der drehsymmetrischen und der antimetrischen Belastung mitn=0 undn=1 gelten die Gleichungen der Schnitt- und Verschiebungsgrößen für eine beliebige Meridianfunktion=(). Die Betrachtungen der Randbedingungen offener Schalen bei harmonischer Belastung geben über die infinitesimalen Deformationen einer drehsymmetrischen Membran mit überall negativer Krümmung Aufschluß.     

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     

Exact Solutions of the Hydrodynamic Equations Derived from Partially Invariant Solutions

The paper proposes a heuristic approach to constructing exact solutions of the hydrodynamic equations based on the specificity of these equations. A number of systems of hydrodynamic equations possess the following structure: they contain a reduced system of n equations and an additional equation for an extra function w. In this case, the reduced system, in which w = 0, admits a Lie group G. Taking a certain partially invariant solution of the reduced system with respect to this group as a seed:rdquo; solution, we can find a solution of the entire system, in which the functional dependence of the invariant part of the seed solution on the invariants of the group G has the previous form. Implementation of the algorithm proposed is exemplified by constructing new exact solutions of the equations of rotationally symmetric motion of an ideal incompressible liquid and the equations of concentrational convection in a plane boundary layer and thermal convection in a rotating layer of a viscous liquid.