The Extended Malkus–Robbins Dynamo as a Perturbed Lorenz System

Recent investigations of some self-exciting Faraday-disk homopolar dynamos [Hide, R. and Moroz, I. M., Physica D 134, 1999, 387–301; Moroz, I. M. and Hide, R., International Journal of Bifurcation and Chaos 2000, 2701–2716; Moroz, I. M., International Journal of Bifurcation and Chaos 13, 2003, 147–161; Moroz, I. M., International Journal of Bifurcation and Chaos, to appear] have yielded the classic Lorenz equations as a special limit when one of the principal bifurcation parameters is zero. In this paper we focus upon one of those models [Moroz, I. M., International Journal of Bifurcation and Chaos 13, 2003, 147–161] and illustrate what happens to some of the lowest order unstable periodic orbits as this parameter is increased from zero.     

Performances of Unknown Input Observers for Chaotic LPV Maps in a Stochastic Context

Jointly modeling chaotic maps as Linear Parameter Varying (LPV) systems and using Unknown Input Observers (UIOs) for retrieving the information in a secure communication scheme has previously been motivated in a deterministic context [Millerioux, G. and Daafouz, J. International Journal of Bifurcation and Chaos 14(4), 2004]. In this paper, some new theoretical results from a control theory point of view, concerning the design in a stochastic and so more realistic context of UIOs for chaotic LPV systems is provided. The design of such observers is expressed in terms of the resolution of a finite set of Matrix Inequalities constraints and guarantees some prescribed performances on the state reconstruction error.     

Chaotification for switched linear systems with controllers

This paper shows that two or more switched linear systems can generate chaotic dynamical behaviors by an appropriate switching rule as they at least consist of a controllable system and an unstable system with the expanding property. According to the results in the reference （Xie, L. L., Zhou, Y., and Zhao, Y. Criterion of chaos for switched linear systems with antrollers. International Journal of Bifurcation and Chaos, 20（12）, 4105-4109 （2010））, a nonlinear feedback gain is needed to generate chaotic dy- namics. A linear feedback control is usually used to approximate the nonlinear one for simulation. In order to obtain the exact control, as a main result of this paper, the con- troller is constructed by Russell＇s result, and a block diagram is included to interpret the realization of the controller. Numerical simulations are given to illustrate the generated chaotic dynamical behavior of the switched linear systems with some parameters and show the effects of the constructed controller.     

Multibody Formulation with Large Bending Finite Elements (LBFE)

The goal of this paper is to present a flexible multibody formulation for Euler-Bernoulli beams involving large displacements. This method is based on a discretisation of internal and kinetic energies. The beam is represented by its line of centroids and each section is oriented by a frame defined by three Euler angles. We apply a finite element formulation to describe the evolution of these angles along the neutral fibre and describe the internal energy. The kinetic energy is approximated as the one of two rigid bars tangent to the neutral fibre at the ends of the beam. We derive the equations of motion from a Lagrange formulation. These equations are solved using the Newmark method or/and the Newton-Raphson technique. We solve some very classic problems taken from the literature as the curved beam presented by Simo [Simo, J. C., ‘A three-dimensional finite-strain rod model. the three-dimensional dynamic problem. Part I’, Comput. Meths. Appl. Mech. Engrg. 49, 1985, 55–70; Simo, J. C. and Vu-Quoc, L., ‘A three-dimensional finite-strain rod model, Part II: Computationals aspects’, Comput. Meths. Appl. Mech. Engrg. 58, 1988, 79–116] and Lee [Lee, Kisu, ‘Analysis of large displacements and large rotations of three-dimensional beams by using small strains and unit vectors’, Commun. Numer. Meth. Engrg. 13, 1997, 987–997] or the rotational rod presented by Avello [Avello, A., Garcia de Jalon, J., and Bayo, E., ‘Dynamics of flexible multibody systems using cartesian co-ordinates and large displacement theory’, Int. J. Num. Methods in Engineering 32, 1991, 1543–1563] and Simo [Simo, J. C. and Vu-Quoc, L., ‘On the dynamics of flexible beams under large overall motions – the planar case. Part I’ Jour. of Appl. Mech. 53, 1986, 849–854; Simo, J. C. and Vu-Quoc, L., ‘On the dynamics of flexible beams under large overall motions – the planar case. Part II’, Jour. of Appl. Mech. 53, 1986, 855–863].     

