Nonlinear analysis of the forced response of a beam with three mode interaction

An analysis is presented for the primary resonance of a clamped-hinged beam, which occurs when the frequency of excitation is near one of the natural frequencies,_n . Three mode interaction (₂ 3₁ and _{3 1} + 2₂) is considered and its influence on the response is studied. The case of two mode interaction (₂ 3₁) is also considered to compare it with the case of three mode interaction. The straight beam experiencing mid-plane stretching is governed by a nonlinear partial differential equation. By using Galerkin's method the governing equation is reduced to a system of nonautonomous ordinary differential equations. The method of multiple scales is applied to solve the system. Steady-state responses and their stability are examined. Results of numerical investigations show that there exists no significant difference between both modal interactions' influences on the responses.     

Oscillatory measurements of silicone oils —Loss and storage modulus master curves

The work describes a way to obtain loss modulus and storage modulus master curves from oscillatory measurements of silicone oils.The loss modulus master curve represents the dependence of the viscous flow behavior on · ₀ ^* and the storage modulus master curve — the dependence of the elastic flow behavior on · ₀ ^* .The relation between the values of the loss modulus and storage modulus master curves (at a certain frequency) is a measurement of the viscoelastic behavior of a system. The G/G-ratio depends on · ₀ ^* which leads to a viscoelastic master curve. The viscoelastic master curve represents the relation between the elastic and viscous oscillatory flow behavior.     

Evaluation of the dynamic-mechanical properties of solids from a cooperative model for stress relaxation

For many solid materials the stress relaxation process obeys the universal relationF = – (d/d lnt)_max = (0.1 ± 0.01) ( ₀ — _i ), regardless of the structure of the material. Here denotes the stress,t the time, ₀ the initial stress of the experiment and _i the internal stress. A cooperative model accounting for the similarity in relaxation behaviour between different materials was developed earlier. Since this model has a spectral character, the concepts of linear viscoelasticity are used here to evaluate the corresponding prediction of the dynamic mechanical properties, i.e. the frequency dependence of the storageE () and lossE () moduli. Useful numerical approximations ofE () andE () are also evaluated. It is noted that the universal relation in stress relaxation had a counterpart in the frequency dependence ofE (). The theoretical prediction of the loss factor for high-density polyethylene is compared with experimental results. The agreement is good.     

Flow in porous media II: The governing equations for immiscible,two-phase flow

The Stokes flow of two immiscible fluids through a rigid porous medium is analyzed using the method of volume averaging. The volume-averaged momentum equations, in terms of averaged quantities and spatial deviations, are identical in form to that obtained for single phase flow; however, the solution of the closure problem gives rise to additional terms not found in the traditional treatment of two-phase flow. Qualitative arguments suggest that the nontraditional terms may be important when / is of order one, and order of magnitude analysis indicates that they may be significant in terms of the motion of a fluid at very low volume fractions. The theory contains features that could give rise to hysteresis effects, but in the present form it is restricted to static contact line phenomena.Roman Letters (, = , , and ) A interfacial area of the- interface contained within the macroscopic system, m² - A _e area of entrances and exits for the -phase contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A ^* 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² - g gravity vector, m²/s - H mean curvature of the- interface, m^–1 - H area average of the mean curvature, m^–1 - H – H , deviation of the mean curvature, m^–1 - I unit tensor - K Darcy's law permeability tensor, m² - K permeability tensor for the-phase, m² - K viscous drag tensor for the-phase equation of motion - K viscous drag tensor for the-phase equation of motion - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the-phase, m - n unit normal vector pointing from the-phase toward the-phase (n = –n ) - p _c p – P , capillary pressure, N/m² - p pressure in the-phase, N/m² - p intrinsic phase average pressure for the-phase, N/m² - p – p , spatial deviation of the pressure in the-phase, N/m² - r ₀ radius of the averaging volume, m - t time, s - v velocity vector for the-phase, m/s - v phase average velocity vector for the-phase, m/s - v intrinsic phase average velocity vector for the-phase, m/s - v – v , spatial deviation of the velocity vector for the-phase, m/s - V averaging volume, m³ - V volume of the-phase contained within the averaging volume, m³ Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - viscosity of the-phase, Nt/m² - surface tension of the- interface, N/m - viscous stress tensor for the-phase, N/m² - / kinematic viscosity, m²/s     

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     

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     

