According to the work-energy theorem, the power delivered to a particle by the net force on it is equal to the time rate of change of its kinetic energy    ,
where is the net force acting on the particle. As a result, if the net force has a component in the direction of motion, the kinetic energy of the particle increases. Conversely, if the net force has a component in the opposite direction of motion, the kinetic energy of the particle decreases. In each case, the normal component of the net force changes the direction of motion of the particle. If the net force is strictly normal to the velocity vector of the particle, it only changes the direction of motion of the particle without affecting its kinetic energy and speed. To this end, it is stated that normal force does not do work    .
While developing a computer simulation, a moving particle was successively impacted perpendicular to its direction of motion to generate a circular path. During the simulation, it was noticed that the speed of the particle increased after each impact even though the impulsive force was perpendicular to the instantaneous velocity of the particle. This indicated that in this case the normal force did in fact do work on the particle. Although this is not a significant discovery by any scale, it warrants some discussion and should be brought to the attention of the computational science community.
In this article, a fairly detailed analysis of the dynamics of a particle subjected to a net force acting normal to its direction of motion is presented. The case of an impulsive force as well as the case of a continuous force, such as the magnetic force on a moving charged particle, or gravitational force on a circularly orbiting planet or satellite, is considered. The question of the change of kinetic energy and speed of the particle in each case is addressed, and it is shown that unless one is careful enough during the simulation, the normal force can do work on the particle and change its kinetic energy and speed.
To prevent the kinetic energy and speed from changing, one can always rescale the speed of the particle after each application of the normal force. This, however, is not a common practice in computational physics and, as a result, the kinetic energy and speed of the particle increase over time and error accumulates  . In this article, a different algorithm is suggested to avoid this problem. The algorithm is simple and more efficient than rescaling the speed in every step of the calculations.
2. The Impulsive Force
Let us begin by considering a particle of mass m moving with a constant velocity v along the x axis. At a certain time, an impulsive force F acts on the particle in the y direction for a short time interval . The velocity of the particle after the action of the force is given by
which simply indicates that the speed of the particle has increased to
This difficulty seems to result from the finite duration of the time interval over which the particle changes its direction of motion slightly, causing the force to gain a component in the direction of motion during this time interval. However, the problem can not be resolved even if the time interval is shrunk to zero and the force is reduced to a Dirac delta function. To see this, consider an impulsive force given by
where here denotes the Dirac delta function, and the force acts on the particle at time . In this equation, p is a constant having the units of momentum (as the delta function has the units of reciprocal time). The velocity of the particle after the action of the force is given by
where, again, the speed of the particle increases. The conclusion is that it is impossible to change the direction of motion of a particle by an impulsive force acting normal to its velocity vector without increasing its speed and kinetic energy. If the speed of the particle is to remain constant and only the direction of motion is to change, the impulsive force must act somewhat “backward”. To see this, consider a particle of mass m which is initially moving along the x axis and is subjected to an impulsive force that acts on it for a short time interval as it passes through the origin (Figure 1). The force acts at an angle backward from normal to the direction of motion as shown. This force causes the x com- ponent of velocity to decrease, but at the same time generates a y component. The velocity of the particle after the action of the force is given by
Now choose so that the speed of the particle remains constant, i.e.,
Therefore, if the particle is acted upon by an impulsive force that is not normal to its velocity vector, but makes a backward angle given by Equation (8), its speed and the kinetic energy remain unchanged. This force changes the direction of motion of the particle only. It is interesting to note that the new direction of motion makes an angle with the initial direction of motion. To see this, note that if we call this angle , then from Equation (6) we have
Figure 1. A particle of mass m moving initially along the x axis is subjected to an impulsive force that acts on it during a short time period as the particle passes through the origin. and are the velocities of the particle before and after the impact, respectively.
Substituting for from Equation (8), reduces Equation (9) to
which proves the statement. Furthermore, this result, together with Equations (2) and (8), reveal that the new direction of motion is the same as what it would have been had the force acted perpendicular to the direction of motion.
3. The Continuous Force
Now consider the case of a continuous force that acts on a particle perpendicular to its velocity vector at all times. Figure 2 shows the path of such a particle during an infinitesimal time interval starting when the particle is at the origin. The x axis is chosen to coincide with the velocity vector of the particle at the beginning of this time interval. The force does not need to have a constant magnitude, however, if we choose small enough then the magnitude of the force can be considered to be constant. The velocity vector of the particle at the end of the time interval is , which makes an angle with the x axis.
Let us use the force at the beginning of the infinitesimal interval to calculate the speed at the end of the interval (Figure 2). This gives
from which we obtain
Figure 2. A normal force acting on a particle during an infinitesimal time interval. The coordinate system is chosen so that at the beginning of the time interval the particle is passing through the origin with velocity in the x direction.
where the last equality is written because is infinitesimal and the higher- order terms are negligible. Since according to Equation (12), is second-order in , integration over time gives a first-order increase in the speed of the particle  . On the other hand, if we use the force at the end of the interval to calculate , we obtain
However, since and are both infinitesimal quantities, we can write
where the last equality can be obtain either from Equation (11), or from Equation (13) by neglecting the second- and higher-order terms in . Therefore, Equation (14) reduces to
where again because is infinitesimal, the higher-order terms are negligible. In this case, integration over time gives a first-order decrease in the speed of the particle. Therefore, using the initial value of the force over each interval slightly increases the speed of the particle, whereas using the final value of the force slightly decreases it.
Consider now the average of the two forces and over the infinitesi- mal interval,
Using this average force, the components of the velocity at the end of the time interval are given by
from which we find
which together give
and substituting for from Equation (15), which is also valid in this case as can be seen from Equation (18), reduces Equation (19) to
Finally, substituting for from Equation (15) gives
where the terms higher than 4th order have been neglected. Since is fourth order in , integration over time gives a change in the speed that is third-order in  . Consequently, the speed of the particle remains constant to a high order.
It is interesting to note that the angle that the average force in the infinitesimal interval makes with the y axis is exactly (Figure 2). To see this, let us denote the angle between the average force and the y axis by . Then from Equation (17) we have
which proves the statement. Furthermore, since is infinitesimal, from Equation (15) we have
which is the same result as Equation (8) for the backward angle in the case of the impulsive force. It should be noted, however, that Equation (8) is exact whereas Equation (25) is an approximation obtained for a continuous force under the assumption of infinitesimal .
To sum up, the motion of the particle under a continuous force that is always perpendicular to its velocity vector can be viewed as consisting of a succession of infinitesimal steps or intervals. At the beginning of each step, the particle is acted upon by a force that is not perpendicular to the direction of motion but slightly pointing backward by a backward angle that is given by Equation (8) or (25). As a result of the action of this force, the speed and the kinetic energy of the particle remain unchanged over the interval, but the direction of motion changes by an angle equal to twice the backward angle.
It is impossible to change the direction of motion of a particle by a normal impulsive force without changing its kinetic energy and speed. If the latter is to remain unchanged, the impulsive force must act at a backward angle from the normal to the direction of motion.
This situation is also realized in the case of a continuous normal force, where the dynamics of the particle can be viewed as a sequence of infinitesimal steps. At the beginning of each step an impulsive force acts on the particle according to the rule described above. This conserves the kinetic energy and speed of the particle to a high degree of accuracy and only changes its direction of motion. This algorithm is particularly useful when numerical calculations of dynamics of a particle under a continuously varying net normal force, such as a charged particle moving in a magnetic field or a particle moving under a central force field, are concerned. Since in all computer simulations and numerical calculations finite time intervals are involved, the above algorithm provides an efficient tool for achieving a high degree of constancy in the kinetic energy and speed of the particle. Simply using a normal force at the beginning or at the end of the finite intervals will result in an increase or decrease of the kinetic energy and speed of the particle, respectively.
Finally, as was pointed out in the Introduction section, one can always let the force act perpendicular to the direction of motion of the particle, then rescale its speed after the action of the force. However, this process is more involved and results in a less efficient algorithm. The goal of this article, however, is not to just find a way to keep the speed and kinetic energy of the particle constant. The objective here, among other things, is to show that when a force normal to the velocity of a particle acts on it in a time interval , it does work on it regardless of how short is, unless .