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.     

The shear viscosity dependence on concentration,molecular weight,and shear rate of polystyrene solutions

The solution viscosity of narrow molecular weight distribution polystyrene samples dissolved in toluene and trans-decalin was investigated. The effect of polymer concentration, molecular weight and shear rate on viscosity was determined. The molecular weights lay between 5 10⁴ and 24 10⁶ and the concentrations covered a range of values below and above the critical valuec ^*, at which the macromolecular coils begin to overlap. Flow curves were generated for the solutions studied by plotting log versus log . Different molecular weights were found to have the same viscosity in the non-Newtonian region of the flow curves and follow a straight line with a slope of – 0.83. A plot of log ₀ versus logM _w for 3 wt-% polystyrene in toluene showed a slope of approximately 3.4 in the high molecular weight regime. Increasing the shear rate resulted in a viscosity that was independent of molecular weight. The sloped (log)/d (logM _w ) was found to be zero for molecular weights at which the corresponding viscosities lay on the straight line in the power-law region.On the basis of a relation between _sp and the dimensionless productc · [], simple three-term equations were developed for polystyrene in toluene andt-decalin to correlate the zero-shear viscosity with the concentration and molecular weight. These are valid over a wide concentration range, but they are restricted to molar masses greater than approximately 20000. In the limit of high molecular weights the exponent ofM _w in the dominant term in the equations for both solvents is close to the value 3.4. That is, the correlation between _sp andc · [] results in a sloped(log _sp)/d(logc · []) of approximately 3.4/a at high values ofc · [] wherea is the Mark-Houwink constant. This slope of 3.4/a is also the power ofc in the plot of ₀ versusc at high concentrations. a Mark-Houwink constant - B ₁,B ₂,B _n constants - c concentration (g · cm^–3) - c ^* critical concentration (g · cm^–3) - K, K constants - K _H Huggins constant - M molecular weight - M _c critical molecular weight - M _n number-average molecular weight - M _w weight-average molecular weight - n sloped(log _sp)/d (logc · []) at highc · [] - PS polystyrene - T temperature (K) - shear rate (s^–1) - critical shear rate (s^–1) - viscosity (Pa · s) - ₀ zero-shear viscosity (Pa · s) - _s solvent viscosity (Pa · s) - _sp specific viscosity - [] intrinsic viscosity (cm³ · g^–1) - dynamic viscosity (Pa · s) - | ^*| complex dynamic viscosity (Pa · s) - angular frequency (rad/s) - density of polymer solution (g · cm^–3) - ₁₂ shear stress (Pa) Dedicated to Prof. Dr. J. Schurz on the occasion of his 60^th birthday.Excerpt from the dissertation of Reinhard Kniewske: Bedeutung der molekularen Parameter von Polymeren auf die viskoelastischen Eigenschaften in wäßrigen und nichtwäßrigen Medien, Technische Universität Braunschweig 1983.     

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.     

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.     

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/.     

Relative viscosities of polymer melts filled with inorganic fillers

Summary Viscosities of polyethylene melts filled with natural calcium carbonate (NC), glass fiber (GF), and carbon fiber and polypropylene melts filled with NC, talc, and GF were examined. Constant relative viscosities _r were obtained on these filled systems if use was made of the definition of _r suggested in an earlier paper. The relationships between _r and the volume fraction of the filler were discussed qualitatively in conjunction with earlier data. The magnitude of _r for the material filled with fiber increases with the increase of the average aspect ratio of the fiber.

---

Zusammenfassung Es werden die Viskositäten von Polyäthylenschmelzen untersucht, die mit natürlichem Kalziumkarbonat (NC), Glasfasern (GF) und Karbonfasern, sowie von Propylenschmelzen, die mit Kalziumkarbonat, Kreide und Glasfasern gefüllt sind. Man findet bei diesen gefüllten Systemen konstante relative Viskositäten _r , wenn man dafür die in einer früheren Untersuchung eingeführte Definition zugrundelegt. Die Beziehungen zwischen _r und dem Volumenanteil des Füllstoffs werden unter Bezugnahme auf frühere Ergebnisse qualitativ diskutiert. Die Größe von _r wächst bei den mit Fasern gefüllten Systemen mit zunehmendem Verhältnis von Länge zu Dicke der Fasern.

