Modeling Experimental Nonlinear Dynamics and Chaotic ScenariosView this Special Issue
Research Article | Open Access
Minvydas Ragulskis, Edita Sakyte, Jesús M. Seoane, Miguel A. F. Sanjuán, "A New Mechanical Model for Particle Transport by Surface Waves and Applications", Mathematical Problems in Engineering, vol. 2009, Article ID 731358, 17 pages, 2009. https://doi.org/10.1155/2009/731358
A New Mechanical Model for Particle Transport by Surface Waves and Applications
We present a study of the behavior of a ball under the influence of gravity on a platform. A propagating surface wave travels on the surface of the platform while the platform remains motionless. This is a modification of the classical bouncing ball problem and describes the transport of particles by surface waves. Phase and velocity maps cannot be expressed in the explicit form due to implicit formulations, and no formal analytical analyses is possible. Numerical analysis shows that the transition to chaos is produced via a period doubling routewhich is a common property for classical bouncers. These numerical analysis have been carried out for the conservative and for the viscous cases and also for elastic and for inelastic collisions. The bouncing process can be sensitive to the initial conditions and can be useful for control techniques which can dramatically increase the effectiveness of particle transport in practical applications. Finally, we also consider the mechanical model of a particle sliding on a surface which is also important because it has important physical implications such as the transportation of thin films in biomedical applications, among others.
A particle falling down, in a constant gravitational field, on a moving platform is called a bouncing ball problem, or a bouncer. This model was suggested more than thirty years ago [1, 2] as an alternative to the Fermi-Ulam model  of cosmic ray acceleration . In the ensuing years many approaches to the bouncer model have been studied theoretically and experimentally [5–8]. It has been proved to be a useful system for experimentally exploring several new nonlinear effects [9, 10]. Moreover, it has been implemented into a number of engineering applications [11, 12].
The bouncer model can be briefly characterized by the following basic statements. (i) Maps derived for the bouncer model can be exactly iterated for any time function describing the moving platform [7, 9] (though usually the platform is assumed to oscillate with a single frequency). (ii) The ball-platform collisions can be characterized by a coefficient of restitution changing from for a perfectly elastic case to for a completely inelastic situation. (iii) The chaotic bouncer can be easily used to relate theoretical predictions to experimental results, [9, 10] what makes it a paradigm model in nonlinear dynamics.
In this paper we assume that a particle is falling down in a constant gravitational field on a stationary platform. A propagating surface wave travels on the surface of the platform while the platform remains motionless. Such a model can be used to describe the transport of particles by propagating surface waves, which is an important problem with numerous applications. Powder transport by piezoelectrically excited ultrasonic surface waves , manipulation of bioparticles using traveling wave electrophoresis [14, 15], and conveyance of submerged buoys in coastal waters  are just a few examples of problems involving the interaction between propagating waves and transported bouncing particles.
This paper is organized as follows. In Section 2 we present a complete description of our model, the bouncer system. Section 3 presents numerical simulations showing that the modified bouncer model possesses such an inherent chaotic dynamics. These results are carried out for both, the conservative and the dissipative cases. The case of a particle sliding on the surface is fully analyzed in Section 4. Conclusions and discussions of the main results of this paper are presented in Section 5.
2. Model Description
We consider the two-dimensional system shown in Figure 1, where the surface of an elastic plate is represented by a solid line which coincides with the -axis in the state of equilibrium. A point of the surface in the state of equilibrium is translated to coordinates when a wave process takes place. This translation is sensitive to time and coordinate : where the functions and determine deflections from the state of equilibrium.
Explicitly, the longitudinal and transverse displacements of the medium at the surface of flat boundary with travelling Rayleigh wave can be expressed like  where and are the longitudinal and transverse displacements, the coordinate of the surface point of the medium before the wave process took place; is a constant, the wavenumber and the density. can be found from the following algebraic equation: and the angular velocity can be found from the following transcendental equation: where , , and , the Poisson’s ratio and the first Lamé constant. It can be noted that the ratio between the amplitudes of transverse and longitudinal deformations depends on . In usual elastic media it is quite normal that the transverse displacement is about 1.5 times larger than the longitudinal displacement . The motion of a point in the medium is an ellipse. Also, the direction of the velocity of the particles at the peaks of the wave is opposite to the direction of wave propagation.
Rayleigh waves are dispersive due to a dependence of the wave’s speed on its wavelength. Typical example is Rayleigh waves in the Earth where waves with a higher frequency travel more slowly than those with a lower frequency. Rayleigh waves thus often appear spread out on seismograms recorded at distant earthquake recording stations . Surface acoustic waves (SAWs) generated by SAW devices on rough anisotropic materials also experience considerable dispersion . On the other hand, film waves generated on a surface of a finite liquid bed  can be characterized by a single frequency wave component. Therefore we concentrate on a one single frequency steady-state Rayleigh wave propagation and disregard dispersion.
Whenever a traveling nondispersive Rayleigh surface wave occurs in a medium, it can be characterized by a retrograde elliptic motion of the particles of that medium: where and are longitudinal and transverse amplitudes of the oscillations; is the angular frequency, and is the wave number. Remind that in a usual elastic medium it is quite normal for the transverse displacement to be about 1.5 times larger than the longitudinal displacement .
The coordinates of the particle are denoted as . Assume that the particle is in contact with the surface at time moment , then the following constrain takes place: where is to be found from the following algebraic equality (where and are given and is the unknown):
In other words, the instantaneous shape of the surface cannot be described by an explicit function. Nevertheless, the tangent to the surface at the point with abscise can be expressed explicitly: where is the angle between the tangent and the -axis. Instantaneous velocities (- and -components) of the point of the surface in contact with the particle can be expressed as and accordingly.
The governing equations of motion of a particle in a free flight mode are where top dots denote full derivative by time, is the mass of the particle, is the coefficient of viscous damping of the media above the surface, and is the free fall acceleration. Initial conditions ; ; ; yield partial solutions:
The free flight stage continues until the particle collides with the surface. Unfortunately, it is impossible to determine the explicit time moment of the collision due to the fact that the instantaneous shape of the surface cannot be expressed by an explicit function. Instead, one has to use iterative numerical techniques in order to determine the exact moment of the bounce.
Localization of the root (the time moment of the collision) is performed using a time marching technique starting from the initial conditions until where is the time step; ; is the step number for which (4.8) is satisfied for the first time, and is the solution of (4.4) at fixed : and , are determined by (4.7). Solution of (4.9) also requires an iterative numerical algorithm.
When the root is localized in the interval , one needs to fine down the value of using an iterative computational algorithm. This iterative algorithm can be a most simple bisection method, though more sophisticated algorithms comprising the golden section rule or Newton's iterations for example can be used instead until the desirable accuracy is achieved. As the collision moment is fined down in every iteration, the coordinate (corresponding to the collision point ) is also made more precise. Initially, ; every iteration helps to reach a better accuracy.
Such iterative method of determination of the collision moment leads us to the important conclusion that phase and velocity maps cannot be expressed in an explicit form, and no formal analytical analysis is possible.
Nevertheless, the geometrical coordinates of the point of collision are and can be reconstructed using computational techniques. Velocities of the particle just before the collision are and . Similarly, instantaneous velocities of the surface in contact with the particle are and .
Projections of the particle’s velocities just before the collision to the normal and to the tangent to the surface at the contact point can be expressed in the following form: where the angle is determined from (4.5) at the point of collision.
Analogously, projections of velocities of the point of the surface in contact with the particle to the normal and to the tangent take the following form:
Then, the velocities of the particle just after the collision (in the normal and tangent directions) are where is the coefficient of restitution for the collision in the normal direction. This constant is a measure of the energy loss at each impact. For elastic collisions , and for inelastic collisions. Coefficient determines the friction between the particle and the surface at the moment of collision. There is no friction between the particle and the surface when . The utmost value represents the situation when the projection of the particle's velocity (immediately after the impact) and the projection of the surface's point velocity to the tangent are equal.
The free flight stage starts over again immediately after the collision, and the initial conditions are
The presented model is a modification of the classic bouncer model which can be derived assuming and . In that case and the model becomes explicit.
3. Complex Dynamics in the Modified Bouncer Model
We will demonstrate that the modified bouncer model possesses such an inherent feature as chaotic dynamics. Moreover, we will show that the sensitivity to initial conditions can be exploited for the control of the process of conveyance. We will show these results for the conservative or nonviscous case () and for the viscous case ().
3.1. Nonviscous Case
We take , for which the media above the surface is non-viscous, (elastic collisions) and (no damping generated by sliding).
The dynamics of a bouncing particle on a surface of a propagating wave is very sensitive to the initial conditions if the dynamics is Hamiltonian. Apparently, it is possible to find such a set of initial conditions which lead to regular and periodic dynamics. This is illustrated in Figure 2 where collision heights are plotted versus initial velocity . successive collisions are used for every discrete value of to produce this diagram. The initial condition produces a period motion at , , , , , , and, . This is illustrated in Figure 3.
We plot the trajectory of the particle in for better visual interpretation. Collision moments are marked as black dots. At every moment of collision we also plot the instantaneous shape of the surface (one can note that the instantaneous shape of the surface is not harmonic). One can clearly see the difference in the complexity of the particle dynamics at (see Figure 3) and (see Figure 4).
3.2. Viscous Case
For this case, we assume that collisions are completely elastic (), and there is no tangential friction between the particle and the surface (), but the media above the surface is viscous fixing the value of as in . We use the reduced impact representation, where the height of the bouncing ball is sampled at each impact with the surface (impact sampling). Since the system is dissipative, we plot the bouncing process after the initial transients cease down (see Figure 5(a)). We skip 1500 successive bounces before starting to plot the collision heights for every discrete value of . Parameter is varied following the rule ; , which helps to expand the cascade of period doubling bifurcations. The control parameter in our case is not the amplitude of the platform’s oscillation but the velocity of the wave propagation; the collision height is used instead of collision velocity for a reduced impact representation. Moreover, the media above the surface of the plate is viscous. It appears that the transition to chaos via a period doubling route observed for a classical bouncer  is observed also for a particle bouncing on a surface of a stationary platform with a propagating wave traveling on its surface.
A phenomenological model could be used to exemplify the bifurcation diagram presented in Figure 5(a). The logistic map  is probably the simplest model ever used to study the transition to chaos via a period doubling route. Simple computational experiments with appropriately chosen parameter values of the logistic map would illustrate the universality of the bifurcation diagram in Figure 5(a).
An important parameter characterizing the effectiveness of the transport is the average longitudinal velocity of the particle . We average it over a long period of time after the initial transients cease down. In order to calculate a nondimensional quantity we divide it from the velocity of the traveling wave . Thus, the average velocity of conveyance is equal to the velocity of the traveling wave if is equal to as shown in Figure 5(b).
It is interesting to observe that the particle is transported with the average velocity of the traveling wave until the period 3 bouncing mode after a cascade of period doubling bifurcations (see Figure 5(a)). The particle’s average transportation velocity drops down only when the period 3 bouncing mode experiences its own cascade of period doubling bifurcations. External damping forces acting to the particle prevent its motion with the average wave’s velocity in the direction of the wave propagation when this velocity becomes large enough (even though the collisions are elastic). Also, the bouncing process is insensitive to initial conditions–eventually it converges to the one and only attractor shown in Figure 5(a) (at fixed ).
Figures 6(a)–6(d) show different dynamical behaviors of the transient processes for the elastic case () once we fixed the parameter values as follows: , , , , and . We obtain both, periodic and chaotic motions depending on the value of the parameter . Figures 6(a)–6(c) show for , and period , period and period processes, respectively. Chaotic bouncing for the value is shown in Figure 6(d).
The situation becomes different when collisions are inelastic, as shown in . In the presence of inelastic collisons, a vanishing bouncing process takes place (complete chattering ) when the particle sets into the state of rest on a slope of the propagating wave. The term complete chattering is used in literature to describe the process when the time interval between inelastic bounces tends to zero and the ball finally “sticks” to the surface of the oscillating platform. As we mentioned previously, a complete description of this phenomenon is given in .
4. The Sliding Particle Model
In this section we thoroughly analyze the case in which the particle is sliding on the surface instead of the case in which it is falling down on it, analyzed previously. Our motivation is the following. Conveyance of particles and bodies by propagating waves is an important scientific and engineering problem with numerous applications. Manipulation of bioparticles and gene expression profiling using traveling wave dielectrophoresis [14, 25, 26], segregation of particles in suspensions subject to ac electric fields , transport of sand particles and oil spills in coastal waters [28, 29], powder transport by piezoelectrically excited ultrasonic waves [13, 22], transportation of thin films in biomedical applications  are just a few examples of problems involving interaction between propagating waves and transported objects.
We now describe, as in the bouncer model, the equations of motions of our sliding particle model.
It is assumed that a mass particle is in contact with the deformed surface at a point at a time moment (Figure 7). A point of the surface in the equilibrium state is translated to coordinates at time moment . This translation is sensitive to time and coordinate : where and are predefined functions.
The condition that the particle is located on the surface leads to the following constraint: where is to be found from the following algebraic equality (in which is given and is unknown):
Though the instantaneous shape of the oscillating surface cannot be described by an explicit function, the tangent to the surface at the point can be expressed as
Instantaneous velocities of the surface's point in the direction of - and -axis can be expressed as follows: where dots denote derivatives by .
When a mass particle slides on the surface, it does not necessarily move in contact with one point of the surface. Therefore is no longer a constant. Thus,
The condition that the mass particle continuously slides on the surface brings another constraint into force (the relative velocity in the normal direction to the surface at the contact point must be zero):
Differentiation of (4.9) yields
Then the relative sliding velocity of the particle on the surface can be expressed as
The condition of dynamic equilibrium leads to the following system of equations:
where is the mass of the particle; is the reaction force at the contact point; is the gravity acceleration; is the friction force between the mass particle and the surface. The system of equations in (4.12) is in force when . Otherwise the particle jumps off the oscillating surface.
It is assumed that the friction force is linear. Thus can be expressed like where is the coefficient of linear friction.
Finally, the governing equation of motion can be derived from (4.12). Elementary transformations and substitutions lead to the following explicit differential equation: where
A major obstacle is eliminated, and direct numerical time marching techniques can be used for integration of (4.14)–computation of and is straightforward if the coordinate is given at time (4.1). Existence of a stability of the dynamic equilibrium can be analyzed explicitly.
But before proceeding with the analysis of dynamic equilibrium the following observation can be done. If kinematic relationships describing a traveling Rayleigh wave are in force, the change of variables transforms (4.14) to the following autonomous form: where
An important conclusion can be done. Dynamics of a particle sliding on the surface of a propagating Rayleigh wave cannot be chaotic. This is due to the fact that the governing equation of motion is a second-order autonomous ordinary differential equation with smooth parameter functions .
Equation (4.8) yields the dynamic equilibrium which represents a motion of the particle on a slope of the propagating wave with the velocity of its propagation: where is a constant. Then, it follows from (4.4) that
The term is bounded, therefore (4.20) will be in force when where is a constant satisfying the equality . Moreover, conditions of existence of the dynamic equilibrium are similar in terms of or :
As mentioned earlier, the explicit governing equation is formulated in terms of , not the coordinate of the contact point . First, coordinates of the unstable saddle point are determined. Then coordinates of the same saddle point are calculated in the frame using the relationship in (4.3). Forward and reverse time marching techniques are used to construct basin boundaries of attractors when partial solutions of (4.14) are sought from the infinitesimal surrounding of the saddle point.
The described computational technique is used to construct basin boundaries of the system’s attractors (Figure 8). Solutions in terms of (forward and reverse) are visualized only. It can be noted that two stable attractors can coexist–a stable equilibrium point and a stable limit cycle. Shaded regions in Figure 8 correspond to a basin (attracting set of initial conditions) of stable equilibrium points, while white region corresponds to a basin of the limit cycle. The phase plane in Figure 8 is periodic by 2 and can be visualized in cylindrical coordinates, but the plane representation is clearer.
A special attention should be paid to dashed line intervals on basin boundaries. Equation (4.14) describes a motion of a particle on the surface of a propagating wave. This governing equation of motion holds until the reaction force in (4.12) is positive. Whenever gets equal or lower than zero, the particle looses a contact with the surface and starts a free fly in a gravitational field until it bounces on the surface again. Therefore, the moment when the particle looses the contact with, the surface is detected, and the trajectory is marked by a dashed line. It can be noted that such motions occur only at relatively high particle velocities (Figure 8).
Conveyance of a particle by a propagating Rayleigh wave is a nonlinear problem, so such effects as the coexistence of stable attractors should not be astonishing. Stable equilibrium point type attractor in Figure 8 corresponds to a surf-type motion on a slope of a propagating wave; stable limit cycle corresponds to a motion with an average velocity much lower than the velocity of the propagating wave. Coexistence of attractors (a stable equilibrium point and a stable limit cycle) enables development of motion control strategies based on a small external impulses which can bring the system from the regime of motion with small average velocity into motion with the propagating wave’s velocity . Such attractor control strategy is illustrated in Figure 9 where the particle first oscillates in the limit cycle, and then a small external impulse kicks it to the basin of attraction of the stable focus point.
It can be noted that the up-mentioned control strategy can be implemented only when the stable equilibrium point and the stable limit cycle coexist. Thus, it would be impossible to transport a sand particle with the velocity of the propagating wave by an acoustic surface Rayleigh wave. Nevertheless, such attractor control strategies could be implemented for transportation of biomedical objects on the surface of an undulation film . The sliding particle model presented in this section also exhibits a very rich dynamics as in the case of the bouncer model. In particular, the sliding particle model should also have the sensitivity to the initial conditions for certain sets of parameter values as occuring in the bouncer model (see ). The sensitivity of transient processes to initial conditions takes place for both, the bouncer model and the sliding particle model.
5. Conclusions and Discussion
Transport of particles by surface waves is an important scientific and engineering problem, with numerous practical applications, including MEMS (micro-electro-mechanical systems) used to manipulate objects like particles or cells. We show that this problem is a modification of the classical bouncer model which is considered as a paradigm model in nonlinear physics. The formulations of our model are implicit, thus phase and velocity maps cannot be expressed in explicit form.
Chaotic dynamics of a conveyed particle is not an unexpected fact due to the complexity of the constitutive model. More surprising is the rich dynamical behavior in models comprising dissipative dynamics, elastic and inelastic collisions. It appears that the transition to chaos via a period doubling route is a universal property for bouncers and is observed in our model of particles transport in both, conservative and viscous media. Moreover, the sensitivity to initial conditions can be useful for control techniques which can dramatically increase the effectiveness of particles transport by surface waves. These results are relevant in the sense that we have also found the sensitivity to the initial conditions for the sliding particle model, which may have important applications in practical implementations as powder transport by piezoelectrically excited ultrasonic waves, transport of sand particles, among others.
Though the numerical analysis was concentrated on the dimensionless system only, theoretical and experimental investigation of dry particle conveyance and its control is a definite object for future research.
This work was supported by the Spanish Ministry of Education and Science under project number FIS2006-08525 and by Universidad Rey Juan Carlos and Comunidad de Madrid under project number URJC-CM-2007-CET-1601.
- G. M. Zaslavsky, “The simplest case of a strange attractor,” Physics Letters A, vol. 69, no. 3, pp. 145–147, 1978.
- L. D. Pustilnikov, “On Ulam's problem,” Theoretical and Mathematical Physics, vol. 57, p. 1035, 1983.
- M. A. Lieberman and A. J. Lichtenberg, “Stochastic and adiabatic behavior of particles accelerated by periodic forces,” Physical Review A, vol. 5, no. 4, pp. 1852–1866, 1972.
- E. Fermi, “On the origin of the cosmic radiation,” Physical Review, vol. 75, no. 8, pp. 1169–1174, 1949.
- B. V. Chirikov, “A universal instability of many-dimensional oscillator systems,” Physics Reports, vol. 52, no. 5, pp. 264–379, 1979.
- A. J. Lichtenberg, M. A. Lieberman, and R. H. Cohen, “Fermi acceleration revisited,” Physica D, vol. 1, no. 3, pp. 291–305, 1980.
- P. Pieranski and J. Malecki, “Noisy precursors and resonant properties of the period-doubling modes in a nonlinear dynamical system,” Physical Review A, vol. 34, no. 1, pp. 582–590, 1986.
- R. M. Everson, “Chaotic dynamics of a bouncing ball,” Physica D, vol. 19, no. 3, pp. 355–383, 1986.
- K. Wiesenfeld and N. B. Tufillaro, “Suppression of period doubling in the dynamics of a bouncing ball,” Physica D, vol. 26, no. 1–3, pp. 321–335, 1987.
- Z. J. Kowalik, M. Franaszek, and P. Pierański, “Self-reanimating chaos in the bouncing-ball system,” Physical Review A, vol. 37, no. 10, pp. 4016–4022, 1988.
- M.-O. Hongler, P. Cartier, and P. Flury, “Numerical study of a model of vibro-transporter,” Physics Letters A, vol. 135, no. 2, pp. 106–112, 1989.
- M.-O. Hongler and J. Figour, “Periodic versus chaotic dynamics in vibratory feeders,” Helvetica Physica Acta, vol. 62, no. 1, pp. 68–81, 1989.
- M. Mracek and J. Wallaschek, “A system for powder transport based on piezoelectrically excited ultrasonic progressive waves,” Materials Chemistry and Physics, vol. 90, no. 2-3, pp. 378–380, 2005.
- M. S. Talary, J. P. H. Burt, J. A. Tame, and R. Pethig, “Electromanipulation and separation of cells using travelling electric fields,” Journal of Physics D, vol. 29, no. 8, pp. 2198–2203, 1996.
- C.-F. Chou, J. O. Tegenfeldt, O. Bakajin et al., “Electrodeless dielectrophoresis of single- and double-stranded DNA,” Biophysical Journal, vol. 83, no. 4, pp. 2170–2179, 2002.
- R. H. Plaut, A. L. Farmer, and M. M. Holland, “Bouncing-ball model of ‘dry’ motions of a tethered buoy,” Journal of Vibration and Acoustics, vol. 123, no. 3, pp. 333–339, 2001.
- J. D. Achenbach, Wave Propagation in Elastic Solids, Elsevier, New York, NY, USA, 1984.
- L. D. Landau and E. M. Lifshitz, Theory of Elasticity, Pergamon Press, Oxford, UK, 1986.
- K. Aki and P. G. Richards, Quantitative Seismology, Freeman, New York, NY, USA, 1980.
- C. M. Flannery and H. Von Kiedrowski, “Dispersion of surface acoustic waves on rough anisotropic materials,” in Proceedings of the IEEE Ultrasonics Symposium, vol. 1, pp. 583–586, 2001.
- M. Ragulskis and K. Koizumi, “Applicability of attractor control techniques for a particle conveyed by a propagating wave,” Journal of Vibration and Control, vol. 10, no. 7, pp. 1057–1070, 2004.
- M. Ragulskis and M. A. F. Sanjuán, “Transport of particles by surface waves: a modification of the classical bouncer model,” New Journal of Physics, vol. 10, Article ID 083017, 2008.
- J. M. Luck and A. Mehta, “Bouncing ball with a finite restitution: chattering, locking, and chaos,” Physical Review E, vol. 48, no. 5, pp. 3988–3997, 1993.
- R. M. May, “Simple mathematical models with very complicated dynamics,” Nature, vol. 261, no. 5560, pp. 459–467, 1976.
- C.-F. Chou, J. O. Tegenfeldt, O. Bakajin et al., “Electrodeless dielectrophoresis of single- and double-stranded DNA,” Biophysical Journal, vol. 83, no. 4, pp. 2170–2179, 2002.
- L. Cui and H. Morgan, “Design and fabrication of travelling wave dielectrophoresis structures,” Journal of Micromechanics and Microengineering, vol. 10, no. 1, pp. 72–79, 2000.
- A. D. Dussaud, B. Khusid, and A. Acrivos, “Particle segregation in suspensions subject to high-gradient ac electric fields,” Journal of Applied Physics, vol. 88, no. 9, pp. 5463–5473, 2000.
- W. N. Hassan and J. S. Ribberink, “Transport processes of uniform and mixed sands in oscillatory sheet flow,” Coastal Engineering, vol. 52, no. 9, pp. 745–770, 2005.
- S. D. Wang, Y. M. Shen, and Y. H. Zheng, “Two-dimensional numerical simulation for transport and fate of oil spills in seas,” Ocean Engineering, vol. 32, no. 13, pp. 1556–1571, 2005.
- R. C. Hilborn, Chaos and Nonlinear Dynamics, Oxford University Press, New York, NY, USA, 1994.
Copyright © 2009 Minvydas Ragulskis et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.