# American Institute of Mathematical Sciences

September  2018, 23(7): 2935-2950. doi: 10.3934/dcdsb.2018112

## Smooth to discontinuous systems: A geometric and numerical method for slow-fast dynamics

 1 School of Mathematics, Georgia Tech, Atlanta, GA 30332, USA 2 Dipartimento di Matematica, Univ. of Bari, I-70100, Bari, Italy

Received  August 2017 Revised  November 2017 Published  March 2018

Fund Project: This work was carried out while the second author was visiting the School of Mathematics of the Georgia Institute of Technology, whose hospitality is gratefully acknowledged

We consider a smooth planar system having slow-fast motion, where the slow motion takes place near a curve γ. We explore the idea of replacing the original smooth system with a system with discontinuous right-hand side (DRHS system for short), whereby the DRHS system coincides with the smooth one away from a neighborhood of γ. After this reformulation, in the region of phase-space where γ is attracting for the DRHS system, we will obtain sliding motion on γ and numerical methods apt at integrating for sliding motion can be applied. Moreover, we further bypass resolving the sliding motion and monitor entries (transversal) and exits (tangential) on the curve γ, a fact that can be done independently of resolving for the motion itself. The end result is a method free from the need to adopt stiff integrators or to worry about resolving sliding motion for the DRHS system. We illustrate the performance of our method on a few problems, highlighting the feasibility of using simple explicit Runge-Kutta schemes, and that we obtain much the same orbits of the original smooth system.

Citation: Luca Dieci, Cinzia Elia. Smooth to discontinuous systems: A geometric and numerical method for slow-fast dynamics. Discrete & Continuous Dynamical Systems - B, 2018, 23 (7) : 2935-2950. doi: 10.3934/dcdsb.2018112
##### References:
 [1] L. Dieci and L. Lopez, A survey of Numerical Methods for IVPs of ODEs with Discontinuous right-hand side, Journal of Computational and Applied Mathematics, 236 (2012), 3967-3991. doi: 10.1016/j.cam.2012.02.011. [2] A. F. Filippov, Differential Equations with Discontinuous Right-Hand Sides, Mathematics and Its Applications, Kluwer Academic, Dordrecht, 1988. [3] J. K. Hale, Ordinary Differential Equations, Krieger Publishing Co, Malabar, 1980. [4] M. Levi, Qualitative analysis of the periodically forced relaxation oscillations, Memoirs AMS, 32 (1981), ⅵ+147 pp. [5] N. Levinson, A second order differential equation with singular solutions, Annals of Mathematics, 50 (1949), 127-153. doi: 10.2307/1969357. [6] S. Natsiavas, Dynamics of piecewise linear oscillators with van der Pol type damping, International Journal of Non Linear Mechanics, 26 (1991), 349-366. doi: 10.1016/0020-7462(91)90065-2. [7] A. Roberts and P. Glendinning, Canard-like phenomena in piecewise-smooth Van der Pol systems, Chaos: An Interdisciplinary Journal of Nonlinear Science, Chaos, 24 (2014), 023138, 11pp. [8] J. Sotomayor and M. A. Teixeira, Regularization of discontinuous vector fields, International Conference on Differential Equations, Lisboa, (1998), 207-223. [9] D. W. Storti and R. H. Rand, A simplified model of coupled relaxation oscillators, Int.l J. Nonlin. Mechanics, 22 (1987), 283-289. doi: 10.1016/0020-7462(87)90020-5. [10] W. Wasow, Asymptotic Expansions for Ordinary Differential Equations, Dover Publications, 1987.

show all references

