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

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.     

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     

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.     

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.     

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.     

Surface waves in deep water due to arbitrary periodic surface pressure

Summary The expressions for the velocity potential and the surface elevation are obtained for the three dimensional motion on surface waves generated by a surface pressure F(x, y)e^{i t}, t>0, on the surface of a deep sea. The wave integrals are asymptotically evaluated for large distances and times by the method of stationary phase. By a direct passage to the limit t, the steady state solution of the problem is evaluated in exact terms. The results are illustrated for a class of physically plansible pressure distributions on the surface. Some features of the wave motions are discussed.

---

Sommario Si ottengono le espressioni per il potenziale della velocità e la elevazione della superficie nel moto tridimensionale di onde superficiali generate da una pressione superficiale F(x, y)e^{i t}, t>0, sulla superficie di un mare profondo. Gli integrali di onde sono calcolati asintoticamente per grandi distanze e tempi con il metodo della fase stazionaria. Con un passaggio diretto al limite t, si calcola in termini esatti la soluzione di regime del problema. I risultati vengono illustrati per una classe di distribuzioni di pressione superficiale fisicamente plausibile, Si discutono alcune caratteristiche dei moti ondosi.

     

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.     

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.     

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.     

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.     

Flow over a body with development of a steady separation zone as re → ∞

This paper discusses formulation of the total problem of flow of an incompressible liquid over a body, with formation of a closed stationary separation zone as Re . The scheme used is based on the method of matched asymptotic expansions [1]. Following [1], it is postulated that the separated zone is developed (i.e., it is not infinitely fragmented and does not vanish as Re ), and the flow inside it has a definite degree of regularity with respect to Re. With these hypotheses we can use the Prandtl-Batchelor theorem [2], which states that, in the limit as Re , a region of circulating flow becomes vortex flow of an inviscid liquid with constant vorticity . Therefore, a basis for constructing matched asymptotic expansions is the vortex-potential problem (the problem of determining a stream function , satisfying the equation = 0 in the region of translational motion and the equation = in a certain region, unknowna priori, of circulating motion). In the general case the solution of the vortex-potential problem depends on two parameters: the total pressure p_o and the vorticity in the separated zone. These parameters appear in the condition for matching the solutions of the first and second boundary-layer approximations (at the boundary of the separated zone for the end Re values) with the corresponding solutions for the inviscid flow. It is shown in the present paper that the conditions for matching the cyclic boundary layer with the external translational flow are the same additional relations which allow us to close the total problem. Thus, in using the method of matched asymptotic expansions to solve the problem of flow over a body with closed stationary separation zones one must simultaneously consider no less than two approximations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 28–37, March–April, 1978.The authors thank G. Yu. Stepanov for discussion of the paper and valuable comments.     

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     

Evolution of Domains of Attraction of a Forced Beam with Two-Mode Interaction

When a dissipative dynamical system has multiple attractors, it is a task to determine and recognize the global domain of attraction of each attractor. In this paper we study the global behavior of a forced hinged-clamped beam with two-mode interaction. The equation of motion of the beam is reduced to four first-order autonomous ordinary differential equations. The system has an internal resonance condition of 2 31, where 1 and 2 denote natural frequencies of the first and second modes, respectively. When the excitation frequency is near 1, the system can have three equilibrium solutions, among which two are asymptotically stable and one is unstable. We examine how the domains of attraction of two stable equilibrium solutions evolve as the forcing frequency is varied across jump points. By using a special plane which contains all equilibrium points, called the principal plane, the global domains of attraction can be discussed more effectively. Results show that knowledge of this evolution helps us better understand the jump phenomenon.     

Concentration dependence of the critical viscoelastic properties of gelatin at the gel point

Gelatin gel properties have been studied through the evolution of the storage [G()] and the loss [G()] moduli during gelation or melting near the gel point at several concentrations. The linear viscoelastic properties at the percolation threshold follow a power-law G()G() and correspond to the behavior described by a rheological constitutive equation known as the Gel Equation. The critical point is characterized by the relation: tan = G/G = cst = tan ( · /2) and it may be precisely located using the variations of tan versus the gelation or melting parameter (time or temperature) at several frequencies. The effect of concentration and of time-temperature gel history on its variations has been studied. On gelation, critical temperatures at each concentration were extrapolated to infinite gel times. On melting, critical temperatures were determined by heating step by step after a controlled period of aging. Phase diagrams [T = f(C)] were obtained for gelation and melting and the corresponding enthalpies were calculated using the Ferry-Eldridge relation. A detailed study of the variations of A with concentration and with gel history was carried out. The values of which were generally in the 0.60–0.72 range but could be as low as 0.20–0.30 in some experimental conditions, were compared with published and theoretical values.     

Taylor–Couette flow with independently rotating end plates

