Effects of Multi Global Bifurcations on Basin Organization,Catastrophes and Final Outcomes in a Driven Nonlinear Oscillator at the 2T-Subharmonic Resonance

The behavior of the escape driven oscillator at the 2T-periodic subharmonic resonance is considered, and the mechanism of generating different fractal patterns of the basins of attraction of coexisting attractors, as well as its effects on the unpredictable asymptotic system behaviors, are the main points of interest. The analysis is based on the numerical study of the sudden qualitative changes of the structure of basin-phase portraits, the changes implied by multi global bifurcations. Attention is focused on two qualitatively different regions of control space: the region prior to the subcritical flip bifurcation, where all three attractors (2T-periodic, T-periodic and the attractor at infinity) coexist, and the region after the bifurcation, where only two attractors (2T-periodic and the attractor at infinity) coexist. In particular, the concept of the global (homoclinic and heteroclinic) bifurcations is extended to the latter region, where the arising flip saddle (instead of the direct saddle) is involved in the events. The possible forms of unpredictable outcomes, which arise in both regions of control parameters, are pointed out.     

Ši’lnikov Chaos in the Generalized Lorenz Canonical Form of Dynamical Systems

This paper studies the generalized Lorenz canonical form of dynamical systems introduced by elikovský and Chen [International Journal of Bifurcation and Chaos 12(8), 2002, 1789]. It proves the existence of a heteroclinic orbit of the canonical form and the convergence of the corresponding series expansion. The ilnikov criterion along with some technical conditions guarantee that the canonical form has Smale horseshoes and horseshoe chaos. As a consequence, it also proves that both the classical Lorenz system and the Chen system have ilnikov chaos. When the system is changed into another ordinary differential equation through a nonsingular one-parameter linear transformation, the exact range of existence of ilnikov chaos with respect to the parameter can be specified. Numerical simulation verifies the theoretical results and analysis.     

Flow through a helically coiled annulus

In this paper the flow is studied of an incompressible viscous fluid through a helically coiled annulus, the torsion of its centre line taken into account. It has been shown that the torsion affects the secondary flow and contributes to the azimuthal component of velocity around the centre line. The symmetry of the secondary flow streamlines in the absence of torsion, is destroyed in its presence. Some stream lines penetrate from the upper half to the lower half, and if is further increased, a complete circulation around the centre line is obtained at low values of for all Reynolds numbers for which the analysis of this paper is valid, being the ratio of the torsion of the centre line to its curvature.Nomenclature A =constant - a outer radius of the annulus - b unit binormal vector to C - C helical centre line of the pipe - D rL - g 1000 - K Dean number=Re² - L 1+r sin - M (L ²+ ² r ²)^1/2 - n unit normal vector to C - P, P pressure and nondimensional pressure - p ₀, p pressures of O(1) and O() - Re Reynolds number=aW ₀/ - (r, , s), (r, , s) coordinates and nondimensional coordinates - nonorthogonal unit vectors along the coordinate directions - r ₀ radius of the projection of C - t unit tangent vector to C - V _r, V , V _s velocity components along the nonorthogonal directions - V_r, V, V _s nondimensional velocity components along - W ₀ average velocity in a straight annulus Greek symbols , curvature and nondimensional curvature of C - U, V, W lowest order terms for small in the velocity components along the orthogonal directions t - _r, , _s first approximations to V _r , V, V _s for small - =/=/ - kinematic viscosity - density of the fluid - , torsion and nondimensional torsion of C - , stream function and nondimensional stream function - nondimensional streamfunction for U, V - a inner radius of the annulus After this paper was accepted for publication, a paper entitled On the low-Reynolds number flow in a helical pipe, by C.Y. Wang, has appeared in J. Fluid. Mech., Vol 108, 1981, pp. 185–194. The results in Wangs paper are particular cases of this paper for =0, and are also contained in [9].     

A Capillary Microstructure of the Wetting Front

This article reports the experimental results of a study of the wetting-front microscale structure formed only by capillary forces in homogeneous and random etched glass capillary models. In the homogeneous model, water propagates through the capillary system, evenly filling the capillaries across the direction of flow. Air is trapped by the pinch-off mechanism inside the pore bodies in the form of individual bubbles. The experiments specified three consecutive steps of the pinch-off mechanism, film flow, snap-off, and interface movement. In the random model, both the bypass and pinch-off, forming bypass/cut-off mechanism, create residual air structure. Bypass traps air inside large capillary-pore aggregates which are bounded by small-diameter capillaries in where pinch-off traps air in the adjacent pores. An analysis of the residual air distribution versus depth below the surface in the homogeneous and random micromodels made it possible to identify three successive zones, namely a transition zone, a transmission zone, and a wetting-and-front zone. In the transition zone, the residual air content increases with depth from zero to the constant value in the transmission zone where it remains practically constant. The capillary processes within the wetting-and-front combined zone govern air replacement with wetting and formation of the transmission zone.     

