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??