Chaos and integrability in a nonlinear wave equation

We consider the parametrized family of equations _tt ,u- _xx u-au+u ₂ ² u=O,x(0,L), with Dirichlet boundary conditions. This equation has finite-dimensional invariant manifolds of solutions. Studying the reduced equation to a four-dimensional manifold, we prove the existence of transversal homoclinic orbits to periodic solutions and of invariant sets with chaotic dynamics, provided that =2, 3, 4,.... For =1 we prove the existence of infinitely many first integrals pairwise in involution.     

The theory of a vibrating-rod viscometer

The paper presents a complete theory for a viscometer based upon the principle of a circular-section rod, immersed in a fluid, performing transverse oscillations perpendicular to its axis. The theory is established as a result of a detailed analysis of the fluid flow around the rod and is subject to a number of criteria which subsequently constrain the design of an instrument. Using water as an example it is shown that a practical instrument can be designed so as to enable viscosity measurement with an accuracy of ±0.1%, although it is noted that many earlier instruments failed to satisfy one or more of the newly-established constraints.Nomenclature A, D constants in equation (46) - A _m , B _m , C _m , D _m constants in equations (50) and (51) - A _j , B _j constants in equation (14) - a _j ⁺ , a _j ^– wavenumbers given by equation (15) - C _f drag coefficient defined in equation (53) - c speed of sound - D _b drag force of fluid b - D ₀ coefficient of internal damping - E extensional modulus - f(z) initial deformation of rod - f(), F _m () functions of defined in equation (41) - F force in the rod - force per unit length near t=0 - F dimensionless force per unit length near t=0 - g _m amplitude of transient force - G modulus of rigidity - h, h* functions defined by equations (71) and (72) - H functions defined by equation (69) and (70) - I second moment of area - I _0,1, J _0,1, K _0,1 modified Bessel functions - k, k functions defined in equations (2) - L half-length of oscillator - Ma Mach number - m _b added mass per unit length of fluid b - m _s mass per unit length of solid - n _j eigenvalue defined in equations (15) and (16) - R radius of rod - R _c radius of container - 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 - y ₀ initial lateral displacement - y ₁, y ₂ successive maximum lateral displacement - z axial coordinate - dimensionless tension - dimensionless mass of fluid - dimensionless drag of fluid - amplification factor - logarithmic decrement in a fluid - _a , _b logarithmic decrement in fluids a and b - ₀ logarithmic decrement in vacuo - _j logarithmic decrement in mode j in a fluid - spatial resolution of amplitude - _v voltage resolution - r, , , _s, , increments in R, , , _s , , - dimensionless amplitude of oscillation - dimensionless axial coordinate - angular coordinate - _f thermal conductivity of fluid - viscosity of fluid - viscosity of fluid calculated on assumption that * - _a , _b viscosity of fluids a and b - _m constants in equation (10) - dimensionless displacement - _j j the component of - density of fluid - _a , _b density of fluids a and b - _s density of tube or rod material - dimensionless radial coordinate - * dimensionless radius of container - dimensionless times - spatial component of defined in equation (11) - _j , _tm jth, mth component of - dimensionless streamfunction - ⁰, ¹ components of in series expansion in powers of - streamfunction - dimensionless frequency (based on ) - angular frequency - ₀ angular frequency in absence of fluid and internal damping - _j angular frequency in mode j in a fluid - _a , _b frequencies in fluids a and b     

Similar solutions and the asymptotics of filtration equations

In this paper we consider the asymptotic behavior of solutions of the quasilinear equation of filtration as t. We prove that similar solutions of the equation u _t = (u )_xx asymptotically represent solutions of the Cauchy problem for the full equation u _t = [(u)]_xx if (u) is close to u for small u.     

Lp-estimates for the nonlinear spatially homogeneous Boltzmann equation