On the effects of uniform high suction on the rotationally symmetric flow of a conducting liquid near a stationary disc in the presence of a transverse magnetic field

In the present paper an attempt has been made to find out effects of uniform high suction in the presence of a transverse magnetic field, on the motion near a stationary plate when the fluid at a large distance above it rotates with a constant angular velocity. Series solutions for velocity components, displacement thickness and momentum thickness are obtained in the descending powers of the suction parameter a. The solutions obtained are valid for small values of the non-dimensional magnetic parameter m (= 4 _e ² H ₀ ² /) and large values of a (a2).Nomenclature a suction parameter - E electric field - E _r , E , E _z radial, azimuthal and axial components of electric field - F, G, H reduced radial, azimuthal and axial velocity components - H magnetic field - H _r , H , H _z radial, azimuthal and axial components of magnetic field - H ₀ uniform magnetic field - H* displacement thickness and momentum thickness ratio, */ - h induced magnetic field - h _r , h , h _z radial, azimuthal and axial components of induced magnetic field - J current density - m nondimensional magnetic parameter - p pressure - P reduced pressure - R Reynolds number - U ₀ representative velocity - V velocity - V _r , V , V _z radial, azimuthal and axial velocity components - w ₀ uniform suction through the disc. - density - electrical conductivity - kinematic viscosity - _e magnetic permeability - a parameter, (/)^1/2 z - a parameter, a - * displacement thickness - momentum thickness - angular velocity     

Numerisch einfache Beziehungen zwischen Verlust- und Speicherkomponente des dynamischen Schermoduls und der dynamischen Nachgiebigkeit

Zusammenfassung Es werden Formeln hergeleitet, die analog den Kramers-Kronig-Beziehungen Verlust- und Speicherkomponente des komplexen dynamischen Schermoduls bzw. der komplexen dynamischen Nachgiebigkeit miteinander in Verbindung bringen, die aber den Vorteil haben, daß kein Integral mehr auftritt und daß sich der Einfluß der jeweiligen Kurvenausläufer für 0 und in Form einer einfach gebauten Summe aus Stützstellenwerten darstellen läßt, wobei als Stützstellen die Punkte 2 ^j ₀ mitj= ±1, ±2, benützt werden.Die Formeln enthalten zusätzlich die logarithmische Ableitung einer der jeweiligen Meßgrößen an der Stelle ₀. Es ist jedoch möglich, Näherungsformeln aufzustellen, welche die Ableitung nicht mehr enthalten. Der relative Fehler, der dabei in Kauf zu nehmen ist, beträgt weniger als 0,5%.Durch Linearkombination der Gleichungen für entsprechende Stützstellen kann der Anteil der unendlichen Summen eliminiert werden, so daß man Beziehungen erhält, welche zur Kontrolle einer Messung beider Komponenten des Moduls beziehungsweise der Nachgiebigkeit dienen können. Um den Preis einer geringeren Genauigkeit, wobei die maximale Abweichung relativ zur Verlustkomponente angegeben wird, können solche Kontrollbeziehungen schon für wenige, mindestens aber drei Stützstellen aufgestellt werden.Darüber hinaus sind aufgrund der Formeln einseitige Schranken für die GrenzwerteG , usw. aus den Meßwerten berechenbar. Die Schranken sind um so schärfer, je größer der gemessene Bereich ist. Bis zu einem gewissen Grade kann auch der Verlauf der Funktionen über die Grenzen des Meßbereichs hinaus abgeschätzt werden.

---

Summary Formulae are derived, which relate the loss modulus to the storage modulus and the loss compliance to the storage compliance as the Kramers-Kronig equations do. The advantage of the formulae given is, that they involve simple sums of functional values instead of integrals. Additionally the logarithmic derivative of for example the storage modulus at ₀ must be known to calculate loss modulus at ₀. But there are also given approximation formulae, which do not involve any derivative. The relative error of all approximations is smaller than 0.5%.By combining the equations at different values of the influence of unaccessible parts of one of the components can be eliminated. In this way relations are obtained, which are very useful for checking measurements. The error of each of these approximate relations can be given relative to the loss component. As a further application, formulae are discussed, which allow the calculation of an upper limit forG and a lower limit for the viscosity, when both components of dynamic modulus or of compliance are measured. The lower or upper limits are the sharper the greater the-range where measurements are available.