##### References:
 [1] L. Dieci and L. Lopez, A survey of Numerical Methods for IVPs of ODEs with Discontinuous right-hand side, Journal of Computational and Applied Mathematics, 236 (2012), 3967-3991. doi: 10.1016/j.cam.2012.02.011. [2] A. F. Filippov, Differential Equations with Discontinuous Right-Hand Sides, Mathematics and Its Applications, Kluwer Academic, Dordrecht, 1988. [3] J. K. Hale, Ordinary Differential Equations, Krieger Publishing Co, Malabar, 1980. [4] M. Levi, Qualitative analysis of the periodically forced relaxation oscillations, Memoirs AMS, 32 (1981), ⅵ+147 pp. [5] N. Levinson, A second order differential equation with singular solutions, Annals of Mathematics, 50 (1949), 127-153. doi: 10.2307/1969357. [6] S. Natsiavas, Dynamics of piecewise linear oscillators with van der Pol type damping, International Journal of Non Linear Mechanics, 26 (1991), 349-366. doi: 10.1016/0020-7462(91)90065-2. [7] A. Roberts and P. Glendinning, Canard-like phenomena in piecewise-smooth Van der Pol systems, Chaos: An Interdisciplinary Journal of Nonlinear Science, Chaos, 24 (2014), 023138, 11pp. [8] J. Sotomayor and M. A. Teixeira, Regularization of discontinuous vector fields, International Conference on Differential Equations, Lisboa, (1998), 207-223. [9] D. W. Storti and R. H. Rand, A simplified model of coupled relaxation oscillators, Int.l J. Nonlin. Mechanics, 22 (1987), 283-289. doi: 10.1016/0020-7462(87)90020-5. [10] W. Wasow, Asymptotic Expansions for Ordinary Differential Equations, Dover Publications, 1987.