This paper studies L^p-estimates for solutions of the nonlinear, spatially homogeneous Boltzmann equation. The molecular forces considered include inverse k^th-power forces with k > 5 and angular cut-off.The main conclusions are the following. Let f be the unique solution of the Boltzmann equation with f(v,t)(1 + ¦v²¦)^(s ₁ ^{+ /p)/2} L¹, when the initial value f ₀ satisfies f ₀(v) 0, f ₀(v) (1 + ¦v¦²)^(s ₁ ^{+ /p)/2} L¹, for some s₁ 2 + /p, and f ₀(v) (1 + ¦v¦²)^s/2 L^p. If s 2/p and 1 < p < , then f(v, t)(1 + ¦v¦²)^{(s s} ₁)^/2 L^p, t > 0. If s >2 and 3/(1+ ) < p < , thenf(v,t) (1 + ¦v¦²)^(s(s ₁ ^{+ 3/p))/2} L^p, t > 0. If s >2 + 2C₀/C₁ and 3/(l + ) < p < , then f(v,t)(1 + ¦v¦²)^s/2 L^p, t > 0. Here 1/p + 1/p = 1, x y = min (x, y), and C₀, C₁, 0 < 1, are positive constants related to the molecular forces under consideration; = (k – 5)/ (k – 1) for k^th-power forces.Some weaker conclusions follow when 1 < p 3/ (1 + ).In the proofs some previously known L-estimates are extended. The results for L^p, 1 < p < , are based on these L-estimates coupled with nonlinear interpolation.     

Nonlinear rheological behavior of a concentrated spherical silica suspension

Linear and nonlinear viscoelastic properties were examined for a 50 wt% suspension of spherical silica particles (with radius of 40 nm) in a viscous medium, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture. The effective volume fraction of the particles evaluated from zero-shear viscosities of the suspension and medium was 0.53. At a quiescent state the particles had a liquid-like, isotropic spatial distribution in the medium. Dynamic moduli G^* obtained for small oscillatory strain (in the linear viscoelastic regime) exhibited a relaxation process that reflected the equilibrium Brownian motion of those particles. In the stress relaxation experiments, the linear relaxation modulus G(t) was obtained for small step strain (0.2) while the nonlinear relaxation modulus G(t, ) characterizing strong stress damping behavior was obtained for large (>0.2). G(t, ) obeyed the time-strain separability at long time scales, and the damping function h() (–G(t, )/G(t)) was determined. Steady flow measurements revealed shear-thinning of the steady state viscosity () for small shear rates (< ^–1; = linear viscoelastic relaxation time) and shear-thickening for larger (>^–1). Corresponding changes were observed also for the viscosity growth and decay functions on start up and cessation of flow, + (t, ) and _– (t, ). In the shear-thinning regime, the and dependence of ₊(t,) and _–(t,) as well as the dependence of () were well described by a BKZ-type constitutive equation using the G(t) and h() data. On the other hand, this equation completely failed in describing the behavior in the shear-thickening regime. These applicabilities of the BKZ equation were utilized to discuss the shearthinning and shear-thickening mechanisms in relation to shear effects on the structure (spatial distribution) and motion of the suspended particles.Dedicated to the memory of Prof. Dale S. Parson     

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.     

Mass transport in flow of sources

Steady and unsteady local concentration has been determined analytically for two- und three-dimensional sources and is presented for various boundary-concentrations, volumetric flows and diffusion coefficients. The steady cases have been evaluated numerically. In addition an unsteady two-dimensional mass transport has been evaluated.

---

Stofftransport in Quellströmungen

Zusammenfassung Es wurden die stationäre und instationäre örtliche Konzentration von einer zwei- und drei-dimensionalen Quellströmung als Funktion verschiedener Randkonzentrationen, verschiedener Stromvolumen und Diffusionskoeffizienten analytisch bestimmt. Die stationären Fälle wurden numerisch ausgewertet. Außerdem wurde ein zwei-dimensionaler instationärer Stofftransport behandelt.

---

Nomenclature a inner radius of circle (2-dimensional case), inner radius of sphere (three-dimensional case) - b } >a outer radius of circle (2-dimensional case), outer radius of sphere (three-dimensional case) - c concentration - c ₁,c ₂ given concentration at the boundariesr=a andb resp - c _i initial concentration at the timet=0 - D diffusion coefficient - I _n +1/2 modified spherical Bessel function - J _v ,Y _v Bessel function ofv-th order and first and second kind resp - k =b/a} > 1 diameter ratio - P _n ^o () Legendre polynomials - ¯ r, polar coordinates - r, , spherical coordinates - t time - u velocity in radial direction - V ₀ volumetric flow - ₀ V/4D flow parameter for two-dimensional flow - ₀ V ₀/8 D flow parameter for three-dimensional flow - _mn eigenvalues - _mn te] ² =n ² + ₀ ² ,=cos =r/a roots of determinant (28)     