Effect of rigid boundary on propagation of torsional surface waves in porous elastic layer

The paper presents the effect of a rigid boundary on the propagation of torsional surface waves in a porous elastic layer over a porous elastic half-space using the mechanics of the medium derived by Cowin and Nunziato (Cowin, S. C. and Nunziato, J. W. Linear elastic materials with voids. Journal of Elasticity, 13(2), 125–147 (1983)). The velocity equation is derived, and the results are discussed. It is observed that there may be two torsional surface wave fronts in the medium whereas three wave fronts of torsional surface waves in the absence of the rigid boundary plane given by Dey et al. (Dey, S., Gupta, S., Gupta, A. K., Kar, S. K., and De, P. K. Propagation of torsional surface waves in an elastic layer with void pores over an elastic half-space with void pores. Tamkang Journal of Science and Engineering, 6(4), 241–249 (2003)). The results also reveal that in the porous layer, the Love wave is also available along with the torsional surface waves. It is remarkable that the phase speed of the Love wave in a porous layer with a rigid surface is different from that in a porous layer with a free surface. The torsional waves are observed to be dispersive in nature, and the velocity decreases as the oscillation frequency increases.     

Symmetry of Mountain Pass Solutions of Some Vector Field Equations

We prove radial symmetry (or axial symmetry) of the mountain pass solution of variational elliptic systems − AΔu(x) + ∇ F(u(x)) = 0 (or − ∇.(A(r) ∇ u(x)) + ∇ F(r,u(x)) = 0,) u(x) = (u ₁(x),...,u _N (x)), where A (or A(r)) is a symmetric positive definite matrix. The solutions are defined in a domain Ω which can be , a ball, an annulus or the exterior of a ball. The boundary conditions are either Dirichlet or Neumann (or any one which is invariant under rotation). The mountain pass solutions studied here are given by constrained minimization on the Nehari manifold. We prove symmetry using the reflection method introduced in Lopes [(1996), J. Diff. Eq. 124, 378–388; (1996), Eletron. J. Diff. Eq. 3, 1–14].     

Fluid transfer between tubes in interacting capillary bundle models

The interacting capillary bundle model proposed by Dong et al. [Dong, M., Dullien, F.A.L., Zhou, J.: Trans. Porous Media 31, 213–237 (1998); Dong, M., Dullien, F.A.L., Dai, L., Li, D.: Trans. Porous Media 59, 1–18 (2005); Dong, M., Dullien, F.A.L., Dai, L., Li, D.: Trans. Porous Media 63, 289–304 (2006)] has simulated correctly various aspects of immiscible displacement in porous media, such as oil production histories at different viscosity ratios, the effects of water injection rate and of the oil–water viscosity ratio on the shape of the displacement front and the independence of relative permeabilities of the viscosity ratio. In the interacting capillary bundle model pressure equilibrium was assumed at any distance x measured along the bundle. Interaction between the capillaries also results in transfer of fluids across the capillaries. In the first part of this paper the process of fluid transfer between two capillaries is analysed and an algebraic expression for this flow is derived. Consistency with the assumption of pressure equilibration requires that all transfer must take place at the positions of the oil/water menisci in the tubes without any pressure drop. It is shown that fluid transfer between the tubes has no effect on the predictions obtained with the model. In the second part of the paper the interacting tube bundle model is made more realistic by assuming fluid transfer between the tubes all along the single phase flow regions across a uniform resistance, resulting in pressure differences throughout the single phase regions between the fluids present in the different tubes. The results of numerical simulations obtained with this improved interacting capillary bundle model show only small differences in the positions of the displacement front as compared with the predictions of the idealized model.     

