Scaling law in saddle-node bifurcations for one-dimensional maps: a complex variable approach

The study of transient dynamical phenomena near bifurcation thresholds has attracted the interest of many researchers due to the relevance of bifurcations in different physical or biological systems. In the context of saddle-node bifurcations, where two or more fixed points collide annihilating each other, it is known that the dynamics can suffer the so-called delayed transition. This phenomenon emerges when the system spends a lot of time before reaching the remaining stable equilibrium, found after the bifurcation, because of the presence of a saddle-remnant in phase space. Some works have analytically tackled this phenomenon, especially in time-continuous dynamical systems, showing that the time delay, τ, scales according to an inverse square-root power law, τ∼(μ−μ _c )^−1/2, as the bifurcation parameter μ, is driven further away from its critical value, μ _c . In this work, we first characterize analytically this scaling law using complex variable techniques for a family of one-dimensional maps, called the normal form for the saddle-node bifurcation. We then apply our general analytic results to a single-species ecological model with harvesting given by a unimodal map, characterizing the delayed transition and the scaling law arising due to the constant of harvesting. For both analyzed systems, we show that the numerical results are in perfect agreement with the analytical solutions we are providing. The procedure presented in this work can be used to characterize the scaling laws of one-dimensional discrete dynamical systems with saddle-node bifurcations.     

An investigation of electromagnetic wave propagation in plasma by shock tube

This paper presents the electromagnetic wave propagation characteristics in plasma and the attenuation coefficients of the microwave in terms of the parameters he, v, w, L, wb. The φ800 mm high temperature shock tube has been used to produce a uniform plasma. In order to get the attenuation of the electromagnetic wave through the plasma behind a shock wave, the microwave transmission has been used to measure the relative change of the wave power. The working frequency is f = (2-35)GHz (ω=2πf, wave length A =15cm-8mm). The electron density in the plasma is ne = (3×10^10-1×10^14) cm^-3. The collision frequency v = (1×10^8-6×10^10) Hz. The thickness of the plasma layer L = (2-80)cm. The electron circular frequency ωb=eBo/me, magnetic flux density B0 = (0-0.84)T. The experimental results show that when the plasma layer is thick (such as L/λ≥10), the correlation between the attenuation coefficients of the electromagnetic waves and the parameters ne,v,ω, L determined from the measurements are in good agreement with the theoretical predictions of electromagnetic wave propagations in the uniform infinite plasma. When the plasma layer is thin (such as when both L and A are of the same order), the theoretical results are only in a qualitative agreement with the experimental observations in the present parameter range, but the formula of the electromagnetic wave propagation theory in an uniform infinite plasma can not be used for quantitative computations of the correlation between the attenuation coefficients and the parameters ne,v,ω, L. In fact, if ω<ωp, v^2<<ω^2, the power attenuations K of the electromagnetic waves obtained from the measurements in the thin-layer plasma are much smaller than those of the theoretical predictions. On the other hand, if ω>ωp, v^2<<ω^2 (just v≈f), the measurements are much larger than the theoretical results. Also, we have measured the electromagnetic wave power attenuation value under the magnetic field and without a magnetic field. The result indicates that the value measured under the magnetic field shows a distinct improvement.     

Heat transfer enhancement in grooved channels due to flow bifurcations