Travelling, non-periodic patterns in nonlinear stability problems

In this short note we present results on the existence of several classes of travelling, non-periodic solutions of the complex Ginzburg-Landau equation. First we give a very short introduction to the G-L equation and show its importance in nonlinear stability theory. We then study the G-L equation with complex coefficients and establish the existence of a 2-parameter family of quasi-periodic solutions and two different types of one-parameter families of heteroclinic orbits; all members of these families travel with a well-defined wave-speed. The heteroclinic solutions correspond to (travelling) soliton-like localized structures which connect different (stable) periodic patterns. Mathematically, these families of travelling solutions (quasi-periodic and heteroclinic) are continuations into the complex case of the stationary solutions of the real G-L equation.     

Properties of Chaotic and Regular Boundary Crisis in Dissipative Driven Nonlinear Oscillators

The phenomenon of the chaotic boundary crisis and the related concept of the chaotic destroyer saddle has become recently a new problem in the studies of the destruction of chaotic attractors in nonlinear oscillators. As it is known, in the case of regular boundary crisis, the homoclinic bifurcation of the destroyer saddle defines the parameters of the annihilation of the chaotic attractor. In contrast, at the chaotic boundary crisis, the outset of the destroyer saddle which branches away from the chaotic attractor is tangled prior to the crisis. In our paper, the main point of interest is the problem of a relation, if any, between the homoclinic tangling of the destroyer saddle and the other properties of the system which may accompany the chaotic as well as the regular boundary crisis. In particular, the question if the phenomena of fractal basin boundary, indeterminate outcome, and a period of the destroyer saddle, are directly implied by the structure of the destroyer saddle invariant manifolds, is examined for some examples of the boundary crisis that occur in the mathematical models of the twin-well and the single-well potential nonlinear oscillators.     

Shear viscosity behavior of highly concentrated suspensions at low and high shear-rates

A method is proposed for determining the shear viscosity behavior of highly concentrated suspensions at low and high shear-rates through the use of a formulation that is a function of three parameters signifying the effects of particle size distribution. These parameters are the intrinsic viscosity [], a parametern that reflects the level of particle association at the initiation of motion and the maximum packing concentration _m. The formulation reduces to the modified Eilers equation withn = 2 for high shear rates. An analytical method was used for the calculation of maximum packing concentration which was subsequently correlated with the experimental values to account for the surface induced interaction of particles with the fluid. The calculated values of viscosities at low and high shear-rates were found to be in good agreement with various experimental data reported in literature. A brief discussion is also offered on the reliability of the methods of measuring the maximum packing concentration. _r = /₀ relative viscosity of the suspension - volumetric concentration of solids - k _n coefficient which characterizes a specific effect of particle interactions - _m maximum packing concentration - _r,0 relative viscosity at low shear-rates - [] intrinsic viscosity - n, n parameter that reflects the level of particle interactions at low and high shear-rates, respectively - _r, relative viscosity at high shear-rates - (_m)_s, (_m)_i, (_m)_l packing factors for small, intermediate and large diameter classes - v _s, v_i, v_l volume fractions of small, intermediate and large diameter classes, respectively - _si, _sl coefficient to be used in relating a smaller to an intermediate and larger particle group, respectively - _is, _il coefficient to be used in relating an intermediate to a smaller and larger particle group, respectively - _ls, _li coefficient to be used in relating a larger to a smaller and intermediate particle group, respectively - _m0 maximum packing concentration for binary mixtures - _m,e measured maximum packing concentration - _m,c calculated maximum packing concentration     

Finite Elastic Deformations of a Tyre Modelled as an Ideal Fibre-Reinforced Shell

Vehicle tyres are anisotropic inhomogeneous fibre-reinforced shells which undergo finite elastic deformations. Calculation of their stress and deformation fields is a difficult task and is normally performed using the finite element technique. In this paper an attempt is made to provide an approximate analysis of the deformation field modelling the tyre as an ideal fibre-reinforced material. Radial-ply tyres are reinforced by a belt of fibres running around the wheel in the circumferential direction under the tread of the tyre. A second set of fibres lies in each radial cross-section, of the tyre and runs from the bead wire which seats against one wheel rim to the bead wire at the other wheel rim. We shall assume each radial cross-section of the tyre is in a state of plane strain and is formed from an arch of fibre-reinforced composite material which is reinforced in the hoop direction. This composite is assumed to be an ideal material which is inextensible in the fibre-direction and is incompressible. The plane-strain deformations of this section are examined and then used to analyse the deformation of the tyre as a whole.     