Comment on the Shiner–Davison–Landsberg Measure

The complexity measure from Shiner et al. [Physical Review E 59, 1999, 1459–1464] (henceforth abbreviated as SDL-measure) has recently been the subject of a fierce debate. We discuss the properties and shortcomings of this measure, from the point of view of our recently constructed fundamental, statistical mechanics-based measures of complexity C_s(γ,β) [Stoop et al., J. Stat. Phys. 114, 2004, 1127–1137]. We show explicitly, what the shortcomings of the SDL-measure are: It is over-universal, and the implemented temperature dependence is trivial. We also show how the original SDL-approach can be modified to rule out these points of critique. Results of this modification are shown for the logistic parabola.     

Group classification of one-dimensional nonisentropic equations of fluids with internal inertia

A systematic application of the group analysis method for modeling fluids with internal inertia is presented. The equations studied include models such as the nonlinear one-velocity model of a bubbly fluid (with incompressible liquid phase) at small volume concentration of gas bubbles (Iordanski Zhurnal Prikladnoj Mekhaniki i Tekhnitheskoj Fiziki 3, 102–111, 1960; Kogarko Dokl. AS USSR 137, 1331–1333, 1961; Wijngaarden J. Fluid Mech. 33, 465–474, 1968), and the dispersive shallow water model (Green and Naghdi J. Fluid Mech. 78, 237–246, 1976; Salmon 1988). These models are obtained for special types of the potential function W(r,[(r)\dot],S){W(\rho,\dot \rho,S)} (Gavrilyuk and Teshukov Continuum Mech. Thermodyn. 13, 365–382, 2001). The main feature of the present paper is the study of the potential functions with W _S  ≠ 0. The group classification separates these models into 73 different classes.     

A thermo-mechanically-coupled large-deformation theory for amorphous polymers in a temperature range which spans their glass transition

Amorphous thermoplastic polymers are important engineering materials; however, their non-linear, strongly temperature- and rate-dependent elastic-viscoplastic behavior is still not very well understood, and is modeled by existing constitutive theories with varying degrees of success. There is no generally agreed upon theory to model the large-deformation, thermo-mechanically-coupled, elastic-viscoplastic response of these materials in a temperature range which spans their glass transition temperature. Such a theory is crucial for the development of a numerical capability for the simulation and design of important polymer processing operations, and also for predicting the relationship between processing methods and the subsequent mechanical properties of polymeric products. In this paper we extend our recently published theory [Anand, L., Ames, N. M., Srivastava, V., Chester, S. A., 2009. A thermo-mechanically-coupled theory for large deformations of amorphous polymers. Part I: formulation. International Journal Plasticity 25, 1474–1494; Ames, N. M., Srivastava, V., Chester, S. A., Anand, L., 2009. A thermo-mechanically coupled theory for large deformations of amorphous polymers. Part II: applications. International Journal of Plasticity 25, 1495–1539] to fill this need.     

On the Implementation of the Eigenvalue Method for Limit Cycle Determination in Nonlinear Systems

In many practical systems, limit cycles can be predicted with suitable precision by frequency domain methods using describing functions. Within such an approach, limit cycles can be predicted using the “eigenvalue method” [Somieski, G., Nonlinear Dynamics 26(1), 2001, 3–22]. This contribution presents a novel and advantageous implementation of this method, using singular value instead of eigenvalue calculations, and enhancing computational efficiency by avoiding a so called “frequency iteration”. An erratum to this article is available at .     

Effect of yaw on the starting point of curvature for reflected diffracted shock wave (subsonic and sonic cases)