The flow bifurcation scenario and heat transfer characteristics in grooved channels, are investigated by direct numerical simulations of the mass, momentum and energy equations, using the spectral element methods for increasing Reynolds numbers in the laminar and transitional regimes. The Eulerian flow characteristics show a transition scenario of two Hopf bifurcations when the flow evolves from a laminar to a time-dependent periodic and then to a quasi-periodic flow. The first Hopf bifurcation occurs to a critical Reynolds number Rec₁ that is significantly lower than the critical Reynolds number for a plane-channel flow. The periodic and quasi-periodic flows are characterized by fundamental frequencies ω₁ and m· ω₁+n·ω₂, respectively, with m and n integers. Friction factor and pumping power evaluations demonstrate that these parameters are much higher than the plane channel values. The time-average mean Nusselt number remains mostly constant in the laminar regime and continuously increases in the transitional regime. The rate of increase of this Nusselt number is higher for a quasi-periodic than for a periodic flow regime. This higher rate originates because better flow mixing develops in quasi-periodic flow regimes. The flow bifurcation scenario occurring in grooved channels is similar to the Ruelle-Takens-Newhouse transition scenario of Eulerian chaos, observed in symmetric and asymmetric wavy channels.     

Rheo-Dielectric behavior of low molecular weight liquid crystals.

Dielectric relaxation behavior was examined for 4-4′-n-pentyl-cyanobiphenyl (5CB) and 4-4′-n-heptyl-cyanobiphenyl (7CB) under flow. In quiescent states at all temperatures examined, both 5CB and 7CB exhibited dispersions in their complex dielectric constant ε*(ω) at characteristic frequencies ω _c above 10⁶ rad s^–1. This dispersion reflected orientational fluctuation of individual 5CB and 7CB molecules having large dipoles parallel to their principal axis (in the direction of C≡N bond). In the isotropic state at high temperatures, these molecules exhibited no detectable changes of ε*(ω) under flow at shear rates . In contrast, in the nematic state at lower temperatures the terminal relaxation intensity of ε*(ω) as well as the static dielectric constant ε′(0) decreased under flow at . This rheo-dielectric change was discussed in relation to the flow effects on the nematic texture (director distribution) and anisotropy in motion of individual molecules with respect to the director. Received: 14 April 1998 Accepted: 29 July 1998     

Revisiting the Drainage Relative Permeability Measurement by Centrifuge Method Using a Forward–backward Modeling Scheme

Measurement of drainage relative permeability by the centrifuge method was first introduced by Hagoort (SPE J. 29(3):139–150, 1980). It has been shown that capillary end effects can cause error in the measurement of relative permeability if a minimum rotational speed is not honoured. To determine the minimum rotational speed that makes the capillary end effect negligible, ω _min, we propose that the value of capillary-gravity number, N _cg, should be of the order of 10⁻² or smaller. This conclusion is based on the use a Forward–backward scheme consisting of a forward numerical simulator developed for centrifuge experiments and applying Hagoort’s method as a backward model. The article presents the use of this Forward–backward scheme as a powerful tool for error analysis such as determining the impact of capillary end effects. By using this loop, we first determine ω _min for specific core and fluid properties. Later, we generalize the ω _min calculations by using the definition of N _cg as a “rule of thumb” for designing relative permeability experiments by centrifuge method. We also demonstrate another use of this loop for controlling the quality of the experimental data.     

Rheo-dielectric behavior of low molecular weight liquid crystal 2. Behavior of 8CB in nematic and smectic states

Rheo-dielectric behavior was examined for 4−4^′−n-octyl-cyanobiphenyl (8CB) having large dipoles parallel to its principal axis (in the direction of the C≡N bond). In the quiescent state at all temperatures (T) examined, orientational fluctuation of the 8CB molecules was observed as dielectric dispersions at characteristic frequencies ω_c>10⁶ s⁻¹. In the isotropic state at high T, no detectable changes of the complex dielectric constant ɛ^*(ω) were found under slow flow at shear rates ˙γ≫ω_c. In the nematic state at intermediate T, the terminal relaxation intensity of ɛ^*(ω) was decreased under such slow flow. In the smectic state at lower T, the flow effect became much less significant. These results were related to the flow-induced changes of the liquid crystalline textures in the nematic and smectic states, and the differences of the rheo-dielectric behavior in these states are discussed in relation to a difference of the symmetry of molecular arrangements in the nematic and smectic textures. Received: 1 October 1998 Accepted: 13 January 1999     

