Spanwise vorticity measurements in a perturbed boundary layer

Hot-wire measurements of the spanwise vorticity fluctuation _z have been carried out in a turbulent boundary layer subjected to concentrated suction, applied through a porous wall strip. The results indicate that, relative to the no-suction case, the rms value of _z is significantly reduced in the near-wall region, with this reduction increasing with the suction rate . The reduction near the wall suggests an alteration in the dynamics of the layer. As the Reynolds number increases, the dynamics of the layer in the near-wall region become more intense and suction becomes less effective.     

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     

Normal forms for random diffeomorphisms

Given a dynamical system (,, ,) and a random diffeomorphism (): ^{d d} with fixed point at x=0. The normal form problem is to construct a smooth near-identity nonlinear random coordinate transformation h() to make the random diffeomorphism ()=h()^–1() h() as simple as possible, preferably linear. The linearization D(, 0)=:A() generates a matrix cocycle for which the multiplicative ergodic theorem holds, providing us with stochastic analogues of eigenvalues (Lyapunov exponents) and eigenspaces. Now the development runs pretty much parallel to the deterministic one, the difference being that the appearance of turns all problems into infinite-dimensional ones. In particular, the range of the homological operator is in general not closed, making the conceptof-normal form necessary. The stochastic versions of resonance and averaging are developed. The case of simple Lyapunov spectrum is treated in detail.     

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.     

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.     

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.     

Generation of vortices in a partially filled,rapidly rotating cylinder

We describe a system in which vortices are shed from a cylindrical free surface approximately centered in a rotating flow. Shedding is controlled by the parameter =2 g/ ² d, where g, , d denote gravity, rotation rate and the diameter of the free surface. We find vortex shedding for >0.162 and no vortex shedding for < 0.0847. The range depends on the aspect ratio L/d, where L is the column length, in a nonmonotonic fashion. These results are independent of viscosity and surface tension for small values of these parameters.Now at Martin Marietta, Orlando Aerospace, PO Box 5837, Mail Point 150, Orlando, FL 32855, USA     

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.     

Harmonic solutions of micropolar elastodynamics

The equations of micropolar elastodynamics are considered for an unbounded continuum subjected to a body force and a body couple. These act harmonically with the same real frequency , but with individual arbitrary spatial distributions. Over a harmonic state, the displacement and microrotation are related to two radiation conditioned harmonic vectors, each acquiring three eigenvalue contributions, assuming a noncritical -frequency. Altogether, four distinct eigenvalues are admissible. If ²<2² ₀, ₀ being a frequency parameter of the continuum, two of these are real while two are purely imaginary. But if ²<2² ₀, then all admissible eigenvalues are real. Each eigenvalue contribution resolves into a series of Hankel and Bessel functions coupled to Hankel type transforms of: (i) spherical integrals which, in turn, can be expanded via spherical harmonics for the 3-dimensional problem, (ii) circular integrals for the 2-dimensional problem. Axisymmetric and spherically symmetric results are deduced in 3-dimensions. Asymptotic solutions are also established; they disclose long-range formation of radially attenuated spherical (or circular) waves propagating with, generally, anisotropic amplitudes but, invariably, isotropic eikonals.If, in the absence of a body couple, a body force acts radially in 3-dimensions with a spherically symmetric strength, then the elastic displacement behaves likewise while the microrotation vanishes identically. Another application is made to a 2-dimensional problem for a 1 × 3 source system of body force plus body couple without longitudinal variation but with magnitudes symmetric about a longitudinal axis.As approaches a certain critical frequency , dependent solely on the continuum, at least two eigenvalues approach the same value. The phenomenon is explored for a continuum consistent with ²<2² ₀ and under the hypothesis ²<2² ₀. All admissible eigenvalues are then real throughout an -neighbourhood of . Here, two associated eigenvalue contributions behave singularly. Nevertheless, their essential singularities cancel out within the relevant combination. Examination of a far-field suggests that critical frequency attainment sets off a slow instability in the 2-dimensional configuration. In the 3-dimensional configuration, however, it preserves stability and eliminates radial attenuation; an exact solution is formulated for this case.     

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.     

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

An analogue study for flame flickering

An analogue experiment is proposed to simulate flame flickering comprising a free ascending column fed on its side with a light gas (helium) emerging from a vertical slot in ambient air. The convective motion of the helium jet is considered to represent the motion of burnt gases of buoyant jet flames. The helium jet is accelerated by buoyancy effects and the flow field is similar to that of burnt gases observed for real buoyant flames. The vertical velocity profile of the steady helium jet is measured at different vertical distances. The unsteady helium jet is also studied by measuring the instability frequency as a function of ambient pressure at different injection flow rates, and by analyzing the tomography images of the helium jet. The instability morphology is the same as that observed on real buoyant flames. We conclude that this type of instability can be approximately characterized by the maximum vertical velocityu _max, and the distance betweenu _max in the helium ascending column andu = o in the ambient air. For this type of instability the local vorticity is proportional to which can be influenced by gravity and ambient pressure. Theoretical prediction of the instability frequency as a function of gravity and ambient pressure has been obtained, and is in good agreement with the experimental results.List of symbols C ₁,C ₂ constants - F instability frequency - F _c critical frequency - F _m the most amplified frequency - F (K, ) function defined in (11) - g gravitational acceleration - g reduced gravity acceleration g(₀-^*)/^* - k real wave number of the disturbance - K reduced wave numberK=2k - K _c reduced wave number of the critical instability mode - K _m nondimensional wavenumber of the most amplified mode - L vertical characteristic length (in x direction) - P ambient pressure - u local vertical buoyant velocity (inx direction) - u _max local maximum vertical velocity - v local velocity component iny direction (horizontal) - V ₀ injection velocity of helium (iny direction) - x vertical distance measured from the leading edge of boundary layer - y horizontal distance measured from the exit plane of the vertical slot - Z(K, ) function defined in equation (11) Greek symbols distance betweenu _max in the helium ascending column andu = o in the ambient air - - wavelength of instability - _c critical wavelength - _m the most amplified wavelength - ^* helium density at slot exit - ₀ ambient air density - ^* helium dynamic viscosity at slot exit - v ^* helium kinematic viscosity at slot exit - complex number presented in disturbancee ^i(kx+t) - _i imaginary part of , representing the amplification rate of disturbance - _r real part of , where ( _r /k) represents the group velocity - reduced complex number of , defined      

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     