---

---

Mit 5 Abbildungen und 6 Tabellen     

Neue kritische Drehzahlen von einfach besetzten Wellen

Zusammenfassung Wird eine mit einer einzelnen Scheibe besetzte, masselose Welle mit einer periodischen Längskraft P ₀üP ₁ cos t oder mit einem periodischen Torsionsmoment W ₀üW ₁ cos t belastet, so ergibt es unter der Annahme, daß bei umlaufenden Wellen die Frequenz mit der Drehzahl der Welle gemäß der Beziehung = gekoppelt ist, auch bei konstantem und um so mehr bei rein harmonisch periodischem Antriebsmoment der Scheibe neben der kritischen Drehgeschwindigkeit erster Art = _k unendlich viele kritische Zustände zweiter Art bei den Drehgeschwindigkeiten = _k/¦ n±1¦ (n=1, 2, 3, ...), falls die Scheibe etwas exzentrisch auf die Welle gekeilt ist. Die Gefährlichkeit der kritischen Drehzahlen zweiter Art ist jedoch viel geringer als die der kritischen Drehzahl erster Art.Außer den genannten kritischen Drehzahlen ergeben sich ohne Dämpfung noch unendlich viele kritische Drehzahlbereiche in der Nähe des 2/m v-fachen der kritischen Drehzahl _k, in denen die mit der Scheibe besetzte Welle exponentiell mit der Zeit anwachsende Biegeschwingungen ausführt. Die Breite dieser Bereiche und die Ausweichgeschwindigkeit der Scheibe nehmen mit wachsendem P ₁ bzw. W ₁ zu. Da die Ausweichgeschwindigkeit exponentiell mit der Zeit anwächst, unterscheiden sich diese neuen kritischen Drehzahlen dritter Art wesentlich von den kritischen Drehzahlen erster und zweiter Art, bei denen die Scheibe ja mit konstanter Geschwindigkeit ausweicht. Dieser Unterschied kommt besonders dann zur Geltung, wenn man noch den Einfluß einer geschwindigkeitsproportionalen Dämpfung berücksichtigt. Eine solche begrenzt bei den kritischen Zuständen erster und zweiter Art das seitliche Ausweichen der Scheibe, während bei den neu gefundenen kritischen Drehzahlen dritter Art die Amplitude der Biegeschwingung nach wie vor unbegrenzt anwächst. Allerdings fallen infolge der Dämpfung die Bereiche mit höheren Ordnungszahlen m weg, und die übrigen werden um so kleiner, je größer die dazugehörende Ordnungszahl m ist. In der Praxis müssen — wie an zwei Zahlenbeispielen gezeigt wurde — ein erster Linie die Instabilitätsbereiche beachtet werden, die in der Nähe von 2 _k/ und _k/ liegen.     

Representation of the rheological properties of polymer melts in terms of complex fluidity

In dynamic rheological experiments melt behavior is usually expressed in terms of complex viscosity ^* () or complex modulusG ^* (). In contrast, we attempted to use the complex fluidity ^* () = 1/µ ^* () to represent this behavior. The main interest is to simplify the complex-plane diagram and to simplify the determination of fundamental parameters such as the Newtonian viscosity or the parameter of relaxation-time distribution when a Cole-Cole type distribution can be applied. ^* () complex shear viscosity - () real part of the complex viscosity - () imaginary part of the complex viscosity - G ^* () complex shear modulus - G() storage modulus in shear - G() loss modulus in shear - J ^* () complex shear compliance - J() storage compliance in shear - J() loss compliance in shear - shear strain - rate of strain - angular frequency (rad/s) - shear stress - loss angle - ^* () complex shear fluidity - () real part of the complex fluidity - () imaginary part of the complex fluidity - ₀ zero-viscosity - ₀ average relaxation time - h parameter of relaxation-time distribution     