One approximate solution of the Nekrasov problem

An approximate solution ω = A[ω, μ] of the nonlinear integral Nekrasov equation is obtained by successive replacement of the kernel of the integral operator by a close one. The solution is sought not directly at the bifurcation point μ₁ = 3 of the linearized equation ω = μL[ω] but at the point μ = 1 at which operator A[ω, μ], remaining nonlinear in ω, is linear in μ. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 6, pp. 50–56, November–December, 2007.     

On the completeness of oscillation spaces

The oscillation spaces introduced by Jaffard are a variation on the definition of Besov spaces for either s ≥ 0 or s ≤ −d/p. On the contrary, the spaces for −d/p < s < 0 cannot be sharply imbedded between Besov spaces with almost the same exponents, and, thus, they are new spaces of really different nature. Their norms take into account correlations between the positions of large wavelet coefficients through the scales. Several numerical studies uncovered such correlations in several settings including turbulence, image processing, traffic, finance, etc. These spaces allow one to capture oscillatory behaviors that are left undetected by Sobolev or Besov spaces. Unlike Sobolev spaces (respectively, Besov spaces B _p ^s,q (ℝ^d)), which are expressed by simple conditions on wavelet coefficients as ℓ^p norms (respectively, mixed ℓ^p − ℓ^q norms), oscillation spaces are written as ℓ^p averages of local C ^s ′ norms. In this paper, we prove the completeness of oscillation spaces in spite of such a mixture of two norms of different kinds. Published in Neliniini Kolyvannya, Vol. 8, No. 4, pp. 435–443, October–December, 2005.     

New Explicit Expression of Barnett-Lothe Tensors for Anisotropic Linear Elastic Materials

The three Barnett-Lothe tensors H, L, S appear often in the Stroth formalism of two-dimensional deformations of anisotropic elastic materials [1–3]. They also appear in certain three-dimensional problems [4, 5]. The algebraic representation of H, L, S requires computation of the eigenvalues p^v(v=1,2,3) and the normalized eigenvectors (a, b). The integral representation of H, L, S circumvents the need for computing p ^v(v=1,2,3) and (a, b), but it is not simple to integrate the integrals except for special materials. Ting and Lee [6] have recently obtained an explicit expression of H for general anisotropic materials. We present here the remaining tensors L, S using the algebraic representation. They key to our success is the obtaining of the normalization factor for (a, b) in a simple form. The derivation of L and S then makes use of (a, b) but the final result does not require computation of (a, b), which makes the result attractive to numerical computation. Even though the tensor H given in [6] is in terms of the elastic stiffnesses C^μ v while the tensors L, S presented here are in terms of the reduced elastic compliances s′ _μv , the structure of L, S is similar to that of H. Following the derivation of H, we also present alternate expressions of L, S that remain valid for the degenerate cases p ₁ p ₂ and p₁=p₂ = p ₃. One may want to compute H, L, S using either C _μv or s′ _μv v, but not both. We show how an expression in C^μ v can be converted to an expression in s′ _μv v, and vice versa. As an application of the conversion, we present explicit expressions of the extic equation for p in C^μ v and s′ _μv v. This revised version was published online in August 2006 with corrections to the Cover Date.     

Wave Packets, Signaling and Resonances in a Homogeneous Waveguide

We solve the initial-boundary-value linear stability problem for small localised disturbances in a homogeneous elastic waveguide formally by applying a combined Laplace – Fourier transform. An asymptotic evaluation of the solution, expressed as an inverse Laplace – Fourier integral, is carried out by means of the mathematical formalism of absolute and convective instabilities. Wave packets, triggered by perturbations localised in space and finite in time, as well as responses to sources localised in space, with the time dependence satisfying e^−iωt + O(e^−ɛt ), for t → ∞, where Im ω₀ = 0 and ω > 0 , that is, the signaling problem, are treated. For this purpose, we analyse the dispersion relation of the problem analytically, and by solving numerically the eigenvalue stability problem. It is shown that due to double roots in a wavenumber k of the dispersion relation function D(k, ω), for real frequencies ω, that satisfy a collision criterion, wave packets with an algebraic temporal decay and signaling with an algebraic temporal growth, that is, temporal resonances, are present in a neutrally stable homogeneous waveguide. Moreover, for any admissible combination of the physical parameters, a homogeneous waveguide possesses a countable set of temporally resonant frequencies. Consequences of these results for modelling in seismology are discussed. This revised version was published online in August 2006 with corrections to the Cover Date.     

