LBM D3Q19 LID DRIVEN CAVITY INSTABILITY

Hello,
i’m developing a 3d lid driven cavity case, but i’ve some problems with the code, the density diverges. So where could be the problem? I can’t understand it…( I imposed omega=0.3158 Re=4.5 xDim=yDim=zDim=61 uMax=0.1 nu=0.8888 tau=3.1666 Ma=0.1732)

Thanks for the help.