Impacts of local dispersion and first-order decay on solute transport in randomly heterogeneous porous media

Stochastic subsurface transport theories either disregard local dispersion or take it to be constant. We offer an alternative Eulerian-Lagrangian formalism to account for both local dispersion and first-order mass removal (due to radioactive decay or biodegradation). It rests on a decomposition of the velocityv into a field-scale componentv , which is defined on the scale of measurement support, and a zero mean sub-field-scale componentv _s , which fluctuates randomly on scales smaller than. Without loss of generality, we work formally with unconditional statistics ofv _s and conditional statistics ofv . We then require that, within this (or other selected) working framework,v _s andv be mutually uncorrelated. This holds whenever the correlation scale ofv is large in comparison to that ofv _s . The formalism leads to an integro-differential equation for the conditional mean total concentration c which includes two dispersion terms, one field-scale and one sub-field-scale. It also leads to explicit expressions for conditional second moments of concentration cc. We solve the former, and evaluate the latter, for mildly fluctuatingv by means of an analytical-numerical method developed earlier by Zhang and Neuman. We present results in two-dimensional flow fields of unconditional (prior) mean uniformv . These show that the relative effect of local dispersion on first and second moments of concentration dies out locally as the corresponding dispersion tensor tends to zero. The effect also diminishes with time and source size. Our results thus do not support claims in the literature that local dispersion must always be accounted for, no matter how small it is. First-order decay reduces dispersion. This effect increases with time. However, these concentration moments c and cc of total concentrationc, which are associated with the scale below, cannot be used to estimate the field-scale concentrationc directly. To do so, a spatial average over the field measurement scale is needed. Nevertheless, our numerical results show that differences between the ensemble moments ofc and those ofc are negligible, especially for nonpoint sources, because the ensemble moments ofc are already smooth enough.     

The theory of a vibrating-rod densimeter

The paper presents a theory of a device for the accurate determination of the density of fluids over a wide range of thermodynamic states. The instrument is based upon the measurement of the characteristics of the resonance of a circular section tube, or rod, performing steady, transverse oscillations in the fluid. The theory developed accounts for the fluid motion external to the rod as well as the mechanical motion of the rod and is valid over a defined range of conditions. A complete set of working equations and corrections is obtained for the instrument which, together with the limits of the validity of the theory, prescribe the parameters of a practical design capable of high accuracy.Nomenclature A, B, C, D constants in equation (60) - A _j , B _j constants in equation (18) - a _j ⁺ , a _j ^– wavenumbers given by equation (19) - C _f drag coefficient defined in equation (64) - C _f ^/0 , C _f ^/1 components of C _f in series expansion in powers of - c speed of sound - D _b drag force of fluid b - D ₀ coefficient of internal damping - E extensional modulus - force per unit length - F _j ⁺ , F _j ^– constants in equation (24) - f, g functions of defined in equations (56) - G modulus of rigidity - I second moment of area - K constant in equation (90) - k, k constants defined in equations (9) - L half-length of oscillator - Ma Mach number - m _a mass per unit length of fluid a - m _b added mass per unit length of fluid b - m _s mass per unit length of solid - n _j eigenvalue defined in equation (17) - P power (energy per cycle) - P _a , P _b power in fluids a and b - p pressure - R radius of rod or outer radius of tube - R _c radius of container - R _i inner radius of tube - r radial coordinate - T tension - T _visc temperature rise due to heat generation by viscous dissipation - t time - v _r , v radial and angular velocity components - y lateral displacement - z axial coordinate - dimensionless tension - _a dimensionless mass of fluid a - _b dimensionless added mass of fluid b - _b dimensionless drag of fluid b - dimensionless parameter associated with - ₀ dimensionless coefficient of internal damping - dimensionless half-width of resonance curve - dimensionless frequency difference defined in equation (87) - spatial resolution of amplitude - R, , , _s , increments in R, , , _s , - dimensionless amplitude of oscillation - dimensionless axial coordinate - ratio of to - _a , _b ratios of to for fluids a and b - angular coordinate - parameter arising from distortion of initially plane cross-sections - _f thermal conductivity of fluid - dimensionless parameter associated with - viscosity of fluid - _a , _b viscosity of fluids a and b - dimensionless displacement - _j jth component of - density of fluid - _a , _b density of fluids a and b - _s density of tube or rod material - density of fluid calculated on assumption that * - dimensionless radial coordinate - * dimensionless radius of container - dimensionless times - _{rr rr}, _r radial normal and shear stress components - spatial component of defined in equation (13) - _j jth component of - dimensionless streamfunction - ⁰, ¹ components of in series expansion in powers of - phase angle - _r phase difference - _ra , _rb phase difference for fluids a and b - streamfunction - _j jth component defined in equation (22) - dimensionless frequency (based on ) - _a , _b dimensionless frequency in fluids a and b - _s dimensionless frequency (based on _s ) - angular frequency - ₀ resonant frequency in absence of fluid and internal damping - _r resonant frequency in absence of internal fluid - _ra , _rb resonant frequencies in fluids a and b - dimensionless frequency - dimensionless frequency when _a vanishes - dimensionless frequencies when _a vanishes in fluids a and b - dimensionless resonant frequency when _a , _b, _b and ₀ vanish - dimensionless resonant frequency when _a , _b and _b vanish - dimensionless resonant frequency when _b and _b vanish - dimensionless frequencies at which amplitude is half that at resonance     