---

---

With 6 figures and 3 tables     

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.     

Saint-Venant's principle in linear two-dimensional elasticity for non-striplike domains

Let be a simply connected domain in the x ₁-x ₂ plane which lies within the strip 0<x ₂, is a simple closed piecewise smooth curve. Let l= [(x ₁, x ₂): (x ₁, x ₂) and x ₁>0], _l = [(x ₁ x ₂): (x ₁ ,x ₂) and x ₁>1>0].Suppose that a two-dimensional homogeneous isotropic elastic body occupies , that a self-equilibrated stress loading is applied to - l, and that l is stress-free. Knowles [2] and Flavin [6] showed that the elastic energy in _l decays exponentially with respect to l with an exponential decay constant of the form k/b, where k is a universal constant. It is shown here that a decay constant of the form c/ may be obtained where c is a universal constant and is a characteristic dimension of , which is more appropriate than b for general non-striplike domains. In addition, an appropriate decay theorem is obtained for coil-like domains.     

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.     

Development characteristics of drag-reducing surfactant solution flow in a duct

Development characteristics of dilute cationic surfactant solution flow have been studied through the measurements of the time characteristics of surfactant solution by birefringence experiments and of the streamwise mean velocity profiles of surfactant solution duct flow by a laser Doppler velocimetry system. For both experiments, the concentration of cationic surfactant (oleylbishydroxymethylethylammonium chloride: Ethoquad O/12) was kept constant at 1000 ppm and the molar ratio of the counter ion of sodium salicylate to the surfactants was at 1.5. From the birefringence experiments, dilute surfactant solution shows very long retardation time corresponding to micellar shear induced structure formation. This causes very slow flow development of surfactant solution in a duct. Even at the end of the test section with the distance of 112 times of hydraulic diameter form the inlet, the flow is not fully developed but still has the developing boundary layer characteristics on the duct wall. From the time characteristics and the boundary layer development, it is concluded that the entry length of 1000 to 2000 times hydraulic diameter is required for fully developed surfactant solution flow.List of abbreviations and symbols A₁, A₂ Coefficients for time constant fitting [-] - B Breadth of the test duct [m] - C₁, C₂ Coefficients for time constant fitting [-] - D Pipe diameter [m] - D_H Hydraulic diameter [m] - g Impulse response function [Pa] - H Width of the test duct [m] - n Index of Bird-Carreau model [-] - Re Reynolds number (=U_mD_H/) - Re_D Pipe Reynolds number (=U_mD/) - Re_x Streamwise distance Reynolds number (=U₀x/) - T Absolute temperature [K] - t Time [s] - t_a Retardation time [s] - t_b Build-up time [s] - t_x Relaxation time [s] - t_x1, t_x2 Relaxation time for double time constant fitting [s] - t Time constant in Bird-Carreau model [s] - U Time mean velocity [m/s] - U_m Bulk mean velocity [m/s] - U_max Maximum velocity in a pipe [m/s] - U₀ Main flow velocity [m/s] - u Friction velocity [m/s] - x, y Coordinates [m] - Shear rate [s^–1] - Mean shear rate [s^–1] - n Birefringence [-] - 99% boundary layer thickness [m] - Solution viscosity [Pa·s] - _P, _S Surfactant and solvent viscosity [Pa·s] - ₀, Zero and infinite viscosity of Bird-Carreau model [Pa·s] - Characteristic time in Maxwell model [s] - Water kinematic viscosity [m²/s] - Density [kg/m³] - Solution shear stress [Pa] - _P, _S Surfactant and solvent shear stress [Pa] - Time in convolution [s]     

Bäcklund transformations for the equation ∂2u/∂x1∂x1 + ∂2u/∂x2∂x2−f (u)