Elliptic Two-Dimensional Invariant Tori for the Planetary Three-Body Problem

The spatial planetary three-body problem (i.e., one star and two planets, modelled by three massive points, interacting through gravity in a three dimensional space) is considered. It is proved that, near the limiting stable solutions given by the two planets revolving around the star on Keplerian ellipses with small eccentricity and small non-zero mutual inclination, the system affords two-dimensional, elliptic, quasi-periodic solutions, provided the masses of the planets are small enough compared to the mass of the star and provided the osculating Keplerian major semi-axes belong to a two-dimensional set of density close to one.     

Bifurcation and Chaos in an Apparent-Type Gyrostat Satellite

Attitude dynamics of an asymmetrical apparent gyrostat satellite has been considered. Hamiltonian approach and Routhian are used to prove that the dynamics of the system consists of two separate parts, an integrable and a non-integrable. The integrable part shows torque free motion of gyrostat, while the non-integrable part shows the effect of rotation about the earth and asphericity of the satellites inertia ellipsoid. Using these results, theoretically when the non-integrable part is eliminated, we are able to design a satellite with exactly regular motion. But from the engineering point of view the remaining errors of manufacturing process of the mechanical parts cause that the non-integrable part can not be eliminated, completely. So this case can not be achieved practically. Using Serret–Andoyer canonical variable the Hamiltonian transformed to a more appropriate form. In this new form the effect of the gravity, asphericity, rotational motion and spin of the rotor are explicitly distinguished. The results lead us to another way of control of chaos. To suppress the chaotic zones in the phase space, higher rotational kinetic energy can be used. Increasing the parameter related to the spin of the rotor causes the systems phase space to pass through a heteroclinic bifurcation process and for the sufficiently large magnitude of the parameter the heteroclinic structure can be eliminated. Local bifurcation of the phase space of the integrable part and global heteroclinic bifurcation of whole systems phase space are presented. The results are examined by the second order Poincaré surface of section method as a qualitative, and the Lyapunov characteristic exponents as a quantitative criterion.     

Heteroclinic orbits in singular systems: A unifying approach

We consider singularly perturbed systems , such that=f(, _o, 0). _o ^m , has a heteroclinic orbitu(t). We construct a bifurcation functionG(, ) such that the singular system has a heteroclinic orbit if and only ifG(, )=0 has a solution=(). We also apply this result to recover some theorems that have been proved using different approaches.     

A Direct Method for the Identification of the Permeability Field Based on Flux Assimilation by a Discrete Analog of Darcy's Law

Inverse models to determine the permeability are generally based on existing forward models for the pressure. The permeabilities are adapted in such a way that the calculated pressures match the specified pressures in a number of points. To assimilate a priori knowledge about the flux, we introduce the flux assimilation method, which is based on the vector potential–pressure formulation of Darcy's law. Thanks to an unconventional discretization technique – the edge-based face element method – not only the specified pressures, but also specified information about the flux density can easily be assimilated. A relatively simple, but insightful analytical example illustrates the potential of this method.     

Experimental Investigation of the Nonlinear Response of a Hanging Cable. Part II: Global Analysis

The nonlinear dynamics of the same experimental model of an internally resonant hanging elastic cable considered in Part I [1] are addressed here from the point of view of the global system behaviour in the control parameter space. Synthetic results of systematic response measurements, made at different amplitudes of the support motion in frequency ranges including meaningful external resonance conditions, are reported and discussed. Attention is devoted to the detection of the most robust classes of motion. Quite complicated overall pictures of regular response regions with variable contributions from different planar and nonplanar cable modes are observed, as well as several regions of quasiperiodic and chaotic responses. Sample quantitative characterizations of nonregular motions are presented. Some experimental results are also observed against the background of the nonlinear dynamic phenomena exhibited by a theoretical model of a continuous cable with four-degrees-of-freedom.     

Permeability Effects of Turbulent Flow Through a Porous Insert in a Backward-Facing-Step Channel

