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.     

On the measurement of lateral velocity derivatives in turbulent flows

Direct numerical simulation data for the lateral velocity derivative u/y at the centreline of a fully developed turbulent channel flow provide reasonable support for Wyngaard's analysis of the error involved in measuring this quantity using parallel hot wires. Numerical data in the wall region of the channel flow also provide a useful indication of how to select the separation between the wires. Justification for this choice is obtained by comparing several measured statistics of u/y with the corresponding numerical data.     

Limitation and improvement of PIV

In this second part of the paper, the Particle Image Distortion (PID) technique is described. It is proposed to overcome the limitations of conventional PIV due to the local deformations u/x, u/y, v/x and v/y in two-dimensional flows. Both simulation and experiment demonstrate that high accuracy and high spatial resolution are possible with this technique. The large time required to compute the cross-correlations, however, limits its wide applications at present.     

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.     

An application of nonsymmetric Lax-Milgram lemma to nonassociated plasticity

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

Three-dimensional numerical simulation of free convection in a plane-parallel flow

The results are given of numerical simulation of convection in Couette and Poiseuille flow for stationary (T*/t = 0) and uniformly increasing (T*/t = = const) mean temperature of the convective layer. Numerical experiments were made for air (Prandtl number Pr = 0.7) in the range of Rayleigh numbers 2 000 Ra 44 000. The results confirm the conclusion drawn in earlier studies that at slightly supercritical Ra the dominant forms in the convection are cylindrical rolls with helical circulation. The rolls are oriented along the direction of the flow and are stationary formations. When a certain value of Ra, which depends on the vertical distribution of the temperature and velocity, is reached, the roll structure is deformed by transverse perturbations. All the considered flow forms have a stabilizing influence on the transverse modes, occurring at larger Rayleigh numbers than is the case for convection in a fluid at rest. The perturbations are displaced at a phase velocity close to the mean velocity of the undisturbed flow. In the considered range of Rayleigh numbers, a shear flow does not have an appreciable influence on the heat transfer, although there is a certain tendency for the Nusselt number to be larger in a shear flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 127–135, May–June, 1979.We thank N. M. Sazanovich and L. I. Derevich for assistance in the calculations.     

Filament stretching equation in Lagrangian coordinates

The general integral of the very simple equation ^21/n/() was found to describe the cross sectional area of filaments of isothermal power law fluids while in transient stretching where is time and is the initial location of fluid molecules at time = 0 given as the distance from a reference point fixed in space. Any such stretching transient given as a solution of the above equation is physically realizable subject to the restrictions > 0 and/ < 0.     

Measurement of instantaneous flow rates in periodically operating injection systems

Research and development work carried out to provide a method to measure accurately instantaneous flow rates in periodically operating injection systems is summarized. The instantaneous flow rate is reconstructed from axial velocity time series measured by a laser Doppler anemometer on the center-line of a capilary pipe flow. The theoretical background, on which the evaluation of the instantaneous flow rate is based is provided. It is shown that the axial mean velocity is sufficient to reconstruct the periodically varying flow rate and the pressure gradient. The application of the proposed method is described and an instrument is suggested that can be employed in many fields where fast, periodically varying flow rates occur and instantaneous information is needed.List of Symbols C _j complex Fourier coefficients - C.C. complex conjugate of different variables - D inner diameter of pipe - f ₀ repetition frequency of valve motion - i complex unit - J ₀ Bessel function of zeroth order - J _l Bessel function of first order - L distance between valve and measuring volume - m integrated fluid mass within one period - P/z mass flow rate - P/z pressure gradient along the pipe axis - p ₀,p _n normalized amplitude of the complex pressure gradient - r radial position - R pipe radius - Re Reynolds number - t time - T time of the period of the repetition valve opening - Ta Taylor number - U axial velocity - V integrated fluid volume within one period - instantaneous volumetric flow rate - viscosity - kinematic viscosity - fluid density - normalized open valve time - angular velocity of Fourier modes     

LDA measurements in low Reynolds number turbulent boundary layer

LDA measurements of the mean velocity in a low Reynolds number turbulent boundary layer allow a direct estimate of the friction velocity U from the value of /y at the wall. The trend of the Reynolds number dependence of / is similar to the direct numerical simulations of Spalart (1988).     

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     

Two-phase flow in heterogeneous porous media II: Numerical experiments for flow perpendicular to a stratified system

