[cojua8]

Hello,

I'm working with basicKinematicCollidingCloud and I'm trying to implement a new particle force: the Basset force (along with drag, lift and other forces).

For this purpose I will use this model https://arxiv.org/pdf/1008.0833 (or another one which needs practically the same input parameters. [Computation of the particle Basset force with a fractional-derivative approach])

So, I'm having the following difficulties:

1. I need to store all the past relative accelerations of the particle for every particle: .
2. The tracking algorithm of OpenFOAM modifies the in every iteration, so the models above becomes invalid or at least not accurate, since they used a fixed . Is there any way to fix the of the particle? or a way to modify the models (this requires more numerics knowledge obviously)