A unified approach to closed-form solutions of moving heat-source problems

A closed-form model for the computation of temperature distribution in an infinitely extended isotropic body with a time-dependent moving-heat sources is discussed. The temperature solutions are presented for the sources of the forms: (i) ₀₁(t)=₀ exp(−λt), (ii) ₀₂(t) =₀(t/t ^*)exp(−λt), and ₀₃(t)=₀[1+a cos(ωt)], where λ and ω are real parameters and t ^* characterizes the limiting time. The reduced (or dimensionless) temperature solutions are presented in terms of the generalized representation of an incomplete gamma function Γ(α,x;b) and its decomposition C _Γ and S _Γ. The solutions are presented for moving, -point, -line, and -plane heat sources. It is also demonstrated that the present analysis covers the classical temperature solutions of a constant strength source under quasi-steady state situations. Received on 13 June 1997     

The Evolution of Stress Intensity Factors and the Propagation of Cracks in Elastic Media

When a crack Γ _s propagates in an elastic medium the stress intensity factors evolve with the tip x(s) of Γ _s . In this paper we derive formulae which describe the evolution of these stress intensity factors for a homogeneous isotropic elastic medium under plane strain conditions. Denoting by ψ=ψ(x,s) the stress potential (ψ is biharmonic and has zero traction along the crack Γ _s ) and by κ(s) the curvature of the crack at the tip x(s), we prove that the stress intensity factors A ₁(s), A ₂(s), as functions of s, satisfy:

Fully developed flow of a viscoelastic film down a vertical cylindrical or planar wall

The one-dimensional, gravity-driven film flow of a linear (l) or exponential (e) Phan-Thien and Tanner (PTT) liquid, flowing either on the outer or on the inner surface of a vertical cylinder or over a planar wall, is analyzed. Numerical solution of the governing equations is generally possible. Analytical solutions are derived only for: (1) l-PTT model in cylindrical and planar geometries in the absence of solvent, b o [(h)\tilde]_s/([(h)\tilde]_s +[(h)\tilde]_p)=0\beta\equiv {\tilde{\eta}_s}/\left({\tilde{\eta}_s +\tilde{\eta}_p}\right)=0, where [(h)\tilde]_p\widetilde{\eta}_p and [(h)\tilde]_s\widetilde{\eta}_s are the zero-shear polymer and solvent viscosities, respectively, and the affinity parameter set at ξ = 0; (2) l-PTT or e-PTT model in a planar geometry when β = 0 and x 1 0\xi \ne 0; (3) e-PTT model in planar geometry when β = 0 and ξ = 0. The effect of fluid properties, cylinder radius, [(R)\tilde]\tilde{R}, and flow rate on the velocity profile, the stress components, and the film thickness, [(H)\tilde]\tilde{H}, is determined. On the other hand, the relevant dimensionless numbers, which are the Deborah, De=[(l)\tilde][(U)\tilde]/[(H)\tilde]De={\tilde{\lambda}\tilde{U}}/{\tilde{H}}, and Stokes, St=[(r)\tilde][(g)\tilde][(H)\tilde]²/([(h)\tilde]_p +[(h)\tilde]_s )[(U)\tilde]St=\tilde{\rho}\tilde{g}\tilde{\rm H}^{2}/\left({\tilde{\eta}_p +\tilde{\eta}_s} \right)\tilde{U}, numbers, depend on [(H)\tilde]\tilde{H} and the average film velocity, [(U)\tilde]\widetilde{U}. This makes necessary a trial and error procedure to obtain [(H)\tilde]\tilde{H} a posteriori. We find that increasing De, ξ, or the extensibility parameter ε increases shear thinning resulting in a smaller St. The Stokes number decreases as [(R)\tilde]/[(H)\tilde]{\tilde{R}}/{\tilde{H}} decreases down to zero for a film on the outer cylindrical surface, while it asymptotes to very large values when [(R)\tilde]/[(H)\tilde]{\tilde{R}}/{\tilde{H}} decreases down to unity for a film on the inner surface. When x 1 0\xi \ne 0, an upper limit in De exists above which a solution cannot be computed. This critical value increases with ε and decreases with ξ.     