Lighthill (Proc. R. Soc. A 198, 454–470, 1949) considered the diffraction of a normal shock wave passing over a small bend. The bend being small Lighthill was able to linearize the flow equations and solved the problem through several mathematical techniques. Following Lighthill (Proc. R. Soc. A 198, 454–470, 1949), Srivastava and Chopra (J. Fluid Mech. 40, 821–831, 1970) extended the work to the diffraction of oblique shock waves. Srivastava (AIAAJ 33, 2230–2231, 1995) considered the problem of starting point of curvature and extended the work to yawed wedges (Srivastava in Proceedings of the 14th International Mach reflection symposium Sun Marina Hotel, Yonezawa, Japan, 1–5 October 2000, pp. 225–249, 2002). Srivastava (Shock waves 13, 323–326, 2003) considered the problem for starting point of curvature when the relative outflow behind reflected shock before diffraction has been subsonic and sonic. The present work is an extension of the work published in Srivastava (Shock waves 13, 323–326, 2003) when the wedge has been yawed through an angle. The results have been obtained for two angles χ = 60° and χ = 40° (χ is the angle of yaw).      

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

Three-dimensional Turbulent Boundary Layer in a Shrouded Rotating System

A thre-dimensional direct numerical simulation is combined with a laboratory study to describe the turbulent flow in an enclosed annular rotor-stator cavity characterized by a large aspect ratio G = (b − a)/h = 18.32 and a small radius ratio a/b = 0.152, where a and b are the inner and outer radii of the rotating disk and h is the interdisk spacing. The rotation rate Ω considered is equivalent to the rotational Reynolds number Re = Ωb ²/ν= 9 .5 × 10⁴ (ν the kinematic viscosity of water). This corresponds to a value at which experiment has revealed that the stator boundary layer is turbulent, whereas the rotor boundary layer is still laminar. Comparisons of the computed solution with velocity measurements have given good agreement for the mean and turbulent fields. The results enhance evidence of weak turbulence by comparing the turbulence properties with available data in the literature (Lygren and Andersson, J Fluid Mech 426:297–326, 2001). An approximately self-similar boundary layer behavior is observed along the stator. The wall-normal variations of the structural parameter and of characteristic angles confirm that this boundary layer is three-dimensional. A quadrant analysis (Kang et al., Phys Fluids 10:2315–2322, 1998) of conditionally averaged velocities shows that the asymmetries obtained are dominated by Reynolds stress-producing events in the stator boundary layer. Moreover, Case 1 vortices (with a positive wall induced velocity) are found to be the major source of generation of special strong events, in agreement with the conclusions of Lygren and Andersson (J Fluid Mech 426:297–326, 2001).     

Asymptotic Stability of Riemann Solutions for a Class of Multidimensional Systems of Conservation Laws with Viscosity

We prove the asymptotic stability of two-state nonplanar Riemann solutions for a class of multidimensional hyperbolic systems of conservation laws when the initial data are perturbed and viscosity is added. The class considered here is those systems whose flux functions in different directions share a common complete system of Riemann invariants, the level surfaces of which are hyperplanes. In particular, we obtain the uniqueness of the self-similar L^∞ entropy solution of the two-state nonplanar Riemann problem. The asymptotic stability to which the main result refers is in the sense of the convergence as t→∞ in L_loc¹ of the space of directions ξ = x/t. That is, the solution u(t, x) of the perturbed problem satisfies u(t, tξ)→R(ξ) as t→∞, in L_loc¹(ℝⁿ), where R(ξ) is the self-similar entropy solution of the corresponding two-state nonplanar Riemann problem.     

Special Fast Diffusion with Slow Asymptotics: Entropy Method and Flow on a Riemannian Manifold