Results are presented from a combined numerical and experimental study of steady bifurcation phenomena in a modified Taylor–Couette geometry where the end plates of the flow domain are allowed to rotate independently of the inner cylinder. The ends rotate synchronously and the ratio between the rate of rotation of the ends _e and the inner cylinder _i defines a control parameter :=_e/_i. Stationary ends favour inward motion along the end walls whereas rotating walls promote outward flow. We study the exchange between such states and focus on two-cell flows, which are found in the parameter range between =0 and =1 for =2. Hence is used as an unfolding parameter. A cusp bifurcation is uncovered as the organizing centre for the stability exchange between the two states. Symmetry breaking bifurcations, which lead to flows that break the mid-plane symmetry are also revealed. Overall, excellent agreement is found between numerical and experimental results. PACS 47.20, 47.11, 47.54     

Dynamic melt viscoelastic behavior of random copolymers

Summary The viscoelastic properties of 65/35 styrenen-butyl methacrylate random copolymers were determined using the Eccentric Rotating Disks device of the Rheometrics Mechanical Spectrometer. Similar to the behavior observed in homopolymers, an increase in the molecular weight of the copolymer resulted in extension of the rubbery plateau and in a reduction in the terminal region. The dynamic complex viscosity showed onset of non-Newtonian flow at higher frequencies, with the non-Newtonian region increasing with increasing molecular weight.The elastic modulus,G, was dependent upon the frequency,, asG ^1.5 in the terminal region, rather than asG ² observed for polystyrene. The viscous modulus,G, was proportional to the frequency,, asG , similar to what is observed for polystyrene. The dynamic viscosity | ^*| at high frequencies showed a region independent of molecular weight where a power law of | ^*| ^0.9 is applicable, consistent with entanglement models. Thy dynamic viscosity at low frequencies in the Newtonian region is related to molecular weight as |^*| . Using WLF equations, the coefficient of expansion, _f , was obtained that, together with glass transition, showed a negative deviation from the Fox-Flory relationship.

---

Zusammenfassung Die viskoelastischen Eigenschaften von statistischen 65/35-Styrol/n-Butyl-Methacrylat-Kopolymeren wurden mit Hilfe einer Maxwell-Rheometer-Anordnung in Verbindung mit dem Mechanischen Spektrometer der Fa. Rheometrics bestimmt. Ähnlich dem bei Homopolymeren beobachteten Verhalten ergab sich auch hier mit wachsendem Molekulargewicht eine Verbreiterung des Kautschuk-Plateaus und eine Verkleinerung des Endbereichs. Die komplexe Viskosität zeigte erst bei höheren Frequenzen das Einsetzen nicht-newtonschen Fließens an, wobei der nichtnewtonsche Bereich mit steigendem Molekulargewicht größer wurde.Der SpeichermodulG ergab sich im Endbereich als proportional zu ^1,5, im Unterschied zu der bei Polystyrol beobachteten Proportionalität mit ². Dagegen war der VerlustmodulG der Frequenz direkt proportional, ähnlich wie es auch bei Polystyrol beobachtet worden war. Die dynamische Viskosität | ^*| zeigte unabhängig vom Molekulargewicht bei hohen Frequenzen einen Bereich, in dem eine Potenz-Beziehung | ^*| ~ ^0,9 herrschte, was auf die Wirkung von Verzweigungen hindeutet. Dagegen galt bei den niedrigen Frequenzen des newtonschen Bereichs|^*| ~ . Mit Hilfe der WLF-Gleichung wurde der Ausdehnungskoeffizient _f bestimmt, der ebenso wie der Glasübergang eine negative Abweichung von der Fox-Flory-Beziehung zeigte.

---

---

With 10 figures and 1 table     

Symplectic approach to the theory of quantized fields. II

We prove that the set D of vector fields on the configuration space B of a field whose 1-parameter groups locally associated are groups of fibre-preserving transformations of B that leave invariant that field in the sense of variational theory, is a Lie algebra with respect to ordinary addition, multiplication by real numbers and Lie brackets. We see that this Lie algebra structure can be carried over to the corresponding set of Noether invariants, which then becomes a Lie algebra in a natural way.Further, we define the n-form of Poincaré-Cartan of a field, and we use it to generalize the Lie algebras D and in a reasonable way. The algebras D and are subalgebras of the new Lie algebras D and introduced. A main result in this connection is the following: the differential d of the n-form of Poincaré-Cartan is –(d+f), where (, d+f) are the field equations on the vertical bundle B.The symplectic manifold of solutions associated with a field is introduced in a formal way and the former Lie algebras D, , D, are interpreted on this manifold. In imitation of the case of analytical dynamics, the main results in this direction are: a) Every vector field of the Lie algebra D defines, in a canonical way, a vector field on the manifold of solutions such that its polar 1-form with respect to the symplectic metric ₂ is the differential of its corresponding Noether invariant, and b) the Lie bracket [, ] of two Noether invariants , is the Noether invariant given by ₂(D, D), where D, D are the vector fields on the manifold of solutions defined, in the sense a), by two infinitesimal generators of , , respectively. This will allow us to regard the Lie algebra as the analogous object in field theory to the Poisson algebra of analytic dynamics.We apply the general formalism to the relativistic theory of non-linear scalar fields, and we compare our results with the formalism developed by I. Segal for this case.     

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.     

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.