The wedge subjected to tractions: a paradox resolved

The classical two-dimensional solution provided by Lévy for the stress distribution in an elastic wedge, loaded by a uniform pressure on one face, becomes infinite when the opening angle 2 of the wedge satisfies the equation tan 2 = 2. Such pathological behavior prompted the investigation in this paper of the stresses and displacements that are induced by tractions of O(r ^–) as r0. The key point is to choose an Airy stress function which generates stresses capable of accommodating unrestricted loading. Fortunately conditions can be derived which pre-determine the form of the necessary Airy stress function. The results show that inhomogeneous boundary conditions can induce stresses of O(r ^–), O(r ^– ln r), or O(r ^– ln² r) as r0, depending on which conditions are satisfied. The stress function used by Williams is sufficient only if the induced stress and displacement behavior is of the power type. The wedge loaded by uniform antisymmetric shear tractions is shown in this paper to exhibit stresses of O(ln r) as r0 for the half-plane or crack geometry. At the critical opening angle 2, uniform antisymmetric normal and symmetric shear tractions also induce the above type of stress singularity. No anticipating such stresses, Lévy used an insufficiently general Airy stress function that led to the observed pathological behavior at 2.     

Generalized Rouse model for a dilute solution of clustered polymers

A theory analogue to tha of Rouse is given, to describe the rheological behavior of dilute solutions consisting of clusters of crosslinked polymers. The frequency-dependent behavior of the dynamic moduli of these fluids differs substantially from that of the well-known Rouse-like fluid (GG^1/2). In our case the storage modulus becomes proportional to ^3/2, while the loss modulus is proportional to . The loss modulus dominates the dynamic behavior for frequencies smaller than the largest normal frequency of the clusters.     

The investigation on the properties and structures of starting vortex flow past a backward-facing step by WBIV technique

An experimental investigation of a starting vortex flow around a backward-facing step was conducted in a water channel. The properties and structures of the flow were investigated by qualitative flow visualization using the hydrogen bubble method and by quantitative velocity and vorticity measurements using White-light Bubble Image Velocimetry (WBIV) — a newly developed PIV method. Some invariant properties and 4-stage structures of starting vortex flow were observed.List of symbols a flow acceleration during starting stage - h height of backward-facing step - d _v dimensionless vortex size - t time - t dimensionless time - U free uniform velocity - u, v streamwise and spanwise velocity components respectively - Re Reynolds number based on a and h - x, y streamwise and spanwise coordinates respectively in flow field - x _c , y _c dimensionless vortex center position - vorticity - ov dimensionless vorticity - _max maximum vorticity - ov _max dimensionless maximum vorticity - circulation - dimensionless circulation - kinematic viscosity This work was supported by the CNSF Grant 1939 100-1-3     

Vorticity characteristics of the turbulent intermediate wake

