arxiv: v1 [math.ds] 20 Sep 2018

Transcription

1 ON NON SMOOTH SLOW FAST SYSTEMS JAIME R. DE MORAES 1 AND PAULO R. DA SILVA 2 arxiv: v1 [math.ds] 20 Sep 2018 Abstract. We deal with non smooth differential systems ż = X(z), z R n, with discontinuity occurring in a codimension one smooth surface Σ. A regularization of X is a 1 parameter family of smooth vector fields X δ, δ > 0, satisfying that X δ converges pointwise to X in R n \ Σ, when δ 0. We work with two known regularizations: the classical one proposed by Sotomayor and Teixeira and its generalization, using non-monotonic transition functions. Using the techniques of geometric singular perturbation theory we study minimal sets of regularized systems. Moreover, non-smooth slow fast systems are studied and the persistence of the sliding region by singular perturbations is analyzed. 1. Introduction One finds in real life and in various branches of science distinguished phenomena whose mathematical models are expressed by piecewise smooth systems and deserve a systematic analysis, see for instance [10, 12, 18]. However sometimes the treatment of such objects is far from the usual techniques or methodologies found in the smooth universe. In fact, for such systems, everything we know from the qualitative theory of dynamical systems has its own versions, starting with the concept of solution. Consider two smooth vector fields X +, X defined in R n. A piecewisesmooth system is ẋ = X(x) with (1) X = 1 [ (1 + sgn(h))x + + (1 sgn(h))x ) ], 2 h : R n R smooth and 0 a regular value of h. The set Σ = {x R n : h(x) = 0} is called switching manifold. In order to define what a solution is, it is necessary, first of all, to agree on what happens in Σ. The points in Σ are classified as regular (if X + and X are transversal to Σ) or singular (if X + or X is tangent to Σ). Moreover the regular points are classified according to Filippov s terminology [7] ( 1 ) (i) Σ w = {x Σ : (X + h.x h)(x) > 0} is the sewing region; (ii) Σ s = {x Σ : (X + h.x )h(x) < 0} is the sliding region Mathematics Subject Classification. 34C25, 34C45, 34D15, 70K70. Key words and phrases. periodic solutions, invariant manifolds, singular perturbations, slow and fast motions. 1 As usual, we denote Xf = f.x. 1

2 2 JAIME R. DE MORAES AND PAULO R. DA SILVA To be more precise we subdivide Σ s in attracting sliding Σ s a (if X + h < 0 and X h > 0) and repelling sliding or escape Σ s r (if X + h > 0 and X h < 0). The orbits of X by Σ w are naturally concatenated. On Σ s is defined the sliding vector field X Σ as a linear convex combination of X + and X which is tangent to Σ. The orbits by Σ s follow the flow of X Σ, a linear convex combination of X + and X tangent to Σ, that is (2) X Σ = (X+.h)X (X.h)X + (X + X. ).h The vector field X Σ is called sliding vector field. While Filippov said how the flow of a piecewise smooth vector field behaves when finding the set of discontinuity, Sotomayor and Teixeira (see [24]) addressed the problem by seeking smooth approximations which were called regularization. A regularization of X is a family of smooth vector fields X δ depending on a parameter δ > 0 such that X δ converges simply to X in R n \ Σ when δ goes to zero. The Sotomayor-Teixeira regularization (ST-regularization) is the one parameter family X δ given by ( 1 + ϕ(h/δ) ) ( 1 ϕ(h/δ) ) (3) X δ = X + + X 2 2 where ϕ : R [ 1, 1] is a smooth function satisfying that ϕ(t) = 1 for t 1, ϕ(t) = 1 for t 1 and ϕ (t) > 0 for t ( 1, 1). The regularization is smooth for δ > 0 and satisfies that X δ = X + on {h δ} and X δ = X on {h δ}. The flow of the regularized vector field proposed by Sotomayor Teixeira, after the limit process, is exactly the flow idealized by Filippov. In 2005 Silva, Teixeira and Buzzi, strongly inspired by Freddy Dumortier, wrote the article [3]. They proved that the regularization proposed by Sotomayor and Teixeira generates a singular perturbation problem. This provides a very important application of GSP theory (geometric singular perturbation theory). Joint with Llibre they published [13, 15, 16, 17]. They studied regularization problems in R n, the double regularization in the case in which the discontinuity has codimension 1 (intersection of two planes) and the regularization in more degenerate surfaces. Bonet-Revés, Larrosa, M-Seara, Kristiansen, Uldall and Hogan also studied the singular perturbation problem arising from regularization. More precisely, they analyzed the regularization of fold-fold singularities where bifurcation and canard boundary cycles can occur. In [1, 2] the authors use of asymptotic analysis and the extension of the critical manifold to non normally hyperbolic points. In [11] the authors strongly use the blow up techniques developed by Dumortier Roussarie [5] to deal with the same problem. Our contribution to the general theory is to investigate the following problems.