In the present work, a k– model, based on the work of Lee and Howell (Proceedings of the ASME-JSME Thermal Engineering Hawaii, 1987), is rigorously derived based on time average of spatially averaged Navier–Stokes equations. The model is then employed to solve for a flow in a backward-facing step channel with a porous insert. The numerical solver is modified from the STREAM code (Lien and Leschziner, Comput. Meth. Appl. Mech. Eng. 114 (1994a) 123–148), and it has been validated against the experimental data of Seegmiller and Driver (AIAA Journal 23 (1985) 163–171). The code is then used to perform simulation for cases with a porous insert. The resistance of the porous insert can be altered by changing its permeability (), Forchheimers constant (F), or thickness (b). The goal is to examine the influence of each parameter on the resulting flow and turbulent kinetic energy (k) distributions. It is discovered that, by increasing the resistance of the insert, flow eventually enters a transitional regime towards relaminarization. This is due to the contribution of Darcys and Forchheimers terms in the governing equations, and modifying these two terms changes the levels of P_k and, hence, k and . Generally speaking, lowering or raising F results in a greater suppression of P_k than , causing the flow to relaminarize. Meanwhile, if the pore size is reasonably large to sustain turbulence within the porous media, increasing b reduces but does not eliminate the turbulent activity in the porous insert.     

Inhomogeneous ‘longitudinal’ plane waves in a deformed elastic material

By definition, a homogeneous isotropic compressible Hadamard material has the property that an infinitesimal longitudinal homogeneous plane wave may propagate in every direction when the material is maintained in a state of arbitrary finite static homogeneous deformation. Here, as regards the wave, homogeneous means that the direction of propagation of the wave is parallel to the direction of eventual attenuation; and longitudinal means that the wave is linearly polarized in a direction parallel to the direction of propagation. In other words, the displacement is of the form u = ncos k(n · x – ct), where n is a real vector. It is seen that the Hadamard material is the most general one for which a longitudinal inhomogeneous plane wave may also propagate in any direction of a predeformed body. Here, inhomogeneous means that the wave is attenuated, in a direction distinct from the direction of propagation; and longitudinal means that the wave is elliptically polarized in the plane containing these two directions, and that the ellipse of polarization is similar and similarly situated to the ellipse for which the real and imaginary parts of the complex wave vector are conjugate semi-diameters. In other words, the displacement is of the form u = {S exp i(S · x – ct)}, where S is a complex vector (or bivector). Then a Generalized Hadamard material is introduced. It is the most general homogeneous isotropic compressible material which allows the propagation of infinitesimal longitudinal inhomogeneous plane circularly polarized waves for all choices of the isotropic directional bivector. Finally, the most general forms of response functions are found for homogeneously deformed isotropic elastic materials in which longitudinal inhomogeneous plane waves may propagate with a circular polarization in each of the two planes of central circular section of the ⁿ -ellipsoid, where is the left Cauchy-Green strain tensor corresponding to the primary pure homogeneous deformation.     

Critique of the computational Preston tube method

The general validity of the Computational Preston Tube Method (CPM) proposed by Nitsche et al. is disputed. The method involves the determination of both the skin-friction coefficient c _f and the von Karman constant from fitting two near-wall mean velocity measurements to a generalized van Driest family of velocity profiles. It is demonstrated by a detailed examination of measured mean velocity profiles that for at least one flow situation for which success with the CPM has been claimed, the laminar-turbulent transition of a flat-plate boundary layer, the method must either fail completely or lead to non-unique results. It is concluded that the applicability of the CPM is restricted to flows for which the normal law-of-the-wall is valid.     

Alloy separation of a binary mixture in a stressed elastic sphere

We consider the static state of a spherical isotropic binary elastic solid mixture whose boundary is given a uniform radial displacement. The elastic volumetric strain energy is given by the classical quadratic form from linear elasticity theory,% MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqef0uAJj3BZ9Mz0bYu% H52CGmvzYLMzaerbd9wDYLwzYbItLDharqqr1ngBPrgifHhDYfgasa% acOqpw0xe9v8qqaqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbba9q8Wq% Ffea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qqQ8frFve9Fve9Ff0dme% GabaqaaiGacaGaamqadaabaeaafiaakeaacaqGxbGaeyypa0ZaaSaa% aeaacaGGXaaabaGaaiOmaaaacaGG7bGaaCOUdiaacIcacaqGJbGaam% ykaiaadIcacaWG0bGaamOCaiaadwgacaWGPaWaaWbaaSqabeaacaGG% YaaaaOGaey4kaSIaaiOmaiabeY7aTjaacIcacaqGJbGaaiykaiaacY% hacaWGLbGaeyOeI0IaaiiiamaalaaabaGaamymaaqaaiaadodaaaGa% aiikaiaadshacaWGYbGaamyzaiaacMcacaWGXaGaaiiFamaaCaaale% qabaGaaiOmaaaakiaac2hacaGGUaaaaa!63E0!\[{\text{W}} = \frac{1}{2}\{ {\mathbf{\kappa }}({\text{c}})(tre)^2 + 2\mu ({\text{c}})|e - \frac{1}{3}(tre)1|^2 \} .\]Here, e is the infinitesimal strain tensor, c[0, 1] is the volumetric concentration of the mixture, and (·) and (·) are the (positive) bulk and shear material moduli, respectively, which are given functions of the concentration. As a function of c and e, the strain energy function is generally nonconvex. Thus, we consider the nonconvex problem of minimizing the potential energy of the body, among all spatial concentration and displacement fields, subject to a given boundary displacement and a fixed amount of component materials. Assuming spherical symmetry, we find that the two component materials must be separated in the optimal state of minimum potential energy. The harder material forms the central core of the sphere, and the softer material is segregated into a surrounding shell. This behavior is remindful of a general notion in metallurgy that in the casting of materials the harder material tends to migrate toward the center.Partial support of the NSF under grant MSS-9024637 and Alliant Techsystems Inc. is gratefully acknowledged.Professor R. Bartelletti of the Università di Pisa is gratefully acknowledged.     

