# 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  September 2018 Early access  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:

show all references

##### References:
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, 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] Chunhua Shan. Slow-fast dynamics and nonlinear oscillations in transmission of mosquito-borne diseases. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021097 [5] 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 [6] 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 [7] Antonella Zanna. Symplectic P-stable additive Runge—Kutta methods. Journal of Computational Dynamics, 2021  doi: 10.3934/jcd.2021030 [8] 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 [9] 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 [10] Renato Huzak. Cyclicity of the origin in slow-fast codimension 3 saddle and elliptic bifurcations. Discrete & Continuous Dynamical Systems, 2016, 36 (1) : 171-215. doi: 10.3934/dcds.2016.36.171 [11] 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 [12] Chungang Shi, Wei Wang, Dafeng Chen. Weak time discretization for slow-fast stochastic reaction-diffusion equations. Discrete & Continuous Dynamical Systems - B, 2021, 26 (12) : 6285-6310. doi: 10.3934/dcdsb.2021019 [13] 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 [14] 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 [15] Elisa Giesecke, Axel Kröner. Classification with Runge-Kutta networks and feature space augmentation. Journal of Computational Dynamics, 2021, 8 (4) : 495-520. doi: 10.3934/jcd.2021018 [16] 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 [17] Manuel F. Rañada. Quasi-bi-Hamiltonian structures and superintegrability: Study of a Kepler-related family of systems endowed with generalized Runge-Lenz integrals of motion. Journal of Geometric Mechanics, 2021, 13 (2) : 195-208. doi: 10.3934/jgm.2021003 [18] 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 [19] 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 [20] 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

2020 Impact Factor: 1.327