High Oscillating Frequency in LBM

Hi,

I wish to model moving boundary problem e.g. a 2D flapping airfoil at higher frequency (say >40Hz). I tried with BGK model but every time solution diverged despite increasing mesh size and decreasing time step. I am doing in my own code,

Any help regarding the problem??

Regards
Ali