On the Navier�CStokes equations with rotating effect and prescribed outflow velocity

We consider the equations of Navier–Stokes modeling viscous fluid flow past a moving or rotating obstacle in \mathbb R^d{\mathbb R^d} subject to a prescribed velocity condition at infinity. In contrast to previously known results, where the prescribed velocity vector is assumed to be parallel to the axis of rotation, in this paper we are interested in a general outflow velocity. In order to use L ^p -techniques we introduce a new coordinate system, in which we obtain a non-autonomous partial differential equation with an unbounded drift term. We prove that the linearized problem in \mathbb R^d{\mathbb R^d} is solved by an evolution system on L^p_s(\mathbb R^d){L^p_{\sigma}(\mathbb R^d)} for 1 < p < ∞. For this we use results about time-dependent Ornstein–Uhlenbeck operators. Finally, we prove, for p ≥ d and initial data u₀ ? L^p_s(\mathbb R^d){u_0\in L^p_{\sigma}(\mathbb R^d)}, the existence of a unique mild solution to the full Navier–Stokes system.     

Existence of the minimal positive solution of some nonlinear elliptic systems when the nonlinearity is the sum of a sublinear and a superlinear term

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｔｈｉｓｐａｐｅｒ,ｗｅｃｏｎｓｉｄｅｒｔｈｅｅｌｌｉｐｔｉｃｓｙｓｔｅｍ(1λ)　-Δｕ=ｆ(λ,ｘ,ｕ)-ｖ　　(ｉｎΩ),-Δｖ=δｕ-γｖ(ｉｎΩ),ｕ=ｖ=0(ｏｎΩ),ｗｈｅｒｅΩｉｓａｓｍｏｏｔｈｂｏｕｎｄｅｄｄｏｍａｉｎｉｎＲＮ(Ｎ≥2)ａｎｄλｉｓａｒｅａｌｐａｒａｍｅｔｅｒ.Ｔｈｅｓｏｌｕｔｉｏｎｓ(ｕ,ｖ)ｏｆｔｈｉｓｓｙｓｔｅｍｒｅｐｒｅｓｅｎｔｓｔｅａｄｙｓｔａｔｅｓｏｌｕｔｉｏｎｓｏｆｒｅａｃｔｉｏｎｄｉｆｆｕｓｉｏｎｓｙｓｔｅｍｓｄｅｒｉｖｅｄｆｒｏｍｓｅｖｅｒａｌａｐ…     

Generalized Hopf Bifurcation for Planar Vector Fields via the Inverse Integrating Factor