On nonlinear problems of mixed type: A qualitative theory using infinite-dimensional center manifolds

We consider the equation a(y)u_xx+div_y(b(y)_yu)+c(y)u=g(y, u) in the cylinder (–l,l)×, being elliptic where b(y)>0 and hyperbolic where b(y)<0. We construct self-adjoint realizations in L₂() of the operatorAu= (1/a) div_y(b_yu)+(c/a) in the case ofb changing sign. This leads to the abstract problem u_xx+Au=g(u), whereA has a spectrum extending to + as well as to –. For l= it is shown that all sufficiently small solutions lie on an infinite-dimensional center manifold and behave like those of a hyperbolic problem. Anx-independent cross-sectional integral E=E(u, u_x) is derived showing that all solutions on the center manifold remain bounded forx ±. For finitel, all small solutionsu are close to a solution on the center manifold such that u(x)-(x) Ce ^-(1-|x|) for allx, whereC and are independent ofu. Hence, the solutions are dominated by hyperbolic properties, except close to the terminal ends {±1}×, where boundary layers of elliptic type appear.     

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.     

Existence and regularity for solutions of the Korteweg-de Vries equation

The construction suggested by an inverse-scattering analysis establishes the existence of solutions u(x, t) of the Korteweg-de Vries equation subject to an initial condition u(x, 0)=U(x), where U has certain regularity and decay properties. It is assumed that UC³(), that U is piecewise of class C ⁴, and that U ^(j) decays at an algebraic rate for j4. The faster the decay of U ^(j) the smoother the solution will be for t0. If U and its first four derivatives decay faster than ¦x¦^–n for all n, then the solution will be infinitely differentiable for t0. For t>0, the decay rate of u(x, t) as x + increases with the decay rate of U; but the decay rate as x - depends on the regularity of U. A solution u ₁ of the Korteweg-de Vries equation such that u ₁(·, 0)C() may fail to remain in class C for all time if u ₁(x, 0) does not decay fast enough as ¦x¦.This research was performed in part as a Visiting Member of the Courant Institute of Mathematical Science.     

Analysis of melt spinning transients in Lagrangian coordinates

Transients in melt spinning of isothermal power law and Newtonian fluids were found to be governed by an extremely simple partial differential equation ² ( ^1/n )/() = 0 in Lagrangian coordinates where is the cross-sectional area,n the power law exponent, the time and the the time at which a fluid molecule constituting the spinline left the spinneret. The general integral ^1/n =f() +g () of the above governing equation containing two arbitrary functions represents physically attainable spinline transients. Hitherto unknown analytical transient solutions of the above governing equation were obtained for the response of isothermal constant tension spinlines to a stepwise change in tension, spinneret hole area, extrusion speed or extrusion viscosity and for the starting transient in gravitational spinning. Linearized perturbation solutions and the stability limit of the spinline derived from the above new found nonlinear solutions were in agreement with previous findings and the above nonlinear response of the spinline to a step increase in the spinneret hole area was found to be equivalent to Orowan's tandem cylinder model of dent growth in filament stretching.     

Exact Solutions of the Hydrodynamic Equations Derived from Partially Invariant Solutions

The paper proposes a heuristic approach to constructing exact solutions of the hydrodynamic equations based on the specificity of these equations. A number of systems of hydrodynamic equations possess the following structure: they contain a reduced system of n equations and an additional equation for an extra function w. In this case, the reduced system, in which w = 0, admits a Lie group G. Taking a certain partially invariant solution of the reduced system with respect to this group as a seed:rdquo; solution, we can find a solution of the entire system, in which the functional dependence of the invariant part of the seed solution on the invariants of the group G has the previous form. Implementation of the algorithm proposed is exemplified by constructing new exact solutions of the equations of rotationally symmetric motion of an ideal incompressible liquid and the equations of concentrational convection in a plane boundary layer and thermal convection in a rotating layer of a viscous liquid.     

The Thermodynamic Significance of the Local Volume Averaged Temperature

