# 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] Weiwei Liu, Jinliang Wang, Yuming Chen. Threshold dynamics of a delayed nonlocal reaction-diffusion cholera model. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020316 [2] Manil T. Mohan. First order necessary conditions of optimality for the two dimensional tidal dynamics system. Mathematical Control & Related Fields, 2020  doi: 10.3934/mcrf.2020045 [3] Cuicui Li, Lin Zhou, Zhidong Teng, Buyu Wen. The threshold dynamics of a discrete-time echinococcosis transmission model. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020339 [4] Shao-Xia Qiao, Li-Jun Du. Propagation dynamics of nonlocal dispersal equations with inhomogeneous bistable nonlinearity. Electronic Research Archive, , () : -. doi: 10.3934/era.2020116 [5] Ebraheem O. Alzahrani, Muhammad Altaf Khan. Androgen driven evolutionary population dynamics in prostate cancer growth. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020426 [6] Vieri Benci, Marco Cococcioni. The algorithmic numbers in non-archimedean numerical computing environments. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020449 [7] Héctor Barge. Čech cohomology, homoclinic trajectories and robustness of non-saddle sets. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020381 [8] Ying Lin, Qi Ye. Support vector machine classifiers by non-Euclidean margins. Mathematical Foundations of Computing, 2020, 3 (4) : 279-300. doi: 10.3934/mfc.2020018 [9] Sergey Rashkovskiy. Hamilton-Jacobi theory for Hamiltonian and non-Hamiltonian systems. Journal of Geometric Mechanics, 2020, 12 (4) : 563-583. doi: 10.3934/jgm.2020024 [10] A. M. Elaiw, N. H. AlShamrani, A. Abdel-Aty, H. Dutta. Stability analysis of a general HIV dynamics model with multi-stages of infected cells and two routes of infection. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020441 [11] José Luis López. A quantum approach to Keller-Segel dynamics via a dissipative nonlinear Schrödinger equation. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020376 [12] Yangrong Li, Shuang Yang, Qiangheng Zhang. Odd random attractors for stochastic non-autonomous Kuramoto-Sivashinsky equations without dissipation. Electronic Research Archive, 2020, 28 (4) : 1529-1544. doi: 10.3934/era.2020080 [13] Pengyu Chen. Non-autonomous stochastic evolution equations with nonlinear noise and nonlocal conditions governed by noncompact evolution families. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020383 [14] Lin Shi, Xuemin Wang, Dingshi Li. Limiting behavior of non-autonomous stochastic reaction-diffusion equations with colored noise on unbounded thin domains. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5367-5386. doi: 10.3934/cpaa.2020242 [15] Dan Zhu, Rosemary A. Renaut, Hongwei Li, Tianyou Liu. Fast non-convex low-rank matrix decomposition for separation of potential field data using minimal memory. Inverse Problems & Imaging, , () : -. doi: 10.3934/ipi.2020076 [16] Bernard Bonnard, Jérémy Rouot. Geometric optimal techniques to control the muscular force response to functional electrical stimulation using a non-isometric force-fatigue model. Journal of Geometric Mechanics, 2020  doi: 10.3934/jgm.2020032

2019 Impact Factor: 0.649