Existence and Uniqueness of Time-Periodic Physically Reasonable Navier-Stokes Flow Past a Body

Let be a three-dimensional exterior domain of class C^2,, 0<<1. Assume that a Navier-Stokes liquid is moving in under the action of a body force F that is time-periodic of period T, and that the velocity of the liquid is zero at spatial infinity. In this paper we show that, if F satisfies suitable conditions, and its norm, in appropriate function spaces, is sufficiently small, there is at least one time-periodic strong solution. Furthermore, the velocity field v of such a solution decays to zero for large |x| as |x|^–1 and its spatial gradient decays as |x|^–2, both uniformly in time. In addition, the pressure p decays like |x|^–2 and its gradient like |x|^–3, for almost all t[0,T]. In the special case where F is time-independent, these solutions are also time-independent and coincide with Finns physically reasonable solutions [4]. Moreover, we show that our time-periodic solutions are unique in a very large class, namely, the class of time-periodic weak solutions satisfying the energy inequality and with corresponding pressure fields verifying mild summability conditions in ×[0,T].     

Numerische Untersuchung des elastisch-plastischen Verhaltens eines Fadenverbundwerkstoffes mit metallischer Matrix Eine Fallstudie

Übersicht Die Fehlervorhersage eines einachsigen Fadenverbundwerkstoffes mit metallischer Matrix erfordert die Kenntnis des nichtlinearen elastisch-plastischen Verhaltens in mikroskopischer Abmessung. Als Fallstudie wurde ein spezieller Verbundwerkstoff mit einem FEM-Programm untersucht. Im plastischen Bereich wurde das klassische, v. Mises Potential' mit dem neuen Übergangsfließpotential unter Berücksichtigung der plastischen Volumendehnung verglichen. Unter transversaler Normalbelastung zeigte der Verbund deutliche Unterschiede in der Zunahme und der Ausdehnung der lokalen plastischen Zonen. Alle kritischen Verformungsstadien wurden von dem Übergangsfließpotential bei geringerer Belastung erreicht.

---

Numerical investigation of the elastic-plastic behaviour of a fibre-reinforced composite with a metallic matrix

Summary Failure prediction of unidirectional fibre-reinforced composite with a metallic matrix needs knowledge about the nonlinear elastic-plastic behaviour in a microscopic scale. A specific composite was investigated using a FEM-program as a case study. In the plastic range the classical v. Mises Potential was compared with the new Transition Flow Potential, taking into account the plastic volume dilatation. Subjected to transverse normal loading the composite showed evident differences in the increase and the spread of locally plastic regions. All critical deformation states were reached by the Transition Flow Potential at lower loading.

     

Analytic Solution to a Problem of Seepage in a Chequer-Board Porous Massif

A study is made of steady two-dimensional seepage in a porous massif composed by a double-periodic system of white and black chequers of arbitrary conductivity. Rigorous matching of Darcy's flows in zones of different conductivity is accomplished. Using the methods of complex analysis, explicit formulae for specific discharge are derived. Stream lines, travel times, and effective conductivity are evaluated. Deflection of marked particles from the natural direction of imposed gradient and stretching of prescribed composition of these particles enables the elucidation of the phenomena of transversal and longitudinal dispersion. A model of pure advection is related with the classical one-dimensional vective dispersion equation by selection of dispersivity which minimizes the difference between the breakthrough curves calculated from the two models.