We consider the asymptotic behaviour of positive solutions u(t, x) of the fast diffusion equation ${u_t=\Delta (u^{m}/m)= {\rm div}\,(u^{m-1} \nabla u)}We consider the asymptotic behaviour of positive solutions u(t, x) of the fast diffusion equation u_t=D(u^m/m) = div (u^m-1 ?u){u_t=\Delta (u^{m}/m)= {\rm div}\,(u^{m-1} \nabla u)} posed for x ? \mathbb R^d{x\in\mathbb R^d}, t > 0, with a precise value for the exponent m = (d − 4)/(d − 2). The space dimension is d ≧ 3 so that m < 1, and even m = −1 for d = 3. This case had been left open in the general study (Blanchet et al. in Arch Rat Mech Anal 191:347–385, 2009) since it requires quite different functional analytic methods, due in particular to the absence of a spectral gap for the operator generating the linearized evolution. The linearization of this flow is interpreted here as the heat flow of the Laplace– Beltrami operator of a suitable Riemannian Manifold (\mathbb R^d,g){(\mathbb R^d,{\bf g})}, with a metric g which is conformal to the standard \mathbb R^d{\mathbb R^d} metric. Studying the pointwise heat kernel behaviour allows to prove suitable Gagliardo–Nirenberg inequalities associated with the generator. Such inequalities in turn allow one to study the nonlinear evolution as well, and to determine its asymptotics, which is identical to the one satisfied by the linearization. In terms of the rescaled representation, which is a nonlinear Fokker–Planck equation, the convergence rate turns out to be polynomial in time. This result is in contrast with the known exponential decay of such representation for all other values of m.     

Dynamics of a Non-Autonomous ODE System Occurring in Coagulation Theory

We consider a constant coefficient coagulation equation with Becker–D?ring type interactions and power law input of monomers J ₁(t) = α t ^ω, with α > 0 and . For this infinite dimensional system we prove solutions converge to similarity profiles as t and j converge to infinity in a similarity way, namely with either or constants, where is a function of t only. This work generalizes to the non-autonomous case a recent result of da Costa et al. (2004). Markov Processes Relat. Fields 12, 367–398. and provides a rigorous derivation of formal results obtained by Wattis J. Phys. A: Math. Gen. 37, 7823–7841. The main part of the approach is the analysis of a bidimensional non-autonomous system obtained through an appropriate change of variables; this is achieved by the use of differential inequalities and qualitative theory methods. The results about rate of convergence of solutions of the bidimensional system thus obtained are fed into an integral formula representation for the solutions of the infinite dimensional system which is then estimated by an adaptation of methods used by da Costa et al. (2004). Markov Processes Relat. Fields 12, 367–398.      

Deduction of Saint-Venant"s Principle for the Asymptotic Residue of Elastic Rods

Let u(ε) be a rescaled 3-dimensional displacement field solution of the linear elastic model for a free prismatic rod Ω^ε having cross section with diameter of order ε, and let u ⁽⁰⁾ –Bernoulli–Navier displacement – and u ⁽²⁾ be the two first terms derived from the asymptotic method. We analyze the residue r(ε) = u(ε) − (u ⁽⁰⁾ + ε² u ⁽²⁾) and if the cross section is star-shaped, we prove such residue presents a Saint-Venant"s phenomenon near the ends of the rod. This revised version was published online in August 2006 with corrections to the Cover Date.     

Eigenvalues of Self-Similar Solutions of the Dafermos Regularization of a System of Conservation Laws via Geometric Singular Perturbation Theory

The Dafermos regularization of a system of n conservation laws in one space dimension admits smooth self-similar solutions of the form u=u(X/T). In particular, there are such solutions near a Riemann solution consisting of n possibly large Lax shocks. In Lin and Schecter (2004, SIAM. J. Math. Anal. 35, 884–921), eigenvalues and eigenfunctions of the linearized Dafermos operator at such a solution were studied using asymptotic expansions. Here we show that the asymptotic expansions correspond to true eigenvalue–eigenfunction pairs. The proofs use geometric singular perturbation theory, in particular an extension of the Exchange Lemma.     

Convergence of Solutions to the Boltzmann Equation in the Incompressible Euler Limit

 We consider here the problem of deriving rigorously, for well-prepared initial data and without any additional assumption, dissipative or smooth solutions of the incompressible Euler equations from renormalized solutions of the Boltzmann equation. This completes the partial results obtained by Golse [B. Perthame and L. Desvillettes eds., Series in Applied Mathematics 4 (2000), Gauthier-Villars, Paris] and Lions & Masmoudi [Arch. Rational Mech. Anal. 158 (2001), 195–211]. (Accepted June 6, 2002) Published online December 3, 2002 Communicated by Y. BRENIER