In this paper we study the maximum number of limit cycles that can bifurcate from a focus singular point p ₀ of an analytic, autonomous differential system in the real plane under an analytic perturbation. We consider p ₀ being a focus singular point of the following three types: non-degenerate, degenerate without characteristic directions and nilpotent. In a neighborhood of p ₀ the differential system can always be brought, by means of a change to (generalized) polar coordinates (r, θ), to an equation over a cylinder in which the singular point p ₀ corresponds to a limit cycle γ ₀. This equation over the cylinder always has an inverse integrating factor which is smooth and non-flat in r in a neighborhood of γ ₀. We define the notion of vanishing multiplicity of the inverse integrating factor over γ ₀. This vanishing multiplicity determines the maximum number of limit cycles that bifurcate from the singular point p ₀ in the non-degenerate case and a lower bound for the cyclicity otherwise. Moreover, we prove the existence of an inverse integrating factor in a neighborhood of many types of singular points, namely for the three types of focus considered in the previous paragraph and for any isolated singular point with at least one non-zero eigenvalue.     

Three-mode interactions in harmonically excited systems with quadratic nonlinearities

An investigation is presented of the response of a three-degree-of-freedom system with quadratic nonlinearities and the autoparametric resonances ₃2₂ and ₂2₁ to a harmonic excitation of the third mode, where the _m are the linear natural frequencies of the system. The method of multiple scales is used to determine six first-order nonlinear ordinary differential equations that govern the time variation of the amplitudes and phases of the interacting modes. The fixed points of these equations are obtained and their stability is determined. For certain parameter values, the fixed points are found to lose stability due to Hopf bifurcations and consequently the system exhibits amplitude-and phase-modulated motions. Regions where the amplitudes and phases display periodic, quasiperiodic, and chaotic time variations and hence regions where the overall system motion is periodically, quasiperiodically, and chaotically modulated are determined. Using various numerical simulations, we investigated nonperiodic solutions of the modulation equations using the amplitudeF of the excitation as a control parameter. As the excitation amplitudeF is increased, the fixed points of the modulation equations exhibit an instability due to a Hopf bifurcation, leading to limit-cycle solutions of the modulation equations. AsF is increased further, the limit cycle undergoes a period-doubling bifurcation followed by a secondary Hopf bifurcation, resulting in either a two-period quasiperiodic or a phase-locked solution. AsF is increased further, there is a torus breakdown and the solution of the modulation equations becomes chaotic, resulting in a chaotically modulated motion of the system.     

Analysis of regular and chaotic dynamics of the Euler-Bernoulli beams using finite difference and finite element methods

Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions(symmetric and non-symmetric)are studied in this work.Reliability of the obtained results is verified by the finite difference method(FDM)and the finite element method(FEM)with the Bubnov-Galerkin approximation for various boundary conditions and various dynamic regimes(regular and non-regular).The influence of boundary conditions on the Euler-Bernoulli beams dynamics is studied mainly,dynamic behavior vs.control parameters { ωp,q0 } is reported,and scenarios of the system transition into chaos are illustrated.     

An analogue rotated vector field of polynomial system

1　Ｉｎｔｒｏｄｕｃｔｉｏｎ·ＤｅｆｉｎｉｔｉｏｎｏｆＡｎａｌｏｇｕｅＲｏｔａｔｅｄＶｅｃｔｏｒＳｕｐｐｏｓｉｎｇｔｈａｔａｌｉｍｉｔｃｙｃｌｅｉｓｌｏｃａｔｅｄｉｎａｒｏｔａｔｅｄｖｅｃｔｏｒｆｉｅｌｄｏｆｐｏｌｙｎｏｍｉａｌｓｙｓｔｅｍｔｈａｔｄｅｐｅｎｄｓｏｎａｐａｒａｍｅｔｅｒα,ａｎｄｗｈｅｎαｍｏｎｏｔｏｎｏｕｓｌｙｃｈａｎｇｅｓ,ｔｈｉｓｌｉｍｉｔｃｙｃｌｅｗｉｌｌｍｏｎｏｔｏｎｏｕｓｌｙｅｘｐａｎｄ(ｏｒｒｅｄｕｃｅ)ｗｉｔｈｔｈｅα.Ｂｕｔｍｏｒｅｔｈａｎｏｎｅｎｅｉｇｈｂｏｕｒｉｎ…