Since the temperature is not an additive function, the traditional thermodynamic point of view suggests that the volume integral of the temperature has no precise physical meaning. This observation conflicts with the customary analysis of non-isothermal catalytic reactors, heat pipes, driers, geothermal processes, etc., in which the volume averaged temperature plays a crucial role. In this paper we identify the thermodynamic significance of the volume averaged temperature in terms of a simple two-phase heat transfer process. Given the internal energy as a function of the point temperature and the density

Extended Korn's inequalities and the associated best possible constants

We are concerned with the coerciveness of the strain energy E(u) (in linear elasticity) associated with a displacement vector u on the Sobolev space H¹ () or its subspaces, a domain in ⁿ representing an isotropic elastic body—certain specific cases are called Korn's inequalities. Sufficient (and necessary) conditions on the Lamé moduli for E(·) to be coercive (or uniformly positive) on such spaces are given, and the associated best possible constants are obtained for some cases.     

A light transmittance probe for interfacial area measurements of dispersed two phase flows with small particles

An optical probe measuring interfacial area () by light attenuation has been designed with a special emphasis on flows with sub-millimetric particles. It permits measurements in liquid-liquid or gas-liquid dispersions without need of introducing empirical correcting factors for the standard exponential decay law of light intensity while keeping an extended application range. This probe was successfully tested with an air-glass particle flow, the parameters of which were carefully determined basically by hold-up methods. The volume fraction of the dispersed phase was varied between 0.05% and 5%, and the particle size between 10 m and 300 m.List of symbols D diameter of spherical particle - D _S Sauter diameter - E ₀ irradiance on a surface perpendicular to light propagation 226E;=(1/l) averaged density function along y axis - f density function of a dispersion - f ₁, f ₂ focal length of the lenses L ₁, L ₂ - g granulometry function of a powder (probability density) - h granulometry function of a powder (unnormalized) - I ₀, I light beam intensity respectively before and after passing through the dispersion - j volumetric powder flow - K ₁, K ₂, K ₃ dimensionless constants - l optical path length of the beam in the dispersion - L experimental pipe width along x axis - m mass of a sample - n optical index of the continuous phase - p _a, p ₀ respectively slope of _a and ₀ straight line - r distance between particles - S _d scattering cross-section - V volume of dispersion - averaged particles velocity - x, y, z spatial coordinates - interfacial area - _a absolute interfacial area (by unit volume of dispersion) - ₀ interfacial area measured by light attenuation method - _d angle (around the initial direction of light propagation) within which a particle diffracts - _dr detector aperture angle - light wavelength - _d scattering cross section by unit volume of dispersion - light beam diameter - ₁, ₂ L₁, and L₂ lenses diameters - local volumetric fraction of dispersed phase - averaged fraction of dispersed phase along x axis - ₂ averaged fraction of dispersed phase along x and y axis - volumetric mass of particles     

A visualization study of the interaction of a free vortex with the wake behind an airfoil

The two-dimensional interaction of a single vortex with a thin symmetrical airfoil and its vortex wake has been investigated in a low turbulence wind tunnel having velocity of about 2 m/s in the measuring section. The flow Reynolds number based on the airfoil chord length was 4.5 × 10³. The investigation was carried out using a smoke-wire visualization technique with some support of standard hot-wire measurements. The experiment has proved that under certain conditions the vortex-airfoil-wake interaction leads to the formation of new vortices from the part of the wake positioned closely to the vortex. After the formation, the vortices rotate in the direction opposite to that of the incident vortex.List of symbols c test airfoil chord - C vortex generator airfoil chord - TA test airfoil - TE test airfoil trailing edge - TE _G vortex generator airfoil trailing edge - t dimensionless time-interval measured from the vortex passage by the test airfoil trailing edge: gDt=(T-T- _TE)·U/c - T time-interval measured from the start of VGA rotation - U free stream velocity - U vortex induced velocity fluctuation - VGA vortex generator airfoil - y distance in which the vortex passes the test airfoil - Z vortex circulation coefficient: Z=/(U · c/2) - vortex generator airfoil inclination angle - vortex circulation - vortex strength: =/2     

On the Discretization of Distributed-Parameter Systems with Quadratic and Cubic Nonlinearities