Bäcklund transformations for the equation ² u/x¹x ¹+ ²u/x ² x ²=f (u)here f is an arbitrary function) is studied in this paper, using the procedure of Wahlquist and Estabrook (WEP). We conclude that the condition d ² f/du ²=f is sufficient for the existence of Bäcklund transformations for the equation of our interest. A special case of our results leads to the conclusion of Leibbrandt ^[1,2].     

On the prediction of riblet performance with engineering turbulence models

The paper reports the outcome of a numerical study of fully developed flow through a plane channel composed of ribleted surfaces adopting a two-equation turbulence model to describe turbulent mixing. Three families of riblets have been examined: idealized blade-type, V-groove and a novel U-form that, according to computations, achieves a superior performance to that of the commercial V-groove configuration. The maximum drag reduction attained for any particular geometry is broadly in accord with experiment though this optimum occurs for considerably larger riblet heights than measurements indicate. Further explorations bring out a substantial sensitivity in the level of drag reduction to the channel Reynolds number below values of 15 000 as well as to the thickness of the blade riblet. The latter is in accord with the trends of very recent, independent experimental studies.Possible shortcomings in the model of turbulence are discussed particularly with reference to the absence of any turbulence-driven secondary motions when an isotropic turbulent viscosity is adopted. For illustration, results are obtained for the case where a stress transport turbulence model is adopted above the riblet crests, an elaboration that leads to the formation of a plausible secondary motion sweeping high momentum fluid towards the wall close to the riblet and thereby raising momentum transport.Nomenclature c _f Skin friction coefficient - c _f Skin friction coefficient in smooth channel at the same Reynolds number - k Turbulent kinetic energy - K ⁺ k/ _w - h Riblet height - S Riblet width - H Half height of channel - Re Reynolds number = volume flow/unit width/ - Modified turbulent Reynolds number - R _t turbulent Reynolds numberk ²/ - P _k Shear production rate ofk, _t (U _i /x _j + U _j /x _i ) U _i /x _j - dP/dz Streamwise static pressure gradient - U _i Mean velocity vector (tensor notation) - U Friction velocity, _w/ where _w=–H dP/dz - W Mean velocity - W _b Bulk mean velocity through channel - y ⁺ yU /v. Unless otherwise stated, origin is at wall on trough plane of symmetry - Kinematic viscosity - _t Turbulent kinematic viscosity - Turbulence energy dissipation rate - Modified dissipation rate – 2(k ^1/2/x _j )² - Density - _k , Effective turbulent Prandtl numbers for diffusion ofk and      

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.     

On the sensitive dynamical system and the transition from the apparently deterministic process to the completely random process

The detailed analysis of the dynamical process of coin tossing is made. Through calculations, it is illustrated how and why the result is extremely sensitive to the initial conditions. It is also shown that, as the initial height of the mass center of the coin increases, the final configuration, i.e. head or tail, becomes more and more sensitive to the initial parameters (the initial velocity angular velocity, and the initial orientation), the coefficient of the air drag, and the energy absorption factor of the surface on which the coin bounces. If we keep the head upward initially but allow a small range for the change of some other initial parameters, the frequency that the final configuration is head, would be 1 if the initial height h of the mass center is sufficiently small, and would be clo to 1/2 if h is sufficiently large. An interesting question is how this frequency changes continuously from 1 to 1/2 as h increases. Detailed calculations show that such a transition is very similar to the transition from laminar to turbulent flows. A basic difference between the transition stage and the completely random stage is indicated: In the completely random stage, the deterministic process of the individual case is extremely sensitive to the initial conditions and the dynamical parameters, out the statistical properties of the ensemble are insensitive to the small changes of the initial conditions and the dynamical parameters. On the contrary, in the transition stage, both the deterministic process of the individual case and the statistical properties of the ensemble are sensitive to the initial conditions and the dynamical parameters. The mechanism for this feature of the transition stage is the existence of the long-train structure in the parameter space. The illuminations of this analysis on some other random phenomena are discussed.     

Beitrag zur Auswertung von Anlauf-Vorgängen in viskoelastischen Flüssigkeiten

