# American Institute of Mathematical Sciences

March  2020, 12(1): 53-84. doi: 10.3934/jgm.2020004

## Rolling and no-slip bouncing in cylinders

 1 Department of Mathematics and Statistics, Mount Holyoke College, 50 College St, South Hadley, MA 01075, USA 2 Department of Mathematics, Tarleton State University, Box T-0470, Stephenville, TX 76401, USA 3 Department of Mathematics and Statistics, Washington University, Campus Box 1146, St. Louis, MO 63130, USA

* Corresponding author: Renato Feres

Received  February 2019 Revised  June 2019 Published  January 2020

We compare a classical non-holonomic system—a sphere rolling against the inner surface of a vertical cylinder under gravity—with certain discrete dynamical systems called no-slip billiards in similar configurations. A feature of the former is that its height function is bounded and oscillates harmonically up and down. We investigate whether similar bounded behavior is observed in the no-slip billiard counterpart. For circular cylinders in dimension $3$, no-slip billiards indeed have bounded orbits, and very closely approximate rolling motion, for a class of initial conditions we call transversal rolling impact. When this condition does not hold, trajectories undergo vertical oscillations superimposed to overall downward acceleration. Concerning different cross-sections, we show that no-slip billiards between two parallel hyperplanes in arbitrary dimensions are always bounded even under a constant force parallel to the plates; for general cylinders, when the orbit of the transverse system (a concept relying on a factorization of the motion into transversal and longitudinal components) has period two, the motion, under no forces, is generically not bounded. This commonly occurs in planar no-slip billiards.

Citation: Timothy Chumley, Scott Cook, Christopher Cox, Renato Feres. Rolling and no-slip bouncing in cylinders. Journal of Geometric Mechanics, 2020, 12 (1) : 53-84. doi: 10.3934/jgm.2020004
##### References:

show all references