Approximate methods for analyzing the vibrations of an Euler--Bernoulli beam resting on a nonlinear elastic foundation are discussed. The cases of primary resonance ( _n ) and subharmonic resonance of order one-half ( 2 _n ), where is the excitation frequency and _n is the natural frequency of the nth mode of the beam, are investigated. Approximate solutions based on discretization via the Galerkin method are contrasted with direct application of the method of multiple scales to the governing partial-differential equation and boundary conditions. The amplitude and phase modulation equations show that single-mode discretization leads to erroneous qualitative as well as quantitative predictions. Regions of softening (hardening) behavior of the system, the spatial dependence of the response drift, and frequency-response curves are numerically evaluated and compared using both approaches.     

The value of the critical exponent for reaction-diffusion equations in cones

Let D R ^N be a cone with vertex at the origin i.e., D = (0, )x where S ^N–1 and x D if and only if x = (r, ) with r=¦x¦, . We consider the initial boundary value problem: u _t = u+u ^p in D×(0, T), u=0 on Dx(0, T) with u(x, 0)=u ₀(x) 0. Let ₁ denote the smallest Dirichlet eigenvalue for the Laplace-Beltrami operator on and let ₊ denote the positive root of (+N–2) = ₁. Let p ^* = 1 + 2/(N + ₊). If 1 < p < p ^*, no positive global solution exists. If p>p ^*, positive global solutions do exist. Extensions are given to the same problem for u _t=+¦x¦ u ^p .This research was supported in part by the Air Force Office of Scientific Research under Grant # AFOSR 88-0031 and in part by NSF Grant DMS-8 822 788. The United States Government is authorized to reproduce and distribute reprints for governmental purposes not withstanding any copyright notation therein.     

Optical co-ordinates in general relativity 2 — the problem of distance

Summary Let denote the congruence of null geodesics associated with a given optical observer inV ₄. We prove that determines a unique collection of vector fieldsM₍₎ ( =1, 2, 3) and ₍₀₎ overV ₄, satisfying a weak version of Killing's conditions.This allows a natural interpretation of these fields as the infinitesimal generators of spatial rotations and temporal translation relative to the given observer. We prove also that the definition of the fieldsM₍₎ and ₍₀₎ is mathematically equivalent to the choice of a distinguished affine parameter f along the curves of, playing the role of a retarded distance from the observer.The relation between f and other possible definitions of distance is discussed.

---

Sommario Sia la congruenza di geodetiche nulle associata ad un osservatore ottico assegnato nello spazio-tempoV ₄. Dimostriamo che determina un'unica collezione di campi vettorialiM₍₎ ( =1, 2, 3) e ₍₀₎ inV ₄ che soddisfano una versione in forma debole delle equazioni di Killing. Ciò suggerisce una naturale interpretazione di questi campi come generatori infinitesimi di rotazioni spaziali e traslazioni temporali relative all'osservatore assegnato. Dimostriamo anche che la definizione dei campiM₍₎, ₍₀₎ è matematicamente equivalente alla scelta di un parametro affine privilegiato f lungo le curve di, che gioca il ruolo di distanza ritardata dall'osservatore. Successivamente si esaminano i legami tra f ed altre possibili definizioni di distanza in grande.

---

---

Work performed in the sphere of activity of: Gruppo Nazionale per la Fisica Matematica del CNR.     

Impulsive motion of a semi-infinite flat plate in a viscous fluid

The unsteady laminar boundary layer flow is investigated for a semi-infinite flat plate subjected to impulsive motion. An approximate solution is obtained by utilizing Meksyn's method. These results vary smoothly from Rayleigh's unsteady solution to the steady state solution of Blasius. Results are compared to those of Lam and Crocco.Nomenclature A expansion coefficient, see eq. (13) - a expansion coefficient, see eq. (10) - B expansion coefficient, see eq. (14) - b expansion coefficient, see eq. (12) - G function defined by eq. (6) - U free stream velocity - u velocity in x direction - v velocity in y direction - x coordinate along plate - y coordinate normal to plate Greek symbols (l, ) incomplete gamma function - function defined by eq. (15) - y(U/x) ^1/2 - kinematic viscosity - x/Ut - (Uvx)^1/2 f(, )