In this work, we make use of numerical experiments to explore our original theoretical analysis of two-phase flow in heterogeneous porous media (Quintard and Whitaker, 1988). The calculations were carried out with a two-region model of a stratified system, and the parameters were chosen be consistent with practical problems associated with groundwater flows and petroleum reservoir recovery processes. The comparison between theory (the large-scaled averaged equations) and experiment (numerical solution of the local volume averaged equations) has allowed us to identify conditions for which the quasi-static theory is acceptable and conditions for which a dynamic theory must be used. Byquasi-static we mean the following: (1) The local capillary pressure,everywhere in the averaging volume, can be set equal to the large-scale capillary pressure evaluated at the centroid of the averaging volume and (2) the large-scale capillary pressure is given by the difference between the large-scale pressures in the two immiscible phases, and is therefore independent of gravitational effects, flow effects and transient effects. Bydynamic, we simply mean a significant departure from the quasi-static condition, thus dynamic effects can be associated with gravitational effects, flow effects and transient effects. To be more precise about the quasi-static condition we need to refer to the relation between the local capillary pressure and the large-scale capillary pressure derived in Part I (Quintard and Whitaker, 1990). Herep _c ¦_y represents the local capillary pressure evaluated at a positiony relative to the centroid of the large-scale averaging volume, and {p _c }¦_x represents the large-scale capillary pressure evaluated at the centroid.In addition to{p _c } ^c being evaluated at the centroid, all averaged terms on the right-hand side of Equation (1) are evaluated at the centroid. We can now write the equations describing the quasi-static condition as , , This means that the fluids within an averaging volume are distributed according to the capillary pressure-saturation relationwith the capillary pressure held constant. It also means that the large-scale capillary pressure is devoid of any dynamic effects. Both of these conditions represent approximations (see Section 6 in Part I) and one of our main objectives in this paper is to learn something about the efficacy of these approximations. As a secondary objective we want to explore the influence of dynamic effects in terms of our original theory. In that development only the first four terms on the right hand side of Equation (1) appeared in the representation for the local capillary pressure. However, those terms will provide an indication of the influence of dynamic effects on the large-scale capillary pressure and the large-scale permeability tensor, and that information provides valuable guidance for future studies based on the theory presented in Part I.Roman Letters A scalar that maps {}*/t onto - A scalar that maps {}*/t onto - A interfacial area between the -region and the -region contained within, m² - A interfacial area between the -region and the -region contained within, m² - A interfacial area between the -region and the -region contained within, m² - a vector that maps ({}*/t) onto , m - a vector that maps ({}*/t) onto , m - b vector that maps ({p}– g) onto , m - b vector that maps ({p}– g) onto , m - B second order tensor that maps ({p}– g) onto , m² - B second order tensor that maps ({p}– g) onto , m² - c vector that maps ({}*/t) onto , m - c vector that maps ({}*/t) onto , m - C second order tensor that maps ({}*/t) onto , m² - C second order tensor that maps ({}*/t) onto . m² - D third order tensor that maps ( ) onto , m - D third order tensor that maps ( ) onto , m - D second order tensor that maps ( ) onto , m² - D second order tensor that maps ( ) onto , m² - E third order tensor that maps () onto , m - E third order tensor that maps () onto , m - E second order tensor that maps () onto - E second order tensor that maps () onto - p _c =(), capillary pressure relationship in the-region - p _c =(), capillary pressure relationship in the-region - g gravitational vector, m/s² - largest of either or - - - i unit base vector in thex-direction - I unit tensor - K local volume-averaged-phase permeability, m² - K local volume-averaged-phase permeability in the-region, m² - K local volume-averaged-phase permeability in the-region, m² - {K } large-scale intrinsic phase average permeability for the-phase, m² - K –{K }, large-scale spatial deviation for the-phase permeability, m² - K –{K }, large-scale spatial deviation for the-phase permeability in the-region, m² - K –{K }, large-scale spatial deviation for the-phase permeability in the-region, m² - K ^* large-scale permeability for the-phase, m² - L characteristic length associated with local volume-averaged quantities, m - characteristic length associated with large-scale averaged quantities, m - I _i i = 1, 2, 3, lattice vectors for a unit cell, m - l characteristic length associated with the-region, m - ; characteristic length associated with the-region, m - l _H characteristic length associated with a local heterogeneity, m - - n unit normal vector pointing from the-region toward the-region (n =–n ) - n unit normal vector pointing from the-region toward the-region (n =–n ) - p pressure in the-phase, N/m² - p local volume-averaged intrinsic phase average pressure in the-phase, N/m² - {p } large-scale intrinsic phase average pressure in the capillary region of the-phase, N/m² - p local volume-averaged intrinsic phase average pressure for the-phase in the-region, N/m² - p local volume-averaged intrinsic phase average pressure for the-phase in the-region, N/m² - p –{p }, large scale spatial deviation for the-phase pressure, N/m² - p –{p }, large scale spatial deviation for the-phase pressure in the-region, N/m² - p –{p }, large scale spatial deviation for the-phase pressure in the-region, N/m² - P _c p –{p }, capillary pressure, N/m² - {p_c}^c large-scale capillary pressure, N/m² - r ₀ radius of the local averaging volume, m - R ₀ radius of the large-scale averaging volume, m - r position vector, m - , m - S /, local volume-averaged saturation for the-phase - S ^* {}*{}*, large-scale average saturation for the-phaset time, s - t time, s - u , m - U , m² - v -phase velocity vector, m/s - v local volume-averaged phase average velocity for the-phase in the-region, m/s - v local volume-averaged phase average velocity for the-phase in the-region, m/s - {v } large-scale intrinsic phase average velocity for the-phase in the capillary region of the-phase, m/s - {v } large-scale phase average velocity for the-phase in the capillary region of the-phase, m/s - v –{v }, large-scale spatial deviation for the-phase velocity, m/s - v –{v }, large-scale spatial deviation for the-phase velocity in the-region, m/s - v –{v }, large-scale spatial deviation for the-phase velocity in the-region, m/s - V local averaging volume, m³ - V volume of the-phase in, m³ - V large-scale averaging volume, m³ - V capillary region for the-phase within, m³ - V capillary region for the-phase within, m³ - V _c intersection of m³ - V volume of the-region within, m³ - V volume of the-region within, m³ - V () capillary region for the-phase within the-region, m³ - V () capillary region for the-phase within the-region, m³ - V () , region in which the-phase is trapped at the irreducible saturation, m³ - y position vector relative to the centroid of the large-scale averaging volume, m Greek Letters local volume-averaged porosity - local volume-averaged volume fraction for the-phase - local volume-averaged volume fraction for the-phase in the-region - local volume-averaged volume fraction for the-phase in the-region - local volume-averaged volume fraction for the-phase in the-region (This is directly related to the irreducible saturation.) - {} large-scale intrinsic phase average volume fraction for the-phase - {} large-scale phase average volume fraction for the-phase - {}* large-scale spatial average volume fraction for the-phase - –{}, large-scale spatial deviation for the-phase volume fraction - –{}, large-scale spatial deviation for the-phase volume fraction in the-region - –{}, large-scale spatial deviation for the-phase volume fraction in the-region - a generic local volume-averaged quantity associated with the-phase - mass density of the-phase, kg/m³ - mass density of the-phase, kg/m³ - viscosity of the-phase, N s/m² - viscosity of the-phase, N s/m² - interfacial tension of the - phase system, N/m - , N/m - , volume fraction of the-phase capillary (active) region - , volume fraction of the-phase capillary (active) region - , volume fraction of the-region ( + =1) - , volume fraction of the-region ( + =1) - {p }– g, N/m³ - {p }– g, N/m³     

