Scaling properties of the regular dynamics for a dissipative bouncing ball model

Some scaling properties of the regular dynamics for a dissipative version of the one-dimensional Fermi accelerator model are studied. The dynamics of the model is given in terms of a two-dimensional nonlinear area contracting map. Our results show that the velocities of saddle fixed points (saddle velocities) can be described using scaling arguments for different values of the control parameter. (c) 2007 Elsevier B.V. All rights reserved.

The Feigenbaum's delta for a high dissipative bouncing ball model

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

A bouncing ball model with two nonlinearities: a prototype for Fermi acceleration

Some dynamical properties of a bouncing ball model under the presence of an external force modelled by two nonlinear terms are studied. The description of the model is made by the use of a two-dimensional nonlinear measure-preserving map on the variable's velocity of the particle and time. We show that raising the straight of a control parameter which controls one of the nonlinearities, the positive Lyapunov exponent decreases in the average and suffers abrupt changes. We also show that for a specific range of control parameters, the model exhibits the phenomenon of Fermi acceleration. The explanation of both behaviours is given in terms of the shape of the external force and due to a discontinuity of the moving wall's velocity.

Consequences of quadratic frictional force on the one dimensional bouncing ball model

Some dynamical properties of the one dimensional Fermi accelerator model, under the presence of frictional force are studied. The frictional force is assumed as being proportional to the square particle's velocity. The problem is described by use of a two dimensional non linear mapping, therefore obtained via the solution of differential equations. We confirm that the model experiences contraction of the phase space area and in special, we characterized the behavior of the particle approaching an attracting fixed point. © 2007 American Institute of Physics.

The feigenbaumes δ for a high dissipative bouncing ball model

We have studied a dissipative version of a one-dimensional Fermi accelerator model. The dynamics of the model is described in terms of a two-dimensional, nonlinear area-contracting map. The dissipation is introduced via inelastic collisions of the particle with the walls and we consider the dynamics in the regime of high dissipation. For such a regime, the model exhibits a route to chaos known as period doubling and we obtain a constant along the bifurcations so called the Feigenbaum's number 8.

Some dynamical properties of a classical dissipative bouncing ball model with two nonlinearities

Some dynamical properties for a bouncing ball model are studied. We show that when dissipation is introduced the structure of the phase space is changed and attractors appear. Increasing the amount of dissipation, the edges of the basins of attraction of an attracting fixed point touch the chaotic attractor. Consequently the chaotic attractor and its basin of attraction are destroyed given place to a transient described by a power law with exponent -2. The parameter-space is also studied and we show that it presents a rich structure with infinite self-similar structures of shrimp-shape. © 2013 Elsevier B.V. All rights reserved.

Transport and dynamical properties for a bouncing ball model with regular and stochastic perturbations

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Global ballistic acceleration in a bouncing-ball model

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Feedback control of impact dynamics: The bouncing ball revisited

We study the the design of a tracking controller for the popular bouncing ball model: the continuous-time actuation of a table is used to control the impacts of the table with a bouncing ball. The proposed control law uses the impact times as the sole feedback information. We show that the acceleration of the table at impact plays no role in the stability analysis but is an important parameter for the robustness of the feedback system to model uncertainty, in particular to the uncertainty on the coefficient of restitution. © 2006 IEEE.

A theoretical characterization of scaling properties in a bouncing ball system

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Parameter space for a dissipative Fermi-Ulam model

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Breaking down the Fermi acceleration with inelastic collisions

The phenomenon of Fermi acceleration is addressed for a dissipative bouncing ball model with external stochastic perturbation. It is shown that the introduction of energy dissipation (inelastic collisions of the particle with the moving wall) is a sufficient condition to break down the process of Fermi acceleration. The phase transition from bounded to unbounded energy growth in the limit of vanishing dissipation is characterized.

Enclosing the behavior of a hybrid automaton up to and beyond a Zeno point

Even simple hybrid automata like the classic bouncing ball can exhibit Zeno behavior. The existence of this type of behavior has so far forced a large class of simulators to either ignore some events or risk looping indefinitely. This in turn forces modelers to either insert ad-hoc restrictions to circumvent Zeno behavior or to abandon hybrid automata. To address this problem, we take a fresh look at event detection and localization. A key insight that emerges from this investigation is that an enclosure for a given time interval can be valid independent of the occurrence of a given event. Such an event can then even occur an unbounded number of times. This insight makes it possible to handle some types of Zeno behavior. If the post-Zeno state is defined explicitly in the given model of the hybrid automaton, the computed enclosure covers the corresponding trajectory that starts from the Zeno point through a restarted evolution.

Particle fluid mechanics in shear flows, acoustic waves, and shock waves