Measurements of the lateral components _j (j=2 and 3) of the vorticity fluctuation vector have been made, using a vorticity probe consisting of two X-wires, in the intermediate wake of a circular cylinder. The effect of the spatial resolution of the probe on the measurement of _j has been studied. As the spatial resolution impairs, the variance and flatness factor of _j decrease whereas the skewness of _j increases. Reasonably accurate values of _j ² can be obtained by applying spectral corrections for the spatial resolution effect.Near the beginning of the intermediate wake, the variance of ₂ is larger than that of ₃ due to the significant contribution from ribs which connect consecutive spanwise roll vortices. This difference decreases with downstream distance. Also, the presence of the rolls is reflected by a local extremum in the skewness of ₃ on each side of the wake centerline. The magnitude of the extremum decreases with downstream distance.The support of the Australian Research Council is gratefully acknowledged.     

Das Dispersionsmodell

Zusammenfassung Ein Vergleich im Frequenzbereich zeigt, daß bei der Berechnung der Verweilzeitverteilung mit dem Dispersionsmodell das endlich begrenzte System für Péclet-Zahlen Pe > 10 mit guter Näherung durch ein einseitig unbegrenztes System und für Pe > 50 durch ein beidseitig unbegrenztes System ersetzt werden kann.

---

The dispersion model. A comparison of approximations

A comparison in the frequency domain shows that for the determination of the residence time distribution with the dispersion model the finitely restricted system may be substituted with good approximation for Peclet numbers Pe > 10 by a one-side unrestricted system and for Pe > 50 by a both-side unrestricted system.

---

Bezeichnungen A Rohrquerschnitt - A=A mittlerer Strömungsquerschnitt in der Schüttschicht - Konzentration (Partialdichte) der Bezugskomponente i - Bezugskonzentration nach Gl. (5) - c_i Konzentration (Dichte) der reinen Bezugskomponente i - D axialer Dispersionskoeffizient - E Fehler im Frequenzbereich nach Gl. (36) - G(,) Übertragungsfunktion - G(,i) Frequenzgang - L Länge der Schüttschicht - m Masse - Massenstrom - Péclet-Zahl - s Laplace-Variable - t Zeit - t Impulsbreite - v Strömungsgeschwindigkeit im leeren Rohr - mittlere axiale Strömungsgeschwin digkeit in der Schüttschicht - V=AL Zwischenraumvolumen der Schüttschicht - x Ortskoordinate - (t) Dirac-Stoss - Porosität - dimensionslose Zeit - dimensionslose Konzentration - Laplace-Transformierte der Konzentration - Fourier-Transformierte der Konzentration - dimensionslose Ortskoordinate - =s dimensionslose Laplace-Variable - mittlere Verweilzeit - Kreisfrequenz - = dimensionslose Kreisfrequenz Indices A Ausgang - D Dispersion - E Eingang - i Bezugskomponente - K Konvektion Mitteilung Nr. 44 des Institutes für Mess-und Regel-technik der Eidgenössischen Technischen Hochschule Zürich (Vorsteher: Prof. Dr. P. Profos)     

Approximate First Integrals of a Galaxy Model

First-order approximate first integrals (conserved quantities)of a Hamiltonian dynamical system with two degrees of freedomwhich arises in the modeling of central part of a deformed galaxy [1] havebeen obtained based on the approximate Noether symmetries for resonances₁=₂, ₁=2₂ and 2₁=3₂. Furthermore,KAM curves have been obtained analytically and they have been compared with thenumerical ones on the Poincaré surface of section.     

The extension of Beltrami's ellipsoid to anisotropic bodies

