Some characteristics of the M function methodology to describe the reinjection process in chaotic intermittency

doi:10.15406/oajmtp.2018.01.00029

Open Access Journal of

eISSN: 2641-9335

Mathematical and Theoretical Physics

Research Article Volume 1 Issue 4

Some characteristics of the M function methodology to describe the reinjection process in chaotic intermittency

Ezequiel del R o,² Luis F Guti rrez Marcantoni,¹ Sergio Elaskar¹

Verify Captcha

Regret for the inconvenience: we are taking measures to prevent fraudulent form submissions by extractors and page crawlers. Please type the correct Captcha word to see email ID.

¹Departamento de Aeronautica, Universidad Nacional de Cordoba and CONICET, Argentina
²E.T.S.I. Aeronautica y Espacio, Universidad Politecnica de Madrid, Spain

Correspondence: Sergio Elaskar, Departamento de Aeronáutica, FCEFyN, Instituto de Estudios Avanzados en Ingenierııía y Tecnologıía, IDIT, Universidad Nacional de Córdoba and CONICET, Córdoba, Argentina, Tel +5435 1535 3800

Received: July 29, 2018 | Published: August 27, 2018

Citation: Elaskari S, del Río E, Gutiérrez Marcantoni L. Some characteristics of the M function methodology to describe the reinjection process in chaotic intermittency. Open Acc J Math Theor Phy. 2018;1(4):168-173 DOI: 10.15406/oajmtp.2018.01.00029

Download PDF

Abstract

The M function methodology allows calculating the statistical properties of intermittency in a broad class of one-dimensional maps. It has shown to be very accurate in type I, II, III and V intermittencies; and it also includes the uniform reinjection as a particular case. This paper studies some properties of theM function methodology. We establish the conditions that a reinjection probability density function must verify to obtain $m = 1 / 2,$ $m = 1 / 2,$ and we describe new pathological cases where the reinjection is not uniform, but the characteristic relation is the same that for uniform reinjection.

Keywords: m function methodology, reinjection process, chaotic intermittency, one-dimensional maps

Introduction

Intermittency is a traditional route to chaos, where a system evolves from regular behavior to chaotic behavior and trajectories alternate chaotic bursts and regular phases. The regular or laminar phases correspond to regions of pseudo-equilibrium or pseudo-periodic solutions, while the bursts are regions with chaotic evolution. In the early eighties, intermittency was classified, according to the system Floquet multipliers or the local Poincaré map eigenvalues, into three different types known as I, II and III.^1‒4 Nevertheless, more recent advances have included other types such as V, X, on-off, in-out, ring and eyelet.^5‒10In engineering, biology, physics and chemistry there are several systems in which chaotic intermittency has been observed.^11‒21 Furthermore, intermittency has been found in economics and medicine systems.^22‒24 Consequently, a more accurate intermittency description might help to increase the knowledge about all these phenomena.

Intermittency behavior is defined by both the local map around the unstable or vanished fixed point and the reinjection mechanism.^1,2,4 The reinjection mechanism maps trajectories from the chaotic zone to the laminar one, which is described by the reinjection probability density function (RPD). Accordingly, the accurate evaluation of the RPD function has a strong influence to describe the intermittency phenomenon correctly. Note that, the calculation of the RPD function from data series (experimental or numerical), is not a simple activity owing in no small amount of data needed and the statistical fluctuations involved in the numerical computations or experimental measurements. In consequence, many approaches have been utilized to describe the RPD function, where the most common one was a constant RPD (uniform reinjection).

A more general methodology to achieve the reinjection probability density function has been elaborated in the last decade, which is named the M function methodology. It includes the uniform reinjection as only a particular case. This methodology has shown to be very accurate for a wide class of maps showing type I, II, III and V intermittencies.^4,25‒38 In this paper, we analyze some characteristics of this methodology. We have shown the M function methodology works very accurately for several maps, but there are pathological cases where it can only partially describe the intermittency reinjection process.

Evaluation of the RPD function