Three different categories of flow problems of a fluid containing small particles are being considered here. They are: (i) a fluid containing small, non-reacting particles (Parts I and II); (ii) a fluid containing reacting particles (Parts III and IV); and (iii) a fluid containing particles of two distinct sizes with collisions between two groups of particles (Part V).

Part I

A numerical solution is obtained for a fluid containing small particles flowing over an infinite disc rotating at a constant angular velocity. It is a boundary layer type flow, and the boundary layer thickness for the mixture is estimated. For large Reynolds number, the solution suggests the boundary layer approximation of a fluid-particle mixture by assuming W = W_p. The error introduced is consistent with the Prandtl’s boundary layer approximation. Outside the boundary layer, the flow field has to satisfy the “inviscid equation” in which the viscous stress terms are absent while the drag force between the particle cloud and the fluid is still important. Increase of particle concentration reduces the boundary layer thickness and the amount of mixture being transported outwardly is reduced. A new parameter, β = 1/Ω τ_v, is introduced which is also proportional to μ. The secondary flow of the particle cloud depends very much on β. For small values of β, the particle cloud velocity attains its maximum value on the surface of the disc, and for infinitely large values of β, both the radial and axial particle velocity components vanish on the surface of the disc.

Part II

The “inviscid” equation for a gas-particle mixture is linearized to describe the flow over a wavy wall. Corresponding to the Prandtl-Glauert equation for pure gas, a fourth order partial differential equation in terms of the velocity potential ϕ is obtained for the mixture. The solution is obtained for the flow over a periodic wavy wall. For equilibrium flows where λ_v and λ_T approach zero and frozen flows in which λ_v and λ_T become infinitely large, the flow problem is basically similar to that obtained by Ackeret for a pure gas. For finite values of λ_v and λ_T, all quantities except v are not in phase with the wavy wall. Thus the drag coefficient C_D is present even in the subsonic case, and similarly, all quantities decay exponentially for supersonic flows. The phase shift and the attenuation factor increase for increasing particle concentration.

Part III

Using the boundary layer approximation, the initial development of the combustion zone between the laminar mixing of two parallel streams of oxidizing agent and small, solid, combustible particles suspended in an inert gas is investigated. For the special case when the two streams are moving at the same speed, a Green’s function exists for the differential equations describing first order gas temperature and oxidizer concentration. Solutions in terms of error functions and exponential integrals are obtained. Reactions occur within a relatively thin region of the order of λ_D. Thus, it seems advantageous in the general study of two-dimensional laminar flame problems to introduce a chemical boundary layer of thickness λ_D within which reactions take place. Outside this chemical boundary layer, the flow field corresponds to the ordinary fluid dynamics without chemical reaction.

Part IV

The shock wave structure in a condensing medium of small liquid droplets suspended in a homogeneous gas-vapor mixture consists of the conventional compressive wave followed by a relaxation region in which the particle cloud and gas mixture attain momentum and thermal equilibrium. Immediately following the compressive wave, the partial pressure corresponding to the vapor concentration in the gas mixture is higher than the vapor pressure of the liquid droplets and condensation sets in. Farther downstream of the shock, evaporation appears when the particle temperature is raised by the hot surrounding gas mixture. The thickness of the condensation region depends very much on the latent heat. For relatively high latent heat, the condensation zone is small compared with Ʌ_D.

For solid particles suspended initially in an inert gas, the relaxation zone immediately following the compression wave consists of a region where the particle temperature is first being raised to its melting point. When the particles are totally melted as the particle temperature is further increased, evaporation of the particles also plays a role.

The equilibrium condition downstream of the shock can be calculated and is independent of the model of the particle-gas mixture interaction.

Part V

For a gas containing particles of two distinct sizes and satisfying certain conditions, momentum transfer due to collisions between the two groups of particles can be taken into consideration using the classical elastic spherical ball model. Both in the relatively simple problem of normal shock wave and the perturbation solutions for the nozzle flow, the transfer of momentum due to collisions which decreases the velocity difference between the two groups of particles is clearly demonstrated. The difference in temperature as compared with the collisionless case is quite negligible.

An Energy Conserving Finite Difference Scheme for Simulation of Collisions

Nonlinear phenomena play an essential role in the sound production process of many musical instruments. A common source of these effects is object collision, the numerical simulation of which is known to give rise to stability
issues. This paper presents a method to construct numerical schemes that conserve the total energy in simulations of one-mass systems involving collisions, with no conditions imposed on any of the physical or numerical parameters.
This facilitates the adaptation of numerical models to experimental data, and allows a more free parameter adjustment in sound synthesis explorations. The energy preservedness of the proposed method is tested and demonstrated though several examples, including a bouncing ball and a non-linear oscillator, and implications regarding the wider applicability are discussed.