Summary The spectral decomposition of the compliance, stiffness, and failure tensors for transversely isotropic materials was studied and their characteristic values were calculated using the components of these fourth-rank tensors in a Cartesian frame defining the principal material directions. The spectrally decomposed compliance and stiffness or failure tensors for a transversely isotropic body (fiber-reinforced composite), and the eigenvalues derived from them define in a simple and efficient way the respective elastic eigenstates of the loading of the material. It has been shown that, for the general orthotropic or transversely isotropic body, these eigenstates consist of two double components, ₁ and ₂ which are shears (₂ being a simple shear and ₁, a superposition of simple and pure shears), and that they are associated with distortional components of energy. The remaining two eigenstates, with stress components ₃, and ₄, are the orthogonal supplements to the shear subspace of ₁ and ₂ and consist of an equilateral stress in the plane of isotropy, on which is superimposed a prescribed tension or compression along the symmetry axis of the material. The relationship between these superimposed loading modes is governed by another eigenquantity, the eigenangle .The spectral type of decomposition of the elastic stiffness or compliance tensors in elementary fourth-rank tensors thus serves as a means for the energy-orthogonal decomposition of the energy function. The advantage of this type of decomposition is that the elementary idempotent tensors to which the fourth-rank tensors are decomposed have the interesting property of defining energy-orthogonal stress states. That is, the stress-idempotent tensors are mutually orthogonal and at the same time collinear with their respective strain tensors, and therefore correspond to energy-orthogonal stress states, which are therefore independent of each other. Since the failure tensor is the limiting case for the respective _x, which are eigenstates of the compliance tensor S, this tensor also possesses the same remarkable property.An interesting geometric interpretation arises for the energy-orthogonal stress states if we consider the projections of _x in the principal₃D stress space. Then, the characteristic state ₂ vanishes, whereas stress states ₁, ₃ and ₄ are represented by three mutually orthogonal vectors, oriented as follows: The ₃ and ₄ lie on the principal diagonal plane (₃₁₂) with subtending angles equaling (–/2) and (-), respectively. On the positive principal ₃-axis, is the eigenangle of the orthotropic material, whereas the ₁-vector is normal to the (₃₁₂)-plane and lies on the deviatoric -plane. Vector ₂ is equal to zero.It was additionally conclusively proved that the four eigenvalues of the compliance, stiffness, and failure tensors for a transversely isotropic body, together with value of the eigenangle , constitute the five necessary and simplest parameters with which invariantly to describe either the elastic or the failure behavior of the body. The expressions for the _x-vector thus established represent an ellipsoid centered at the origin of the Cartesian frame, whose principal axes are the directions of the ₁-, ₃- and ₄-vectors. This ellipsoid is a generalization of the Beltrami ellipsoid for isotropic materials.Furthermore, in combination with extensive experimental evidence, this theory indicates that the eigenangle alone monoparametrically characterizes the degree of anisotropy for each transversely isotropic material. Thus, while the angle for isotropic materials is always equal to _i = 125.26° and constitutes a minimum, the angle || progressively increases within the interval 90–180° as the anisotropy of the material is increased. The anisotropy of the various materials, exemplified by their ratiosE _L/2G_L of the longitudinal elastic modulus to the double of the longitudinal shear modulus, increases rapidly tending asymptotically to very high values as the angle approaches its limits of 90 or 180°.     

Flow visualization and analysis of aerodynamically detuned supersonic rotors

The effect of aerodynamic detuning on the supersonic steady and unsteady blade passage flow field is experimentally investigated on a free surface water table by means of color Schlieren and shadowgraph flow visualization techniques. Two aerodynamic detuning mechanisms are considered: (1) alternate circumferential spacing of adjacent airfoils; (2) the replacement of alternate airfoils with splitters. The steady flow visualization demonstrates the significant effect of aerodynamic detuning on the passage flow field and, in particular, the shock wave-airfoil surface intersection locations. The unsteady flow visualization studies show the importance of the interblade phase angle. A mathematical model is also described and utilized to demonstrate the enhanced aeroelastic stability associated with the altered cascade passage shock wave structure due to these aerodynamic detuning mechanisms.List of symbols a dimensionless perturbation sonic velocity - C airfoil chord - I _x mass moment of inertia - k reduced frequency, k = C/2 - K spring constant - M _R dimensionless unsteady aerodynamic moment - M Mach number - P _split split splitter circumferential spacing - P _start splitter leading edge location - R reference full chord airfoil - R _s reference splitter - Sp _chord splitter chord length - u dimensionless perturbation chordwise velocity - v dimensionless perturbation normal velocity - amplitude of oscillation - interblade phase angle - level of aerodynamic detuning - undamped natural torsional frequency - ₀ reference frequency - flutter frequency