Let us consider a general one-dimensional map: $x n + 1 = F (x n),$ $x n + 1 = F (x n),$ which shows chaotic intermittency. To describe the reinjection process, we should evaluate the RPD function, $φ (x),$ $φ (x),$ which determines the probability density that trajectories are reinjected into a point x inside the laminar interval.^1,2,4

A new theoretical scheme to evaluate the intermittency reinjection process is the M function methodology, in which the RPD is indirectly obtained. First, the $M (x)$ $M (x)$ function must be evaluated:^4,25‒39

$M (x) = ⎧ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ \begin{matrix} \frac{\int_{ˆ x}^{x} τ ϕ (τ) d τ}{\int_{ˆ x}^{x} ϕ (τ) d τ}, i f x \int ˆ x ϕ (τ) d τ \neq 0, 0, o t h e r w i s e, \end{matrix}$ $M (x) = {\begin{matrix} \frac{\int_{\hat{x}}^{x} τ ϕ (τ) d τ}{\int_{\hat{x}}^{x} ϕ (τ) d τ}, i f \int_{\hat{x}}^{x} ϕ (τ) d τ \neq 0, \\ 0, o t h e r w i s e, \end{matrix}$ (1)

Where $ˆ x$ $\hat{x}$ is the lower boundary of reinjection. From the data series the $ˆ x$ $\hat{x}$ calculation is straightforward:

$M (x) ≅ \frac{1}{N} \sum_{J = 1}^{N} x_{j}$ $M (x) ≅ \frac{1}{N} \sum_{J = 1}^{N} x_{j}$ (2)

Where the reinjection points ${x, j}_{j = 1}^{N}$ ${x, j}_{j = 1}^{N}$ must be sorted from the lowest to the highest,^25‒33 i.e. $x_{j} \leq x_{j + 1} .$ $x_{j} \leq x_{j + 1} .$

Previous studies have shown a linear $M (x)$ $M (x)$ for a wide class of maps with type I, II, III and V intermittencies:^4,25‒39

$M (x) = {\begin{matrix} m (x - ˆ x) + x, i f ˆ x \leq x \leq c, 0, o t h e r w h i s e, \end{matrix}$ $M (x) = {\begin{matrix} m (x - \hat{x}) + x, i f \hat{x} \leq x \leq c, \\ 0, o t h e r w h i s e, \end{matrix}$ (3)

c is the limit of the laminar interval.

In the previous equation the main parameter is $m \in (0, 1)$ $m \in (0, 1)$ (the $M (x)$ $M (x)$ function slope) which is determined by the nonlinear map. When $M (x)$ $M (x)$ is a linear function, it has been shown the RPD function is a power law:^{4,25-26,28,31,38,39}

$ϕ (x) = b (α) {(x - ˆ x)}^{α}, with α = \frac{2 m - 1}{1 - m}$ $ϕ (x) = b (α) {(x - \hat{x})}^{α}, with α = \frac{2 m - 1}{1 - m}$ (4)

Where $b (α) = \frac{α + 1}{{(c - ˆ x)}^{(α + 1)}}$ $b (α) = \frac{α + 1}{{(c - \hat{x})}^{(α + 1)}}$ is the normalization parameter.

From Eq. (4) two parameters are only needed to describe the RPD function: m and $ˆ x,$ $\hat{x},$ which are calculated from $M (x) .$ $M (x) .$ m is the slope of $M (x),$ $M (x),$ and it verifies $M (ˆ x) = ˆ x,$ $M (\hat{x}) = \hat{x},$ i.e. it allows to calculate $ˆ x .$ $\hat{x} .$ Then, the $M (x)$ $M (x)$ function stores the map nonlinear information.^4,35,36

Finally, note that uniform reinjection $(α = 0)$ $(α = 0)$ is only a particular case of the new theoretical formulation when $m = 1 / 2.$ $m = 1 / 2.$

Some characteristic of the M(x) function

In this section some properties of the $M (x)$ $M (x)$ function are studied. We have special interest to describe the behavior of $ϕ (x)$ $ϕ (x)$ and $M (x)$ $M (x)$ functions, mainly when the $M (x)$ $M (x)$ function slope is $m = 1 / 2.$ $m = 1 / 2.$

In the following Theorem, will expose two conditions that the RPD must verify to obtain $m = 1 / 2$ $m = 1 / 2$ in Eq. (3).

Theorem:

If the reinjection probability density function, $φ (x),$ $φ (x),$ satisfies the following conditions:

$ϕ (ˆ x) \neq 0$ $ϕ (\hat{x}) \neq 0$
${\frac{d ϕ (x)}{d x} ∣ ∣}_{x = ˆ x} i s b o u n d e d$ ${\frac{d ϕ (x)}{d x} |}_{x = \hat{x}} i s b o u n d e d$

Then the $M (x)$ $M (x)$ function can be approximated as $M (x) = x / 2$ $M (x) = x / 2$ for points close enough to $x = ˆ x .$ $x = \hat{x} .$

Proof:

To prove this hypothesis, the slope of the $M (x)$ $M (x)$ function for x close to $ˆ x$ $\hat{x}$ can be derived from Eq. (1):

$lim x \to ˆ x \frac{d M (x)}{d x} = lim x \to ˆ x \frac{x ϕ (x) \int_{ˆ x}^{x} ϕ (τ) d τ - ϕ (x) \int_{ˆ x}^{x} τ ϕ (τ) d τ}{{⎛ ⎝ \int_{ˆ x}^{x} ϕ (τ) d τ ⎞ ⎠}^{2}}$ $\lim_{x \to \hat{x}} \frac{d M (x)}{d x} = \lim_{x \to \hat{x}} \frac{x ϕ (x) \int_{\hat{x}}^{x} ϕ (τ) d τ - ϕ (x) \int_{\hat{x}}^{x} τ ϕ (τ) d τ}{{(\int_{\hat{x}}^{x} ϕ (τ) d τ)}^{2}}$ (5)

Which by the L’Hôpital rule, gives

$\lim_{x \to \hat{x}} \frac{d M (x)}{d x} = \lim_{x \to \hat{x}} \frac{ϕ (x) \int_{\hat{x}}^{x} ϕ (τ) d τ + x ϕ^{'} (x) \int_{\hat{x}}^{x} ϕ (τ) d τ - ϕ^{'} (x) \int_{\hat{x}}^{x} τ ϕ (τ) d τ}{2 ϕ (x) \int_{\hat{x}}^{x} ϕ (τ) d τ}$ (6)

The last expression results:

$\lim_{x \to \hat{x}} \frac{d M (x)}{d x} = \lim_{x \to \hat{x}} [\frac{1}{2} + \frac{1}{2} \frac{x ϕ^{'} (x)}{ϕ (x)} - \frac{1}{2} \frac{ϕ^{'} (x)}{ϕ (x)} \frac{\int_{\hat{x}}^{x} τ ϕ (τ) d τ}{\int_{\hat{x}}^{x} ϕ (τ) d τ}]$ (7)

Applying the L’Hôpital rule only on the last term, we can obtain:

$\lim_{x \to \hat{x}} \frac{d M (x)}{d x} = \lim_{x \to \hat{x}} [\frac{1}{2} + \frac{1}{2} \frac{x ϕ^{'} (x)}{ϕ (x)} - \frac{1}{2} \frac{x ϕ^{'} (x)}{ϕ (x)}]$ (8)

Because $ϕ (x) \neq 0$ and ${\frac{d ϕ (x)}{d x} |}_{x = \hat{x}}$ is bounded, we obtain 1/2 for the limit given in Eq. 6.

This Theorem generalizes previous developments in which $\hat{x} = x_{0} = 0,$ where $x_{0} = 0$ is the vanished or unstable fixed point.^4,26 For the Theorem, the LBR $(\hat{x})$ can be any point inside the laminar interval, being the vanished or unstable fixed point $x_{0}$ only a particular case.

If the conditions (1) and (2) are verified, then for a small laminar interval $(c \to 0),$ the RPD can be uniform, and the classical theory is valid. However, there could be RPDs satisfying the two conditions and not having uniform reinjection; we call these RPDs as pathological cases.

Pathological case for type II intermittency

Type II intermittency starts from a subcritical Hopf bifurcation.^2,4 Then, two complex-conjugate eigenvalues of the system leave the unit circle. When a control parameter $(ε)$ is less than a critical value, the map has a stable fixed point. On the other hand, if the control parameter rises above this threshold, a bifurcation leads to a new chaotic attractor. The old attractor remains as a subset of the new attractor.

In several previous studies the Theorem of Section 3 has been verified for type I, II, III and V intermittencies with and without lower boundary of reinjection.^{4,^25‒38}

In this section, we study the intermittency process for a map showing type II intermittency. The map is an extension of those studied^25,40 and it is given by:

$F (x) = {\begin{matrix} F_{1} (x) = (1 - ε) x - a x^{3}, x < x_{r} \\ F_{2} (x) = w {(\frac{x - x_{r}}{1 - x_{r}})}^{γ_{2}}, x > x_{r} \end{matrix}$ (9)

$x_{r}$ is obtained from $F_{1} (x_{r}) = 1,$ and $0 \leq w \leq 1.$ The map has a fixed point $x_{0} = 0$ which is stable for $- 2 < ε < 0.$ It becomes unstable for $ε > 0$ and type II intermittency can happen. The laminar interval is $[0, c] .$ In Eq. (9), $F_{1} (x)$ function is the typical local map for type II intermittency and the nonlinear function, $F_{2} (x),$ generates the reinjection process. Therefore, the exponents $γ_{1}$ and $γ_{2}$ will drive the reinjection mechanism.

To numerically analyze type II intermittency phenomenon for this map, we consider the following parameters: $γ_{1} = 2,$ $γ_{1} = 0.87,$ $w = 1 / 2,$ $ε = 0.001$ and $N = 100000,$ and we apply the M function methodology to calculate the RPD function. Figure 1 displays the results. Red points represent the numerical data and blue lines the theoretical results calculated using Eqs. (2) and (4).

Figure 1 Type II intermittency.

$M (x)$ and

$ϕ (x)$ functions for

$γ_{1} = 2,$

$γ_{2} = 0.87,$ w=0.5,

$ε = 0.001$ and N=100000. Red points represent the numerical data, and the blue line the theoretical results calculated using Eqs. 2 and 4.

From Figure 1(A) we can observe the theoretical RPD cannot approach accurately the numerical RPD. Beside this, the numerical RPD verifies the conditions $ϕ (\hat{x}) \neq 0$ and bounded ${\frac{d ϕ}{d x} |}_{x = \hat{x}};$ therefore following the previous theorem the $M (x)$ function slope must be $m = 1 / 2.$ Using the date shown in Figure 1(B) we can calculate $m ≅ 0.506 ≅ 0.5$ in agreement with the theory.

Following Section 2, if $m ≅ 0.5,$ the RPD must be constant (uniform reinjection). However, the numerical data in Figure 1(B) (red points) do not correspond to uniform reinjection. Therefore, it is a pathological case for the M function methodology; then the theoretical RPD (calculated using the M function methodology) does not show high accuracy with the numerical data.

To analyze the intermittency behavior for this pathological test, we study the map at points previous to reinjection (pre-reinjection points).³⁹ If the reinjection points are $x_{n + 1} = F (x_{n}),$ we study the map derivative at $x_{n}$ points; because they have influence on the RPD function.35 The derivative of map (9) can be written as:

$\frac{d F (x)}{d x} = {\begin{cases} \frac{d F_{1} (x)}{d x} = (1 + ε) + 3 a x^{2}, x \leq x_{r} \\ \frac{d F_{2} (x)}{d x} = \frac{w}{{(1 - x_{r})}^{γ_{1}}} γ_{1} {(x - x_{r})}^{γ_{1} - 1} + \frac{1 - w}{{(1 - x_{r})}^{γ_{2}}} γ_{2} {(x - x_{r})}^{γ_{2} - 1}, x > x_{r} \end{cases}$ (10)

Figure 3 shows the Eq. (10) (red line) inside the interval $[x_{r} ≅ 0.682, F_{2}^{- 1} (c) ≅ 0.726473]$ which contains all the pre-reinjection points, i.e $x_{n}$ points previous to reinjection. Moreover, this figure shows the derivative of the map used in^25,40 with $γ = 0.4$ (blue line) and $γ = 1$ (green line):

$F_{1} (x) = {\begin{cases} F_{1} (x) = (1 + ε) x + a x^{3}, x \leq x_{r} \\ \frac{d F_{2} (x)}{d x} = {(\frac{x - x_{r}}{1 - x_{r}})}^{γ}, x > x_{r} \end{cases}$ (11)

$x_{r}$ verifies $F_{1} (x_{r}) = 1,$ and $x_{0} = 0$ is a fixed point of the map, which is stable for $- 2 < ε < 0.$

Note some differences between the three curves. The derivative of Eq. (9) (red line) has lower values regarding the derivative of Eq. (11) (blue line) –which is important for points close to $x_{r}$ – and the second derivative of map (9) goes from negative to positive values inside the interval [0.682, 0.726473]. This derivative behavior generates a different reinjection process between maps (9) and (11).

The RPD for $γ = 0.4$ is given in Figure 2(B), note as $x \to \hat{x},$ the derivative $d F (x) / d x \to \infty,$ which implies $ϕ (x) \to 0.$ Also, for $γ = 0.4,$ the derivative $d F (x) / d x$ decreases inside the laminar interval.^41‒47 On the other hand, for Eq. (9), the derivative ${F (x) / d x |}_{x \to \hat{x}}$ does not tend to infinite, but it takes higher values than in other points in the laminar interval. This result might seem contradictory with the analytical limit of Eq. (10) which tends to infinite when $x \to \hat{x} .$ This contradiction can be explained because Figure 3 is obtained by a discrete process where there are a large but finite number of reinjected points, and not by a continuous process with infinite number of reinjected points (like as of Eq. (10)).

Figure 2 Type II intermittency.

$M (x)$ and

$ϕ (x)$ functions for

$γ = 0.4,$

$ε = 0.0001,$ c = 0.1 and N = 20000. Red points represent the numerical data, and the blue line the theoretical results calculated using Eqs. (2) and (4).

As the derivative ${F (x) / d x |}_{x \to \hat{x}}$ reaches higher values than in other points inside the laminar interval for map (9), then the RPD acquires lower and non-zero values close to $\hat{x}$ (Figure 1(B)). As x increases, $d F (x) / d x$ decreases and $ϕ (x)$ increases until $d F (x) / d x = 0$ (inlection point), from this point $d F (x) / d x$ grows very softly and $ϕ (x)$ decreases very softly too. The inflection point is $x_{i n} ≅ 0.707,$ it can be calculated by:

$x_{i n} = x_{r} + (1 - x_{r}) {(\frac{1 - w}{w} \frac{γ_{1} (1 - γ_{2})}{γ_{2} (γ_{1} - 1)})}^{\frac{1}{γ_{1} - γ_{2}}}$ (12)

Finally, we highlight that for uniform reinjection $d F (x) / d x$ is constant (green line in Figure 3).

Figure 3 Type II intermittency. The derivative of maps (9) and (11). The parameters for map 9 are

$γ_{1} = 2,$

$γ_{2} = 0.87,$ w=0.5,

$ε = 0.001$ and N=100000. The parameters for map 11 are

$γ = 0.4$ (blue) and

$γ = 1$ (green),

$ε = 0.001.$ Red points correspond to map 9. Blue and green points correspond to map (11)

$γ = 0.4$ with and

$γ = 1$ respectively. Figure (B) is a zoom of Figure (A).

The characteristic relation is an important expression to determine the intermittency behavior. Using the M function methodology, this relation for type II intermittency is:⁴

$\bar{l} α ε^{- β}$ (13)

Where the critical exponent $β$ is given:²⁵

$β = \frac{3 - α - 2}{3 - 1} = \frac{3 (1 - m) - 1}{(3 - 1) (1 - m)} t h e s a m e t h i n g$ (14)

and $\hat{l}$ is the average laminar length, which can be calculated using:

$\bar{I} = \int_{x_{0} - c}^{x_{0} + c} ϕ (x) l (x, c) d x$ (15)

$l (x, c)$ is the laminar length for each point in the laminar interval, which measures the number of iterations inside the laminar interval.^2,4

To study the characteristic relation, we have carried out several numerical tests. Table 1 shows m and $α$ for different $ε .$ From the table, we can observe that $α$ does not depend on the control parameter, $ε .$

From Table 1, we have found $α ≅ 0$ and $m ≅ 1 / 2;$ then the theoretical evaluation for the characteristic relation exponent, given by Eq. (14), is $β ≅ 1 / 2.$ To verify this exponent, we calculate $\bar{l}$ for the tests included in Table 1. The results are exposed in Figure 4, which shows the numerical characteristic relation for the map (9). Red points are the numerical data, and the blue line is its linear interpolation. The slope of this straight line is approximately -0.495, which is very close to the theoretical exponent $β = 1 / 2$ predicted using the M function methodology.^4,25 Another consequence of the previous Theorem is that the traditional values of the characteristic relation exponent, $β = 1 / 2,$ can be obtained not only for uniform reinjection but also for any RPD holding the conditions (1) and (2). Figure 4 shows this behavior. Therefore, the test studied in this subsection does not have a uniform RPD, but its characteristic relation is equal to those obtained with uniform reinjection.

$ε$	m	$α$
0.05	0.508	0.035
0.01	0.507	0.032
0.005	0.504	0.0189
0.001	0.506	0.024
0.0005	0.504	0.0188
0.0001	0.506	0.024
0.00001	0.505	0.0229
0.000001	0.505	0.0229

Table 1 Values of m, for different Parameters: w=0.5 and c=0.1

Figure 4 Type II intermittency. Characteristic relation for

$γ_{1} = 2,$

$γ_{1} = 0.87,$ and w=0.5. Red points are the numerical average laminar length, and the blue line the linear interpolation.

Analysis and conclusion

In this paper, we have analyzed a theorem describing some properties of the $M (x)$ and $ϕ (x)$ functions. We have been able to identify numerically two reinjection mechanisms satisfying this theorem. One of them corresponds to uniform reinjection process, and the M function methodology captures and describes it very well. Pathological cases like Eq. (9) give the other one. For these cases, the $M (x)$ function slope verifies $m = 1 / 2,$ but there is not uniform reinjection, and the M unction methodology would seem to introduce errors in the RPD evaluation, but it calculates the characteristic relation correctly.

The map derivative for pathological cases at points previous to reinjection, $x_{n} - x_{n + 1} = F (x_{n})$ where $x_{n + 1}$ are the reinjected points–shows a different behavior regarding to the non-pathological cases (Figure 3). For pathological behavior, the map derivative does not tend to infinity close to $\hat{x}$ which is produced by the discrete process and the finite number of reinjected points. Another aspect to highlight is that inside the laminar interval the second derivative goes from negative to positive values.

TheM function methodology has theoretically established that the characteristic relation exponent, $β,$ depends on $α$ using Eqs. (13) and (14). Here, we have found numerically that the exponent, $β,$ depends on $α$ in the same way for both pathological cases and uniform reinjection processes. The characteristic relation for all cases with $α ≅ 0.5$ results:

$\bar{l} α ε^{- β}$ (16)

It is a valuable property because the characteristic relation is used to describe the intermittency process. Therefore, when we obtain from experimental or numerical data a characteristic relation like Eq. (16), we have to study other parameters to determine the intermittency behavior, because it can be a process with uniform reinjection or a pathological case.

In some previous papers, we have verified that the M function methodology works accurately in several maps with different intermittency types. Besides, here we have shown there are pathological cases where this methodology partially describes the intermittency reinjection process. For these cases, theM function methodology captures very well de characteristic relation, but it could introduce errors in the RPD function evaluation. This phenomenon occurs because a discrete reinjection process around $x_{0} = 0$ is studied; this process has a large but finite number of reinjected points. In this case none of the two exponents $γ$ in Eq. (9)( $γ_{1} = 2$ and $γ_{2} = 0.87$ ) dominate the reinjection mechanism. However, if we study a process with an infinite number of reinjected points, the RPD will tend to zero for $x_{0} = 0$ as shown Eq. (10), the Theorem of Section 4 would not be applied, and the M function methodology would work correctly.