I am developing a solver to simulate particle resuspension. The particles are initially deposited on a wall and then exposed to a flow field. I wish to calculate a coefficient x for each particle as a criterion. This coefficient needs to use the fluid velocity at the particle. I have tried to declare and define the variable x in my atom.h and atom.c solver files, how should I call the fluid velocity at the particle location at each time step?
I need to update this coefficient x in a time loop and return it to the liggghts, does anyone have a good suggestion? Thanks