Numerical Investigation of the Natural Convection of Water in the Neighborhood of the Density Inversion Point For Grashof Numbers up to 106

The natural convection of fresh water in a square cell is considered at a temperature close to the density inversion temperature for Grashof numbers 2.9 · 10⁴ Gr 10⁶. As a result of the numerical investigation, one steady-state and three self-oscillating regimes are found in addition to the three steady-state flows previously detected earlier and described for low Grashof numbers ( 0 Gr 2 · 10⁵). The basic characteristics of the unsteady flows are analyzed by means of the Fourier method, the fundamental oscillation frequencies are found, and the flow evolution and the variation of the oscillation characteristics with increase in the Gr number are considered.     

Modeling pressure-strain rate correlations in turbulence theory

On the basis of a spectral representation of the rapid part _ij,2 of the correlation tensor p(u _i /x _j ) using Cramer's theorem the inequality _ij,2(U _j /x _i )0 is obtained. As distinct from the realizability conditions, it can serve as a direct and very rigorous test of the adequacy of model expressions for _ij,2. In particular, it is shown that the best known of such expressions do not satisfy this test.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.2, pp. 42–46, March–April, 1992.     

Stability criterion for steady-state burning of powders

A model of a powder with a variable surface temperature is considered. A stability criterion for the steady-state burning of powders at constant pressure is found on the assumption that all the processes in the gas phase and the reaction layer of the condensed phase are inertialess. It is shown that the stability region is determined by two parameters: K=(T₁ – t₀) ( ln m/T₀)p and r=(T₁/T₀)p (here T₁ is the surface temperature, T₁ is the initial temperature of the powder, m is the mass burning rate, p is pressure). Burning is always stable if < 1. If > 1, the burning process is stable only when r > ( – 1)²/ / + 1.The author expresses his appreciation to A. S. Kompaneits, O. I. Leipunskii, A. G. Istratov, V. B. Librovich, and S. S. Novikov for discussing his work.     