3 NON SMOOTH SLOW FAST SYSTEMS 3 Applying techniques of the geometric singular perturbation theory we study the limit periodic sets (equilibrium points and periodic orbits contained in the switching manifold) obtained as limit of orbits of X δ when δ 0 and we give alternative proofs of some results of [23]. See Section 2. The effect of breaking the monotonicity condition of the transition function ϕ used in the regularization process (3). We review the concepts of sliding and sewing and their dependence on the regularization process considered. What is the relation between such sets and those idealized by Filippov? See Section 3. Non smooth slow fast systems with sliding points in the critical manifold. We generalize the results of [4] considering the new concept of sliding and sewing. See Sections 4 and 5. In Section 2 we show how the GSP-theory can be used to get information about the regularized vector field. From our previous papers we know that the trajectories of a piecewise smooth vector field are obtained solving a slow fast system with critical manifold being the graphic of a smooth function defined in Σ s. Moreover, the projection of the reduced flow, in Σ, is the sliding flow. We prove that hyperbolic equilibrium points p (respectively periodic orbits γ) of the sliding vector field (2) are limit of sequences of hyperbolic equilibrium points p δ (respectively periodic orbits γ δ ) of the regularized vector field (3). Besides, the dimensions of stable and unstable manifolds of p δ (respectively γ δ ) are determined. See Theorem 1. In Section 3 we propose a more general definition of sewing and sliding points (called r-sewing and r-sliding points, respectively). Sewing and sliding points are defined depending on the choice of the regularization. Roughly speaking, a point p is a sewing point for a regularization r if around p the flow of r is transversal to Σ. A point p is a sliding point for r if there exists a sequence of invariant manifolds of r tending to a neighborhood of p in Σ. We prove that the sewing region Σ w contains the r-sewing region Σ w r and the r-sliding region Σ s r contains the sliding region Σ s. Besides the sliding vector field on Σ s can be smoothly extended on Σ s r. See Figure 1 and Theorem 2. In Section 4 we study non smooth slow fast systems. Let F, G : R n+1 [0, + ) R n, H : R n+1 [0, + ) R and h : R n+1 [0, + ) R be smooth. We consider systems of the kind { F (x, y, ε), if h(x, y, ε) 0 ẋ =, εẏ = H(x, y, ε). G(x, y, ε), if h(x, y, ε) 0 We prove that the r-sliding region in the critical manifold H(x, y, 0) = 0 persists for small ε > 0. See Theorem 3. In Section 5 we give one more definition of sewing and sliding points using any continuous combination (not necessary convex) of X + and X. We state and prove that sliding regions obtained via continuous combination

4 4 JAIME R. DE MORAES AND PAULO R. DA SILVA (A) (B) (C) Figure 1. Figure (A) exhibits the flow of a piecewise smooth vector fields. Figure (B) is obtained using the ST regularization. The sliding region is given in red. Figure (C) is obtained using a r regularization. Note that Σ s Σ s r. of X + and X also are persistent by singular perturbation. See Theorem Fenichel s Theory and Sliding Vector Fields. In order to simplify our explanation we take local local coordinates such that R n = R n 1 R, Σ = {(x, y) R n 1 R : y = 0}. Let X +, X : R n R n be smooth vector fields. Denote X + = (f 1, g 1 ) and X = (f 2, g 2 ). The piecewise smooth vector field which we consider in this section is (4) X = 1 [ (1 + sgn(y))x + + (1 sgn(y))x ) ]. 2 The sliding vector field (2) becomes (5) X Σ = (X+.y)X (X.y)X + (X + X ).y = ( ) f2 g 1 f 1 g 2, 0. g 1 g 2 The trajectories of the ST-regularized vector field X δ given by (3) are the the solutions of the differential system { ẋ = (f1 + f (6) P δ : 2 )/2 + ϕ(y/δ)(f 1 f 2 )/2 ẏ = (g 1 + g 2 )/2 + ϕ(y/δ)(g 1 g 2 )/2. Consider the polar blow up (x, y, δ) = Φ(x, θ, r) = (x, r cos θ, r sin θ). We get the system { ẋ = (f1 + f 2 ) /2 + ϕ (cot θ) (f 1 f 2 ) /2, P r : r θ = sin θ [(g 1 + g 2 ) /2 + ϕ (cot θ) (g 1 g 2 ) /2].

5 NON SMOOTH SLOW FAST SYSTEMS 5 ȳ 1 x θ 1 Figure 2. Correspondence between the polar coordinates (r, θ) and the directional coordinates (x, ȳ). The phase portrait on the semi cylinder r = 0, (θ, x) [0, π] R n 1 is the central projection of the phase portrait of system (8) with δ = 0, ȳ [ 1, 1]. For r = 0, Pr has two limit problems, the reduced (7) and the layer, as we will see below. (7) P0 reduced : { ẋ = (f1 + f 2 ) /2 + ϕ (cot θ) (f 1 f 2 ) /2, 0 = sin θ [(g 1 + g 2 ) /2 + ϕ (cot θ) (g 1 g 2 ) /2]. P 0 layer : { x = 0, θ = sin θ [(g 1 + g 2 ) /2 + ϕ (cot θ) (g 1 g 2 ) /2]. The polar blow up is not the most suitable for calculations. For this reason, we perform the directional blow-up. The directional blow up consists in the following change of coordinates (x, y, δ) = Γ(x, y, δ) = (x, δy, δ). We observe that the direcional blow up and the polar blow up are essentially the same. In fact, if we consider the map G(x, θ, r) = (x, cot θ, r sin θ) then Γ G = Φ. The direcional blow up applied in system (6) gives (8) Pδ : ẋ = α(x, ȳ, δ), δ ȳ = β(x, ȳ, δ) with α = (f 1 + f 2 )/2 + ϕ(ȳ)(f 1 f 2 )/2, β = (g 1 + g 2 )/2 + ϕ(ȳ)(g 1 g 2 )/2 and the f 1, f 2, g 1, g 2 evaluated at (x, δȳ) GSP-theory. Systems as (8) are known in the literature as slow fast systems. General slow fast systems are systems of the kind (9) ẋ = α(x, y, δ), δẏ = β(x, y, δ),

6 6 JAIME R. DE MORAES AND PAULO R. DA SILVA where x R n, y R k, δ 0 and α and β are smooth functions. Taking δ = 0 in (9) we obtain the reduced system ẋ = α(x, y, 0), β(x, y, 0) = 0. The set S = {β(x, y, 0) = 0} is called critical (or slow) manifold. The time scale τ = t/δ transforms system (9) in the fast system (10) ẋ = δα(x, y, δ), ẏ = β(x, y, δ). Taking δ = 0 in (10) we get the layer system. We say that a point (x 0, y 0 ) S is normally hyperbolic if the real parts of the eigenvalues of D y β(x 0, y 0, 0) are nonzero. Let N S be a compact normally hyperbolic set. Consider (x, y) N and suppose that D y β(x, y, 0) has k s eigenvalues with negative real parts and k u eigenvalues with positive real parts. The following result ensures the persistence of normally hyperbolic sets in S as invariant manifolds of system (9), for small values of δ > 0. Proposition 1 (Fenichel, [6]). Let N S be a j dimensional compact normally hyperbolic manifold with a (j + j s ) dimensional local stable manifold W s and a (j + j u ) dimensional local unstable manifold W u. Then there exists a family N δ such that the following statements hold. (a) N 0 = N. (b) N δ is an invariant manifold of (9) with a (j + j s + k s ) dimensional local stable manifold Nδ s and a (j + ju + k u ) dimensional local unstable manifold Nδ u. Proposition 2. Let X be the piecewise smooth vector field (4) and X δ its ST-regularization (3). If Σ = Σ sl then all points on the critical manifold S of system (8) are normally hyperbolic. In particular S is a graphic of a smooth function ȳ = h(x) with (x, 0) Σ sl. Moreover, the projection of the reduced flow, on Σ, is the sliding flow of (5). Proof. For completeness of the text we rewrite the proof originally presented in [14]. The critical manifold S is defined by the equation (g 1 +g 2 )+ϕ(ȳ)(g 1 g 2 ) = 0, with g 1, g 2 evaluated at (x, 0). If (x, 0) Σ s then g 1 g 2 < 0 and thus (x, ȳ) S ϕ(ȳ) = g 1 + g 2 g 1 g 2. Since g 1+g 2 g 1 g 2 1 and g 1+g 2 g 1 g 2 = 1 if and only if g 1 g 2 = 0 it follows that Dȳβ(x, ȳ, 0) = ϕ (ȳ)(g 1 g 2 )(x, 0) 0. In fact, for (x, 0) Σ s we have g 1+g 2 g 1 g 2 < 1, ȳ ( 1, 1) and thus ϕ (ȳ) > 0. It concludes the proof of the assertion about the critical manifold. To see the relation between the reduced flow and the sliding vector vector field observe that taking ϕ(ȳ) = g 1+g 2 g 1 g 2 the system (8) becomes ẋ = f 2g 1 f 1 g 2 g 1 g 2

7 NON SMOOTH SLOW FAST SYSTEMS 7 which is exactly the same equation of the trajectories of X Σ. Theorem 1. Let X be the piecewise smooth vector field (4). If Q is a l dimensional compact invariant manifold of X Σ given by (5) with a l + l s dimensional local stable manifold, a l + l u dimensional local unstable manifold and l s +l u = n 1. If Q Σ s a then there are a neighborhood V of Q in R n and δ 0 > 0 such that for 0 < δ < δ 0, X δ has a invariant manifold Q δ V with (l + l s + 1) dimensional stable manifold and (l + l u ) dimensional unstable manifold. Proof. Assume l = 0. Thus Q = (x 0, 0) is a hyperbolic equilibrium point of X Σ. The trajectories of the ST-regularization X δ are the solutions of system (6). Performing the directional blow up we get system (8). Moreover, Q is an equilibrium of ẋ = f 2g 1 f 1 g 2 g 1 g 2, ȳ = h(x) with S : ȳ = h(x) defined implictly by (g 1 + g 2 )/2 + ϕ(ȳ)(g 1 g 2 )/2 = 0. The hypothesis g 1 (Q) < 0 and g 2 (Q) > 0 implies that S satisfies the attractiveness condition: all points in S are normally hyperbolic; k s = 1 and k u = 0. Since the reduced system and the sliding system have the same equations we conclude that (x 0, y 0 ), with y 0 = h(x 0 ), is an equilibrium of the reduced system with j s = l s and j u = l u. Then, applying Proposition 1 we conclude that there exists Q δ, an equilibrium point of system (6), with (l s + 1)- dimensional stable manifold and (l u )-dimensional unstable manifolds. Example 1. Let X : R 2 R 2 be the piecewise smooth vector field (4) with X + = (0, 1) and X = (x, y + 1). The sliding vector field (5) is ( x ) X Σ = 2, 0. Q = (0, 0) Σ s a is an repelling equilibrium point of X Σ, that is, l s = 0 and l u = 1. Thus, Theorem 1 says that for small δ, X δ has an equilibrium Q δ of the kind saddle. Note that this can be verified directly with simple calculations. In fact, X δ = ( x 2 ( 1 ϕ ( y δ )), y 2 + ϕ ( y δ ) ( y 2 2 and ( y0 ) Q δ = (0, y 0 ), with ϕ = y 0 δ y 0 2. The existence of Q δ is guaranteed by the fact that the graphics of ϕ(y/δ) y and intersect at some value of y ( 1, 1). Moreover the eigenvalues y 2 of the linearized system at (0, y 0 ) are λ 1 = 1 2, λ 2 = 1 ( y y ) 0 ϕ. 1 y 0 2. δ δ 2 Since y 0 ( 1, 1) we have λ 2 < 0. Soon Q δ is a saddle. ))

8 8 JAIME R. DE MORAES AND PAULO R. DA SILVA Example 2. Let X : R 3 R 3 be the piecewise smooth vector field (4) with ( ) X + = 0, 2x 1 + 2x 2 ( x 21 + x22 1), 1 and X = ( ) 2x 2 + 2x 1 ( x 21 + x22 1), 0, 1. The sliding region is Σ = (x 1, x 2 ) plane and the sliding vector field is ( ) X Σ = x 2 + x 1 ( x x2 2 1), x 1 + x 2 ( x 21 + x22 1), 0. The equilibrium point p 0 = (0, 0, 0) and the limit cycle γ 0 : x x2 2 = 1, of X Σ satisfy the hypotheses of Theorem 1 with l = 0 and l = 1 respectively. 3. Non smooth Systems and Regularization Consider piecewise smooth system (1) defined in an open set U R n with X +, X : U R n, h : U R smooth and assume that 0 is a regular value of h. The following regularization will be refered as r regularization: (11) X δ r = (1 + ψ (p, h/δ)) X (1 ψ (p, h/δ)) X, 2 where ψ : Σ R [ 1, 1] is a more general smooth transition function satisfying that ψ(p, t) = 1 for t 1 and ψ(p, t) = 1 for t 1. The definitions of r-sewing and r-sliding points were introduced in [21] and depend of the regularization r considered. More precisely, p Σ is a r-sewing point if there exist an open neighborhood U R n 1 R of p and local coordinates (x, y) defined in U such that: (a) Σ = {y = 0}; (b) for each sufficiently small δ > 0, the vector field v(x, y) = (0, 1) is a generator of (11) in U. p Σ is a r-sliding point if there exist an open neighborhood U R n of p and a family of smooth manifolds S δ U satisfying: (a) S δ is invariant for r; (b) for each compact K U, the sequence S δ K converges to Σ K, as δ goes to zero according to Hausdorff distance. We denote Σ w r and Σ s r the r-sewing and r-sliding regions, respectively and we assume local coordinates (x, y) R n 1 R such that h(x, y) = y, X + = (f 1, g 1 ), X = (f 2, g 2 ) with (f i (x, y), g i (x, y)) R n 1 R, i = 1, 2.

9 NON SMOOTH SLOW FAST SYSTEMS 9 Σ p S δ Figure 3. A r-sliding point p defined by r-regularization. Theorem 2. Consider a non smooth system (1) and a r-regularization (11). Suppose that ψ t (x, t) 0, for all x Σs and t < 1. Then the following statements hold. (a) Σ w r Σ w and Σ s Σ s r. (b) If g 1 g 2 in Σ s r \Σ s, then the sliding vector field X Σ can be smoothly extended from Σ s to Σ s r. Proof. Consider x 0 Σ w r. Since v(x, y) = (0, 1) is a generator of (11), X (x 0, 0) and X + (x 0, 0) point to the same hand side. So (X +.y)(x.y)(x 0, 0) > 0 and we conclude that x 0 Σ w. Therefore Σ w r Σ w. Now, consider x 0 Σ s. Taking the directional blow up y = δy, the r-regularization (11) becomes the slow fast system (12) ẋ = α(x, y, δ), δẏ = β(x, y, δ) where and α(x, y, δ) = 1/2 ( (1 + ψ(x, y))f 1 + (1 ψ(x, y)) ) f 2 β(x, y, δ) = 1/2 ( (1 + ψ(x, y))g 1 + (1 ψ(x, y)) ) g 2, with the functions f i and g i, i = 1, 2, evaluated at (x, δy). The reduced system (13) ẋ = α(x, y, 0), is defined for (x, y) in the critical manifold S = {β(x, y, 0) = 0}. Consider y 0 < 1 such that (x 0, y 0 ) S. Note that (g 1 g 2 )(x 0, 0) 0 because g 1 g 2 (x 0, 0) < 0. So, β y (x 0, y 0, 0) = ψ y (x 0, y 0 ).(g 1 g 2 )(x 0, 0) 0

10 10 JAIME R. DE MORAES AND PAULO R. DA SILVA and thus (x 0, y 0 ) is a normally hyperbolic point. The Fenichel s result ensures the existence of an invariant manifold S δ of (12) such that S δ S 0 as δ 0, according to Hausdorff s distance. Thus x 0 Σ s r and it concludes the proof of item (a). Since (g 1 g 2 )(x, 0) 0 for x Σ s r, the equation β(x, y, 0) = 0 provides (14) ψ (x, y) = g 1 + g 2 g 1 g 2. The dynamics at (x, y) obtained from (13) and (14) is given by ẋ = f 2g 1 f 1 g 2 g 1 g 2. Hence X Σ can be smoothly extended on Σ s r. It concludes the proof of item (b). The following example illustrates Theorem 2. Example 3. Consider the non smooth system { (x2 1, 1), if x (15) ẋ = (ẋ 1, ẋ 2 ) = 1 0, (x 2, 1), if x 1 0. The sliding region is Σ s =]0, 1[ and the sliding vector field (2) is X Σ = (0, 2x 2 + 1). The r-regularization of (15) is the 1 parameter family x 1 = 1/2 ( 2x 2 1 ψ (x 1 /δ, x 2 ) ), x 2 = ψ (x 1 /δ, x 2 ). Suppose that partial derivative ψ t (t, x 2) vanishes only at (t, x 2 ) = (a 0, b 0 ) with 1 < a 0 < 1 and b 0 > 1. Applying the directional blow up x 1 = δx 1 we obtain the slow fast system (16) δẋ 1 = 1/2 ( 2x 2 1 ψ (x 1, x 2 ) ), ẋ 2 = ψ(x 1, x 2 ). The critical manifold S is given implicitly by x 2 = ψ(x 1, x 2 ) and it is a curve connecting (x 1, x 2 ) = ( 1, 0) to (x 1, x 2 ) = (1, 1). All points in S 0 = S \ (a 0, b 0 ) are normally hyperbolic. So Fenichel s result ensures the existence of an invariant manifold S δ of (16) converging to S 0 according to Hausdorff s distance. Thus Σ s r = ]0, b 0 [. The reduced system in S 0 is ẋ 2 = 2x Thus X Σ can be smoothly extended on Σ s r = ]0, b 0 [. Note that the r-sliding region Σ s r obtained contains the sliding region Σ s = ]0, 1[ defined by Filippov. See Figure 4.

11 NON SMOOTH SLOW FAST SYSTEMS 11 (a 0, b 0) b Figure 4. Slow and fast dynamics of Example (3). The black region is the r-sewing region, the red region is the sliding region according to Filippov convention. The sliding region in blue appears considering the r-regularization of system (15). 4. Non smooth Slow Fast Systems A non smooth slow fast system is { F (x, y, ε), if h(x, y, ε) 0 (17) ẋ =, εẏ = H(x, y, ε) G(x, y, ε), if h(x, y, ε) 0 where x R n, y R and ε > 0 is a small parameter. For each ε 0 let Σ ε be the switching manifold of (17), i.e. Σ ε = {h(x, y, ε) = 0}. We assume that Σ 0 and S = {H(x, y, 0) = 0} are transversal and H y (x, y, 0) 0 on S. Taking ε = 0 in (17) we have the reduced system { F (x) if h(x) 0, (18) ẋ = H(x) = 0, G(x) if h(x) 0, where F (x) = F (x, y(x), 0), G(x) = G(x, y(x), 0), h(x) = h(x, y(x), 0) and H(x) = H(x, y(x), 0). System (18) is defined in the critical manifold S = { H(x) = 0}. We denote Σ w r,0, Σs r,0 and Σw r,ε, Σ s r,ε the r-sewing and r-sliding regions of systems (17) and (18) respectively. The r-regularizations of (17) and (18) are and ẋ = 1/2 ( (1 + ψ (h/δ, p)) F + (1 ψ (h/δ, p)) G ), εẏ = H ẋ = 1/2 ( (1 + ψ( h/δ, p)) F + (1 ψ( h/δ, p)) G ), H = 0, respectively, where ψ : R Σ [ 1, 1] is a more general smooth transition function satisfying that ψ(t, p) = 1 for t 1 and ψ(t, 1) = 1 for t 1.

12 12 JAIME R. DE MORAES AND PAULO R. DA SILVA The next theorem states that r-sliding points p 0 Σ s r,0 persist under effect of singular perturbations with the additional assumptions: (19) H y (p 0, 0) 0, h x (p 0, 0) 0, h y 0, (20) F.h(p0, 0) G.h(p ψ 0, 0), t (t, p 0) 0 for 1 < t < 1. The assumptions are necessary for applying the change of coordinates described in the proof and for ensuring the normal hyperbolicity of p 0. In [22] Sieber and Kowalczyk show that stable periodic motion with sliding is not robust under effect of singular perturbations. Fridman ([8, 9]) also studies periodic motion considering the last assumption of (19). In [4] the authors provide examples showing that sliding regions are not persistent with respect singular perturbations if this assumption is not considered. Theorem 3. Consider a non smooth slow fast system (17) and p 0 Σ s r,0 satisfying the assumptions (19) and (20). Then the following statements hold. (a) There exist sufficiently small ε 0 > 0 and a family of r-sliding points {p ε : ε (0, ε 0 )} of system (17) such that p ε p 0 as ε 0, according to Hausdorff distance. (b) If p 0 is an equilibrium point (or periodic orbit) of the sliding vector field associated to reduced system (18) then there exist sufficiently small ε 1 > 0 and a family of equilibrium points (or periodic orbits) {p ε : ε (0, ε 1 )} of the sliding vector field associated to system (17) such that p ε p 0 as ε 0, according to Hausdorff distance. Example 4. Consider the non smooth slow fast system { (x2 1, 1 + ε), if x (21) (ẋ 1, ẋ 2 ) = 1 0, (x 1 + x 2 + ε, x ε), if x 1 0, The corresponding reduced system { (x2 1, 1), if x (22) (ẋ 1, ẋ 2 ) = 1 0, (x 1 + x 2, 1), if x 1 0, εẏ = y. is defined on the plane {y = 0}. Note that (0, 0, 0) and (0, 1, 0) are fold points and Σ s = ]0, 1[. See Figure 5 (A). The r-regularization of (21) is ẋ 1 = 1/2 ( 1 + x 1 + 2x 2 + ε (x 1 + ε + 1)ψ(x 1 /δ, x 2, y) ), ẋ 2 = 1/2 ( x 2 (2 + x 2 2ε)ψ(x 1 /δ, x 2, y) ), εẏ = y. Assume that the partial derivative ψ t (t, x 2, 0) vanishes only at (t, x 2, 0) = (a 0, b 0, 0) with a 0 = 0 and b 0 = 2. Applying the directional blow up

13 (A) y NON SMOOTH SLOW FAST SYSTEMS 13 (B) y (0, 1, 0) (1, 1, 0) (0, 0, 0) x 1 x 1 ( 1, 0, 0) (0, 2, 0) x 2 x 2 Figure 5. In (A) we have the folds points (0, 0, 0) and (0, 1, 0). The sliding region is Σ s = ]0, 1[. In (B) we obtain after the blow up a smooth curve (critical manifold) connecting the fold points. The point (0, 2, 0) is a non normal hyperbolic point. x 1 = δx 1 the previous system becomes (23) δẋ 1 = 1/2 ( 1 + δ x 1 + 2x 2 + ε (δ x 1 + ε + 1)ψ(x 1, x 2, y) ), ẋ 2 = 1/2 ( x 2 (2 + x 2 2ε)ψ(x 1, x 2, y) ), εẏ = y. The reduced system (ε = δ = 0) associated to (23) is 0 = 1 + 2x 2 ψ (x 1, x 2, 0), ẋ 2 = 1/2 ( x 2 (2 + x 2 )ψ(x 1, x 2, 0) ), 0 = y, and it is a r regularization of (22). The slow manifold S = {(x 1, (1 + ψ(x 1, x 2, 0))/2, 0)} is a curve connecting the points ( 1, 0, 0) and (1, 1, 0). All points in S are normally hyperbolic (parameter δ), except the point (0, 2, 0). According to Fenichel s result S 0 = S {(0, 2, 0)} persists for the system δẋ 1 = 1/2 ( 1 + δ x 1 + 2x 2 (δ x 1 + 1)ψ(x 1, x 2, y) ), ẋ 2 = 1 2( x2 (2 + x 2 )ψ (x 1, x 2, y) ), 0 = y, i.e, there exists an invariant manifold S 0,δ of the previous system converging to S 0, as δ 0. So, Σ s r,0 = ] 2, 1[. For each small δ > 0, S 0,δ is normally hyperbolic (parameter ε) for the previous system. Thus, there exists an invariant manifold S ε,δ of (23) such that ε 0 δ 0 S ε,δ S 0,δ S 0, that is, there exists a r-sliding region Σ s r,ε of (21) converging to Σ s r,0. The sliding vector field of (22) is (24) x 1 = 0, ẋ 2 = 1 x 2 x 2 2, y = 0.

14 14 JAIME R. DE MORAES AND PAULO R. DA SILVA It has the equilibrium points p ± 0 (0, = ( 1 ± ) 5)/2, 0. The sliding vector field of (21) is the slow fast system x 1 = 0, ẋ 2 = x εx 2 x 2 + ε 2 2ε + 1, εẏ = y. ε + 1 Note that its reduced system coincides with (24). It has the equilibrium points ( p ± ε = 0, ( 1 + 2ε ± ) ) 5 12ε + 8ε 2 /2, 0, converging respectively to p + 0 and p 0 as ε goes to zero. Proposition 3. Consider the non smooth slow fast system (17) and p 0 Σ s r,0 satisfying the last two assumptions of (19). Then there exist local coordinates around (p, ε) = (p 0, 0) such that h(x, y, ε) = x 1. Proof. Without lost of generality we can suppose that h/ x 1 (p 0, 0) 0. Applying the change of coordinates x 1 = h(x, y, ε), x i = x i and y = y, for i = 2,..., n, we obtain ẋ 1 = h x 1 ẋ h x n ẋ n + h y ẏ + h ε ε = h x 1 ẋ h ẋ n + h x n y H ε + h ε ε. Since h/ y 0 and ε = 0, the previous expression becomes ẋ 1 = h x 1 ẋ h x n ẋ n. Since h/ x 1 0 the determinant of the change of coordinates matrix is nonzero and system (17) becomes { F (x, y, ε), if x1 0 ẋ =, εẏ = H(x, y, ε). G(x, y, ε), if x 1 0 where x = (x 1,..., x n ). Proof of Theorem 3. Suppose that p 0 Σ s r,0. According to Proposition 3 we can assume that h(x, y, ε) = x 1 in (17) around (p, ε) = (p 0, 0). Denote p = (x 2,..., x n, y). The switching manifold becomes Σ ε = {(0, p)} for each ε 0. The r-regularization of (17) is the 2 parameters (ε and δ) family ẋ = 1/2 ( (1 + ψ(x 1 /δ, p))f + (1 ψ(x 1 /δ, p))g ), εẏ = H. Note that it is a slow fast system (parameter ε) and its reduced system is a r-regularization of reduced system (18). Applying the directional blow up

15 NON SMOOTH SLOW FAST SYSTEMS 15 x 1 = δx 1 we obtain the three time scale singular perturbation problem (25) δẋ 1 = α 1 (x 1, p, ε, δ), ẋ i = α i (x 1, p, ε, δ), εẏ = H(x 1, p, ε, δ), for i = 2,..., n, where α i (x 1, p, ε, δ) = 1/2 ( (1 + ψ(x 1, p))f i + (1 ψ(x 1, p))g i ), for i = 1,..., n, with the function F i and G i evaluated at (δx 1, p, ε) and H(x 1, p, ε, δ) = H(δ x 1, p, ε). The reduced system associated to (25) (ε = δ = 0) is the following (26) (27) 0 = α 1 (x 1, p, 0, 0), ẋ i = α i (x 1, p, 0, 0), 0 = H(x 1, p, 0, 0), for i = 2,..., n. Let S be the critical manifold given by equations (26) and (27) and (t 0, p 0 ) S. Since H y (t 0, p 0, 0, 0) = H y (0, p 0, 0) 0, Proposition 1 says that the point (t 0, p 0 ) persists for the system (28) 0 = α 1 (x 1, p, ε, 0), ẋ i = α i (x 1, p, ε, 0), εẏ = H(x 1, p, ε, 0), for i = 2,..., n. Indeed there exists a compact set S 0 S containing (t 0, p 0 ) and a family of invariant manifolds S ε of (28), such that S ε S 0 as ε 0, according to Hausdorff distance. See Figure 6. Since α 1 (x 1, p, 0, 0) = ψ (x 1, p).(f 1 G 1 )(0, p, 0) 0, x 1 x 1 for (x 1, p) S 0, by continuity we have that α 1 x 1 (x 1, p, 0, ε) 0, for (x 1, p) S ε and sufficiently small ε > 0. Thus Proposition 1 ensures the persistence of S ε for system (25). More specifically, there exists a family of invariant manifolds S ε,δ of (25) such that S ε,δ S ε as δ 0. Therefore there exists a family of r-sliding points p ε of (17) satisfying that p ε p 0 as ε 0, according to Hausdorff s distance. It concludes the proof of statement (a). Note that the sliding vector field associated to (17) is the slow fast system ẋ i = F 1G i G 1 F i F 1 G 1 (0, p, ε), εẏ = H(0, p, ε),

16 16 JAIME R. DE MORAES AND PAULO R. DA SILVA y α 1 (x, p, ε, 0) = 0 S ε H(x, p, 0, 0) = 0 S 0 x 1 x 2 Figure 6. Geometric situation in the blow up process described in the proof of Theorem 3. for i = 2,..., n, and its reduced system ẋ i = F 1G i G 1 F i F 1 G 1 (0, p, 0), H(0, p, 0) = 0, has the dynamics of the sliding vector field associated to (18). So the proof of statement (b) follows directly from Proposition Continuous Combinations of Non Smooth Systems Now we consider another way to define sliding points. Instead of considering a convex combination of vectors X + (p) and X (p) we consider a continuous combination. This convention is given in [19] and [20]. Consider a non smooth system (1). A continuous combination of X + and X is a 1 parameter family of smooth vector fields X(λ, p) with (λ, p) [ 1, 1] Σ satisfying that X(1, p) = X + (p) and X( 1, p) = X (p). We denote [X +, X ] c = { X(λ, p), λ [ 1, 1]}. Consider and a regular point p Σ. (i) We say that p is a c-sewing point and denote p Σ w c if X.h(λ, p) 0 for all λ ( 1, 1). (ii) We say that p is a c-sliding point and denote p Σ s c if there exists λ ( 1, 1), such that X.h(λ, p) = 0. We say that X(λ(p), p) is a c-sliding vector field if for each p Σ s c there exists λ(p) ( 1, 1) such that X.h(λ(p), p) = 0.

17 NON SMOOTH SLOW FAST SYSTEMS 17 There may be more than one possible sliding on p, see Figure 7. In [20] the authors prove that Σ w c Σ w and Σ s Σ s c. Here we study the persistence of c-sliding points via slow fast systems considering the continuous combination Ỹ (λ, x, y, ε) = ( X(λ, x, y, ε), H(x, y, ε)/ε), where X(λ, x, y, ε) is a continuous combination of X + (x, y, ε) and X (x, y, ε). We denote Σ w c,0, Σs c,0 and Σw c,ε, Σ s c,ε the c-sewing and c-sliding regions of systems (17) and (18) respectively. The next theorem provides results like the ones given in Theorem 3 however for c-sliding points. We consider the assumption (19) and the following one Q (29) λ (λ, p 0, 0) 0, where Q(λ, p, ε) = Ỹ. h(λ, p, ε) and λ satisfies the equation Q(λ, p 0, 0) = 0 for p 0 Σ s c,0. Theorem 4. Consider a non smooth slow fast system (17) and p 0 Σ s r,0 satisfying the assumptions (19) and (29).Then the following statements hold. (a) There exist sufficiently small ε 0 > 0 and a family of c-sliding points {p ε : ε (0, ε 0 )} of system (17) such that p ε p 0 as ε 0, according to Hausdorff distance. (b) If p 0 is an equilibrium point (or periodic orbit) of the c-sliding vector field associated to reduced system (18) then there exist sufficiently small ε 1 > 0 and a family of equilibrium points (or periodic orbits) {p ε : ε (0, ε 1 )} of the c-sliding vector field associated to system (17) such that p ε p 0 as ε 0, according to Hausdorff distance. Example 5. Consider the non smooth slow fast system { (x2 1 + ε, 1 + ε), if x (30) ẋ = 1 0 (x 2 + ε, 1 + ε), if x 1 0, εẏ = y and the continuous combination Ỹ (λ, x 1, x 2, y, ε) = ( λ 2 (x 2 + ε) (λ + 1)/2, λ 2 λ 1 + ε, y/ε ). The equation λ 2 (x 2 + ε) (λ + 1)/2 = 0 provides λ ε 1 = x 2 + 8ε 4(x 2 + ε), λ ε 2 = x 2 + 8ε. 4(x 2 + ε) Replacing λ ε 1 and λε 2 in Ỹ we obtain two c-sliding vector fields ( ( X λε 8x2 + 8ε ) 2 ) 1 8x2 + 8ε ε = 0, 16(x 2 + ε) 2 + ε 1, y 4(x 2 + ε) ε

18 18 JAIME R. DE MORAES AND PAULO R. DA SILVA and X λε 2 ε = ( ( 8x2 + 8ε ) 2 ) 8x2 + 8ε , 16(x 2 + ε) ε 1, y. 4(x 2 + ε) ε Note that x 2 ( 1 8ε)/8. For λ 1 < 1 and λ 2 < 1 we obtain Σ s c,ε = ] ε, 1 ε[ ]1 ε, + [ R. In ]1 ε, + [ R are defined two c sliding vector fields (X λ 1 ε and X λ 2 ε ) and in ] ε, 1 ε[ R it is only defined X λ 2 ε. The dynamics of X λ 1 ε and X λ 2 ε are governed respectively by slow fast systems (ẋ 1, ẋ 2, ẏ) = X λε 1 ε, (ẋ 1, ẋ 2, ẏ) = X λε 2 ε. The reduced systems are given respectively by (0, (2x 2 1) ( 4x 2 + 8x ) ) (31) (ẋ 1, ẋ 2 ) = and (32) (ẋ 1, ẋ 2 ) = 8x 2 2 (0, (2x 2 1) ( 4x 2 + 8x ) ) 8x 2 2 Now, consider the reduced system of (30) { (x2 1, 1), if x ẋ = 1 0 (x 2, 1), if x 1 0, defined in the plane y = 0 and the continuous combination X(λ, x, y) = ( λ 2 x 2 (λ + 1)/2, λ 2 λ 1). Using the same process to get Σ s c,ε we obtain that, y = 0,, y = 0. λ 1,2 = 1 ± 1 + 8x 2, 4x 2 and Σ s c,0 = ]0, 1[ ]1, + [. In ]1, + [ are defined two c sliding vector fields X λ 1 0 and X λ 2 0 and they have the dynamics described respectively by (31) and (32). In ]0, 1[ it is only defined X λ 2 0. The sliding vector fields Xλ 1 0 and X λ 2 0 have a normally hyperbolic equilibrium point p 0 = 1/2. The c sliding vector fields X λ 1 ε and X λ 2 ε have the same equilibrium points satisfying that p ± ε p 0 as ε 0. p ± ε = 4ε3 + 8ε 2 5ε + 2 ± 5ε 2 4ε 3 4 (ε 2, 2ε + 1) Proof of Theorem 4. Consider the non smooth system (17) and p 0 Σ s c,0. According to Proposition 3 we can assume that h(x, y, ε) = x 1 around (p, ε) = (p 0, 0). Denote p = (x 2,..., x n, y) and X(λ, x 1, p) = (s 1 (λ, x 1, p),..., s n (λ, x 1, p)) a continuous combinations of F and G.

19 NON SMOOTH SLOW FAST SYSTEMS 19 x 2 X (0, 3) X λ2 0 (0, 3) X λ1 0 (0, 3) X + (0, 3) 0 x 1 Figure 7. The curve in the continuous combinations of the vectors (in red) X (0, 3) and X + (0, 3) in Example 5. The vectors in red represent the c sliding vector fields X λ 1 0 and X λ 2 0 at (0, 3). The c-sliding regions of (17) and (18) become Σ s c,ε = {(0, x 2,..., x n, y) : λ ( 1, 1), s 1 (λ, x, y, ε) = 0}, Σ s c,0 = {(0, x 2,..., x n, y) : λ ( 1, 1), s 1 (λ, x, y, 0) = 0}. Consider λ such that s 1 (λ, p 0, 0) = 0. Assumption (29) ensures the existence of a neighborhood V of (p 0, 0) such that λ = λ(q, ε), λ(p 0, 0) = λ and s 1 (λ(q, ε), q, ε) = 0 for all (q, ε) V, in particular for q V Σ s c,ε. Therefore there exists p ε Σ s c,ε such that p ε p 0 as ε 0 and statement (a) is proved. The dynamics of the c-sliding vector field Ỹ is given by (33) ẋ i = s i (λ, x, y, ε), εẏ = H(x, y, ε), i = 2,..., n, with λ satisfying the equation s 1 (λ, x, y, ε) = 0. Note that system (33) is a singular perturbation problem. The dynamic of the c-sliding vector field associated to system (18) is given by (34) ẋ i = s i (λ, x, y(x), 0), i = 2,..., n, with λ satisfying s 1 (λ, x, y(x), 0) = 0. Since λ = λ(q, ε) the reduced system associated to system (33) is given by (35) ẋ i = s i (λ(q, 0), x, y(x), 0), i = 2,..., n,

20 20 JAIME R. DE MORAES AND PAULO R. DA SILVA Therefore system (34) coincides with system (35). Now we apply the Fenichel s result for concluding the proof of item (b). Remark. In [19] and [20] is defined a regularization of (1) called nonlinear regularization as (36) ẋ = X (ϕ(h/δ), p), where ϕ is a monotonic transition function. Next lemma says that Σ s c is the sliding region linked to nonlinear regularization (36). Thus Theorem 4 can be proved such as Theorem 3 but considering the nonlinear regularization of systems (17) and (18) respectively. ẋ = Ỹ (ϕ(h/ε), x, y, ε 2), ẋ = X(ϕ( h/ε), x, y(x), 0), H(x) = 0, Lemma 5. Consider a non smooth system (1) with K(λ, p) = X.h(λ, p) and X a continuous combination of X + and X. Suppose that there exists λ ( 1, 1) such that K (37) λ (λ, p 0 ) 0, for p 0 Σ. Thus p 0 is a c-sliding point if and only if p 0 is a sliding point for nonlinear regularization (36). Proof. Consider a system like (1) and a continuous combination X of X + and X. Take local coordinates (x 1, p) = (x 1,..., x n ) such that h(x) = x 1 and X(λ, x) = (s 1 (λ, x 1, p),..., s n (λ, x 1, p)). The switching manifold and the c-sliding region become Σ = {(0, p)}, Σ s c = {(0, p) : λ ( 1, 1), s 1 (λ, 0, p) = 0}, respectively. Consider p 0 Σ s c. Thus there exists λ satisfying s 1 (λ, 0, p) = 0. The nonlinear regularization (36) becomes (38) ẋ 1 = s 1 (ϕ(x 1 /δ), x 1, p), ẋ i = s n (ϕ(h/δ), x 1, p), for i = 2,..., n. Taking the blow up x 1 = δx 1 the previous system becomes the slow fast system δẋ 1 = s 1 (ϕ(x 1 ), δx 1, p), ẋ i = s n (ϕ(x 1 ), δx 1, p), for i = 2,..., n. Define λ = ϕ(x 1 ). Since ϕ (t) > 0 in ( 1, 1) there exists x 1 such that λ = ϕ(x 1 ) with (x 1, p 0) M, where M = {s 1 (λ, 0, p) = 0} is the slow manifold. We claim that (x 1, p 0) is a normally hyperbolic point. In fact, s 1 (λ, 0, p 0 ) = ϕ (x x 1) s 1 1 λ (λ, 0, p 0 ) 0.

21 NON SMOOTH SLOW FAST SYSTEMS 21 Thus the Fenichel s result ensures the existence of an invariant manifold M λ of (38) converging to a compact manifold M 0 M containing p 0. Therefore p 0 is a sliding point for the nonlinear regularization (36). Conversely, if p 0 is a sliding point for the nonlinear regularization (36) satisfying assumption (37) then p 0 is a c-sliding point. In fact, note that the slow manifold M and Σ s c are defined by same equation and ϕ is increasing in ( 1, 1). So, there exist neighborhoods U M and V Σ s c of (λ, p 0 ) and p 0, respectively, and a diffeomorphism ξ : U V satisfying ξ(λ, p 0 ) = p Acknowledgments Jaime R. de Moraes is partially supported by FUNDECT 219/2016. Paulo R. da Silva is partially supported by CAPES and FAPESP. References [1] Bonet-Revés C. and M-Seara T. (2016). Regularization of sliding global bifurcations derived from the local fold singularity of Filippov systems, Discrete Contin. Dynam. Systems 36-7, [2] Bonet-Revés C., Larrosa, J. and M-Seara T. (2018). Regularization around a generic codimension one fold-fold singularity, J.Diff. Equation 265, [3] C.A. Buzzi, P.R. da Silva and M.A. Teixeira, A Singular approach to discontinuous vector fields on the plane, J.Diff. Equations 231 (2006), [4] P.T. Cardin, P.R. da Silva and M.A. Teixeira, On singularly perturbed Filippov systems, European Journal of Applied Mathematics 24 (2013), [5] Dumortier, F. and Roussarie, R. (1996). Canard cycles and center manifolds, Memoirs Amer. Mat. Soc [6] N. Fenichel, Geometric singular perturbation theory for ordinary differential equations, J. Diff. Equations 31 (1979), [7] A.F. Filippov, Differential equations with discontinuous right hand sides, Mathematics and its Applications (Soviet Series), Kluwer Academic Publishers, Dordrecht, [8] L. M. Fridman, Singularly perturbed analysis of chattering in relay control systems, IEEE Transactions on Automatic Control 47 (2002), [9] L. M. Fridman, Slow periodic motions with internal sliding modes in variable structure systems, Int. J. Control 75 (2002), [10] M. Guardia, T.M. Seara and M.A. Teixeira, Generic bifurcations of low codimension of planar Filippov systems, J. Diff. Equations 250 (2011), [11] Kristiansen, K. Uldall and Hogan, S. J. (2015). Regularizations of two-fold bifurcations in planar piecewise smooth systems using blow up, SIAM J. Appl. Dyn. Syst. 14-4, [12] Y.A. Kuznetsov, S. Rinaldi and A. Gragnani, One parameter bifurcations in planar filippov systems, Int. J. Bifurc. Chaos 13 (2003), [13] J. Llibre, P.R. da Silva and M.A. Teixeira, Regularization of discontinuous vector fields via singular perturbation, J. Dynam. Differential Equation 19 (2007), [14] J. Llibre, P.R. da Silva and M.A. Teixeira, Sliding vector fields via slow fast systems, Bulletin of the Belgian Mathematical Society Simon Stevin 15 (2008), [15] J. Llibre, P.R. da Silva M.A. Teixeira, Study of Singularities in non smooth dynamical systems via singular perturbation, SIAM Journal on Applied Dynamical Systems 8 (2009),

22 22 JAIME R. DE MORAES AND PAULO R. DA SILVA [16] J. Llibre, P.R. da Silva and M.A. Teixeira, Sliding vector fields for non-smooth dynamical systems having intersecting switching manifolds, Nonlinearity 28 (2015), [17] Llibre, J., Silva, P.R. and Teixeira, M.A. (2008). Sliding vector fields via slow fast systems, Bulletin of the Belgian Mathematical Society Simon Stevin 15, [18] J. Medrado and J. Torregrosa, Uniqueness of limit cycles for sewing planar piecewise linear systems, Journal of Mathematical Analysis and Applications 431 (2015), [19] D.D. Novaes and M. Jeffrey, Regularization of hidden dynamics in piecewise smooth flows, J. Diff. Equations 259 (2015), [20] D.D. Novaes, I.S. Meza Sarmiento and P.R. da Silva, Nonlinear Regularization of Discontinuous Vector Fields and Singular Perturbation, [21] D. Panazzolo and P.R. da Silva, Regularization of discontinuous foliations: Blowing up and sliding conditions via Fenichel theory, J. Diff. Equations 263 (2017), [22] J. Sieber and P. Kowalczyk, Small scale instabilities in dynamical systems with sliding, Phisica D. Nonlinear Phenomena 239 (2009), [23] Sotomayor, J. and Machado, A.L. (2002). Structurally stable discontinuous vector fields on the plane, Qual. Theory of Dynamical Systems 3, [24] J. Sotomayor and M.A. Teixeira, Regularization of discontinuous vector fields, International Conference on Differential Equations, Lisboa, Equadiff 95 (1996), Curso de Matemática UEMS, Rodovia Dourados Itaum Km 12, CEP Dourados, Mato Grosso do Sul, Brazil. address: jaime@uems.br 2 Departamento de Matemática Instituto de Biociências Letras e Ciências Exatas, UNESP Univ Estadual Paulista, Rua C. Colombo, 2265, CEP São José do Rio Preto, São Paulo, Brazil address: paulo.r.silva@unesp.br