-
N-body linear force law allowing analytic solutions
Authors:
Joseph West,
Sean P. Bartz
Abstract:
We present a pair-wise force law in a system of N particles that produces analytic solutions for arbitrary number of particles, masses, and initial conditions. Each pair of particles interacts via a force that is proportional to the product of their masses and their separation distance, with the force directed radially. We show that, despite the N-body interaction, each particle behaves as if it i…
▽ More
We present a pair-wise force law in a system of N particles that produces analytic solutions for arbitrary number of particles, masses, and initial conditions. Each pair of particles interacts via a force that is proportional to the product of their masses and their separation distance, with the force directed radially. We show that, despite the N-body interaction, each particle behaves as if it interacts only with the center of mass of the system. This effective two-body interaction behaves as Hooke's Law with a common frequency for all particles, with the familiar analytic solutions for the trajectories. With these analytic solutions, it is possible to efficiently simulate a collection of these particles and incorporate other external forces. As an example, we simulate the particles within an adiabatically expanding container and calculate pressure and temperature in both the attractive and repulsive cases.
△ Less
Submitted 2 July, 2024;
originally announced July 2024.
-
Numerical simulation of non-central collisions of spherical magnets
Authors:
Sean P. Bartz,
Jacob Shaw
Abstract:
We present a computational model of non-central collisions of two spherical neodymium-iron-boron magnets, suggested as a demonstration of angular momentum conservation. Our program uses an attractive dipole-dipole force and a repulsive contact force to solve the Newtonian equations of motion for the magnets. We confirm the conservation of angular momentum and study the changes in energy throughout…
▽ More
We present a computational model of non-central collisions of two spherical neodymium-iron-boron magnets, suggested as a demonstration of angular momentum conservation. Our program uses an attractive dipole-dipole force and a repulsive contact force to solve the Newtonian equations of motion for the magnets. We confirm the conservation of angular momentum and study the changes in energy throughout the interaction. Using the exact expression for the dipole-dipole force, including non-central terms, we correctly model the final rotational frequencies, which is not possible with a simple power-law approximation.
△ Less
Submitted 1 December, 2022;
originally announced December 2022.
-
Coefficient of restitution of a linear dashpot on a rigid surface
Authors:
Sean P. Bartz
Abstract:
The linear dashpot model is applied to a single ball bouncing on a rigid surface. It is shown that when gravity is included the coefficient of restitution depends on impact velocity, in contrast to previous work that ignored the effects of gravity. This velocity dependence is most pronounced at low impact velocities and high damping. Previous work has considered the ball to be in contact with the…
▽ More
The linear dashpot model is applied to a single ball bouncing on a rigid surface. It is shown that when gravity is included the coefficient of restitution depends on impact velocity, in contrast to previous work that ignored the effects of gravity. This velocity dependence is most pronounced at low impact velocities and high damping. Previous work has considered the ball to be in contact with the floor when the compression is nonzero, while other analysis terminates the collision earlier, to prevent an attractive force. We compare these models and propose a hybrid between the two. The hybrid model is successful in reproducing experimental results for a cart bouncing repeatedly on a spring.
△ Less
Submitted 22 April, 2022;
originally announced April 2022.
-
Delayed Rebounds in the Two-Ball Bounce Problem
Authors:
Sean P. Bartz
Abstract:
In the classroom demonstration of a tennis ball dropped on top of a basketball, the surprisingly high bounce of the tennis ball is typically explained using momentum conservation for elastic collisions, with the basketball-floor collision treated as independent from the collision between the two balls. This textbook explanation is extended to inelastic collisions by including a coefficient of rest…
▽ More
In the classroom demonstration of a tennis ball dropped on top of a basketball, the surprisingly high bounce of the tennis ball is typically explained using momentum conservation for elastic collisions, with the basketball-floor collision treated as independent from the collision between the two balls. This textbook explanation is extended to inelastic collisions by including a coefficient of restitution. This independent contact model (ICM), as reviewed in this paper, is accurate for a wide variety of cases, even when the collisions are not truly independent. However, it is easy to explore situations that are not explained by the ICM, such as swapping the tennis ball for a ping-pong ball. In this paper, we study the conditions that lead to a "delayed rebound effect," a small first bounce followed by a higher second bounce, using techniques accessible to an undergraduate student. The dynamical model is based on the familiar solution of the damped harmonic oscillator. We focus on making the equations of motion dimensionless for numerical simulation, and reducing the number of parameters and initial conditions to emphasize universal behavior. The delayed rebound effect is found for a range of parameters, most commonly in cases where the first bounce is lower than the ICM prediction.
△ Less
Submitted 7 March, 2022; v1 submitted 29 July, 2020;
originally announced July 2020.