Critical layer formation in a stratified medium with mean shear flows

The problem of the excitation of internal waves with a given wave number k and frequency in a stratified medium with shear flows is considered. The internal wave field of the form v(z)exp(–it+ikx) established as t in a medium without dissipation has a singular point at the level z=z₀ (critical level), at which the flow velocity U(z) coincides with the phase velocity /k. Dissipative effects (viscosity and heat conduction) smooth out this singularity. An exact solution of the model equation describing as t and zz₀ the field excited by oscillating sources activated at t=0 is constructed with allowance for dissipation. This makes it possible to describe the limiting steady-state field, determine the critical layer as the neighborhood of the critical level in which dissipation effects are important, and to estimate its width and the rate of convergence to the limiting steady-state regime. The asymptotic behavior of the fields is examined for Ri1, where Ri is the Richardson number. It is shown that when the well-known Miles stability condition Ri>1/4 is satisfied there are no natural oscillations with a critical level.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 82–93, May–June, 1990.     

Global Behavior of a Nonlinear Quasiperiodic Mathieu Equation

In this paper, we investigate the interaction of subharmonicresonances in the nonlinear quasiperiodic Mathieu equation,x + [ + (cos ₁ t + cos ₂ t)] x + x³ = 0.We assume that 1 and that the coefficient of the nonlinearterm, , is positive but not necessarily small.We utilize Lie transform perturbation theory with elliptic functions –rather than the usual trigonometric functions – to study subharmonic resonances associated with orbits in 2m:1 resonance with a respective driver. In particular, we derive analytic expressions that place conditions on (, , ₁, ₂) at which subharmonic resonance bands in a Poincaré section of action space begin to overlap. These results are used in combination with Chirikov's overlap criterion to obtain an overview of the O() global behavior of equation (1) as a function of and ₂ with ₁, , and fixed.     

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

Hydrodynamic models as applied to the investigation of magnetic plasma stability

In the present paper magnetohydrodynamic models are employed to investigate the stability of an inhomogeneous magnetic plasma with respect to perturbations in which the electric field may be regarded as a potential field (rot E 0). A hydrodynamic model, actually an extension of the well-known Chew-Goldberg er-Low model [1], is used to investigate motions transverse to a strong magnetic field in a collisionless plasma. The total viscous stress tensor is given; this includes, together with magnetic viscosity, the so-called inertial viscosity.Ordinary two-fluid hydrodynamics is used in the case of strong collisions=. It is shown that the collisional viscosity leads to flute-type instability in the case when, collisions being neglected, the flute mode is stabilized by a finite Larmor radius. A treatment is also given of the case when epithermal high-frequency oscillations (not leading immediately to anomalous diffusion) cause instability in the low-frequency (drift) oscillations in a manner similar to the collisional electron viscosity, leading to anomalous diffusion.Notation f particle distribution function - E electric field component - H₀ magnetic field - density - V particle velocity - e charge - m, M electron and ion mass - _i, _e ion and electron cyclotron frequencies - viscous stress tensor - P pressure - r_i Larmor radius - P pressure tensor - t time - frequency - T temperature - collision frequency - collision time - j current density - _i, _e ion and electron drift frequencies - k_x, k_y, k_z wave-vector components - n₀ particle density - g acceleration due to gravity. The authors are grateful to A. A. Galeev for valuable discussion.     

Free damped oscillations of the linear viscoelastic material

This work shows that the complex dynamic viscosity ^* () can be evaluated from experimental data of free damped oscillations (FDO) when upper and lower angular displacements measured in a plate-plate geometry are analyzed by Fourier series. We deduce the equation required to evaluate experimentally the dynamic functions of linear viscoelastic materials in FDO. Two loss tangents can be defined in FDO, which are not necessarily equal to tan () of the forced oscillation test.     

Thin-Walled Beams: the Case of the Rectangular Cross-Section

In this paper we present an asymptotic analysis of the three-dimensional problem for a thin linearly elastic cantilever =×(0,l) with rectangular cross-section of sides and ², as goes to zero. Under suitable assumptions on the given loads, we show that the three-dimensional problem converges in a variational sense to the classical one-dimensional model for extension, flexure and torsion of thin-walled beams. Mathematics Subject Classifications (2000) 474K20, 74B10, 49J45.