Zusammenfassung In diesem Bericht wird die Kraftübertragung zwischen zwei unbegrenzt ausgedehnten und durch eineMaxwell-Flüssigkeit gekoppelten Platten bei beliebigem Zeitgesetz der Bewegung der oberen auf die elastisch aufgehängte untere Platte unter Einbeziehung der Trägheit von unterer Platte undMaxwell-Flüssigkeit berechnet.Diese Fragestellung führt mathematisch auf ein Randwertproblem, dessen eine Randbedingung die Form ³ u/x ² y + ₁ u + ₁ u/y – ₂ ² u/xy = 0 fürx = 0 annimmt, also gemischte Ableitungen nach 2 unabhängigen Veränderlichen enthält.Es wird gezeigt, daß sich die Lösung als unendliche Reihe darstellen läßt, deren 1. Glied aus abklingenden Exponentialfunktionen, alle weiteren hingegen aus gedämpften Schwingungen bestehen.Als spezielle Zeitgesetze für die obere Platte werden die impulsförmige Erregung, die zeitproportionale und die periodische Bewegung untersucht.Es wird abgeleitet, wie sich die Kenngrößen undG (Viskosität und Schubmodul) aus Eigenfrequenz und Dämpfung bzw. aus den beiden Exponenten dere-Funktionen in den beiden ersten Fällen berechnen lassen. Bei periodischer Beanspruchung werden undG aus der Resonanzfrequenz und der Halbwertbreite bzw. dem Phasenwinkel ermittelt.Abschließend wird gezeigt, daß und wie ein Rotationsviskosimeter auf den behandelten Fall ebener Platten zurückgeführt werden kann.

---

Summary In this report the force transmission between two infinitely extended plates coupled by aMaxwell-Fluid is calculated. In this calculation the inertia of the fluid and the influence of the elastically suspended lower plate is included, while the upper plate can perform any time-dependent movement.Mathematically this question leads to a boundery value problem; one boundery condition has the form ³ u/ ² y + ₁ u + ₁ u/y – ₂ ² u/xy = 0 (x = 0) ( ₁, ₂, = constant), that means it contains mixed differential-quotients of two independent variables.The solution is demonstrated as an infinite series, the form of the first term is of two negative exponential functions, all the other terms are damped oscillations.The excitation from a single pulse, the uniform timeproportional and the periodic movement as timedependent laws for the upper plate are analysed.The values of andG (viscosity and shear modulus) are derived from eigenfrequency and damping-constant respectively from the values of the two exponentials. In the periodic case andG are calculated from the resonance frequency and halfwidth or phase-angle respectively.Finally the conditions of the rotation-viscometer are reduced to the case of parallel-plate-viscometer.

---

---

Vorgetragen auf der Jahrestagung der Deutschen Rheologen vom 7.–9. Mai 1973 in Berlin.

---

Mit 6 Abbildungen und 1 Tabelle     

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

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

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).     

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.     

Anwendung des Übertragungsverfahrens bei stark abklingenden Lösungsfunktionen

Übersicht Bei stark abklingenden Funktionen wird die Übertragungsmatrix U() aufgespalten in die Anteilc U ₁() e^– und U ₂() e. Der zweite Term spielt am Rand = 0 keinc Rolle. Die unbekannten Anfangswerte sind über die Matrix U ₁(0) an die bekannten gebunden und eindeutig bestimmbar.

---

Summary For strongly decaying solution functions the transfer matrix U() is splitted into the parts U ₁() e^– and U ₂() e. The second term does not influence at the boundary = 0. The unknown initial values are related by the matrix U ₁(0) to the known values and they can be uniquely determined.

     

Qualitative behavior of dissipative wave equations on bounded domains

The qualitative behavior of solutions of the mixed problem u_tt = u-a(x)u_t in IR x , u=0 on IR x , is studied in the case when a>0 and IRⁿ is bounded. Roughly speaking, if aa_min>0, then solutions decay at least as fast as exp t( –1/2a_min), with the possible exception of a finite dimensional set of smooth solutions whose existence is associated with a phenomenon of overdamping. If a_max is sufficiently small, depending on , then no overdamping occurs.Partially supported by NSF grant NSF GP 34260.This work was partially supported by the National Science Foundation under Grant No. GP 34260