On the Ellipticity of the Middle Surface of a Shell and its Application to the Asymptotic Analysis of ‘Membrane’ Shells

We consider a surface S = (), where ² is a bounded, connected, open set with a smooth boundary and : ³ is a smooth map; let () denote the components of the two-dimensional linearized strain tensor of S and let 0 with length 0 > 0. We assume the the norm ,|| ()||0, in the space V0() = { H¹() × H¹() × L²(); = 0 on 0 } is equivalent to the usual product norm on this space. We then establish that this assumption implies that the surface S is uniformly elliptic and that we necessarily have 0 = .     

Experimental investigation of the resistance to the flow of a conducting fluid in flat insulated ducts in the presence of a transverse magnetic field with allowance for fringe effects and wall roughness

The distribution of pressure, velocity, and electrical potential has been investigated for a mercury flow in insulated rectangular ducts with a large side ratio (Hartmann-type flow). The ranges of variation of the Reynolds, Hartmann, and Stewart numbers were 7·10²R5·10⁵, 0H490, and 0N24, respectively. Special attention is given to the sections of the channel where the flow enters and leaves the magnetic field. In these zones the pressure is sharply nonuniform and the velocity profiles in a plane perpendicular to the field acquire an M shape. A relation is established between the length of the entrance section, where the flow is three-dimensional, and the MHD similarity criteria. It is shown that ducts which are hydraulically smooth in the absence of a magnetic field become increasingly rough as the field grows stronger. Data are obtained on the resistance coefficient for a stabilized flow measured in a magnetic field and on the dependence of the critical Reynolds number on the Hartmann number.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 10–21, July–August, 1971.     

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.     

Stress-strain behavior of segmented polyurethaneureas under pure shear deformation

Elastic properties of segmented polyurethaneureas (SPUs) under pure shear deformation were investigated. Data were analyzed by using strain energy density function (W). The values of the derivatives of W with respect to the invariants of the strain tensor (I _i ; i = 1, 2, 3 ) at zero strain limit were also estimated theoretically and were compared with those estimated by experiment. The limiting value of the derivative with respect to I ₁ (W/I ₁) was shown theoretically to be 5 G/8, while the derivatives with respect to I ₂ and I ₃ (W/I ₂ and W/I ₃ were respectively-G/8 and -3G/8. The theoretical prediction could explain the asymptotic behavior of the derivatives of SPUs as well as isoprene rubber (IR) reported by Kawabata et al. at small I ₁ limit.     

Neck propagation as a shock wave

Neck propagation in the stretching of elastic solid filaments having a yield point was analyzed using the space one-dimensional thin filament governing equations developed previously by the authors and other researchers. Constitutive model for the filament was assumed to be expressible as engineering tensile stress(X) (tensile force) given as a function of elongational strain with the(X) curve having a yield point maxima followed by a minima and a breaking point greater than the yield point maxima. Also incorporated into the model is the hysteresis of irreversible plastic deformation. When inertia is taken into consideration, the thin filament equations were found to reduce to the nonlinear wave equation ² (X)/ ² =C ₁ ² X/ ² where is Lagrangean space coordinate, is time, andC ₁ is inertia coefficient. The above nonlinear wave equation yields a solutionX(, ) having a stepwise discontinuity inX which propagates along the axis. The zero speed limit of the step wave solution was found to describe the above neck propagation occurring in solid filaments. Furthermore, it was recognized that the nonlinear wave equation was known for many years to also govern the plastic shock wave which propagates axially within a metal rod subjected to a very strong impact on its end. The one-dimensional atmospheric shock wave also was known to be governed by the nonlinear wave equation upon making certain simplifying assumptions. The above and other evidences lead to the conclusion that neck propagation occurring in the extension of solid filament obeying the above(X) function can be formally described as a shock wave.     

The rheology of powders

This paper studies the slow flow of powders. It is argued that since powders can flow like liquids, there must be equations similar to those of liquids. The phenomenon of a variable density, dilatancy, is described by an analogue of temperature called the compactivity X. Whereas, in thermal physicsT = E/S, powders are controlled byX = V/S. The equations for, v, T of a liquid are replaced by, v, X. An analogy for free energy is described, and the solution to some simple problems of packing and mixing are offered. As an example of rheology, it is shown that the simplest flow equations produce a transition to plug flow in appropriate circumstances.Delivered as a Gold Medal Lecture at the Golden Jubilee Conference of the British Society of Rheology and Third European Rheology Conference, Edinburgh, 3–7 September, 1990.