Replacing the cubic with piecewise linear
Cubic vs. piecewise linear limit cycles
Tubular (left) and Shift (right) neighborhoods $\mathcal{C}(\rho)$, $\rho = 0.1$
Vector fields inside $\mathcal{C}(\rho)$ for both tubular (left) and shift (right) neighborhoods
Example 3.1. Van der Pol limit cycle: $\beta = 10$
Example 3.1. Modified vector fields for the DRHS reformulation: tubular (left) and shift (right) neighborhoods, $\rho = 0.01$
Example 3.1. Method 2.1: DRHS reformulation and sliding motion, tubular neighborhood
Example 3.1. Method 2.2: DRHS reformulation and no sliding motion. Tubular neighborhood for $C(\rho)$. Stepsize $\tau = 10^{-2}$, and 57 steps are needed for the periodic trajectory
Example 3.1. Method 2.2: DRHS reformulation and no sliding motion. Shift neighborhood for $C(\rho)$. Stepsize $\tau = 10^{-2}$, and 93 steps are needed for the periodic trajectory
Example 3.1. Enlargement of entry on $\gamma$ within $\mathcal{C}(\rho)$. Tubular neighborhood, $\rho = 0.01$, $\beta = 10$. Stepsize $\tau = 10^{-4}$
"Corrugated" Van der Pol limit cycle: $\beta = 20$
Modified vector fields, $\rho = 0.01$ (left) and $\rho = 0.001$ (right)
Exit points, case of $\rho = 0.01$
Method 2.2 on (3.1)
Limit cycle of (3.2), and slow manifold
Modified vector fields for (3.2)
Method 2.2 for (3.2)
 [1] Angela Alberico, Teresa Alberico, Carlo Sbordone. Planar quasilinear elliptic equations with right-hand side in $L(\log L)^{\delta}$. Discrete & Continuous Dynamical Systems - A, 2011, 31 (4) : 1053-1067. doi: 10.3934/dcds.2011.31.1053 [2] Yong Xu, Bin Pei, Rong Guo. Stochastic averaging for slow-fast dynamical systems with fractional Brownian motion. Discrete & Continuous Dynamical Systems - B, 2015, 20 (7) : 2257-2267. doi: 10.3934/dcdsb.2015.20.2257 [3] Sihong Shao, Huazhong Tang. Higher-order accurate Runge-Kutta discontinuous Galerkin methods for a nonlinear Dirac model. Discrete & Continuous Dynamical Systems - B, 2006, 6 (3) : 623-640. doi: 10.3934/dcdsb.2006.6.623 [4] Ilya Schurov. Duck farming on the two-torus: Multiple canard cycles in generic slow-fast systems. Conference Publications, 2011, 2011 (Special) : 1289-1298. doi: 10.3934/proc.2011.2011.1289 [5] Anatoly Neishtadt, Carles Simó, Dmitry Treschev, Alexei Vasiliev. Periodic orbits and stability islands in chaotic seas created by separatrix crossings in slow-fast systems. Discrete & Continuous Dynamical Systems - B, 2008, 10 (2&3, September) : 621-650. doi: 10.3934/dcdsb.2008.10.621 [6] Jie Xu, Yu Miao, Jicheng Liu. Strong averaging principle for slow-fast SPDEs with Poisson random measures. Discrete & Continuous Dynamical Systems - B, 2015, 20 (7) : 2233-2256. doi: 10.3934/dcdsb.2015.20.2233 [7] Alexandre Vidal. Periodic orbits of tritrophic slow-fast system and double homoclinic bifurcations. Conference Publications, 2007, 2007 (Special) : 1021-1030. doi: 10.3934/proc.2007.2007.1021 [8] Renato Huzak. Cyclicity of the origin in slow-fast codimension 3 saddle and elliptic bifurcations. Discrete & Continuous Dynamical Systems - A, 2016, 36 (1) : 171-215. doi: 10.3934/dcds.2016.36.171 [9] Renato Huzak, P. De Maesschalck, Freddy Dumortier. Primary birth of canard cycles in slow-fast codimension 3 elliptic bifurcations. Communications on Pure & Applied Analysis, 2014, 13 (6) : 2641-2673. doi: 10.3934/cpaa.2014.13.2641 [10] D. J. W. Simpson, R. Kuske. Stochastically perturbed sliding motion in piecewise-smooth systems. Discrete & Continuous Dynamical Systems - B, 2014, 19 (9) : 2889-2913. doi: 10.3934/dcdsb.2014.19.2889 [11] Shin-Ichiro Ei, Hirofumi Izuhara, Masayasu Mimura. Infinite dimensional relaxation oscillation in aggregation-growth systems. Discrete & Continuous Dynamical Systems - B, 2012, 17 (6) : 1859-1887. doi: 10.3934/dcdsb.2012.17.1859 [12] Alexandre Caboussat, Allison Leonard. Numerical solution and fast-slow decomposition of a population of weakly coupled systems. Conference Publications, 2009, 2009 (Special) : 123-132. doi: 10.3934/proc.2009.2009.123 [13] Antonia Katzouraki, Tania Stathaki. Intelligent traffic control on internet-like topologies - integration of graph principles to the classic Runge--Kutta method. Conference Publications, 2009, 2009 (Special) : 404-415. doi: 10.3934/proc.2009.2009.404 [14] Da Xu. Numerical solutions of viscoelastic bending wave equations with two term time kernels by Runge-Kutta convolution quadrature. Discrete & Continuous Dynamical Systems - B, 2017, 22 (6) : 2389-2416. doi: 10.3934/dcdsb.2017122 [15] Wenjuan Zhai, Bingzhen Chen. A fourth order implicit symmetric and symplectic exponentially fitted Runge-Kutta-Nyström method for solving oscillatory problems. Numerical Algebra, Control & Optimization, 2019, 9 (1) : 71-84. doi: 10.3934/naco.2019006 [16] Andrea Giorgini. On the Swift-Hohenberg equation with slow and fast dynamics: well-posedness and long-time behavior. Communications on Pure & Applied Analysis, 2016, 15 (1) : 219-241. doi: 10.3934/cpaa.2016.15.219 [17] Fabrice Bethuel, Didier Smets. Slow motion for equal depth multiple-well gradient systems: The degenerate case. Discrete & Continuous Dynamical Systems - A, 2013, 33 (1) : 67-87. doi: 10.3934/dcds.2013.33.67 [18] María del Mar González, Regis Monneau. Slow motion of particle systems as a limit of a reaction-diffusion equation with half-Laplacian in dimension one. Discrete & Continuous Dynamical Systems - A, 2012, 32 (4) : 1255-1286. doi: 10.3934/dcds.2012.32.1255 [19] Hooton Edward, Balanov Zalman, Krawcewicz Wieslaw, Rachinskii Dmitrii. Sliding Hopf bifurcation in interval systems. Discrete & Continuous Dynamical Systems - A, 2017, 37 (7) : 3545-3566. doi: 10.3934/dcds.2017152 [20] C. Connell Mccluskey. Lyapunov functions for tuberculosis models with fast and slow progression. Mathematical Biosciences & Engineering, 2006, 3 (4) : 603-614. doi: 10.3934/mbe.2006.3.603

2017 Impact Factor: 0.972