##### References:
No-slip billiard system between two parallel plates under gravity have bounded orbits in all dimensions. Far left: a simple periodic orbit with gravity turned off. Orbits shown from left to right are under the influence of increasing force
Conditions for a transversal period $2$ orbit in dimension $3$. The projection to $\mathbb{R}^2$ of the velocity $u$ of the center of mass and the angular velocity $\dot{\theta}$ are related by $|\dot{\theta}| = (m r/\mathcal{I})|u\sin\phi|$, where $m$ is the projected disc's mass, $r$ is its radius, and $\mathcal{I}$ is its moment of inertia for the projected (or marginal) mass distribution
Cross-section of the stadium cylinder and the rolling sphere
Height of center of mass of rolling particle in a cylinder with stadium cross-section
When the transversal rolling impact condition does not hold (see Definition 2.8), the particle acquires an overall acceleration downward. The apparent increase in thickness of the height function graph and of the particle's path is due to a small scale vertical zig-zag motion of increasing amplitude. See also Figure 8
Near grazing paths of no-slip billiard particle. On the left-hand side, the transversal rolling impact condition holds for the initial bounce, but on the right-hand side a small deviation of this condition is introduced. Notice the characteristic zig-zag nature of the curve and that it does not return all the way to the initial height. The starting center of mass position is indicated by the small black dot
Comparison of the height functions for the rolling motion in a circular cylinder under a constant downward force (solid line) and the corresponding no-slip billiard motion satisfying the transverse rolling impact condition (small circles). Initial conditions are chosen so that the two processes rotate around the cylinder at the same rate
These graphs correspond to a fixed initial transversal rolling defect $-r\omega\cdot e/u\cdot\tau_a = 1.15.$ (When the transversal rolling impact initial condition holds, this value is $1$.) Here $\omega\cdot e$ is the longitudinal component of the initial angular velocity vector and $u\cdot \tau_a$ is the tangential (to the boundary of the billiard domain) component of the center of mass velocity. As the intercollision flight becomes shorter and motion grazes the cylinder more and more closely, the height function becomes smooth but the falling rate remains essentially unchanged
Notation for the proof of Corollary 1
Cross-sectional projection of initial velocity $\bar{u}_0$ and definition of $\theta$ and $\delta$
Transition from regular to chaotic motion. The moving particle begins from the middle of the lower flat side with linear velocity pointing up and a small angular velocity that causes it to move right after the first collision with the upper flat side. For small values of the angular velocity trajectories never touch the curved sides of the boundary, and the motion along the axis of the cylinder is bounded. If the initial angular velocity is large enough, trajectories move beyond the ends of the flat sides and eventually becomes unstable
Transition from regular to chaotic motion for the transverse dynamics of the stadium cylinder no-slip billiard system, as viewed in the velocity phase portrait. The full velocity space is a disc of radius $1$ as shown in the far right. Initial conditions for the depicted orbits roughly compare to those of Figure 11
Height function for the stadium cylinder for a trajectory in the chaotic regime, corresponding to the short orbit segment on the right of Figure 11 and of Figure 12. To give a sense of the scales involved, the diameter of the stadium is $8$ and velocities are of order $1$. The number of time steps is $4\times 10^4$
 [1] Alberto Bressan, Ke Han, Franco Rampazzo. On the control of non holonomic systems by active constraints. Discrete & Continuous Dynamical Systems - A, 2013, 33 (8) : 3329-3353. doi: 10.3934/dcds.2013.33.3329 [2] Hildeberto E. Cabral, Zhihong Xia. Subharmonic solutions in the restricted three-body problem. Discrete & Continuous Dynamical Systems - A, 1995, 1 (4) : 463-474. doi: 10.3934/dcds.1995.1.463 [3] Vakhtang Putkaradze, Stuart Rogers. Numerical simulations of a rolling ball robot actuated by internal point masses. Numerical Algebra, Control & Optimization, 2021, 11 (2) : 143-207. doi: 10.3934/naco.2020021 [4] Jan Prüss, Laurent Pujo-Menjouet, G.F. Webb, Rico Zacher. Analysis of a model for the dynamics of prions. Discrete & Continuous Dynamical Systems - B, 2006, 6 (1) : 225-235. doi: 10.3934/dcdsb.2006.6.225 [5] Simone Calogero, Juan Calvo, Óscar Sánchez, Juan Soler. Dispersive behavior in galactic dynamics. Discrete & Continuous Dynamical Systems - B, 2010, 14 (1) : 1-16. doi: 10.3934/dcdsb.2010.14.1 [6] Rongchang Liu, Jiangyuan Li, Duokui Yan. New periodic orbits in the planar equal-mass three-body problem. Discrete & Continuous Dynamical Systems - A, 2018, 38 (4) : 2187-2206. doi: 10.3934/dcds.2018090 [7] Juan Manuel Pastor, Javier García-Algarra, Javier Galeano, José María Iriondo, José J. Ramasco. A simple and bounded model of population dynamics for mutualistic networks. Networks & Heterogeneous Media, 2015, 10 (1) : 53-70. doi: 10.3934/nhm.2015.10.53 [8] Paula A. González-Parra, Sunmi Lee, Leticia Velázquez, Carlos Castillo-Chavez. A note on the use of optimal control on a discrete time model of influenza dynamics. Mathematical Biosciences & Engineering, 2011, 8 (1) : 183-197. doi: 10.3934/mbe.2011.8.183 [9] Guirong Jiang, Qishao Lu. The dynamics of a Prey-Predator model with impulsive state feedback control. Discrete & Continuous Dynamical Systems - B, 2006, 6 (6) : 1301-1320. doi: 10.3934/dcdsb.2006.6.1301 [10] Mingxin Wang, Qianying Zhang. Dynamics for the diffusive Leslie-Gower model with double free boundaries. Discrete & Continuous Dynamical Systems - A, 2018, 38 (5) : 2591-2607. doi: 10.3934/dcds.2018109 [11] Xinyuan Liao, Caidi Zhao, Shengfan Zhou. Compact uniform attractors for dissipative non-autonomous lattice dynamical systems. Communications on Pure & Applied Analysis, 2007, 6 (4) : 1087-1111. doi: 10.3934/cpaa.2007.6.1087 [12] Pascal Noble, Sebastien Travadel. Non-persistence of roll-waves under viscous perturbations. Discrete & Continuous Dynamical Systems - B, 2001, 1 (1) : 61-70. doi: 10.3934/dcdsb.2001.1.61 [13] Vieri Benci, Marco Cococcioni. The algorithmic numbers in non-archimedean numerical computing environments. Discrete & Continuous Dynamical Systems - S, 2021, 14 (5) : 1673-1692. doi: 10.3934/dcdss.2020449 [14] Liqin Qian, Xiwang Cao. Character sums over a non-chain ring and their applications. Advances in Mathematics of Communications, 2021  doi: 10.3934/amc.2020134 [15] Zhiming Guo, Zhi-Chun Yang, Xingfu Zou. Existence and uniqueness of positive solution to a non-local differential equation with homogeneous Dirichlet boundary condition---A non-monotone case. Communications on Pure & Applied Analysis, 2012, 11 (5) : 1825-1838. doi: 10.3934/cpaa.2012.11.1825 [16] Wen-Bin Yang, Yan-Ling Li, Jianhua Wu, Hai-Xia Li. Dynamics of a food chain model with ratio-dependent and modified Leslie-Gower functional responses. Discrete & Continuous Dynamical Systems - B, 2015, 20 (7) : 2269-2290. doi: 10.3934/dcdsb.2015.20.2269 [17] Mats Gyllenberg, Jifa Jiang, Lei Niu, Ping Yan. On the classification of generalized competitive Atkinson-Allen models via the dynamics on the boundary of the carrying simplex. Discrete & Continuous Dynamical Systems - A, 2018, 38 (2) : 615-650. doi: 10.3934/dcds.2018027 [18] Davide La Torre, Simone Marsiglio, Franklin Mendivil, Fabio Privileggi. Public debt dynamics under ambiguity by means of iterated function systems on density functions. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021070 [19] Fumihiko Nakamura. Asymptotic behavior of non-expanding piecewise linear maps in the presence of random noise. Discrete & Continuous Dynamical Systems - B, 2018, 23 (6) : 2457-2473. doi: 10.3934/dcdsb.2018055 [20] Hyeong-Ohk Bae, Hyoungsuk So, Yeonghun Youn. Interior regularity to the steady incompressible shear thinning fluids with non-Standard growth. Networks & Heterogeneous Media, 2018, 13 (3) : 479-491. doi: 10.3934/nhm.2018021

2019 Impact Factor: 0.649