
The trajectory equations, and any auxiliary equations describing heat or mass transfer to/from the particle, are solved by stepwise integration over discrete time steps. Integration of time in Equation 22.21 yields the velocity of the particle at each point along the trajectory, with the trajectory itself predicted by
Equations 22.21 and 22.151 are a set of coupled ordinary differential equations. Equation 22.21 can be cast into the following general form
where the term includes accelerations due to all other forces except drag force.
This set can be solved for constant , and by analytical integration. For the particle velocity at the new location we get
The new location can be computed from a similar relationship.
In these equations and represent particle velocities and fluid velocities at the old location. Equations 22.153 and 22.154 are applied when using the analytic discretization scheme.
The set of Equations 22.21 and 22.151 can also be solved using numerical discretization schemes. When applying the Euler implicit discretization to Equation 22.152 we get
When applying a trapezoidal discretization to Equation 22.152 the variables and on the right hand side are taken as averages, while accelerations, , due to other forces are held constant. We get
(22.156) 
The averages and are computed from
The particle velocity at the new location is computed by
For the implicit and the trapezoidal schemes the new particle location is always computed by a trapezoidal discretization of Equation 22.151.
Equations 22.152 and 22.151 can also be computed using a RungeKutta scheme which was published by Cash and Karp [ 50]. The ordinary differential equations can be considered as vectors, where the left hand side is the derivative and the right hand side is an arbitrary function .
We get
with
The coefficients , , and are taken from Cash and Karp [ 50]
This scheme provides an embedded error control, which is switched off, when no Accuracy Control is enabled.
For rotating reference frames, the integration is carried out in the rotating frame with the extra terms described in Equations 22.212 and 22.213, thus accounting for system rotation. Using the mechanisms available for accuracy control, the trajectory integration will be done accurately in time.
The analytic scheme is very efficient. It can become inaccurate for large steps and in situations where the particles are not in hydrodynamic equilibrium with the continuous flow. The numerical schemes implicit and trapezoidal, in combination with Automated Tracking Scheme Selection, consider most of the changes in the forces acting on the particles and are chosen as default schemes. The rungekutta scheme is recommended of nondrag force changes along a particle integration step.
The integration step size of the higherorder schemes, trapezoidal and rungekutta, is limited to a stable range. Therefore it is recommended to use them in combination with Automated Tracking Scheme Selection.